Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,222,596	0	G	A	41.9%	59.4 / 47.1	62	F13F (TTC→TTT)	ymgD	uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (28/8);  new base A (0/26);  total (28/34)
Fisher's exact test for biased strand distribution p-value = 1.42e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACGCTGGCATTACCCATTACCGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCG‑TTTTTAATA‑TTGTAGTCCTGCGCCCA‑GGTCATTGGTGCGAAAACAAACA‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGAACGCCAGTGAC  >  NC_000913/1222449‑1222739                                                                                                                                                       |                                                                                                                                                  aaCGCTGGCATTACCCATTACCGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCG‑TTTTTAATA‑TTGTAGTCCTGCGCCCA‑GGTCATTGGTGCGaaa                                                                                                                                                <  2:922613/151‑1 (MQ=255) aaCGATGGCATTACCCATTACCGCAACAATTTCTGCAAACTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCG‑TTTTTAATA‑TTGTAGTCCTGCGCCCA‑GGTCATTGGTGCGaa                                                                                                                                                 >  2:950622/1‑150 (MQ=255)  aCGCTGGCATTACCCATTACCGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCG‑TTTTTAATA‑TTGTAGTCCTGCGCCCA‑GGTCATTGGTGCGaaaa                                                                                                                                               <  2:67147/151‑1 (MQ=255)     cTGGCATTACCCATTACCGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCG‑TTTTTAATA‑TTGTAGTCCTGCGCCCA‑GGTCATTGGTGCGAaaaca                                                                                                                                             >  2:287344/1‑150 (MQ=255)      tGGCATTACCCATTACCGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCG‑TTTTTAATA‑TTGTAGTCCTGCGCCCA‑GGTCATTGGTGCGAaaacaa                                                                                                                                            >  2:287599/1‑150 (MQ=255)          attacccattaccGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCG‑TTTTTAATA‑TTGTAGTCCTGCGCCCA‑GGTCATTGGTGCGAAAACAAACA‑‑gg                                                                                                                                     >  1:202067/1‑151 (MQ=255)            tacccattaccGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCG‑TTTTTAATA‑TTGTAGTCCTGCGCCCA‑GGTCATTGGTGCGAAAACAAACA‑‑GGcc                                                                                                                                   >  2:622326/1‑151 (MQ=255)              cccattaccGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCG‑TTTTTAATA‑TTGTAGTCCTGCGCCCA‑GGTCATTGGTGCGAAAACAAACA‑‑GGccc                                                                                                                                  <  1:855968/150‑1 (MQ=255)              cccattaccGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCG‑TTTTTAATA‑TTGTAGTCCTGCGCCCA‑GGTCATTGGTGCGAAAACAAACA‑‑GGCCCg                                                                                                                                 >  1:456755/1‑151 (MQ=255)                 attaccGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCG‑TTTTTAATA‑TTGTAGTCCTGCGCCCA‑GGTCATTGGTGCGAAAACAAACA‑‑GGCCCGCCa                                                                                                                              <  2:224459/151‑1 (MQ=255)                        caacaaTTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCATTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCG‑TTTTTAATA‑TTGTAGTCCTGCGCCCA‑GGTCATTGGTGCGAAAACAAACA‑‑GGCCCGCCAGTAAAGc                                                                                                                       >  2:850984/1‑151 (MQ=255)                         aacaaTTTCTGCAACCGGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGAAAGTAATATAGTTTTCGGACGGCAAACCG‑TTTTAAATA‑TTGTAGTCCTGCCCCCA‑GGTCATTGGGGCAAAAACAAACA‑‑GGCCCGCCAGTAAAGCa                                                                                                                      <  1:312619/151‑1 (MQ=255)                          acaaTTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCG‑TTTTTAATA‑TTGTAGTCCTGCGCCCA‑GGTCATTGGTGCGAAAACAAACA‑‑GGCCCGCCAGTAAAGCa                                                                                                                      >  2:1019376/1‑150 (MQ=255)                                 ctccacccTGTGCGTTGTCATTTTTTGCCATTTCGTTGGTTTCTGCGAAAGAAATATAGTTTTCTGACGGCAAACCT‑TTTTAAATT‑TTGTAGTCCGGCGCCAA‑GGTCATGGGGGCAAAAAAAAACG‑‑GGCCCGCCAGTAAAGCAAAtttt                                                                                                                <  2:527347/143‑1 (MQ=255)                                          gtgCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCG‑TTTTTAATA‑TTGTAGTCCTGCGCCCA‑GGTCATTGGTGCGAAAACAAACA‑‑GGCCCGCCAGTAAAGCAAATTTTTTcatcatcat                                                                                                     >  2:688070/1‑151 (MQ=255)                                          gtgCGCGGTAATTTTTTCCCTTTTCGTTGGTTTCGGGACAAGTAAAAAAGGTTTCTGACGGCAACCCT‑TTTTTAATA‑TTGTAGCCCTGCGCCAA‑GGCCTTTGGTACAAAAACAAAAA‑‑GCCCCGCCAGTAAAGCAAATTTTTtcatca                                                                                                         <  2:1026746/147‑1 (MQ=255)                                             cgcTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCG‑TTTTTAATA‑TTGTAGTCCTGCGCCCA‑GGTCATTGGTGCGAAAACAAACA‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTcc                                                                                                  >  2:645475/1‑151 (MQ=255)                                                tGTCATTTTTTGCCATTTCGTTGGTTTCTGCGAAAGTAATATAGTTTTCTGAGGGCAACCCT‑TTTTTAATA‑TTGTAGTCCGGCGCCCA‑GGTTATTGGTGCAAAAAAAAACA‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCtt                                                                                                <  2:1005856/150‑1 (MQ=255)                                                tGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCG‑TTTTTAATA‑TTGTAGTCCTGCGCCCA‑GGTCATTGGTGCGAAAACAAACA‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCtt                                                                                                >  1:790070/1‑150 (MQ=255)                                                  tcggTTTTTGCCTTTCCGTTGGCTTCTGCGCAGGTAATTTAGGTTTCGGACGGCAGCCTT‑TTTTTAATA‑TTGTAGCCCGGCGCCAA‑GGTCATGGTTCCGAAAACAAACA‑‑TGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTa                                                                                               <  1:161605/146‑1 (MQ=255)                                                  tCATTTTTTGCCATTTCGTTGTCTTCTGCGAAAGAAATATAGTTTTCGGACGGAAACCCT‑TTTTTAATA‑TTGTGGCCCTGCCCCCG‑GTTCATTGGGGAAAAAAAAAACA‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTc                                                                                                   <  2:578824/145‑1 (MQ=255)                                                  tCATTTTTTGCCATTTCGTTGGTTTCTGCGCAAGAAATATAGTTTTCTGACGACAACCCT‑TTTTTAATA‑TTGTAGTCCAGCGCCCA‑GGTTATTGGTGCAAAAACAAACA‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCAtt                                                                                                    <  2:250369/144‑1 (MQ=255)                                                      tttttGCCTTTTCTTGGTTTTCTGCGAAAGAAATATAGGTTTCTGACGGCAACCCG‑TTTTTAATA‑TTGTAGTCCTGCCCCCA‑GGTCATGGGTGCAAAAACAAACA‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTAt                                                                                              <  1:174260/146‑1 (MQ=255)                                                      tttttGCATTTTTGTGGGTTTCTGCGAAAGTAATATGGGTTTCGGACGG‑AAACCGTTTTTTAATA‑TTGTAGTCCTGCGCCCG‑GGTCTTGGGGGAGAAAAAAAACG‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTc                                                                                           <  1:692307/149‑1 (MQ=255)                                                        tttGCTTTTTCGTTGGTTTCTGCGAAAGTAATATAGGTTTCGGACGGCAAACCT‑TTTTTAATA‑TTGTAGTCCTGCGCCCA‑GGTCTTGGGTGCAAAAACAAACA‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCAtt                                                                                        <  1:329775/150‑1 (MQ=255)                                                            ccATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCG‑TTTTTAATA‑TTGTAGTCCTGCGCCCA‑GGTCATTGGTGCGAAAACAAACA‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTAcc                                                                                   >  1:752822/1‑151 (MQ=255)                                                               tttCGTTGGCTTCTGGAAAAGTAATATAGTTTTCTGACGGCAAACCG‑TTTTTAATA‑TGGTAGTCCGGCGCCCA‑GGTAATTGGTGAAAAAAAAAACA‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTAc                                                                                    <  1:980817/147‑1 (MQ=255)                                                                ttCTTTGGCTTCTGCCAAAGTAATATAGGTTTCGGCCGGCAACCCT‑TTTTTAATA‑TTGTAGTCCTCCGCCCA‑GGTAATGGGTGCAAAAAAAAACA‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAg                                                                                <  1:262547/150‑1 (MQ=255)                                                                   gTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCG‑TTTTTAATA‑TTGTAGTCCTGCGCCCA‑GGTCATTGGTGCGAAAACAAACA‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAAtt                                                                            >  2:642380/1‑151 (MQ=255)                                                                    ttgttttCTGCGAAAGTAATATAGGTTTCTGACGGCAACCTT‑TTTTTAATA‑TTGTACTCCAGCGCCAG‑GTTTATTGGTGAAAAAAAAAACG‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAAtt                                                                            <  2:566710/145‑1 (MQ=255)                                                                    ttGGTTTCTGCGAAAGTAATATAGTTTTCTGACGGCAACCCT‑TTTTTAATA‑TTGTAGCCCTGCGCCCA‑GGTCATGGGTGCAAAAACAAACA‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAAt                                                                             <  1:731378/149‑1 (MQ=255)                                                                    ttGGTTTCGGCGAAAGTAATATAGTTTTCTGACGGCAACCCT‑TTTTAATTT‑TTGGACCCCTGCGCCCA‑GTTCTTTGGCGAAAAAAAAAACA‑‑GCCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAAt                                                                             <  2:1017575/149‑1 (MQ=255)                                                                         ttCTGAGAAAGAAATAGGGGTTTCTGCGGGCAACCCT‑TTTTTAAGT‑TTGTAGTCCCGCCCCCA‑GTTTTTTGGTGGAAAAAAAAACA‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGcacc                                                                       <  2:965769/150‑1 (MQ=255)                                                                                cAAGTAATATAGGTTTCTGACGGCAAACCG‑TTTTTAATA‑TTGTAGTCCTGCGCCCA‑GGTCATTGGTGCGAAAACAAACA‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTAc                                                               >  2:457957/1‑151 (MQ=255)                                                                                 aattAATATTGTTTTCTGACGGAAACCCT‑TTTTAAATT‑TTGTAGCCCTGCCCCAA‑GTTAATTGGTGGAAAAAAAAACG‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCg                                                                  <  1:859588/144‑1 (MQ=255)                                                                                 aaatAATATAGGTTTCTGACGGCAACCCT‑TTTTTAATA‑TTGTAGTCCTGCCCCCA‑GGTCATTGGTGCAAAAACAAACA‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCg                                                                  <  1:336981/144‑1 (MQ=255)                                                                                 aaGTAATATGGTTTTCTGGCGG‑AAACCC‑TTTTTAAAATTTGTAGCCCGGCGCCCA‑GGTTATTGGTGCAAAACAAAACA‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCAccc                                                                   <  1:268986/146‑1 (MQ=255)                                                                                    tAATATAGGTTTCTGACGGCAAACCG‑TTTTTAATA‑TTGTAGTCCTGCGCCCA‑GGTCATTGGTGCGAAAACAAACA‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGcc                                                           >  1:10485/1‑151 (MQ=255)                                                                                     aaTATGGGTTTCTGACGGAAACCCT‑TTTTAAATT‑TTGTAGCCCTGCGCCCA‑GTTCTTTGGTGCAAAAACAACCG‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGcc                                                           <  1:454915/150‑1 (MQ=255)                                                                                       tataGGTTTCTGACGGCAAACCG‑TTTTTAATA‑TTGTAGTCCTGCGCCCA‑GGTCATTGGTGCGAAAACAAACA‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAAt                                                        >  2:345102/1‑151 (MQ=255)                                                                                         taGGTTTCTGACGGCAAACCG‑TTTTTAATA‑TTGTAGTCCTGCGCCCA‑GGTCATTGGTGCGAAAACAAACA‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATc                                                       >  1:170115/1‑150 (MQ=255)                                                                                           gttttCTGCCGAAAAACTT‑TTTTTAATT‑TTGTACCCCTCCCCCGA‑GGTAATGGGTGAAAAAAAAAACA‑‑GCCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCaa                                                         <  1:42327/144‑1 (MQ=255)                                                                                            gTTTCTGACGGCAAACCG‑TTTTTAATA‑TTGTAGTCCTGCGCCCA‑GGTCATTGGTGCGAAAACAAACA‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCAc                                                     >  2:1030546/1‑149 (MQ=255)                                                                                              ttCTGACGGCAAACCG‑TTTTTAATA‑TTGTAGTCCTGCGCCCA‑GGTCATTGGTGCGAAAACAAACA‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGc                                                  >  2:1048937/1‑150 (MQ=255)                                                                                              ttCTGACGGCAAACCG‑TTTTTAATA‑TTGTAGTCCTGCGCCCA‑GGTCATTGGTGCGAAAACAAACA‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGc                                                  >  2:84271/1‑150 (MQ=255)                                                                                                cTGACGGCAAACCT‑TTTTTAAAG‑TTGAACCCCGGCGCCGG‑GGTAATGGGTGGAAAAAAAAACA‑‑GCCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGc                                                  <  1:285052/148‑1 (MQ=255)                                                                                                 tGACGGCAAACCG‑TTTTTAATA‑TTGTAGTCCTGCGCCCA‑GGTCATTGGTGCGAAAACAAACA‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCt                                              >  2:220710/1‑151 (MQ=255)                                                                                                  gACGGCAAACCG‑TTTTTAATA‑TTGTAGTCCTGCGCCCA‑GGTCATTGGTGCGAAAACAAACA‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCATTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCtt                                             >  1:167055/1‑151 (MQ=255)                                                                                                    cGGCAACCCT‑TTTTTATTT‑TGTTAGCCCGGCCCCCA‑GTTCTTTGGTGCAAAAACAAACA‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTa                                           <  1:628162/151‑1 (MQ=255)                                                                                                    agggaaacctttttttAATA‑TTGTAGTCCTGCGCCCA‑GTTCTTGGGTGAAAAAACAAACA‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCa                                                 <  1:853459/135‑1 (MQ=255)                                                                                                     gacaacctg‑TTTTTATTT‑TTGTAGT‑CGGCCCCCAGGGTTTTTGGGGAAAAAAAAAACG‑‑GACCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCttt                                            <  2:709251/141‑1 (MQ=255)                                                                                                      gCAAACCG‑TTTTTAATA‑TTGTAGTCCTGCGCCCA‑GGTCATTGGTGCGAAAACAAACA‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAAt                                         >  2:420361/1‑151 (MQ=255)                                                                                                         aatctttttatATTA‑TTTTAGCCCCGCGCCAA‑GTTCA‑TGGT‑TGAAAACAAACACTGTCCCGCCAGTAAAGCAAATTTTTTCATCACCATTCCATATTTCATTTTACCCAGAATTGCACCCCCCGTACCGCCAATCACGGCACCTTTAATc                                        <  1:342809/144‑1 (MQ=255)                                                                                                          aCCG‑TTTTTAATA‑TTGTAGTCCTGCGCCCA‑GGTCATTGGTGCGAAAACAAACA‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGcc                                     >  2:828258/1‑151 (MQ=255)                                                                                                                         ttgttgTCCTGCCCCCG‑GTTTATTGGCGCGAAAACAAACG‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGa                              <  1:283072/145‑1 (MQ=255)                                                                                                                                    cccccA‑GTTCTTTGGTGCAAAAACAAACG‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGa             <  2:1062286/149‑1 (MQ=255)                                                                                                                                      cccA‑GGTGTTGGGGGCAAAAACAAACA‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGaa            <  1:719762/150‑1 (MQ=255)                                                                                                                                       ccaaggtttttGGGGCAAAAAAAAACG‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCAt                        <  2:907374/129‑1 (MQ=255)                                                                                                                                        cA‑GGTCATTGGTGCGAAAACAAACA‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGAACGc         >  1:667950/1‑151 (MQ=255)                                                                                                                                            gTCATTGGTGCGAAAACAAACA‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGAACGcc        >  1:720862/1‑149 (MQ=255)                                                                                                                                            gTCATTGGTGCGAAAACAAACA‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGAACGcc        >  1:153176/1‑149 (MQ=255)                                                                                                                                                  ggTGCGAAAACAAACA‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGAACGCCAGTGAc  >  2:550765/1‑149 (MQ=255)                                                                                                                                                       |                                                                                                                                                  AACGCTGGCATTACCCATTACCGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCG‑TTTTTAATA‑TTGTAGTCCTGCGCCCA‑GGTCATTGGTGCGAAAACAAACA‑‑GGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGAACGCCAGTGAC  >  NC_000913/1222449‑1222739

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