BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,828,980	A→C	11.0%	Q12H (CAA→CAC)	xanP →	xanthine permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,828,980	0	A	C	11.0%	303.6 / 2.2	130	Q12H (CAA→CAC)	xanP	xanthine permease
Reads supporting (aligned to +/- strand): ref base A (59/53); new base C (2/16); total (61/69)
Fisher's exact test for biased strand distribution p-value = 8.59e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 134 total aligned reads displayed.

ATAAAAAAGTGCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTGCCAGGCGCTGG  >  NC_000913/3828833‑3829128

ataAAAAAGTGCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAAcc                                                                                                                                                    <  2:2448815/150‑1 (MQ=255)

ataAAAAAGTGCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAAcc                                                                                                                                                    <  1:1200998/150‑1 (MQ=255)

ataAAAAAGTGCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAAcc                                                                                                                                                    >  2:1671893/1‑150 (MQ=255)

ataAAAAAGTGCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAAcc                                                                                                                                                    <  1:2091455/150‑1 (MQ=255)

ataAAAAAGTGCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAAcc                                                                                                                                                    <  1:2091385/150‑1 (MQ=255)

ataAAAAAGTGCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTGTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAAcc                                                                                                                                                    <  2:2448739/150‑1 (MQ=255)

 taAAAAAGTGCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGCGTCAGAAAATGCGCAACCg                                                                                                                                                   >  2:2356027/1‑150 (MQ=255)

  aaaaaaGTGCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCgg                                                                                                                                                  <  2:2711901/150‑1 (MQ=255)

  aaaaaaGTGCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCgg                                                                                                                                                  <  2:2107070/150‑1 (MQ=255)

   aaaaaGTGCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGt                                                                                                                                                 >  1:1749483/1‑150 (MQ=255)

      aaGTGCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGGTgc                                                                                                                                              >  1:1213885/1‑150 (MQ=255)

      aaGTGCATTTATCTAATTTAAATATTTAGATGAGGTTGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTgc                                                                                                                                              <  1:2760361/150‑1 (MQ=255)

           cATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGa                                                                                                                                         >  2:1318256/1‑150 (MQ=255)

            aTTTATCTCGTTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCattt                                                                                                                                        <  2:2240935/150‑4 (MQ=255)

             tttATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACt                                                                                                                                       >  2:1979414/1‑150 (MQ=255)

                aTCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCaa                                                                                                                                    <  2:693663/150‑1 (MQ=255)

                 tctcATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCaaa                                                                                                                                   >  2:2837088/1‑150 (MQ=255)

                   tcATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAAc                                                                                                                                 >  2:2529730/1‑150 (MQ=255)

                      tttAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGTTTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGc                                                                                                                              <  2:1498451/150‑1 (MQ=255)

                      tttAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGc                                                                                                                              <  2:2841740/150‑1 (MQ=255)

                      tttAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGc                                                                                                                              <  2:835857/150‑1 (MQ=255)

                       ttAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTGTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCg                                                                                                                             >  2:1660327/1‑150 (MQ=255)

                              tttAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCCACCGGTTGCGCAGACACAAAAAAGCGAACTGAt                                                                                                                      >  1:1045084/1‑150 (MQ=255)

                              tttAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCCGTTGCGCAGACTCAAAACAGCGAACAGAt                                                                                                                      >  1:379050/1‑150 (MQ=255)

                               ttAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCATAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGAtt                                                                                                                     <  2:1055409/150‑1 (MQ=255)

                               ttAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGAtt                                                                                                                     >  1:95609/1‑150 (MQ=255)

                               ttAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATAGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCACACTCAAAACAGCGAACTGAtt                                                                                                                     >  2:278725/1‑150 (MQ=255)

                                    tGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTAACg                                                                                                                >  1:2899658/1‑150 (MQ=255)

                                       ggATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCt                                                                                                             >  2:1105416/1‑150 (MQ=255)

                                           gtgGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGTGTCAGAAATGGCGCAACGGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGaa                                                                                                         <  2:1326321/150‑1 (MQ=255)

                                           gtgGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGaa                                                                                                         <  2:576922/150‑1 (MQ=255)

                                             ggCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCACCCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGa                                                                                                       >  2:2370982/1‑150 (MQ=255)

                                             ggCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGa                                                                                                       >  2:2370977/1‑150 (MQ=255)

                                             ggCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCGCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGa                                                                                                       <  2:498713/150‑1 (MQ=255)

                                              gCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCCGTTGCGCAGACACAAAACAACGAAAAGATTTAACGTCTTGAAGAt                                                                                                      >  1:1684816/1‑150 (MQ=255)

                                                     tttGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTccgc                                                                                               >  2:382588/1‑150 (MQ=255)

                                                     tttGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTccgc                                                                                               <  2:2107265/150‑1 (MQ=255)

                                                     tttGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTccgc                                                                                               <  2:2107397/150‑1 (MQ=255)

                                                     ttggCTTTTCTCTTTCGGCCTGCTAAACGCCCCGCTTTGCTACCCGGGTTTTGTTTTCTTTGTCTGTTCCCACCGTCGAGCCAGAAAAGGCGCACCCGGTTGCCCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTccgc                                                                                               <  1:384326/147‑1 (MQ=255)

                                                       tGCTTTTCTCTGTCGGTCAGATAAAATGCCCGCTTTGCTACCCCGGGTTTGTTTTCTTTGTCTGTTACCACCCTCGAGTCATAAATTGCGCAACCGGGTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTccgccg                                                                                             <  1:2317267/150‑1 (MQ=255)

                                                       ggCTTATCTCTTCCGGACCGCTAAACGGCCGGCTTTGCTCCCAGGGGATGTTTATCGAGGGCTGTTTCCACCCGCGCGACAGAAAATGCGCACCCGGTGGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTccgccg                                                                                             <  1:2317448/149‑1 (MQ=255)

                                                         cTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCt                                                                                           >  2:1866152/1‑150 (MQ=255)

                                                          ttttCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACGCAAAACAGCGAACTGATTTACCCTCTTGAAGATCGTCCGCCGCtt                                                                                          >  1:1206454/1‑150 (MQ=255)

                                                                 tttCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCaaa                                                                                   <  2:1388660/150‑1 (MQ=255)

                                                                 tttCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCaaa                                                                                   <  2:796355/150‑1 (MQ=255)

                                                                 ttgcGGCCAGCTAAACGGCCCGCTGGGCTACCGGGGGTTTGTTTTGGTTGTGTGTTCCCACCCTCGAGTCAGAACAGGCGCACCGGGGGGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCaaa                                                                                   <  1:498925/147‑1 (MQ=255)

                                                                 tgtCGGTCAGCTAATATGCCGGGTTTGCTCCCACGGGATTGTTTTGGCTGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCACCCGGTTGCGCAGACTCAAAACAGCGAACTGAGTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCaaa                                                                                   <  2:7720/148‑1 (MQ=255)

                                                                     ggTCAGCTAAATTGCCCGCTTTGCTACCACGGTTTTTTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTGGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCt                                                                               <  1:55179/150‑1 (MQ=255)

                                                                     ggTCAGATAAATGGCCCGCTTTGCTACCAGTTTTTTGTTTTGGTTGTCTTTTTCCACCCCCTCGTCATATATTGCGCACCGGGTGGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCt                                                                               <  1:155307/150‑1 (MQ=255)

                                                                      gTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGGGGAACCGGTTGCGCAGACTAAAAAAAGCGAACTGATTTTACGTCTTGTATATAGTCCGCCGCTTGCTCAAAACCTg                                                                              >  1:2269179/1‑150 (MQ=255)

                                                                      gTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAAAAGAGAACTGTTTTACCGTCTTGAAGATCGTACGCCGCTTCCTCAAACCCTg                                                                              >  2:1119208/1‑150 (MQ=255)

                                                                       tCAGATAAACGGCCCTCTTTGATACCAGGGGTTTGTTTTCTTTGTTTGTTTCCACCCTCGAGTAAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGt                                                                             <  1:447675/150‑1 (MQ=255)

                                                                        ccgCTAACCTGCCCGCTTGGCTACCAGGGGTTTGTTGTCGATGTCTGTTTCCCCCCGCGCGGCCGAAACTGCGCACCGGGTGGCGCAGACTCAACACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGtt                                                                            <  1:1390939/148‑1 (MQ=255)

                                                                        cAGCTAACATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGtt                                                                            <  2:1528960/150‑1 (MQ=255)

                                                                         aGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAAACCGTTGCGCAGACTCAAAACAGAGAACTGATTTAAAGTCTTTAAGATATTACGCCGCTTCCACAAAACCTGttt                                                                           >  1:1981151/1‑150 (MQ=255)

                                                                           cTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTgc                                                                         <  2:2417396/150‑1 (MQ=255)

                                                                           cTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAAAAGCGAACTGATTTAACGTCTTGAAGAGAGTCCGCCGCTTCCTCAAACCCGGTTTgc                                                                         >  1:7003/1‑150 (MQ=255)

                                                                           cTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGGTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCGTCCTCAAACCCTGGTTgc                                                                         >  2:332103/1‑150 (MQ=255)

                                                                           cTAAAATGCCCGCTTTGCTACCACGGGATTGTGTTCGATGTCTGTTTCCACCCTCGAGCCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTgc                                                                         <  2:1025676/150‑1 (MQ=255)

                                                                               caggCCCGCTTTGCTACCACGGGATTTTTGTCTTTGTCTGTTCCCCCCCTCGAGTCAGAAAATGCGCAACGGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTgccgcc                                                                     <  1:2473572/147‑1 (MQ=255)

                                                                               atggCCCGCTTTGCTACCACGGGATTGTTTTCGTTGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCACCCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTgccgcc                                                                     <  1:2315750/147‑1 (MQ=255)

                                                                                  gCCCGGTTTGCTACCTGGGTTTTGTTTTGTTTGTGTTTTTCCACCCTCGAGTCAGAAATTGCGCACCGGGTTGCCCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGt                                                                  <  1:2529730/150‑1 (MQ=255)

                                                                                   cccGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAAACGGTTGCGCAGACTCAAAACAGCGAACTGATTTAACGTCTTGAAGATAGTGCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTc                                                                 >  1:1148954/1‑150 (MQ=255)

                                                                                    ccGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAAATGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCa                                                                >  2:376894/1‑150 (MQ=255)

                                                                                     cGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCCACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCCTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCAGTTTGCCGCCAGTCAg                                                               >  1:713598/1‑150 (MQ=255)

                                                                                     cGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGGCTTGAAGATCGTCCGCCGCTTCCTCAAACCCAGTTTGCCGCCTGTCAg                                                               >  1:408518/1‑150 (MQ=255)

                                                                                      gCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGc                                                              >  2:2745118/1‑150 (MQ=255)

                                                                                            cTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCACGAGTCAGAAAATGCGCAACCGGTTTCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTACCTCAAAACCTGTTTGCCGCCTGTCAGCATctgc                                                        >  2:1854843/1‑150 (MQ=255)

                                                                                             tACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATctgct                                                       >  2:1920981/1‑150 (MQ=255)

                                                                                                 aCGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCACCCGGTTGCGCAGACTCAAAAAAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCg                                                   >  2:1069766/1‑150 (MQ=255)

                                                                                                  cGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGa                                                  >  2:2865044/1‑150 (MQ=255)

                                                                                                   gggATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGGAACCCGTTGCGCAGACTCAAAAAAGCGAACAGATTTAACGTCTTTTAGATAGTCCGCCGCGTCCTCAAACCCAGTTTGCCGCCGGTCAGCAACTGCGGGCGAt                                                 >  1:2785954/1‑150 (MQ=255)

                                                                                                    ggtttGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATg                                                <  1:2007742/147‑1 (MQ=255)

                                                                                                    ggATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATg                                                <  2:95609/150‑1 (MQ=255)

                                                                                                    ggATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAAAAGAGAACAGATTTACAGTCTTTAAGATCGTCCGCCCCGTCCTCAAACCCAGGTTGGCGCCTGTCAGCATCTGCTGGCGATg                                                >  1:949635/1‑150 (MQ=255)

                                                                                                          ttttCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGGGCCACCCGTTGCGCAGAATCAAAAACCCGGACGGATTTAACCTCGTTCATATAGTCACCCCCGTACTCAAACCCAGTGTGGCGCCCGTCATCCTATTCTGGGGCTGTTCGtt                                          >  1:742840/1‑150 (MQ=255)

                                                                                                          ttttCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCCGTTGCGCAGACTCAAAACAGCGAACTGATTTAACGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTCGCGAAGTTCGtt                                          >  2:1357379/1‑150 (MQ=255)

                                                                                                            ttCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGc                                        >  2:2178055/1‑150 (MQ=255)

                                                                                                             tCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCg                                       <  2:701356/150‑1 (MQ=255)

                                                                                                               gATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGt                                     <  2:695277/150‑1 (MQ=255)

                                                                                                                  gTCTGTTTCCACCCTCGAGTCAGAACAGGCGCACCGGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGAt                                  <  1:1806111/150‑1 (MQ=255)

                                                                                                                   tCTGTGTCCCCCCTCGAGTCAGAAAATGCGCCACCGGTGGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATc                                 <  2:2557612/150‑1 (MQ=255)

                                                                                                                   tCTGGTGCCACCCTCGAGTCAGAACATGCGCACCCGGTGGCGCAGACTCAACACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCGTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATc                                 <  2:2653535/150‑1 (MQ=255)

                                                                                                                    cTGTGTCCACCCTCGTGTCCGAAACGGCGCCACGGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCa                                <  1:1458315/150‑1 (MQ=255)

                                                                                                                          ccACCCTCGAGTCAGAAAATGCGCACCCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAg                          <  1:1979414/150‑1 (MQ=255)

                                                                                                                          ccACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAg                          >  2:1316712/1‑150 (MQ=255)

                                                                                                                          ccACCCGCGAGCCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAg                          <  2:751392/150‑1 (MQ=255)

                                                                                                                           cACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAgc                         >  2:497762/1‑150 (MQ=255)

                                                                                                                           cACCCTCGAGTCAGAAAATGCGCAACCGGGTGCGCAGACTCAAAACAGCGCACTGAGTTCCCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCGGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCCgc                         >  2:1032125/1‑150 (MQ=255)

                                                                                                                             cccTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAgcgc                       <  1:2583086/150‑1 (MQ=255)

                                                                                                                             cccTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAgcgc                       <  1:163293/150‑1 (MQ=255)

                                                                                                                                  gAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTa                  <  1:947294/150‑1 (MQ=255)

                                                                                                                                    gTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAAt                >  2:658088/1‑150 (MQ=255)

                                                                                                                                       aGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTg             <  2:633265/150‑1 (MQ=255)

                                                                                                                                       aGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTg             >  2:1861803/1‑150 (MQ=255)

                                                                                                                                        gAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTGc            <  2:670172/150‑1 (MQ=255)

                                                                                                                                          aaaTGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTGCCa          <  2:1641114/150‑1 (MQ=255)

                                                                                                                                               gccAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTGCCAGgcgc     >  2:148861/3‑150 (MQ=255)

                                                                                                                                                  aaCCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTGCCAGGCGCTgg  >  2:1706286/1‑150 (MQ=255)

                                                                                                                                                  aaCCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGGTCGTTGCGGTGATCACGCCAGCGCTATTAATCTGCCAGGCGCTgg  <  2:1923986/150‑1 (MQ=255)

ATAAAAAAGTGCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTGCCAGGCGCTGG  >  NC_000913/3828833‑3829128

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