Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,378,802	0	A	C	42.5%	141.6 / 11.0	73	V48G (GTA→GGA)	ampC	penicillin‑binding protein; beta‑lactamase, intrinsically weak
Reads supporting (aligned to +/- strand):  ref base A (40/2);  new base C (0/31);  total (40/33)
Fisher's exact test for biased strand distribution p-value = 1.36e-18
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCCTCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCA‑GG‑T‑AAAGTAATAAGG‑TTTACCC‑TGATAAA‑TTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCT  >  NC_000913/4378654‑4378951                                                                                                                                                           |                                                                                                                                                      cccTCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCA‑GG‑T‑AAAGTAATAAGG‑TTTACAC‑TGATAAA‑TTacc                                                                                                                                                     <  2:510754/151‑1 (MQ=255)  ccTCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCA‑GG‑T‑AAAGTAATAAGG‑TTTACCC‑TGATAAA‑TTaccg                                                                                                                                                    <  1:409294/151‑1 (MQ=255)     cGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCA‑GG‑T‑AAAGTAATAAGG‑TTTACCC‑TGATAAA‑TTaccgcca                                                                                                                                                 <  2:721650/151‑1 (MQ=255)      gAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCA‑GG‑T‑AAAGTAATAAGG‑TTTACCC‑TGATAAA‑TTaccgccac                                                                                                                                                >  1:1073578/1‑151 (MQ=255)         cAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCA‑GG‑T‑AAAGTAATAAGG‑TTTACCC‑TGATAAA‑TTaccgccaccgc                                                                                                                                             >  2:599631/1‑151 (MQ=255)          aaTAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCA‑GG‑T‑AAAGTAATAAGG‑TTTACCC‑TGATAAA‑TTaccgccaccgcc                                                                                                                                            >  1:323680/1‑151 (MQ=255)            tAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCA‑GG‑T‑AAAGTAATAAGG‑TTTACCC‑TGATAAA‑TTACCGCCACCGCCAt                                                                                                                                          >  1:581467/1‑151 (MQ=255)                     caccaAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCA‑GG‑T‑AAAGTAATAAGG‑TTTACCC‑TGATAAA‑TTACCGCCACCGCCATACCCGGGAt                                                                                                                                 >  2:1071761/1‑151 (MQ=255)                     caccaAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCA‑GG‑T‑AAAGTAATAAGG‑TTTACCC‑TGATAAA‑TTACCGCCACCGCCATACCCGGGAt                                                                                                                                 >  2:421092/1‑151 (MQ=255)                       ccaAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCA‑GG‑T‑AAAGTAATAAGG‑TTTACCC‑TGATAAA‑TTACCGCCACCGCCATACCCGGGATCt                                                                                                                               >  1:966592/1‑151 (MQ=255)                         aaccccGCCAGAAAATTTTTTGCTGCCCAACCCAAACTAAAACAACGTTTGTGGTGGGAGGGGCTGTTTTTGGGGGTTGCCCGCATACCCCCA‑GG‑A‑AAAGAAAAAAGT‑TTCCCCT‑TAAAAAA‑TCCCCGCCACCGCCATACCCGGGATCt                                                                                                                               <  2:744204/144‑1 (MQ=255)                            cccGCCATAAATTGTTTGTCTGCCCGAACCAAACTAAAAAAAGGTTGGTGGTGTGGCGGGTTTTTTTTTGGAGATGTCCGCACAGCCACA‑GG‑A‑AAAGAAAAAAGT‑TTTACCC‑TAAAAAT‑TCACCGCCACCGCCATACCCGGGATCtttt                                                                                                                            <  1:59913/147‑1 (MQ=255)                             aCGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCA‑GG‑T‑AAAGTAATAAGG‑TTTACCC‑TGATAAA‑TTACCGCCACCGCCATACCCGGGATCTTttgtt                                                                                                                         >  1:1173888/1‑151 (MQ=255)                              cGCCAGTAAATGTTTTGCGAACCGACCCAAACTAAAACACCTTTTGTGGTGTGAGGGGTTGTTTTTTGGCGATGCCCGCATACCCCCA‑GG‑A‑AAAGAAAAAAGT‑TTTCCCC‑TAAAAAT‑TTCCCCCCACCGCCATACCCGGGATCTTttgt                                                                                                                          <  2:1192692/149‑1 (MQ=255)                               gCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCA‑GG‑T‑AAAGTAATAAGG‑TTTACCC‑TGATAAA‑TTACCGCCACCGCCATACCCGGGATCTTTTGTTGc                                                                                                                       >  1:785790/1‑151 (MQ=255)                                caagaaatttttttGCTGCCCGACCCAAACACAAACACCGTTTGCTGTGTGGGGGGCTTTTTTTGGGCGATGCCCGCACACCCCCG‑GG‑A‑AAAGAAATAGGT‑TTTCCCT‑TAAAAAA‑TTCCCGCCACCGCCATACCCGGGATCTTttgtt                                                                                                                         <  1:379126/138‑1 (MQ=255)                                     aagtGTTTGGCTGCCCAAACCAAACTAAACCAGCTTGTGCTGTGGGGCGGGCTGTTTTTGGGGGTTGCCCTCATACCCCCG‑GG‑G‑AAAGTATTAAGT‑TTTCCCC‑GGATAAA‑TTCCCGCCCCCGCCATACCCGGGATCTTTTGTTGct                                                                                                                      <  1:1152019/143‑1 (MQ=255)                                     aaaTTTTTTGCTGACAGACCCTAACTAAAACACCTTTTGCTGTGTGAGGGGTTGTTTTTTGGCGATGCCCGCATGCCCCCA‑GA‑A‑AAAGTAATAAGT‑TTTCCCC‑TAATAAT‑TTCCCCCCACCGCCATACCCGGGATCTTTTGTTGctc                                                                                                                     <  2:879048/147‑1 (MQ=255)                                     aaaTGTTTTGCGGACCAACCCTACCAAAAACAACTTTTGCGGTGTGACGGCCTGTTTTTTGGCGTTGTCCGCATAGCCCAA‑GG‑A‑AAAGTAATAAGG‑TTTCCCC‑TAAAAAA‑TTCCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCtata                                                                                                                 <  2:298047/151‑1 (MQ=255)                                          ttttGCTGACCAACCCAAACCAAAACAGCTTTTGCTGTGTGAGGGTCTTTTTTTTGGCTATGCCCGCATCCCCCCA‑GG‑A‑AAATTAAAAAGT‑TTTCCCT‑TAATAAT‑TCCCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCg                                                                                                             <  1:1338386/150‑1 (MQ=255)                                          ttttGCTGACCAACCAAAACTAAAACAACTTTTGCTGTGTGAGGGGCTGTTTTGTGGGGATTTCCGCATGCCCCCA‑GG‑T‑AAAAAAATAAGT‑TTTACCT‑TGAAAAC‑TTCCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAg                                                                                                               <  2:1051414/148‑1 (MQ=255)                                           tttGCGGACCGATCCTAACAAAAACACCTTTTGTGGTGGGAGGGGCTGTTTTTTGGCGATGTCCGCA‑GGCCCCA‑GGTT‑AAAGTAATAGGG‑TTTCCCC‑TGATAAA‑TCCCCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCgg                                                                                                            <  2:158224/150‑1 (MQ=255)                                               cggcccGAACCAAACCAAAACAATTTTTGGGGTGGGAGGGGTTGTTTTTGGGCGTTGCCCGAATACCCCCA‑GG‑T‑AAAGTAATAAGG‑TTTCCCC‑TAAAAAT‑TTCCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTaa                                                                                                       <  2:1059534/147‑1 (MQ=255)                                               cTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCA‑GG‑T‑AAAGTAATAAGG‑TTTACCC‑TGATAAA‑TTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTa                                                                                                        >  2:1029158/1‑150 (MQ=255)                                               cTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCA‑GG‑T‑AAAGTAATAAGG‑TTTACCC‑TGATAAA‑TTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTa                                                                                                        >  2:1010492/1‑150 (MQ=255)                                                   ccaacccaaacaaaaaCAACTTTTGTTGTGTGGGGGGTTTTTTTTGGGCGTTGCCCGCATACCCCAA‑GG‑A‑AAAGTAAAAAGG‑TTTCCCA‑TAAAAAA‑TTCCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTa                                                                                                        <  2:35201/133‑1 (MQ=255)                                                   ccaacccTAACAAAAACAACTTTTGCTGTGTGCCGGGCTGTTTTTTGGCGATGCCCGCATACCCCGA‑GG‑A‑AAAGTAATAAGT‑TTTCCCT‑TAAAAAA‑TTCCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTa                                                                                                        <  2:1105171/141‑1 (MQ=255)                                                   ccGAACCAAACAAAAAAAACGTTTGCTGTGTGAGGGGTTGTTTTTTGGCGATGTCCGCATAGCCCAA‑GA‑A‑AAAGAAAAAAGG‑TTTCCCC‑TAATAAT‑TTCCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTaa                                                                                                       <  2:1190892/147‑1 (MQ=255)                                                              tCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCA‑GG‑T‑AAAGTAATAAGG‑TTTACCC‑TGATAAA‑TTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGcaca                                                                                         >  2:575922/1‑150 (MQ=255)                                                                 aacaacGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCA‑GG‑T‑AAAGTAATAAGG‑TTTACCC‑TGATAAA‑TTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAatat                                                                                     >  2:1370146/1‑151 (MQ=255)                                                                 aacaacGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCA‑GG‑T‑AAAGTAATAAGG‑TTTACCC‑TGATAAA‑TTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAatat                                                                                     >  2:594886/1‑151 (MQ=255)                                                                   caacTTGTGGTGTGTGAGGGGCTTTTTTTTGGCGATCTCCGCATCCCCCCA‑GG‑T‑AAAGAAATAGTT‑TTTCCTT‑TTAAAAA‑TGCCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAatat                                                                                     <  1:298239/149‑1 (MQ=255)                                                                   caacGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCA‑GG‑T‑AAAGTAATAAGG‑TTTACCC‑TGATAAA‑TTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCg                                                                                   >  1:573707/1‑151 (MQ=255)                                                                      cGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCA‑GG‑T‑AAAGTAATAAGG‑TTTACCC‑TGATAAA‑TTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTg                                                                                >  2:1417328/1‑151 (MQ=255)                                                                      cGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCA‑GG‑T‑AAAGTAATAAGG‑TTTACCC‑TGATAAA‑TTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTg                                                                                >  2:1221546/1‑151 (MQ=255)                                                                       gTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCA‑GG‑T‑AAAGTAATAAGG‑TTTACCC‑TGATAAA‑TTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTg                                                                                >  2:746627/1‑150 (MQ=255)                                                                       gTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCA‑GG‑T‑AAAGTAATAAGG‑TTTACCC‑TGATAAA‑TTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTg                                                                                >  2:746594/1‑150 (MQ=255)                                                                        tttGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCA‑GG‑T‑AAAGTAATAAGG‑TTTACCC‑TGATAAA‑TTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGAt                                                                              >  2:1144171/1‑151 (MQ=255)                                                                         ttgttgTGTGGGGGCCTTTTTTTTGGCGTTCCCCGCAT‑CCCCCA‑GG‑G‑AAAATTAAAAGGTTTTCCCC‑AGATAAA‑TCCCCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGtt                                                                                 <  2:86022/147‑1 (MQ=255)                                                                             tgtgtgACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCA‑GG‑T‑AAAGTAATAAGG‑TTTACCC‑TGATAAA‑TTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATttgtt                                                                         >  1:407839/1‑151 (MQ=255)                                                                              gtgtgACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCA‑GG‑T‑AAAGTAATAAGG‑TTTACCC‑TGATAAA‑TTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATttgttg                                                                        >  2:1058628/1‑151 (MQ=255)                                                                              gtgggAGGGGCTTTTTTTGGGCGGTGCCCG‑AAAGCCCCC‑AG‑A‑AAAAGAAAAAGG‑TTTTCCCTTAATAAT‑TCCCCCCCTCCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGAtt                                                                             <  1:1396552/146‑1 (MQ=255)                                                                                 tgttggattttttttgGGAGATGCCCGC‑TACCCCCGGGG‑T‑AAAGAAATATGT‑TTCCCCC‑AGAAAAA‑TTCCCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGa                                                                       <  2:1078492/138‑1 (MQ=255)                                                                                     gggCTGTTTTTGGGCGTGGTCCGTATAGCCCAA‑GG‑A‑AAAGAAATAAGT‑TTTCCCC‑TAATAAT‑TCCCCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATttgttg                                                                        <  2:698793/144‑1 (MQ=255)                                                                                     gggCTGCTTTTTGGCGATGTCCGCATAGCCCCA‑GG‑T‑AAAGTAATAAGG‑TTTACCC‑TGATAAA‑TTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGgc                                                                  >  2:801072/1‑150 (MQ=255)                                                                                           ctTTTTGGCGATGTCCGCATAGCCCCA‑GG‑T‑AAAGTAATAAGG‑TTTACCC‑TGATAAA‑TTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAAt                                                           >  2:313412/1‑151 (MQ=255)                                                                                               ttGGTGAGGTCCGCATAGCCCGA‑GG‑A‑AAAGGAAGAAGG‑TTTACCC‑TGATAAA‑TTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTgg                                                       >  1:220717/1‑151 (MQ=255)                                                                                                 ggCGATGTCCGCATAGCCCCA‑GG‑T‑AAAGTAATAAGG‑TTTACCC‑TGATAAA‑TTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTgg                                                       >  1:1000874/1‑149 (MQ=255)                                                                                                 ggCGATGTCCGCATAGCCCCA‑GG‑T‑AAAGTAATAAGG‑TTTACCC‑TGATAAA‑TTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTgg                                                       >  2:658518/1‑149 (MQ=255)                                                                                                  gggATGCCCGCATACCCCCA‑GG‑T‑AAAG‑AAAAAGG‑TTTTCCCTTAAAAAA‑TTCCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGc                                                    <  2:609175/149‑1 (MQ=255)                                                                                                   gtaTGCCCCCATAGCCCAA‑GG‑A‑AAAGTAATAAGT‑TTTCCCT‑AAATAAA‑TACCCGCAACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGc                                                    <  1:1316143/148‑1 (MQ=255)                                                                                                      tGTCCGCATAGCCCCA‑GG‑T‑AAAGTAATAAGG‑TTTACCC‑TGATAAA‑TTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCa                                                   >  2:1222009/1‑148 (MQ=255)                                                                                                         ccctagaGCCACA‑GG‑A‑AAAGTAATAGTT‑TTTCCCC‑TAAAAAT‑TTCCCCCCTCCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAagag                                               <  1:354289/143‑1 (MQ=255)                                                                                                         ccGCATACCCCCA‑GA‑A‑AAAATAATAAGT‑TTTCCCC‑TGAAAAT‑TTCCCGCCCCCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAaga                                                <  2:323730/148‑1 (MQ=255)                                                                                                          ccgcatcccccA‑GG‑AAAAATTAATAGGT‑TTTCCCC‑TAAAAAA‑TCCCCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAagag                                               <  2:83572/142‑1 (MQ=255)                                                                                                          cccaaacccccA‑AG‑A‑AAAAGAAAAAGGTTTTCCCC‑TAATAAA‑TTCCCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAagag                                               <  2:378914/142‑1 (MQ=255)                                                                                                          cGCATAGCCCCA‑GG‑T‑AAAGTAATAAGG‑TTTACCC‑TGATAAA‑TTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCg                                            >  2:231819/1‑151 (MQ=255)                                                                                                           gCATAGCCCCA‑GG‑T‑AAAGTAATAAGG‑TTTACCC‑TGATAAA‑TTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCgg                                           >  1:312896/1‑151 (MQ=255)                                                                                                            cATAGCCCCA‑GG‑T‑AAAGTAATAAGG‑TTTACCC‑TGATAAA‑TTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGt                                          >  2:504863/1‑151 (MQ=255)                                                                                                                ataccccgggcaaaaGAAATAAGG‑TTTACAC‑TGATAAA‑TGCCCCCCCCCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGta                                         <  2:7663/139‑1 (MQ=255)                                                                                                                 ccccG‑GT‑T‑AAAGTAATAAGT‑TTTCCTC‑TAATAAA‑TTCCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGtaat                                       <  2:672450/149‑1 (MQ=255)                                                                                                                 ccccA‑GG‑T‑AAATAAATAAGT‑TTTCCCC‑TGATAAT‑TTCCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGtaatt                                      <  1:1035836/150‑1 (MQ=255)                                                                                                                 ccccA‑GA‑A‑AAAGTAAT‑AGG‑TTTTCCC‑TTAAAAATTTCCCGCCCCCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGtaa                                        <  1:546956/148‑1 (MQ=255)                                                                                                                        g‑T‑AAAGTAATAAGG‑TTTACCC‑TGATAAA‑TTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAgg                               >  2:105379/1‑151 (MQ=255)                                                                                                                         ggaaaataaatagGG‑TTTCCCT‑TAAAAAA‑TTCCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATtaataa                                 <  2:726263/142‑1 (MQ=255)                                                                                                                              aGTAATAAGG‑TTTACCC‑TGATAAA‑TTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGgcgc                            >  2:713892/1‑150 (MQ=255)                                                                                                                              aGTAATAAGG‑TTTACCC‑TGATAAA‑TTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGgcg                             >  1:985549/1‑149 (MQ=255)                                                                                                                                  ataaGG‑TTTACTG‑TGATAAA‑TTAGTGCAACAGACATATCGGGGATGTTTTGTTGATCTATAAGCGGGGTAATTGTGAGATGCAAAATATCGTTGATTTGTTGAGGGGGAGCAAATGTGGAGAAAGAGGCGGTAATTAATAAGGAGGagagg                       >  1:1150806/1‑150 (MQ=255)                                                                                                                                  ataaGG‑TTTACCC‑TGATAAA‑TTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGc                       >  2:1405540/1‑151 (MQ=255)                                                                                                                                  ataaGG‑TTTACCC‑TGATAAA‑TTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGc                       >  2:830011/1‑151 (MQ=255)                                                                                                                                  aaaaGG‑TTTTCCCTTAAAAAA‑TTCCCCCCGCCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCag                          <  2:268103/147‑1 (MQ=255)                                                                                                                                         tttACCC‑TAATAAA‑TTCCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGcgtcgt                  <  2:1041946/150‑1 (MQ=255)                                                                                                                                          ttACCC‑TGATAAA‑TTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGttt                >  1:16551/1‑151 (MQ=255)                                                                                                                                            aCCC‑TGATAAA‑TTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTTATAAGGCGCAGAGCGTCGtttt               >  2:1038875/1‑150 (MQ=255)                                                                                                                                            aCCC‑TGATAAA‑TTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGtttt               >  2:364029/1‑150 (MQ=255)                                                                                                                                             ccc‑TGATAAA‑TTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGtttt               >  1:1145697/1‑149 (MQ=255)                                                                                                                                                  aaTAAA‑AAACCCCCCCCGCCATACCCGGGATCTTTTGTTGCGCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTg              <  1:751335/145‑1 (MQ=255)                                                                                                                                                    tAAA‑TTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACAt         >  1:53513/1‑149 (MQ=255)                                                                                                                                                          tACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCt  >  1:1058162/1‑151 (MQ=255)                                                                                                                                                           accgccaccgccATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAAGCGCAGAGCGTCGTTTTGAACATAGGGt    >  2:804908/1‑148 (MQ=255)                                                                                                                                                           |                                                                                                                                                      CCCTCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCA‑GG‑T‑AAAGTAATAAGG‑TTTACCC‑TGATAAA‑TTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCT  >  NC_000913/4378654‑4378951

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑