Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,717,014	0	T	G	39.3%	167.3 / 14.0	90	G156G (GGT→GGG)	yiaD	multicopy suppressor of bamB; outer membrane lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (20/34);  new base G (34/1);  total (55/35)
Fisher's exact test for biased strand distribution p-value = 1.83e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATATGCCGAACAATGTGACCTTCGACAGCAGCAGCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGTTAACGTGA‑TTGGTTATACCGACAGC‑ACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGCGTTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTACTCAGGGCCTTGGCCCGGCTAACCCAATCGCCAGCAACAGCACCGCAGAAGG  >  NC_000913/3716872‑3717160                                                                                                                                                  |                                                                                                                                                   atatGCCGAACAATGTGACCTTCGACAGCAGCAGCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGTTAACGTGA‑TTGGTTATACCGACAGC‑ACGGGTGGTCACGa                                                                                                                                            <  1:526829/151‑1 (MQ=255)    tGCCGAACCATGTGACCTTCGACAGCAGCAGCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGGTAACGCGA‑TTGGGTATATCCACAAG‑AAGGGGGGGCACGagcc                                                                                                                                         >  2:524442/1‑148 (MQ=255)     gCCGAACAATGTGACCTTCGACAGCAGCAGCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGTTAACGTGA‑TTGGTTATACCGACAGC‑ACGGGTGGTCACGACCt                                                                                                                                         <  1:381797/150‑1 (MQ=255)         aaCAATGTGACCTTCGACAGCAGCAGCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGTTAACGTGA‑TTGGTTATACCGACAGC‑ACGGGTGGTCACGACCTGAACa                                                                                                                                    <  2:1642651/151‑1 (MQ=255)         aaCAATGTGACCTTCGACAGCAGCAGCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGGTAAAGGGG‑TTTGGTTTACCCAAAGC‑ACCGGGGGGCACGGcccg                                                                                                                                        >  2:503614/1‑145 (MQ=255)           cAATGTGACCTTCGACAGCAGCAGCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGTTAACGTGA‑TTGGTTATACCGACAGC‑ACGGGTGGTCACGACCTGAACATg                                                                                                                                  <  1:725165/151‑1 (MQ=255)            aaTGTGACCTTCGACAGCAGCAGCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGCGGGTTAACTGGT‑TTGGTT‑TAACCACAGC‑ACGGGGGGGCACGACCCGaa                                                                                                                                      >  1:499306/1‑146 (MQ=255)            aaTGTGACCTTCGACAGCAGCAGCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGTTAACGTGA‑TTGGTTATACCGACAGC‑ACGGGTGGTCACGACCTGAACATGc                                                                                                                                 <  1:510187/151‑1 (MQ=255)             aTGTGACCTTCGACAGCAGCAGCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGTTAACGTGA‑TTGGTTATACCGACAGC‑ACGGGTGGTCACGACCTGAACATGCg                                                                                                                                <  2:1322467/151‑1 (MQ=255)                 gACCTTCGACAGCAGCAGCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGTTTACGGGA‑TTGGGTATAACCACAGC‑ACGGGGGGGCACGACCTGAAAATGCGt                                                                                                                               >  1:992386/1‑148 (MQ=255)                    cTTCGACAGCAGCAGCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGTTAACGGGA‑TTGGGTATAACGACAAC‑ACGGGGGGGCACCACCCGAACAAGCGTCTCTcc                                                                                                                         >  1:442116/1‑151 (MQ=255)                     ttCGACAGCAGCAGCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGTTTACGTGG‑TTGGGTATACCGACAAC‑ACGGGTGGGCACGGCCCGAACATGCGGCTCTccc                                                                                                                        >  2:1354815/1‑151 (MQ=255)                      tCGACAGCAGCAGCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGGTTAAGTGA‑TTTGGTATAACCCCAGC‑AAGGGGGGGCACCACCTCAAAAATCGTCTCTccc                                                                                                                        >  1:1306866/1‑150 (MQ=255)                       cGACAGCAGCAGCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGTTAACGTGA‑TTGGTTATACCGACAGC‑ACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAg                                                                                                                      <  2:209207/151‑1 (MQ=255)                        gACAGCAGCAGCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGTTAACATGA‑TTGGGTATAACCACAGC‑CCAGGGGGGCACGACCCGGACAAGCGTCTCTcccccc                                                                                                                     >  1:1385562/1‑148 (MQ=255)                        gACAGCAGCAGCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGGTAACGCGA‑TTGGGTTAACCGACAGC‑CCCGCGGGTGACGACCTGAAAAAACGCCTCTccc                                                                                                                        >  2:1522862/1‑148 (MQ=255)                         aCAGCAGCAGCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGTTTACGTGG‑TTGGGTATAACGACAAC‑ACGGGGGGTCACGACCTGAAAATGCCGCTCTCCCCGc                                                                                                                     >  1:1338842/1‑150 (MQ=255)                         aCAGCAGCAGCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGTTAACGTGA‑TTGGTTATACCGACAGC‑ACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCa                                                                                                                    <  2:1606274/151‑1 (MQ=255)                          cagcagcagCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGTTAACGTGA‑TTGGTTATACCGACAGC‑ACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCaa                                                                                                                   <  1:502791/151‑1 (MQ=255)                          cagcagcagCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGTTAACGTGA‑TTGGTTATACCGACAGC‑ACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCa                                                                                                                    <  1:549980/150‑1 (MQ=255)                           agcagcagcGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGTTAACGAGA‑TTGGTTATAACGACAGC‑ACGGGGGGTCCCGCCCCGAAAATGCGGCTCTCCCAGCaa                                                                                                                   >  1:1404521/1‑150 (MQ=255)                             cagcagcGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGTTAACGTGA‑TTGGTTATACCGACAGC‑ACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGt                                                                                                                <  1:1436504/151‑1 (MQ=255)                             cagcagcGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGGTAAACCGG‑TTGGGTTTAACGACAGCAAGGGGGGGCCACGACCAG‑ACAAGCGGCTCTccccacaa                                                                                                                   >  2:30794/1‑143 (MQ=255)                              agcagcGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCCATGGTACTGAAAGAGTATCCGAAAACGGCGGGTGAACG‑GG‑CTGGGTATACAGCCAGA‑ACCGGGGGTCACCAACTGAACATTCGGCTCTccc                                                                                                                        >  2:1366505/1‑142 (MQ=255)                              agcagcGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGTTTACGGGATTTGGTTTTCCCG‑CCGC‑ACGGGGGGGGACGCCCTGTAAATGCGTCTCTcccccccacg                                                                                                                 >  1:1502028/1‑142 (MQ=255)                                cagcGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGTTAACGAGA‑ATGGGTATAACCACCGC‑AACGGTGGGTCCGACCTGCACACTCGTCCCTCCCCGc                                                                                                                     >  2:713690/1‑143 (MQ=255)                                  gcgcGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGTTTAAGTGA‑TTGGGTTTAACGACAGC‑ACGGGGGGGCACGACCTGAACATGCGTCTCTCCCAGCAACGGGCgg                                                                                                            >  1:892299/1‑150 (MQ=255)                                   cgcgACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGTTAACGTGA‑TTGGTTATACCGACAGC‑ACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGAt                                                                                                          <  2:329410/151‑1 (MQ=255)                                      gACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGTTTACCTTT‑TTGGGTATAACGACAGC‑ACGGGGGGGCCCGCACTGAACAATGGGCTCCCCCAGCAACCGGCGGGtttc                                                                                                       >  1:329320/1‑149 (MQ=255)                                      gACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGTTAAAGTGA‑TTGGGTTTAACCACAGC‑ACGGGGGGGCACGACCCGAAAAAGCGTCTCTCCCCGCAAAGGGCGGGTTcc                                                                                                       >  1:870115/1‑151 (MQ=255)                                       aCCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGTTAAAGGGA‑TTGGTTATACCGACAAC‑ACCGGGGGGCACGACCCTAAAATGCGTCTCTCCCCGCAACCTGCGGGtttc                                                                                                       >  1:671066/1‑148 (MQ=255)                                           tGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGTTAACGTGA‑TTGGGTATACCGACAAC‑ACGGGGGGGCACGACCTGAAAATGCGTCTCTCCCCGCAACGTGCGGATTCCGtt                                                                                                    >  1:1492909/1‑149 (MQ=255)                                            gAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGTTAACGTGA‑TTGGTTATACCGACAGC‑ACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTg                                                                                                   <  1:1181956/149‑1 (MQ=255)                                             aaaCCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGTTAACATGA‑ATGGGGATTCCGGACCC‑ACCGGGGGGCACGCCCCCAACACGCGTCTCTCCCTGCAACGtcccggtgccgt                                                                                                     >  2:60189/1‑135 (MQ=255)                                                 cggcggGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGACAGAGTATCCGAAAACGGC‑GGTTCAACCGA‑TTGGTTATAACCACAGC‑ACGGGAGGGCACGCACCGAACCAGCGGCTCTCCCCCCAAAGGGGGGATTTCTTTGCCAgc                                                                                              >  2:367262/1‑149 (MQ=255)                                                 cggcggGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGTTAACGTGA‑TTGGTTTTACCCACAGC‑ACGGGGGGGCACGACCCGAACATGCGTCTCTCCCCGCCAAGTGCGGGTTCCCTTTCCcgcg                                                                                             >  1:224762/1‑150 (MQ=255)                                                   gcgggcgCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGCTAACGTGA‑TTGGGTATAACGACAGC‑ACGGGGGGGCACGGCCTGAACAAGCGTCTCTCCCCGCAACCGGCGGGTTCCGTTTCCa                                                                                                >  2:452975/1‑145 (MQ=255)                                                       gcgcTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGTTAACGTGA‑TTGGGTATACCGACAGC‑ACGGGGGGTCACGACCTGAACAAGCGTCTCTCCCCGCAACGTGCGGATTCCGTTGCCAGCGCGGTg                                                                                        >  1:458843/1‑149 (MQ=255)                                                        cgcTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGTTAACGTGA‑TTGGTTATACCGACAGC‑ACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGCGTTGATc                                                                                     <  1:1044921/151‑1 (MQ=255)                                                        cgcTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGTTAACGTGA‑TTGGTTATACCGACAGC‑ACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGCGTTGATc                                                                                     <  1:569365/151‑1 (MQ=255)                                                         gcTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGTTAAAGTGA‑TTGGGTATACCGACAGC‑ACGGGGGGGCACGACCTGAACATGCGTCTCTCCCCGCAAAGTGCGGATTCCGTTGCCAGCGCGGTGATc                                                                                     >  1:585069/1‑150 (MQ=255)                                                              cacCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGTTAAAGTGG‑TTGGGTATACCGACAGC‑ACGGGGGGGGACGAACTGAAAACGGGGCTCTCCCAGCAAAGTGGGGAGTCCGGTGCCCGCGCGGTGAGCAccc                                                                                 >  2:1178621/1‑149 (MQ=255)                                                              cacCCTGACCGGCGTGGCAATGGTAATGAAAGAGTATCCGAAAACGGC‑GGTTAACGCGA‑TTGGGTTTATCGAAAAC‑ACGGGGGGTGCAGACCTGGAAATTGGACTATCCCCGCAAAGAGCGGATTCCGCTGCCCGCGCGTTGATCAcac                                                                                 >  2:947488/1‑147 (MQ=255)                                                                 ccTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGTTAACGTGA‑TTGGTTATACCGACAGC‑ACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGCGTTGATCACCCAGGGc                                                                            <  2:1617157/151‑1 (MQ=255)                                                                  cTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGTTAACGTGA‑TTGGTTATACCGACAGC‑ACGGGGGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGCGTTGATCACCCAGGGCg                                                                           >  1:1425588/1‑151 (MQ=255)                                                                    gACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGTTAAAGTGA‑TTGGGTATAACGACAGC‑CCGGGGGGTCACGACCCGAACATGCGTCTCTCCCCGCAACGTGCGGATTCCGTTGCCCGCGCGGTGCTCAACCCCGGCGTg                                                                         >  2:1548105/1‑151 (MQ=255)                                                                    gACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGTTAAAGTGA‑TTGGGTATAACCACAGG‑ACGGGGGGTCACCACCCCAACATGGGGCTCTCCCCGCAAAGTGCGGGTTCCGTTTCCCGCGCGGTGGTCACCCAggg                                                                             >  2:1382400/1‑147 (MQ=255)                                                                       cGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGTTAACGTGA‑TTGGTTATACCGACAGC‑ACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGCGTTGATCACCCAGGGCGTGGAc                                                                      <  1:1389401/151‑1 (MQ=255)                                                                       cGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGTTAACGTGA‑TTGGTTATACCGACAGC‑ACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGCGTTGATCACCCAGGGCGTGGAc                                                                      <  2:566654/151‑1 (MQ=255)                                                                        ggCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGTTAACGTGA‑TTGGTTATACCGACAGC‑ACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGCGTTGATCACCCAGGGCGTGGACg                                                                     <  1:1068839/151‑1 (MQ=255)                                                                         gCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGTTAACGTGG‑TTGGGTATACCGACAAC‑ACGGGGGGGCACGAACTGAAAAAGCGTCTCTCCCACCCACCCGCGGAGTCCGTTGGCCGCGCGTTGACCACCCAGGGccggg                                                                        >  2:1021860/1‑143 (MQ=255)                                                                            tGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGTTAACTTGA‑TTGGTTATACCGCCAGC‑ACGGGGGGCCACCACCCGAACATTCGTCTCTCCCAGCAACGCGCGGGTTCCCTTTCCAACGCGTTGATCACCCCGGGCGTGGACGcc                                                                   >  2:1309219/1‑149 (MQ=255)                                                                               cAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGTTAACGTGA‑TTGGTTATACCGACAGC‑ACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGCGTTGATCACCCAGGGCGTGGACGCCAGCCg                                                              <  1:77416/151‑1 (MQ=255)                                                                                aaTGGTACTGAAAGAGTATCCGAAAACTGC‑GGCTAACGGGT‑TTGGCGATACCCGCAAG‑CCAGGTGGCCACGACCTGGACATGCGTCTCTACCATCAACGCGCTGGTTCCCGTGGCAGCCCGGTGGACAACCAGGGGGTGGCCGc                                                                    >  2:1410768/1‑144 (MQ=255)                                                                                aaTGGTACTGAAAGAGTATCCGAAAACGGC‑GGGTAACGGGA‑TTGGGTATAACGACAAC‑ACGGGTGGGGACGAACTGAACATGCGGCTCTCCCAGCAACGTGCGGGTTCCGTTGCCAGCGCGGTGGTCACCCAGGGCGGGGGCGCCAGcc                                                               >  1:781856/1‑149 (MQ=255)                                                                                 aTGGTACTGAAAGAGTATCCGAAAACGGC‑GGTTAACGTGA‑TTGTTTATACCGAAAGC‑ACGGGGGGTCACGACCCGAACAAGCGTCTCTCCCCGCCACGGGCGGAGTCCGGTGCCAGCGCGTTGATCACCCAGGGCGGGGAcgcccgcccc                                                             >  2:1085582/1‑148 (MQ=255)                                                                                  tGGTAATGAAAGAGTATCCGAAAACGGC‑GGTTAAAGGGA‑TTTGTTTTATAGACAGA‑AAGACGGGGCACGAACCGAAACAGCGTTTCTCCCATCCACATGCCGGGTTCATTCCCAGCGCGTGGATCACCCAGAGCGTGGATCCCcgccgc                                                             >  2:1185214/1‑149 (MQ=255)                                                                                        tGAAAGAGTATCCGAAAACGGC‑GGTTAACGTGA‑TTGGTTATACCGACAGC‑ACGGGTGGTCACGCCCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGCGTTGATCCCCCAGGGCGTGGACGGCAGCCGCATCCGt                                                       <  2:1600592/149‑1 (MQ=255)                                                                                        tGAAAGAGTATCCGAAAACGGC‑GGGTAAAGTTA‑TTGGGTATACCGACAAG‑CCGGGGGGGCAAGACCTTACCAAACGTCTTTCCCAGCAAAGTGCGGGATTCGGTGACAGCGAGTTGATACCCCCGCGCGCGGACCCCAGCCGCCACCGTa                                                      >  2:1426347/1‑150 (MQ=255)                                                                                         gAAAGAGTATCCGAAAACGGC‑GGTTAAAGTGA‑TTGGTTATACCGACAGC‑ACGGGGGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGCGTTGATCACCCAGGGCGGGGACGCCAGCCGCATCCGTAc                                                     >  1:1250360/1‑150 (MQ=255)                                                                                              agTATCCGAAAACGGC‑GGTTAACGTGA‑TTGGTTATACCGACAGC‑ACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGCGTTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTACTCAggg                                               <  1:1479531/151‑1 (MQ=255)                                                                                                   ccGAAAACGGC‑GGTTAACGTGA‑TTGGTTATACCGACAGC‑ACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGCGTTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTACTCAGGGCCtt                                           <  2:771753/150‑1 (MQ=255)                                                                                                          cggc‑ggTTAACGTGA‑TTGGTTATACCGACAGC‑CCGGGTGGTCACGACCTGAACATACCTCCCTCCCTGCAACGTGCGGATTCCCTTTCCAGCGCGCTGATCACCCAGGGCGTGGACGCCCGCCGCATCCGTACTCAGGGCCTTGGCCCTGc                                   >  2:1366121/1‑151 (MQ=255)                                                                                                             c‑ggTTAACGTGA‑TTGGTTATACCGACAGC‑ACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGCGTTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTACTCAGGGCCTTGGCCCGGCTaa                                <  2:1076262/151‑1 (MQ=255)                                                                                                              cgcctAACGCGA‑TCGGCTATACCGCCAGC‑ACGGGTAGTCACGACCTGAACATGCGTCTCTCCCGCCAACGTGCGGATTCCGTTGCCAGCGCGTTGATCACCCACGGCGTGGACGCCAGCCGCATCCGTATGCAGGGCCTTGGCCCGGCTa                                 >  2:144116/5‑150 (MQ=255)                                                                                                               ggTTAACGTGA‑TTGGTTATACCGACAGC‑ACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGCGTTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTACTCAGGGCCTTGGCCCGGCTAAc                               >  2:1367695/1‑151 (MQ=255)                                                                                                                 ttAACGTGA‑TTGGTTATACCGACAGC‑ACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCAGTTGCCAGCGCGTTGATCACCCAGGGCGCGGACGCCAGCCGCATCCGTACTCAGGGCCTTGGCCCGGCTAAc                               >  2:1530191/1‑149 (MQ=255)                                                                                                                   aaCGTGA‑TTGGTTATACCGACAGC‑ACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGCGTTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTACTCAGGGCCTTGGCCCGGCTAACCCaa                           <  1:1379959/151‑1 (MQ=255)                                                                                                                      gTGA‑TTGGTTATACCGACAGC‑ACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGCGTTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTACTCAGGGCCTTGGCCCGGCTAACCCAATCg                        <  1:1412244/151‑1 (MQ=255)                                                                                                                       tGA‑TTGGTTATACCGACAGC‑ACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGCGTTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTACTCAGGGCCTTGGCCCGGCTAACCCAATCGc                       >  2:306499/1‑151 (MQ=255)                                                                                                                        gA‑TTGGTTATACCGACAGC‑ACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGCGTTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTACTCAGGGCCTTGGCCCGGCTAACCCAATCGcc                      >  2:858920/1‑151 (MQ=255)                                                                                                                        gA‑TTGGTTATACCGACAGC‑ACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGCGTTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTACTCAGGGCCTTGGCCCGGCTAACCCAATCGcc                      <  1:1355778/151‑1 (MQ=255)                                                                                                                         a‑TTGGTTATACCGGCAGC‑ACGGGGGGTCACGACCGGAGCATGCGGCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGCGGTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTACTCAGGGCCGTGGCCCGGCTAGCCCAATCGCCa                     <  2:958654/151‑1 (MQ=255)                                                                                                                         a‑TTGGTTATACCGACAGC‑ACGGGTGGTCACGACCTGAACATGCGTCTCGCCCAGCAACGTGCGGATTCCGTTGCCAGCGCGTTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTACTCAGGGCCTTGGCCCGGCTAACCCAATCGc                       <  2:1337889/149‑1 (MQ=255)                                                                                                                            tGGTTATACCGACAGC‑ACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGGGCGGATTCCGTTGCCAGCGCGTTGATCACCCAGGGCGTGCACGCCAGCCGCATCCGTACTCAGGGCCTTGGGCCGGCTAACCCAATCGCCAg                    <  2:877/150‑1 (MQ=255)                                                                                                                             ggTTATACCGACAGC‑ACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGCGTTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTACTCAGGGCCTTGGCCCGGCTAACCCAATCGCCAGc                   >  1:139443/1‑150 (MQ=255)                                                                                                                             ggTTATACCGACAGC‑ACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGCGTTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTACTCAGGGCCTTGGCCCGGCTAACCCAATCGCCAGCa                  <  2:348977/151‑1 (MQ=255)                                                                                                                             ggTTATACCGACAGC‑ACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGCGTTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTACTCAGGGCCTTGGCCCGGCTAACCCAATCGCCAGCa                  <  2:6236/151‑1 (MQ=255)                                                                                                                                     cGACAGC‑ACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGCGTTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTACTCAGGGCCTTGGCCCGGCTAACCCAATCGCCAGCAACAGCAcc          >  2:1476134/1‑151 (MQ=255)                                                                                                                                     cGACAGC‑ACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGCGTTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTACTCAGGGCCTTGGCCCGGCTAACCCAATCGCCAGCAACAGCAcc          >  2:98743/1‑151 (MQ=255)                                                                                                                                      gACAGC‑ACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGCGTTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTACTCAGGGCCTTGGCCCGGCTAACCCAATCGCCAGCAACAGCACCg         <  2:869344/151‑1 (MQ=255)                                                                                                                                      gACAGC‑ACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGCGTTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTACTCAGGGCCTTGGCCCGGCTAACCCAATCGCCAGCAACAGCACCg         >  1:1308509/1‑151 (MQ=255)                                                                                                                                        cAGC‑ACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGCGTTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTACTCAGGGCCTTGGCCCGGCTAACCCAATCGCCAGCAACAGCACCGc        >  1:821123/1‑150 (MQ=255)                                                                                                                                         aGC‑ACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGCGTTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTACTCAGGGCCTTGGCCCGGCTAACCCAATCGCCAGCAACAGCACCGc        >  1:499433/1‑149 (MQ=255)                                                                                                                                         aGC‑ACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGCGTTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTACTCAGGGCCTTGGCCCGGCTAACCCAATCGCCAGCAACAGCACCGCAg      <  2:86113/151‑1 (MQ=255)                                                                                                                                          gC‑ACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGCGTTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTACTCAGGGCCTTGGCCCGGCTAACCCAATCGCCAGCAACAGCACCGc        >  1:851632/1‑148 (MQ=255)                                                                                                                                          gC‑ACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGCGTTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTACTCAGGGCCTTGGCCCGGCTAACCCAATCGCCAGCAACAGCACCGc        >  1:927380/1‑148 (MQ=255)                                                                                                                                           c‑ACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGCGTTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTACTCAGGGCCTTGGCCCGGCTAACCCAATCGCCAGCAACAGCACCGCAGaa    <  1:208678/151‑1 (MQ=255)                                                                                                                                           c‑ACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGCGTTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTACTCAGGGCCTTGGCCCGGCTAACCCAATCGCCAGCAACAGCACCGCAGa     >  2:1024543/1‑150 (MQ=255)                                                                                                                                              tgggTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGCGTTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTACTCAGGGCCTTGGCCCGGCTAACCCAATCGCCAGCAACAGCACCGCAGAAgg  >  2:1294184/2‑151 (MQ=255)                                                                                                                                                  |                                                                                                                                                   ATATGCCGAACAATGTGACCTTCGACAGCAGCAGCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGTTAACGTGA‑TTGGTTATACCGACAGC‑ACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGCGTTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTACTCAGGGCCTTGGCCCGGCTAACCCAATCGCCAGCAACAGCACCGCAGAAGG  >  NC_000913/3716872‑3717160

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