BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	36,315	G→A	56.0%	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.0%	83.0 / inf	348	Y79Y (TAC→TAT)	USA300HOU_0026	IS431mec transposase
Reads supporting (aligned to +/- strand): ref base G (67/81); new base A (82/108); total (156/192)
Fisher's exact test for biased strand distribution p-value = 7.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 354 total aligned reads displayed.

TTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTC  >  CP000730/36175‑36455

ttACCAAATTGTTTAATGGGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATg                                                                                                                                              <  2:198336/141‑1 (MQ=17)

  aCCAAGTTGTTTAAGGAGCCGTTTGATAAACGTATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTa                                                                                                                                            <  2:371565/141‑1 (MQ=255)

  aCCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATATATGTCGCTCTGCATCCATGGCACGAGATAAATAGCTCACTTTACGTTTTATTTTCAtata                                                                                                                                            >  2:377235/1‑141 (MQ=11)

        ttGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGtctc                                                                                                                                      <  2:130738/141‑1 (MQ=17)

          gTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTTtcat                                                                                                                                    >  1:100732/1‑141 (MQ=12)

          gTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCtcat                                                                                                                                    >  1:170487/1‑141 (MQ=32)

          gTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCtcat                                                                                                                                    <  1:337707/141‑1 (MQ=32)

           tttAATGAGACGTTTGATAAACGCATCTGATGAATGATTATCTCGTTGCTTAAGCAACCAAATATCTAATGTGTGTTCCTCTGCATACATGGCTCGAAAGAAAAAGATCCATTTTCCTATTCATTTGATTTATGtactcaa                                                                                                                                   >  2:40027/1‑135 (MQ=2)

             tAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGTCCCATTTTCCTTTTATTTTGATGTAGGTCACATCaa                                                                                                                                 >  2:76155/1‑141 (MQ=255)

              aaTGAGACCTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAAt                                                                                                                                <  1:367348/141‑1 (MQ=12)

               aggagaCGTGTGAAAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATa                                                                                                                               <  1:3099/139‑1 (MQ=12)

               aTGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATa                                                                                                                               >  2:56117/1‑141 (MQ=32)

                tGAGTCGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGTTATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATAc                                                                                                                              <  1:219894/141‑1 (MQ=17)

                  agccGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCACATGTGCAGTTTTTGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTGATCAATACGc                                                                                                                            <  1:289686/138‑1 (MQ=11)

                   gaCGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGcc                                                                                                                           <  2:305808/141‑1 (MQ=17)

                    aCGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACTTCTCATCAATACGccc                                                                                                                          >  1:42029/1‑140 (MQ=12)

                      gTTTGATAAACGGATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCAtt                                                                                                                        <  2:246926/141‑1 (MQ=12)

                       tttGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATCCGTCTCATCAATACGCCAtct                                                                                                                       <  2:233521/141‑3 (MQ=255)

                        ttGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTTTGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATTTACGTCTCTTAAATACGTCAattt                                                                                                                      >  1:125635/1‑140 (MQ=11)

                        ttGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGTCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATCCGCCATTTg                                                                                                                      >  1:321443/1‑141 (MQ=17)

                           aTAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGtaa                                                                                                                   <  1:158409/141‑1 (MQ=12)

                              aaCGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGtaataa                                                                                                                <  2:200519/141‑1 (MQ=33)

                              aaCGCATATACTGAATGATTATCTAGATGCTTACTCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCTTTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGtaataa                                                                                                                >  1:292971/1‑141 (MQ=11)

                                cGCATATGCTGAATGATTATCTCGTTTCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCACCAATACGCCGTTTGAAATAAGc                                                                                                              >  1:283701/1‑141 (MQ=255)

                                 gtatatGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCt                                                                                                             <  1:333792/139‑1 (MQ=255)

                                  cATATGCCGAATGATTGTCTCGTTGCTTACGCAACCAAAGATCGAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGtaataaccct                                                                                                            >  2:374246/1‑137 (MQ=11)

                                   atatGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGCCCCTCTGCATCAATGGCACGATATACATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCttt                                                                                                           <  1:28020/141‑1 (MQ=25)

                                           aaTGATTATCTCGTTGCTTACGCAACCAAAGATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTACTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCt                                                                                                   <  2:293304/141‑1 (MQ=255)

                                            aTGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTGATAATCgtttttctttat                                                                                                  >  1:152898/1‑135 (MQ=11)

                                               aTTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCttttt                                                                                               >  2:196778/1‑141 (MQ=255)

                                                 tATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTttct                                                                                             <  2:322024/141‑1 (MQ=255)

                                                  ttcccGTTGCTTGCGTAACCAAAGTTCTAATGTATTTCCCTGTGTATCAATGGTACGATATAAATAGCTCCATTTTCCTTGTATTTTGTGATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTttctt                                                                                            <  2:194900/137‑1 (MQ=11)

                                                   tctcGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTttcttc                                                                                           <  2:102466/141‑1 (MQ=255)

                                                   tctcGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTttcttc                                                                                           <  1:163467/141‑1 (MQ=40)

                                                       ggtactcaCGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCCTTTTCTTCCaaa                                                                                       >  2:103414/8‑141 (MQ=255)

                                                        ttGCTCACGCAACCAAATATCTAATGCATGTCCCCCTGTATCAATGGCACGATATAAATAGCACCATTTTCCTTTTATTTTGATACCCGTCTCATCAATACGCCCTTTGTAATCACCTTCTTGATGCGTTTTCTTCCcaac                                                                                      >  2:239559/1‑140 (MQ=11)

                                                            ttACGCAACCAAATACCTAATGTATTTTCCTCTTCATCAACGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGa                                                                                  <  2:1228/141‑1 (MQ=11)

                                                             tACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTATTTCTTCCAAATTTGAt                                                                                 <  1:364654/141‑1 (MQ=255)

                                                              aCGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCTATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAATTTTGATa                                                                                <  1:251053/141‑1 (MQ=255)

                                                               cGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCTTTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATAc                                                                               <  1:23891/141‑1 (MQ=255)

                                                               cGCAACCAAATATCTAATGGATGTCCCGCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATat                                                                               <  1:63778/141‑2 (MQ=255)

                                                                gCAAACAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATata                                                                              <  2:216560/141‑3 (MQ=255)

                                                                ccAACCAAATATCTAATGTATGTCCCTCTGTATCAATGACACGAAATAAATAGCTCCATTTTCCTTTTATTTTGATCTACGTCTCATCAATACGCCGTTGTTAATAACTCTTTTGTTGCTTTTGCTTCATAAttcatgata                                                                              >  2:267480/2‑134 (MQ=11)

                                                                 cAACCAAATATCTAATGTATGTCCCTCTGTATAAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATTCGTCTCATAAATAGGCCATTTTTTATAAGCATTTTTATGCTTTTTGTTCCAAATTTGATaaaa                                                                             >  1:104592/1‑141 (MQ=11)

                                                                 cAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATaa                                                                             >  2:285797/1‑141 (MQ=38)

                                                                    ccAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAAt                                                                          <  1:340049/141‑1 (MQ=37)

                                                                     cAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAAtt                                                                         >  1:206593/1‑141 (MQ=255)

                                                                     cAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAAtt                                                                         >  1:375413/1‑141 (MQ=255)

                                                                      aaaTATCTAATGTATGTCCCTCTGTATTAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTg                                                                        <  2:254597/141‑1 (MQ=255)

                                                                       aagatCTAATGTATGTGCCTCGGTAGCAATGGCACGATATAAAGAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGGTTTTTTTTGCTTTTTCTTCCAGATTTGATATAAAATTgg                                                                       >  2:333987/4‑141 (MQ=11)

                                                                       aaTATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTgg                                                                       <  1:195483/141‑1 (MQ=255)

                                                                            cTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCat                                                                  <  2:139572/141‑1 (MQ=255)

                                                                            cTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCat                                                                  <  1:56117/141‑1 (MQ=37)

                                                                             tAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAACTTGATATAAAATTGGGGCatc                                                                 >  1:387885/1‑140 (MQ=255)

                                                                                tgtatgtCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTTTAATAAGCTTTTTTATGCTTTTTCTTCAAAATTTGATACAAAATTgacacctatac                                                              >  1:186768/1‑131 (MQ=14)

                                                                                 gtatgtCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCt                                                             >  2:227065/1‑141 (MQ=255)

                                                                                 gtatgtCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCt                                                             <  2:316789/141‑1 (MQ=255)

                                                                                    tttCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGa                                                          <  2:272843/139‑1 (MQ=18)

                                                                                      tCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCGTTTTCTTCCAAACGTGATGCACAATTGGGGGATCTTCTTGCAc                                                        >  1:180757/1‑141 (MQ=14)

                                                                                        ccTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAAccc                                                      >  2:301555/1‑141 (MQ=255)

                                                                                           ctGCATAAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCAAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAAc                                                   <  1:69529/141‑1 (MQ=255)

                                                                                             gCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACgg                                                 <  2:106740/141‑1 (MQ=40)

                                                                                               aTCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTa                                               <  2:154053/141‑1 (MQ=21)

                                                                                                tCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATATGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAACTTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAg                                              >  1:377850/1‑141 (MQ=255)

                                                                                                 cAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTCTTTTCTGCTTTTTCTTCCACATTTGATATAAACTTGGGGTATATTCTTGAACCCATCGGTAgt                                             <  2:35830/141‑2 (MQ=11)

                                                                                                 cAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGa                                             <  1:271157/141‑1 (MQ=21)

                                                                                                 aaaTGGTGCGGTATAAATAGTTGCATTTTTCTTTTATTTTGATGTAGGTCTCATGAATGTGCCTTTTGTAGTAAGCTTTTTTATGCTTTTTCTTTCTAATTTGATTTAAAATTGGGGCGTTTTTTTGATCGCAACGGTAgt                                             >  2:134883/2‑140 (MQ=255)

                                                                                                   aTGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGAcc                                           <  1:109971/141‑1 (MQ=39)

                                                                                                    tGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAAAGGTAGATCg                                          >  1:32571/1‑141 (MQ=255)

                                                                                                     ggCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTTTGCTTTTTCTTCCAAATTTGCTATAAAATTAGGGCATATTATTGAACGCCACgagtgatcgt                                         >  1:10772/1‑132 (MQ=11)

                                                                                                        aCGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCACCGGTAGATCGTTGa                                      >  2:108245/1‑141 (MQ=255)

                                                                                                           atatAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAatg                                   >  1:243346/1‑141 (MQ=255)

                                                                                                            tataAATGGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATTAATACGCCATTTGTAATAAGCTTTTTTATGTTTTTTCTTCCTAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAatga                                  <  2:51290/141‑1 (MQ=11)

                                                                                                            tataAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAatga                                  >  2:149092/1‑141 (MQ=21)

                                                                                                            tataAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATCAAATTGGGGCATATTCTTGAACCCAATGGTAGATCGTTGAatga                                  >  2:4583/1‑141 (MQ=11)

                                                                                                            tataAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAatga                                  >  2:390374/1‑141 (MQ=255)

                                                                                                             ataAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAatgat                                 <  2:349946/141‑1 (MQ=21)

                                                                                                              taAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCAGATTCTTCAAACCACCGGCAGATCTGTGAataatg                                >  1:50953/1‑141 (MQ=11)

                                                                                                                    gCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGGAGACCGTTGAATGATGAACGtt                          >  1:82115/1‑141 (MQ=21)

                                                                                                                    gCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAACATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGtt                          <  1:264108/141‑1 (MQ=255)

                                                                                                                     cTCCATTTTCCTTTTATCTTGATCTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGTGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGttt                         <  2:172714/141‑1 (MQ=11)

                                                                                                                      tCCATTTTTCTTTTATTTTGATATACGTCTCAACAATACGCCATTTGTTATAAGCTTTTTTATGCTTTTTTTTCCAAATTTGATATAAAATTTTTGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTa                        <  2:188482/141‑1 (MQ=11)

                                                                                                                         aTTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTacac                     <  1:228028/141‑1 (MQ=40)

                                                                                                                          ttttCCTTTTATTTTGATATCCGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTTTGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTAcacc                    <  2:311478/141‑1 (MQ=11)

                                                                                                                          ttttCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGGTTGTTGAATGATGAACGTTTAcacc                    <  2:384584/141‑1 (MQ=11)

                                                                                                                           tttCTTTTTATTTTGATATACGTCTGATCAATACGCCATTTGTAATAAGCTTTTTTCTGCTTTTTCTTCCAAATTTGATATAGAATTGGGGCATAGTCTTGACCCCAACGGTAGATCGTTGAATGATGAACGTTTAcacca                   <  2:178242/141‑1 (MQ=11)

                                                                                                                           tttCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTAcacca                   <  1:224928/141‑1 (MQ=18)

                                                                                                                           tttCCTTTTATTTTGATATACGGCTCATCAATACGCCGTTTGTAATAAGCTTTTTTATGCTTTTTCTTCCCAATTTGATATAAAATTGGCGCATGTTCTTGAACCGAACGGTGAATCGTTGAATGGTGAACGTTTAcaccc                   >  2:65178/1‑140 (MQ=11)

                                                                                                                             tCCTTTTATTTCGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACg                 <  1:113942/141‑1 (MQ=21)

                                                                                                                              ccTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGGATGATGAACGTTTACACCACGt                >  1:118701/1‑141 (MQ=255)

                                                                                                                               cTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGtt               >  1:346943/1‑141 (MQ=25)

                                                                                                                               cTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGtt               >  1:137503/1‑141 (MQ=255)

                                                                                                                                tgttATTTTGATATACGTTTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTACCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTc              <  2:84499/139‑1 (MQ=11)

                                                                                                                                   tATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTCTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCt           <  1:303759/141‑1 (MQ=21)

                                                                                                                                    aTTTTGATGTACGTGTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCtt          <  2:169728/141‑1 (MQ=21)

                                                                                                                                    aTTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCtt          <  2:338021/141‑1 (MQ=11)

                                                                                                                                     ttttGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTa         <  2:162852/141‑1 (MQ=21)

                                                                                                                                      tttGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACTACTGAATGATGAACGTTTACACCACGTTCCCTTaa        <  1:192682/141‑1 (MQ=255)

                                                                                                                                            gTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGCATGATGAACGTTTCCACCCCGTTCCCTTAATCtttt  >  1:356129/1‑140 (MQ=255)

                                                                                                                                            gTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTc  >  1:218434/1‑141 (MQ=40)

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: `‑`