BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	36,315	G→A	56.2%	Y79Y (TAC→TAT)	USA300HOU_0026 ←	IS431mec transposase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	36,315	0	G	A	56.2%	231.2 / inf	1187	Y79Y (TAC→TAT)	USA300HOU_0026	IS431mec transposase
Reads supporting (aligned to +/- strand): ref base G (244/264); new base A (298/361); total (556/631)
Fisher's exact test for biased strand distribution p-value = 3.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.47e-01

Only 100 of 1200 total aligned reads displayed.

TTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTC  >  CP000730/36175‑36455

ttACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATg                                                                                                                                              <  2:1101521/141‑1 (MQ=32)

 tACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAATTAGCTCCATTTTCCTTTTATTTTGATGt                                                                                                                                             >  2:736489/1‑141 (MQ=17)

  aCCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTa                                                                                                                                            <  1:309796/141‑1 (MQ=32)

      aaTTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGtc                                                                                                                                        <  2:62015/141‑1 (MQ=12)

        ttGTTTAATTAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTTCGCAACCAATTATCTAATGTATGTTCCTCTCCATCACTGGCACGATATAAATAGCTCCATTTTCCTTTTATTCTGATGTACGtctc                                                                                                                                      <  2:467487/141‑1 (MQ=11)

          gTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCtcat                                                                                                                                    >  1:835750/1‑141 (MQ=17)

             tAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCaa                                                                                                                                 >  1:622761/1‑141 (MQ=17)

             tAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGGTGTACGTCTCATCaa                                                                                                                                 >  2:188214/1‑141 (MQ=17)

             tAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCaa                                                                                                                                 >  1:394981/1‑141 (MQ=32)

             tAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGGCtcctcat                                                                                                                                 >  1:37822/1‑140 (MQ=17)

              actGAGACGTTTGATACACGCATCTGCTGAATGATTATCTCGTTGATTACGCAACCAAATATATAATGTATGCCCCTCTGTATCAATGGCACGATATAAATAGCACCATTTTCCTTTTATTTTGATATTCGTCTCATCAAt                                                                                                                                >  2:316178/3‑141 (MQ=11)

              aaTGAGACGTTTGATAAACGCATATGCTGAATGATCATCCCGTTGCTGACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAAt                                                                                                                                <  1:852298/141‑1 (MQ=12)

               aTGTGACGTTTGATAAACGCATATGTTGAATGATTCTCTCGTTGCTTACGCAACCAAATATCTACTGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATa                                                                                                                               <  2:930022/141‑1 (MQ=11)

               aTGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATa                                                                                                                               <  1:177143/141‑1 (MQ=12)

                 gagaCGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACg                                                                                                                             <  1:949352/141‑1 (MQ=12)

                 gagaCGATTGAAAAACGCATATACTAAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACg                                                                                                                             <  1:963197/141‑1 (MQ=11)

                  ggactgttGAGAAGCGCATATGCTGAATGATTATCTGTTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGGATCAATGGCACGATATAAAGAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGc                                                                                                                            <  2:602737/135‑1 (MQ=11)

                  agaCGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGc                                                                                                                            <  1:892819/141‑1 (MQ=32)

                       tttGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCAttt                                                                                                                       >  2:164307/1‑141 (MQ=32)

                        ttGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTg                                                                                                                      >  2:635207/1‑141 (MQ=17)

                        ttGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTAATTTCGATCTAATGTTTTTCAATACGCCACTTg                                                                                                                      >  2:225405/1‑141 (MQ=11)

                         tGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGt                                                                                                                     >  2:1089376/1‑141 (MQ=17)

                           ggacagcgctatGCGGAATGATTCGCTCGTTGTCTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGtaa                                                                                                                   <  2:959655/132‑1 (MQ=11)

                           aTAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTACATTTTCCTTTTATTTTGATATACGTCTAATCAATACGCCATTTGtaa                                                                                                                   >  2:578137/1‑141 (MQ=255)

                              aaCGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGtaataa                                                                                                                <  1:158538/141‑1 (MQ=18)

                               aCGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGTTTTACGTCTCATCACTACGCCATTTGTAATAAg                                                                                                               >  1:90386/1‑141 (MQ=255)

                                cgtatttGCTGAATGATTATGTTGTTGCTTACGCAACAAAATATCTAAGGTATGTCCCTCAGTATCGCGTGCTGGAGATTTATCGCTCGATTTGCCTTTTATTTTAATATACGTCTCATCAATACGCCATTTGCAATAAGc                                                                                                              <  1:1091942/135‑1 (MQ=2)

                                   atatGCTGAATGATTATCTCGTTGCTTACGCAACCATATATCTAATGCATGTCCCTCTGTATCAATGGCACGAGATTAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAACACGCCATTTGTAAGAAGCttt                                                                                                           <  1:695236/141‑1 (MQ=11)

                                    tatGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCtttt                                                                                                          <  1:352099/141‑1 (MQ=21)

                                    tatGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTTTTTTGTTTTACGTCTACTCAATCTGCCATTTGTACGACGCtttt                                                                                                          >  1:135331/1‑141 (MQ=11)

                                      tGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATTCGTCTCATCAATACGCCATTTGTATTTAACtttttt                                                                                                        >  1:243865/1‑141 (MQ=11)

                                      tGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCTTCAATGGCTCGTTATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGCCTCATCAATACGCCATTTGTAATAAGCtttttt                                                                                                        >  1:78107/1‑141 (MQ=255)

                                       gCTGAATGATTATCTCTTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCTTTTGTAATAAGCttttttt                                                                                                       >  1:551681/1‑140 (MQ=255)

                                       gCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCAATTTGTAATAAGCTTTTTTa                                                                                                       <  1:669237/141‑1 (MQ=38)

                                        cTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGATCCATTTTCCTTTTATTTTGTTTTTCGTCTTCTCAATACGTCTTTTGTTATAAGCtttttttt                                                                                                      >  1:743331/1‑139 (MQ=11)

                                        cTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTAt                                                                                                      <  2:977924/141‑1 (MQ=38)

                                          gcaTGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATTCATGGCCCGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGc                                                                                                    <  2:1060485/139‑1 (MQ=14)

                                          gAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGc                                                                                                    >  2:782287/1‑141 (MQ=255)

                                           aaTGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTTTGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCt                                                                                                   <  2:883362/141‑1 (MQ=255)

                                            aTGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCtt                                                                                                  >  2:1078894/1‑141 (MQ=40)

                                              gATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCtttt                                                                                                >  2:746228/1‑141 (MQ=255)

                                               aTTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCttttt                                                                                               >  2:892978/1‑141 (MQ=40)

                                                ttATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGGTATACGTCTCTTAAATATACCATTTGTAGTTAGATTTTTTATGGtttttt                                                                                              >  1:282703/1‑140 (MQ=11)

                                                ttATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTTATCAATAAGCCATTTGTAATAATCTTTTTTATGCTTTttc                                                                                              >  1:1029352/1‑141 (MQ=14)

                                                ttATCGCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTttc                                                                                              <  1:719676/141‑1 (MQ=255)

                                                  ttctcGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCACTGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCGTTTGTAATAAGCTTTTTTATGCTTTctctt                                                                                            <  2:671721/140‑1 (MQ=14)

                                                  aTCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTttctt                                                                                            <  1:653798/141‑1 (MQ=40)

                                                   tctcGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATGTGTAATAAGCTTTTTTATGCTTTttcttc                                                                                           >  2:864915/1‑141 (MQ=255)

                                                    ctcGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTcc                                                                                          <  1:871182/141‑1 (MQ=255)

                                                         tGCGTAAGCAACGAGATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTCATAAGCTTTTTTATGCTTTTTCTTCCAAAtt                                                                                     >  2:838133/1‑141 (MQ=14)

                                                           cgtaagCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAAtttt                                                                                   >  2:717009/6‑140 (MQ=255)

                                                           cTTACGCAACCAAATATATAATGTATGTCACTCTGCATCAATGCCACGATATAAATAGCACCATTTTCCTTTTATTTTGATGTACGTCTGATCAATACGCCATTTTTAATAAGCTTTTTTATGCTTTTTCTTTCAAttttt                                                                                   >  2:1002335/1‑140 (MQ=14)

                                                              aCGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATa                                                                                >  2:590383/1‑141 (MQ=255)

                                                                  aaCCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCGAAATTTGATATaaa                                                                            <  1:575817/141‑1 (MQ=11)

                                                                   aCCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATCCaaaa                                                                           >  2:39121/1‑141 (MQ=255)

                                                                     cAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCCTTTGTAATAAGCTTTTTTATGTCTTTTCTTCAATTTTTGATACaacctt                                                                         >  2:1031824/1‑137 (MQ=14)

                                                                       aaTTTCCAATGTATGTTCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTgg                                                                       <  2:761138/141‑1 (MQ=11)

                                                                        atatCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTTATGTACGTCTCTTCATTACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTggg                                                                      >  1:422972/1‑141 (MQ=255)

                                                                         tatCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTgggg                                                                     <  1:200023/141‑1 (MQ=255)

                                                                         tatCTAATCTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATCTACGTCTCATCAATACGCCATTTGTAATAAGCTTGTTTATGCTTTTTCTTCCAAATTTTATATACAAttagcg                                                                     >  2:210165/1‑137 (MQ=11)

                                                                           tCTAATGTATGTCCCTCTGTATCACTGGCACGATATCAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTACCAATTTTAATATaaccttgagcta                                                                   >  2:428258/1‑131 (MQ=11)

                                                                           tCTAATGTATCTCCATCTGCATCACTGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCa                                                                   <  1:112689/141‑1 (MQ=255)

                                                                            cTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTTGGGCat                                                                  >  2:441653/1‑141 (MQ=255)

                                                                                tgtatgtCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTc                                                              <  1:220196/141‑1 (MQ=255)

                                                                                   atgtCCCTCTGTATCAAGGGCACGATATAAATAGCTCCCTTTTCCTTTTATTTTGATATACGGCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATAGTCTTg                                                           >  2:251546/1‑141 (MQ=11)

                                                                                        ccTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTTAAccc                                                      >  1:993370/1‑141 (MQ=21)

                                                                                          tctGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATCTAAAATTTGGGCATATTTTTGATACCaa                                                    >  1:480575/1‑141 (MQ=11)

                                                                                           ctGTATCAATGGCACGATATAAATCGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATCCTTTTTCTTTCAAATTTGATATAAAATTGGGGCCTGTTCCTgggatccga                                                   >  1:771372/1‑133 (MQ=11)

                                                                                            tGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAAGAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACg                                                  >  2:799415/1‑141 (MQ=21)

                                                                                            ggtgtCAATGGCACGATATAAAGAGCTCCATTTTCCTTTTATTTTGATAGACGGCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTGTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACg                                                  <  2:555541/137‑1 (MQ=11)

                                                                                             gtaTCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACgg                                                 <  1:922816/139‑1 (MQ=255)

                                                                                                tCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACTTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTTTTCCAATTTTGATATAAAATTGGGGCATTTTCTTTAACCCAAAGGTAg                                              >  2:70008/1‑141 (MQ=11)

                                                                                                       cACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTg                                       <  2:537936/141‑1 (MQ=255)

                                                                                                          gATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAat                                    >  1:256977/1‑141 (MQ=255)

                                                                                                           atatAAATAGCTTCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCGTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAatg                                   <  1:824597/141‑1 (MQ=255)

                                                                                                           atatAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAatg                                   <  2:293196/141‑1 (MQ=18)

                                                                                                           atatAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAatg                                   <  2:825158/141‑1 (MQ=21)

                                                                                                           agataAATAGCTCCATTTTCCTTTTATGTTGATGTACGTCTCATCAATACGCCATGTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAatg                                   <  2:150735/139‑1 (MQ=255)

                                                                                                            tataAATAGCTGCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAatga                                  <  1:459574/141‑1 (MQ=21)

                                                                                                               aaaTAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAAtgatga                               >  1:501211/1‑141 (MQ=255)

                                                                                                                aaTAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTTTTTTTTTTTTTATAAAAAATTGGGGCAAATTTTTGAATCCAGCGGTGGTCCGTTATAtgaggat                              >  2:514763/1‑140 (MQ=11)

                                                                                                                   atcTCCCTTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGt                           <  2:333412/139‑1 (MQ=11)

                                                                                                                   aGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGCATGTTGAACGt                           >  2:478778/1‑141 (MQ=11)

                                                                                                                      tcctttttccttttATTTTTATATACGTCTCATCAATACTCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTTATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTa                        <  1:1055908/141‑1 (MQ=11)

                                                                                                                       cctttttCCTTTTATTTTGATATACGCCTCGTCAATACGACATTTGTAATCAGGTCTTTTCTGCTTTTTCTTCCAAATTTGATAGAAAATTGGGGCATAGTCTTGAACCCGACGGTAGATCGTTGAATGATGAACGTTTac                       <  2:369748/138‑1 (MQ=11)

                                                                                                                        cATTTTCCTTTTATTTCGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGACCGTTTaca                      <  1:338358/141‑1 (MQ=255)

                                                                                                                         aTTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTacac                     <  1:1079110/141‑1 (MQ=255)

                                                                                                                          ttttCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAGTGATGAACGTTTAcacc                    <  1:443985/141‑1 (MQ=16)

                                                                                                                            ttCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTAcaccac                  >  1:933024/1‑141 (MQ=40)

                                                                                                                             tCCTTGTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACg                 <  1:408131/141‑1 (MQ=255)

                                                                                                                              ccTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGGACGTTTACACCACGt                <  2:124817/141‑1 (MQ=18)

                                                                                                                              ccTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGt                <  2:948602/141‑1 (MQ=21)

                                                                                                                               cTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAGTTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGtt               >  1:526894/1‑141 (MQ=25)

                                                                                                                               cTTTTATTTTGATAAACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGtt               >  1:591537/1‑141 (MQ=255)

                                                                                                                                  ttATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTccc            <  2:401856/141‑1 (MQ=21)

                                                                                                                                  ttATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTccc            >  1:1000234/1‑141 (MQ=21)

                                                                                                                                    aTTTTGATGTACGTCTCATCAATACGCCATGTGTAATAAGCTTTCTTCAACATGGGCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCtt          <  1:799415/141‑1 (MQ=11)

                                                                                                                                    aTTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCtt          <  1:150291/141‑1 (MQ=255)

                                                                                                                                       ttGATGTACTTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAat       <  2:487250/141‑1 (MQ=255)

                                                                                                                                           tGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATAttt   >  1:600797/1‑141 (MQ=21)

TTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTC  >  CP000730/36175‑36455

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 12 ≤` `ATCG/ATCG` `< 13 ≤` `ATCG/ATCG` `< 27 ≤` `ATCG/ATCG` `< 41 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`