Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	121,507	0	A	T	44.6%	128.9 / 50.5	92	L15L (CTT→CTA)	aroP	aromatic amino acid transporter
Reads supporting (aligned to +/- strand):  ref base A (11/40);  new base T (41/0);  total (52/40)
Fisher's exact test for biased strand distribution p-value = 2.63e-16
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACGCATGATCAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGATAATCCCTGGCCCTGCGGACTGTATTACGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTAAGGCCGCGCTTTAGCTGCTCGCCGTGCTGTTGACCTTCCATCATGAAACCTCGTGCGGTGGTTGTTTTTTTGATCTACGCAGTGATGCGTGTGTAAGTTTGCAATTCCGTTTGTTGTATTAATTTGTTTACATCAAAGAAGTTTGAATTG  >  NC_000913/121359‑121657                                                                                                                                                     |                                                                                                                                                       gACGCATGATCAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGATAATCCCTGGCCCTGCGGACTGTATTACGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTa                                                                                                                                                        <  2:388600/149‑1 (MQ=255)  aCGCATGATCAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGATAATCCCTGGCCCTGCGGACTGTATTACGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTTAgg                                                                                                                                                     >  1:281995/1‑151 (MQ=255)           cAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGATAATCCCTGGCCCTGCGGACTGTATTACGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTAAGGCCGCGCttt                                                                                                                                            <  1:504959/151‑1 (MQ=255)           cAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGATAATCCCTGGCCCTGCGGACTGTATTACGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTAAGGCCGCGCttt                                                                                                                                            <  2:711589/151‑1 (MQ=255)           cAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGATAATCCCTGACCCTGCGGACTGTATTACGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTAAGGCCGCGCtt                                                                                                                                             <  2:277894/150‑1 (MQ=255)            aGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGATAATCCCTGGCCCTGCGGACTGTATTACGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTTAGGCCCCGCcttt                                                                                                                                           >  2:1380515/1‑150 (MQ=255)            aGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGATAATCCCTGGCCCTGCGGACTGTATTACGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTTAAGCCGCGCttt                                                                                                                                            >  1:1353090/1‑150 (MQ=255)             gAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGATAATCCCTGGCCCTGCGGACTGTATTACGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTAAGGCCGCGCTTTAg                                                                                                                                          <  2:217598/151‑1 (MQ=255)              aaaGGCGATAAAACCAGCAATGGCGTAACCCAGGATAATCCCTGGCCCTGCGGACTGTATTACGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTTAGGCCGCGCttt                                                                                                                                            >  2:1348735/1‑148 (MQ=255)                aGGCGATAAAGCCAGCAACGGCGTAACCCAGGATAATCCCTGGCCCTGCGGACTGTATTACGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGGTTTTAAGGCCGCGCTTTAgct                                                                                                                                        <  2:1113713/150‑1 (MQ=255)                    gATAAAACCAGCAATGGCGTAACCCAGGATAATCCCTGGCCCTGCGGACTGTATTACGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTAAGGCCGCGCTTTAGCTGCTCg                                                                                                                                   <  1:875087/151‑1 (MQ=255)                     aTAAAACCAGCAATGGCGTAACCCAGGATAATCCCTGGCCCTGCGGACTGTATTACGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTTAGGgcgccattt                                                                                                                                            >  1:1243102/1‑136 (MQ=255)                     aTAAAACCAGCAATGGCGTAACCCAGGATAATCCCTGGCCCTGCGGACTGTATTACGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTAAGGCCGCGCTTTAGCTGCTCGc                                                                                                                                  <  2:830359/151‑1 (MQ=255)                       aaaaCCAGCAATGGCGTAACCCAGGATAATCCCTGGCCCTGCGGACTGTATTACGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTAAGGCCGCGCTTTAGCTGCTCGcc                                                                                                                                 <  1:664759/150‑1 (MQ=255)                            cAGCAATGGCGTAACCCAGGATAATCCCTGGCCCTGCGGACTGTATTACGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTTAGGCCGCGCTTTTGCTGCTCGcc                                                                                                                                 >  2:597352/1‑145 (MQ=255)                             aGCAATGGCGTAACCCAGGATAATCCCTGGCCCTGCGGACTGTATTACGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGGTTTTAAAGGCGCCCCTTTGCTGCTccgcggg                                                                                                                              >  2:1112391/1‑141 (MQ=255)                              gCAATGGCGTAACCCAGGATAATCCCTGGCCCTGCGGACTGTATTACGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTAAGGCCGCGCTTTAGCTGCTCGCCGTGCTGtt                                                                                                                         <  2:1055864/151‑1 (MQ=255)                               cAATGGCGTAACCCAGGATAATCCCTGGCCCTGCGGACTGTATTACGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTTAGGCCGCGCTTTTGCTGCCCGCCGGGCTGtt                                                                                                                         >  1:1017851/1‑150 (MQ=255)                               cAATGGCGTAACCCAGGATAATCCCTGGCCCTGCGGACTGTATTACGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTAAGGCCGCGCTTTAGCTGCTCGCCGTGCTGTTg                                                                                                                        <  2:490045/151‑1 (MQ=255)                                aaTGGCGTAACCCAGGATAATCCCTGGCCCTGCGGACTGTATTACGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTAAGGCCGCGCTTTAGCTGCTCGCCGTGCTGTTGa                                                                                                                       <  1:79843/151‑1 (MQ=255)                                  tGGCGTAACCCAGGATAATCCCTGGCCCTGCGGACTGTATTACGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTTAGGCCGCGCTTTTGCTGCTCGCCGTGCTGTTGAcc                                                                                                                     >  2:398294/1‑151 (MQ=255)                                     cGTAACCCAGGATAATCCCTGGCCCTGCGGACTGTATTACGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTTAAGCCCCGCCTTAGCTGCTCGCCCTGCTTTTGACCtt                                                                                                                   >  1:496151/1‑150 (MQ=255)                                               gATAATCCCTGGCCCTGCGGACTGTATTACGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTTAGGCCGCGCTTTTGCTGCTCGCCGGGCTGTTGGCCTTCCCTCATGaa                                                                                                         >  2:662773/1‑150 (MQ=255)                                                aTAATCCCTGGCCCTGCGGACTGTATTACGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTTAGGCCCCGCTTTAGCTGCTCGCCGCGCTGTTGGCCTTCCATCATGaaa                                                                                                        >  1:106403/1‑150 (MQ=255)                                                  aaTCCCTGGCCCTGCGGACTGTATTACGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTTAAGCCGCGCCTTAGCTGCTCGCCGCGCTGTTGACCTTCCCTCATGaaa                                                                                                        >  2:1272641/1‑148 (MQ=255)                                                  aaTCCCTGGCCCTGCGGACTGTATTACGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTAAGGCCGCGCTTTAGCTGCTCGCCGTGCTGTTGACCTTCCATCATGAAAcc                                                                                                      <  2:288524/150‑1 (MQ=255)                                                  aaTCCCTGGCCCTGCGGACTGTATTACGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTAAGGCCGCGCTTTAGCCGGCCCCCGCGCTGTTGGCCCTCCCTCCTGaaa                                                                                                        >  2:531655/1‑148 (MQ=255)                                                   aTCCCTGGCCCTGCGGACTGTATTACGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTTAGGCCGCGCTTTTGCTGCTCGCCGTGCTGTTTACCTTCCATCATGAAAcccc                                                                                                    >  2:564390/1‑149 (MQ=255)                                                   aTCCCTGGCCCTGCGGACTGTATTACGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTAAGGCCGCGCTTTAGCTGCTCGCCGTGCTGTTGACCTTCCATCATGAAACCTc                                                                                                    <  2:102347/151‑1 (MQ=255)                                                    tCCCTGGCCCTGCGGACTGTATTACGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTTAAGCCGCGCTTTAGGTGCTCGCCGTGCTGTTGACCTTCCCTCCTGAAACCTc                                                                                                    >  2:1139730/1‑150 (MQ=255)                                                       cTGGCCCTGCGGACTGTATTACGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTTAAGCCGCGCCTTTACTGCTCCCCCCGGCGTTGGCCTTcccttctgtaaac                                                                                                      >  2:446651/1‑134 (MQ=255)                                                        tGGCCCTGCGGACTGTATTACGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTTAAGGCGCGCTTTTGCTGCTCGCCCGGCTGTTTGCCTTCCATCATGAAACCTCGTg                                                                                                 >  1:608562/1‑149 (MQ=255)                                                         ggCCCTGCGGACTGTATTACGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTAAGGCCGCGCTTTAGCTGCTCGCCGTGCTGTTGACCTTCCATCATGAAACCTCGTGc                                                                                                <  2:990917/149‑1 (MQ=255)                                                          gCCCTGCGGACTGTATTACGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTAAGGCCGCGCTTTAGCTGCTCGCCGTGCTGTTGACCTTCCATCATGAAACCTCGTGCggt                                                                                             <  1:899658/151‑1 (MQ=255)                                                            ccTGCGGACTGTATTACGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTTAGGCCGCGCTTTTGCTGCTCCCCGGGCTGTTGGCCTTCCCTCATGAAACCTCGTGCggggg                                                                                           >  1:1407633/1‑151 (MQ=255)                                                               gCGGACTGTATTACGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTAAGGCCGCGCTTTAGCTGCTCGCCGTGCTGTTGACCTTCCATCATGAAACCTCGTGCggtggt                                                                                          <  2:649150/149‑1 (MQ=255)                                                               gCGGACTGTATTACGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTAAGGCCGCGCTTTAGCTGCTCGCCGTGCTGTTGACCTTCCATCATGAAACCTCGTGCGGTGgtt                                                                                         <  1:798661/150‑1 (MQ=255)                                                                  gACTGTATTACGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTAAGGCCGCGCTTTAGCTGCTCGCCGTGCTGTTGACCTTCCATCATGAAACCTCGTGCGGTGgttgtt                                                                                      <  2:658966/150‑1 (MQ=255)                                                                   aCTGTATTACGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTTAGGGCGCGCTTTAACTGCCCGCCGTGCGGTTGGCCCTTCATCATGGAACCCCGTGCGGGGGGTGtttt                                                                                    >  1:1143685/1‑151 (MQ=255)                                                                   aCTGTATTACGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTTAGGCCGCGCTTTAACTGCTCGCCGTGCTGTTTACCTTTCCTCATGGAACCTCGTGCGGGGGGTGtttt                                                                                    >  2:149356/1‑151 (MQ=255)                                                                   aCTGTATTACGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTAAGGCCGCGCTTTAGCTGCTCGCCGTGCTGTTGACCTTCCATCATGAAACCTCGTGCGGTGGTTGtttt                                                                                    <  2:184115/151‑1 (MQ=255)                                                                      gTATTACGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTAAGGCCGCGCTTTAGCTGCTCGCCGTGCTGTTGACCTTCCATCATGAAACCTCGTGCGGTGGTTGttttttt                                                                                 <  2:1459600/151‑1 (MQ=255)                                                                       tATTACGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTTAAGGCGCGCCTTTGCTGCTCGCCGGGCTGTTGACCTTTCATCATTGAAACTCGGGCGGGGGGTGttttttt                                                                                 >  1:66206/1‑150 (MQ=255)                                                                        aTTACGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTTAGGCCGCGCCTTAGCTGCTCGCCGTGGTGTTGACCTTCCATCATGAAACCTCGGGGGGGGGGTGttttttt                                                                                 >  1:518971/1‑149 (MQ=255)                                                                            cGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTAAGGCCGCGCTTTAGCTGCTCGCCGTGCTGTTGACCTTCCATCATGAAACCTCGTGCGGTGGTTGTTTTTTTGATCTa                                                                           <  2:1243003/151‑1 (MQ=255)                                                                             ggAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTAAAGCCGCGCTTTAGCTGCCCGCCGTGCTGTTGACCTTCCATCATGAAAACTCCTGCGGGGGGTGttttttttgtct                                                                            >  1:362380/1‑144 (MQ=255)                                                                              gAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCCCCCAGCGCGATAAGCTGAATATGGCGGTTTTTAAGGCCGCGCTTTAGCTGCTCGCCGTGCTGTTGACCTTCCATCATGAAACCTCGTGCGGTGGTTGTTTTTTTGATCTACg                                                                         <  2:314148/151‑1 (MQ=255)                                                                              gAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTAAGGCCGCGCTTTAGCTGCTCGCCGTGCTGTTGACCTTCCATCATGAAACCTCGTGCGGTGGTTGTTTTTTTGATCTACg                                                                         <  2:594907/151‑1 (MQ=255)                                                                               aGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTTAGGCCGCGCCTTTGCTGCTCGCCGTGCTGTTGACCTTCCATCATGAAACCTCGTGCGGGGGTTGTTTTTTTTTTCTACg                                                                         >  1:489800/1‑150 (MQ=255)                                                                                  cgcTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTTAGGCCGCGCTTTTGCTGCTCGCCGTGCTGTTGACCTTCCATCATGAAACCTCGTGCGGGGGGTGTTTTTTTTATCTACGccg                                                                      >  1:1126023/1‑148 (MQ=255)                                                                                  cgcTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTTAGGCCGCGCCTTAGCCGCCCGCCGCGCTGGTGACCTTCCCTCCTGAAACCCCGTGCGGGGGTTGGTTTTTTTTTCTACGc                                                                        >  2:451495/1‑148 (MQ=255)                                                                                    cTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTTAGGCCGCGCTTTAGCTGCTCGCCGTGCTGTTGACCCTCCATCATGAAACCTCGTGCGGGGGTTGTTTTTTTTTTCTACGCCGTg                                                                    >  2:1014360/1‑150 (MQ=255)                                                                                     tACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTAAGGGCGCGCTTTTGCTGCTCGCCGCGCTGTTGGGCTTCCCTCATGAAAACTCGTGCGGGGGGTGTTTTTTTGGTCTTCGCAGTg                                                                    >  1:30330/1‑149 (MQ=255)                                                                                       cccAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTTAGGCCGCGCTTTAGCTGCTCGCCGTGCTGTTGACCTTCCCTCATGAAACCTCGTGCGGGGGGTGTTTTTTTTATTTACGCAGTg                                                                    >  1:395612/1‑147 (MQ=255)                                                                                              aTAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGCTTTTAAGGCCGCGCTTTAGCTGCTCGCCGTGCTGTTGACCTTCCATCATGAAACCTCGTGCGGTGGTTGTTTTTTTGATCTACGCAGTGATGCGTGTGTa                                                         <  2:1227116/151‑1 (MQ=255)                                                                                                aaCCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTTAAGGCGCGCTTTTGCTGCTCGCCGTGCTGTTGGCCTTTCCTCCTGAAAACCCCTGCGGGGGGTGTTTTTTTTTTTTTCGCCGTGGTGCGGGGGTaa                                                        >  2:876305/1‑150 (MQ=255)                                                                                                aaCCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTTAAGGCGCGCTTTAGCTGCTCGCCCCGCTGTTGACCTTCCATCATGAAAACTCGTGCGGGGGGTGTTTTTTTGATTTACGCAGTGATGCgtggg                                                           >  2:780370/1‑145 (MQ=255)                                                                                                aaCCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTTAAGCCGCGCTTTAGCTGCTCGCCGGGCTGTTGGCCTTCCCTCATTAAAACTCGTGCGGGGGGTGTTTTTTTTTTCTACGCAGTGATGCgtgtgt                                                          >  2:42427/1‑148 (MQ=255)                                                                                                aaCCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTAAGGCCGCGCTTTAGCTGCTCGCCGTGCTGTTGACCTTCCATCATGAAACCTCGTGCGGTGGTTGTTTTTTTGATCTACGCAGTGATGCGTGTGTaa                                                        <  1:959559/150‑1 (MQ=255)                                                                                                   ccGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTAAGGCCGCGCTTTAGCTGCTCGCCGTGCTGTTGACCTTCCATCATGAAACCTCGTGCGGTGGTTGTTTTTTTGATCTACGCAGTGATGCGTGTGTAAGttt                                                    <  2:1285822/151‑1 (MQ=255)                                                                                                   ccGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTAAGGCCGCGCTTTAGCTGCTCGCCGTGCTGTTGACCTTCCATCATGAAACCTCGTGCGGTGGTTGTTTTTTTGATCTACGCAGTGATGCGTGTGTAAGttt                                                    <  1:1303807/151‑1 (MQ=255)                                                                                                   ccGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTAAGGCCGCGCTTTAGCTGCTCGCCGTGCTGTTGACCTTCCATCATGAAACCTCGTGCGGTGGTTGTTTTTTTGATCTACGCAGTGATGCGTGTGTAAGttt                                                    <  1:1450153/151‑1 (MQ=255)                                                                                                    cGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTTAGGCCGCGCCTTAGCTGCTCGCCGGGCTGTTGAACTTCCATCATGAAACCTCCTGCGGGGGGTGTTTTTTTTTTTTTCGCCGTGGTGGGGGTGTAAGttt                                                    >  2:793658/1‑150 (MQ=255)                                                                                                    cGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTAAGGCCGCGCTTTAGCTGCTCGCCGTGCTGTTGACCTTCCATCATGAAACCTCGTGCGGTGGTTGTTTTTTTGATCTACGCAGTGATGCGTGTGTAAGTTTg                                                   <  1:1046313/151‑1 (MQ=255)                                                                                                     ggTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTTAGGCCGCGCTTTAGCTGCTCGCCGTGCTGTTGGCCCTCCCTCATGAAACCTCGTGCGGGGGGTGTTTTTTTTATTTACGCAGGGATGCGTGTGTAAGttt                                                    >  2:802227/1‑149 (MQ=255)                                                                                                     ggTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTAAGGCCGCGCTTTAGCTGCTCGCCGTGCTGTTGACCTTCCATCATGAAACCTCGTGCGGTGGTTGTTTTTTTGATCTACGCAGTGATGCGTGTGTAAGTTTGc                                                  <  1:147335/151‑1 (MQ=255)                                                                                                           tATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTTAGGCCGCGCTTTTGCTGCTCGCCGTGCTGTTGACCTTCCCTCATGAAAACTCGTGCGGGGGGTGTTTTTTTGATTTACGCAGTGGTGCGTGGGTAAGTTTGCAAtt                                              >  1:884292/1‑149 (MQ=255)                                                                                                             tCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTTAGGCCGCGCTTTTGCTGCTCGCCGCGCTGTTGACCTTCCATCATGAAAACTCGTGCGGGGGGTGTTTTTTTTATCTTCGCAGGGGTGGGTGGGTAAGTTTGCAATTc                                             >  2:1049691/1‑148 (MQ=255)                                                                                                                cgcCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTAAGGCCGCGCTTTAGCTGCTCGCCGTGGTGTTGACCTTCCCTCATGAAAACTCGTGCGGGGGGTGTTTTTTTTTTCTACGCAGGGGTGCGGGTGTAAGTTTGCAATTTCGttt                                        >  1:1000045/1‑150 (MQ=255)                                                                                                                   cACCCAGCGCGATAAGCTGAATATGGCGGTTTTTTAGGCCGCGCTTTAGGTGCTCGCCCGGCTGTTGACCCTCCATCATGAAACCTCGTGCGGGGGGTGTTTTTTTTATCTACCCCGTGATGGGTGTGTAAGGTTTCAATTCCGtttttt                                     >  1:1072202/1‑150 (MQ=255)                                                                                                                        aGCGCGATAAGCTGAATATGGCGGTTTTTTAGGCCGCGCCTTTGCTGCTCGCCGTGCTGTTGGCCTTCCATCATGAAAACCCGTGCGGGGGGTGTTTTTTTTATTTACGCCGGGATGCGGGTGTAAGTTTTTAATTTCGTttgttttttt                                >  2:1217476/1‑145 (MQ=255)                                                                                                                         gcgcGATAAGCTGAATATGGCGGTTTTTTAAGCCCCGCCTTAGCTGCCCCCCGGTCTGTTGGCCTTCCCTCATGAAACCTCGTGCGGGGGGTGTTTTTTTTTTCTACGCAGTGGTGAGTGTGTTAGTTTGCAATTTCGTTTTTTGtttt                                >  2:1283552/1‑149 (MQ=255)                                                                                                                          cgcgATAAGCTGAATATGGCGGTTTTTAAGGCCGCGCTTTAGCTGCTCGCCGTGCTGTTGACCTTCCATCATGAAACCTCGTGCGGTGGTTGTTTTTTTGATCTACGCAGTGATGCGTGTGTAAGTTTGCAATTCCGTTTGTTGTATTAAt                             <  2:1460800/151‑1 (MQ=255)                                                                                                                          cgcgATAAGCTGAATATGGCGGGTTTTTAAGCCGCGCTTTAGCTGGTCGCCCTGGTGTTTGCCTTCCATTATGAAAAATCGGGGGGGGGGTGTTTTTTTTTTCTACGCAGTGGTGCGTGTGGAAGGTTTCAATTTCGTTTTTTGTAtt                                >  2:539453/1‑148 (MQ=255)                                                                                                                           gcgATAAGCTGAATATGGCGGTTTTTTAAGCCGCGCTTTAGCTGCTCGCCGGGCTGTTGGCCTTTCATCATGAAAACTCGTGCGGGGGGTGTTTTTTTTATCTACGCCGTGATGCGTGTGTAAGTTTGCAATTCCGTTTGTTGTATTAAtt                            >  1:756392/1‑151 (MQ=255)                                                                                                                            cgATAAGCTGAATATGGCGGTTTTTAAGGCCGCGCTTTAGCTGCTCGCCGTGCTGTTGACCTTCCATCATGAAACCTCGTGCGGTGGTTGTTTTTTTGATCTACGCAGTGATGCGTGTGTAAGTTTGCAATTCCGTTTGTTGTATTAAttt                           <  2:473663/151‑1 (MQ=255)                                                                                                                              aTAAGCTGAATATGGCGGTTTTTAAGGCCGCGCTTTAGCTGCTCGCCGTGCTGTTGACCTTCCATCATGAAACCTCGTGCGGTGGTTGTTTTTTTGATCTACGCAGTGATGCGTGTGTAAGTTTGCAATTCCGTTTGTTGTATTAATTTGt                         <  1:70503/151‑1 (MQ=255)                                                                                                                                aaGCTGAATATGGCGGTTTTTTAGGCCGCGCTTTAGCTGCTCGCCGGGCTGTTGACCTTCCCTCATTAAACCTCGTGCGGGGGTTGTTTTTTTTTTCTACGCAGTGATGCGTGTGGAAGTTTGCCATTTCGTTTGTTGTATTAATTTGttt                       >  2:297896/1‑151 (MQ=255)                                                                                                                                aaGCTGAATATGGCGGTTTTTAAGGCCGCGCTTTAGCTGCTCGCCGTGCTGTTGACCTTCCATCATGAAACCTCGTGCGGTGGTTGTTTTTTTGATCTACGCAGTGATGCGTGTGTAAGTTTGCAATTCCGTTTGTTGTATTAATTTGtt                        <  2:312517/150‑1 (MQ=255)                                                                                                                                aaGCTGAATATGGCGGTTTTTAAGGCCGCGCTTTAGCTGCTCGCCGTGCTGTTGACCTTCCATCATGAAACCTCGTGCGGTGGTTGTTTTTTTGATCTACGCAGTGATGCGTGTGTAAGTTTGCAATTCCGTTTGTTGTATTAATTTGtt                        <  2:1004652/150‑1 (MQ=255)                                                                                                                                aaGCTGAATATGGCGGTTTTTAAGGCCGCGCTTTAGCTGCTCGCCGTGCTGTTGACCTTCCATCATGAAACCTCGTGCGGTGGTTGTTTTTTTGATCTACGCAGTGATGCGTGTGTAAGTTTGCAATTCCGTTTGTTGTATTAATTTGtt                        <  1:802306/150‑1 (MQ=255)                                                                                                                                  gCTGAATATGGCGGTTTTTTAGGCCGCGCCTTTGCTGCTCGCCGTGGTGTTTACCTTCCCTCATGAAACCTCGTGCGGGGGTTGTTTTTTTTTTTTTCGCCGTGATGCGTGTGTAAGTTTGTAATTCCGTTTGTTGTATTAATtttttt                       >  2:784356/1‑149 (MQ=255)                                                                                                                                  gCTGAATATGGCGGTTTTTTAAGCCGCGCTTTTGCTGCTCGCCGTGCTGTTGACCTTCCCTCATGAAACCTCGTGCGGGGGGTGTTTTTTTTATCTACGCAGTGATGCGTGTGTAAGTTTTCAATTCCGTTTGTTGTATTAATTTGttt                       >  2:639163/1‑149 (MQ=255)                                                                                                                                     gAATATGGCGGTTTTTAAGGCCGCGCTTTAGCTGCTCGCCGTGCTGTTGACCTTCCATCATGAAACCTCGTGCGGTGGTTGTTTTTTTGATCTACGCAGTGATGCGTGTGTAAGTTTGCAATTCCGTTTGTTGTATTAATTTGTTTACa                    <  1:979417/149‑1 (MQ=255)                                                                                                                                     gAATATGGCGGTTTTTAAGGCCGCGCTTTAGCTGCTCGCCGTGCTGTTGACCTTCCATCATGAAACCTCGTGCGGTGGTTGTTTTTTTGATCTACGCAGTGATGCGTGTGTAAGTTTGCAATTCCGTTTGTTGTATTAATTTGTTTACa                    <  1:1380619/149‑1 (MQ=255)                                                                                                                                          tGGCGGTTTTTAAGGCCGCGCTTTAGCTGCTCGCCGTGCTGTTGACCTTCCATCATGAAACCTCGTGCGGTGGTTGTTTTTTTGATCTACGCAGTGATGCGTGTGTAAGTTTGCAATTCCGTTTGTTGTATTAATTTGTTTACATCaa                >  1:756690/1‑148 (MQ=255)                                                                                                                                          tGGCGGTTTTTAAGGCCGCGCTTTAGCTGCTCGCCGTGCTGTTGACCTTCCATCATGAAACCTCGTGCGGTGGTTGTTTTTTTGATCTACGCAGTGATGCGTGTGTAAGTTTGCAATTCCGTTTGTTGTATTAATTTGTTTACATCaa                >  2:1030811/1‑148 (MQ=255)                                                                                                                                           ggCGGTTTTTAAGGCCGCGCTTTAGCTGCTCGCCGTGCTGTTGACCTTCCATCATGAAACCTCGTGCGGTGGTTGTTTTTTTGATCTACGCAGTGATGCGTGTGTAAGTTTGCAATTCCGTTTGTTGTATTAATTTGTTTACATCAaaga             >  2:285441/1‑150 (MQ=255)                                                                                                                                            gCGGTTTTTAAGGCCGCGCTTTAGCTGCTCGCCGTGCTGTTGACCTTCCATCATGAAACCTCGTGCGGTGGTTGTTTTTTTGATCTACGCAGTGATGCGTGTGTAAGTTTGCAATTCCGTTTGTTGTATTAATTTGTTTACATCAaaga             >  2:1173819/1‑149 (MQ=255)                                                                                                                                             cGGTTTTTAAGGCCGCGCTTTAGCTGCTCGCCGTGCTGTTGACCTTCCATCATGAAACCTCGTGCGGTGGTTGTTTTTTTGATCTACGCAGTGATGCGTGTGTAAGTTTGCAATTCCGTTTGTTGTATTAATTTGTTTACATCAAAGAAGt          <  2:324846/151‑1 (MQ=255)                                                                                                                                              ggTTTTTAAGGCCGCGCTTTAGCTGCTCGCCGTGCTGTTGACCTTCCATCATGAAACCTCGTGCGGTGGTTGTTTTTTTGATCTACGCAGTGATGCGTGTGTAAGTTTGCAATTCCGTTTGTTGTATTAATTTGTTTACATCAAATAAGtt         >  1:860291/1‑151 (MQ=255)                                                                                                                                                 ttttAAGGCCGCGCTTTAGCTGCTCGCCGTGCTGTTGACCTTCCATCATGAAACCTCGTGCGGTGGTTGTTTTTTTGATCTACGCAGTGATGCGTGTGTAAGTTTGCAATTCCGTTTGTTGTATTAATTTGTTTACATCAAAGAAGTTTg       >  2:1120106/1‑150 (MQ=255)                                                                                                                                                    tAAGGCCGCGCTTTAGCTGCTCGCCGTGCTGTTGACCTTCCATCATGAAACCTCGTGCGGTGGTTGTTTTTTTGATCTACGCAGTGATGCGTGTGTAAGTTTGCAATTCCGTTTGTTGTATTAATTTGTTTACATCAAAGAAGTTTGAAt    <  1:362329/150‑1 (MQ=255)                                                                                                                                                     aaGGCCGCGCTTTAGCTGCTCGCCGTGCTGTTGACCTTCCATCATGAAACCTCGTGCGGTGGTTGTTTTTTTGATCTACGCAGTGATGCGTGTGTAAGTTTGCAATTCCGTTTGTTGTATTAATTTGTTTACATCAAAGAAGTTTGAATTg  >  1:1467589/1‑151 (MQ=255)                                                                                                                                                     aaGGCCGCGCTTTAGCTGCTCGCCGTGCTGTTGACCTTCCATCATGAAACCTCGTGCGGTGGTTGTTTTTTTGATCTACGCAGTGATGCGTGTGTAAGTTTGCAATTCCGTTTGTTGTATTAATTTGTTTACATCAAAGAAGTTTGAATTg  <  1:1135842/151‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                       GACGCATGATCAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGATAATCCCTGGCCCTGCGGACTGTATTACGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTAAGGCCGCGCTTTAGCTGCTCGCCGTGCTGTTGACCTTCCATCATGAAACCTCGTGCGGTGGTTGTTTTTTTGATCTACGCAGTGATGCGTGTGTAAGTTTGCAATTCCGTTTGTTGTATTAATTTGTTTACATCAAAGAAGTTTGAATTG  >  NC_000913/121359‑121657

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