Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	523,335	0	T	G	38.7%	96.0 / 11.4	64	G25G (GGT→GGG)	rhsD	Rhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor
Reads supporting (aligned to +/- strand):  ref base T (15/23);  new base G (24/0);  total (41/23)
Fisher's exact test for biased strand distribution p-value = 3.13e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGGTCTTGTGTGTCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGG‑TTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGC‑GTGGCGTGCTCGGTGT‑GTCCGGGCGGGATGACTT‑CGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCT  >  NC_000913/523186‑523483                                                                                                                                                       |                                                                                                                                                        cTGGTCTTGTGTGTCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTTCAGGG‑GTTGGCAGgtg                                                                                                                                                        >  1:186677/1‑151 (MQ=255)   gggcTTGTGTGTCATTTTATTTTTTTTTGTTGCTAACACAGAGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGTCGGTCCCATTGTCCAGGG‑TTCGGCAGGTGTa                                                                                                                                                      <  1:799659/148‑1 (MQ=255)        tgtgtgTCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCCCTGTCCAGGG‑GTTGGCAGGTGTaa                                                                                                                                                     >  1:235479/1‑147 (MQ=255)         gtgtgtCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGG‑GTCGGCAGGGGTAAGAAtt                                                                                                                                                >  1:364427/1‑151 (MQ=255)           gtgtCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGG‑TTCGGCAGGTGTAAGAATTg                                                                                                                                               <  1:299991/150‑1 (MQ=255)             gtCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGG‑TTCGGCAGGTGTAAGAATTGgcg                                                                                                                                            <  2:21183/151‑1 (MQ=255)             gtCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGG‑TTCGGCAGGTGTAAGAATTGgc                                                                                                                                             <  1:895923/150‑1 (MQ=255)             gtCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGG‑GTCGGGAGGGGGAAGAATTgggg                                                                                                                                            >  1:1148993/1‑149 (MQ=255)               cATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCCCTCTCCAGGGCTTCGGC‑GGGGTAGGAAATgg                                                                                                                                              >  2:151760/1‑147 (MQ=255)                 ttttattttTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGGCCCATTGTCCAGGG‑GTCGGCGGGGGTTAGAATTGGCGCGcc                                                                                                                                        >  1:16598/1‑151 (MQ=255)                  tttattttTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGG‑GTCGGCAGGGGGAAGAATTGGCGCGccc                                                                                                                                       >  1:1044974/1‑151 (MQ=255)                    tattttTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCCTTGTCCAGGG‑TTCGGCAGGTGGAAGAATTGGGGCGccccc                                                                                                                                     >  1:820818/1‑149 (MQ=255)                     attttTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGG‑GTCGGCAGGGGTAAGAATTGGCGCGCCCAcc                                                                                                                                    >  1:145253/1‑151 (MQ=255)                       ttttttttGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGG‑GTCGGGAGGGGTAAGAATTGGCGCGCCCAcc                                                                                                                                    >  1:127847/1‑149 (MQ=255)                            tttGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCCTTGTTCAGGG‑GTCGGCAGGGGTAAGAATTGGCGCGCCCCCCcgc‑gggg                                                                                                                            >  2:614733/1‑151 (MQ=255)                               gttgCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTTTCCCGGG‑GTCGGCAGGGGGAAGAAATGGCGCGCCCACCcgc‑gtgg                                                                                                                            >  2:276485/1‑148 (MQ=255)                                 tgCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCCTTGTCCCGGG‑TTCGGCGGGGGGAAGAATTGGGGGGcccccc                                                                                                                                    >  2:393235/1‑139 (MQ=255)                                  gCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCCTTGTCCAGGG‑GTCGGCAGGGGGAAGGATTGGCGCGCCCACCcgc‑gtggcgt                                                                                                                         >  1:24680/1‑148 (MQ=255)                                  gCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGACCCATTGGCCCGGG‑GTCGGCAGGGGGAAGAATTGGGGCGCCCAcccccgcgggg                                                                                                                           >  2:130414/1‑138 (MQ=255)                                   cTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGG‑GTCGGCAGGGGTAAGGATTTGGGCGGCCACCggc‑gtggcggg                                                                                                                        >  2:482437/1‑146 (MQ=255)                                    tAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCCGGG‑GTCGGGCGGGGTAAGAAATGGCGCGCCCCCCggc‑ggggcg                                                                                                                          >  2:1108352/1‑145 (MQ=255)                                     aaCACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGG‑GTCGGCAGGTGTAAGAATTGGCGGGCCCACCGGC‑GGGGGGGGGTCgg                                                                                                                   >  1:769582/1‑151 (MQ=255)                                      acacAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGG‑TTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGC‑GTGGCGTGCTCGgt                                                                                                                  <  1:193329/151‑1 (MQ=255)                                       cacaGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCCTTGTCCAGGG‑GTCGGCAGGGGTAAGAATTGGGGCGCCCACCGgc‑gtggcgtggccgg                                                                                                                   >  1:839390/1‑144 (MQ=255)                                        acaGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCCTTCTCCAGGG‑GTCGGGCGGGGTAAGAAATGGCGGGCCCACCggc‑ggcggggg                                                                                                                        >  2:847332/1‑136 (MQ=255)                                         caGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGG‑GTCGGCAGGGGTAAGGATTTGGGCGCCCCCCCGC‑GTGGGGTGGTCCgggt‑g                                                                                                              >  1:396192/1‑151 (MQ=255)                                           gggATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCCTCGTCCCGGG‑CTCGGCATGAGGATGAAATCGCGCGCCCCACGGC‑CGGGTGTGcgccttgtg                                                                                                               >  2:615979/1‑141 (MQ=255)                                                atGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGG‑TTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGC‑GTGGCGTGCTCGGTGT‑GTCcgggc                                                                                                       <  2:811616/151‑1 (MQ=255)                                                atGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGG‑TTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGC‑GTGGCGTGCTCGGTGT‑GTCcggg                                                                                                        <  1:682183/150‑1 (MQ=255)                                                  gAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGG‑TTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGC‑GTGGCGTGCTCGGTGT‑GTCcgggcgg                                                                                                     <  1:163520/151‑1 (MQ=255)                                                   aaCAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGGCCCGGG‑GTCGGCAGGGGGAAGAATTTGGGCGCCCCCCGGC‑GGGGGGTGCTCGGGGT‑GTCccgggggg                                                                                                    >  1:603905/1‑151 (MQ=255)                                                         aaCTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCCGGG‑GTCGGCAGGGGTAAGAAATGGCGCGCCCACCGGC‑GTGGGGTGGTCGGGGT‑GTCcggggggg                                                                                                    >  1:175791/1‑145 (MQ=255)                                                                   gCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGACCCGTTGTCCCCGG‑TTAGGCAGGGGGAAGGATTTGCGGGCCCCCCCGC‑GTGGCGGTCCCGGGGT‑GGACGGGGGGGGGGGATT‑CGGGccgac                                                                                   >  2:962754/1‑147 (MQ=255)                                                                     aCTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTTTTCAGGG‑GTCGGCCGGGGGAAGAATTGGCGCGCCCACCTGC‑GGGGCGGGGTTCGTGT‑GGCCGGGGGGGAGGACCT‑CGGGCAAcc                                                                                   >  1:713431/1‑149 (MQ=255)                                                                                gCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGG‑TTCGGCAGGAGTAAGAATTTGGGCGCCCCCCGGC‑GCGGGGTGGTCGGTGT‑GTCCCGGGGGGGGAACTTCCCGGAAACCCGGaaaaac                                                                          >  1:1069144/1‑146 (MQ=255)                                                                                    aaaaCCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGG‑TTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGC‑GTGGCGTGCTCGGTGT‑GTCCGGGCGGGATGACTT‑CGGGCAACCCGGTAAATCCGctgctg                                                                  <  2:1069887/151‑1 (MQ=255)                                                                                     aaaCCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGG‑TTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGC‑GTGGCGTGCTCGGTGT‑GTCCGGGCGGGATGACTT‑CGGGCAACCCGGTAAATCCGCTGCTgg                                                                 <  2:147138/151‑1 (MQ=255)                                                                                     aaaCCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGG‑TTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGC‑GTGGCGTGCTCGGTGT‑GTCCGGGCGGGATGACTT‑CGGGCAACCCGGTAAATCCGCTGCTgg                                                                 <  1:173985/151‑1 (MQ=255)                                                                                     aaaCCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGG‑GTCGGCAGGGGTAAGAATTTGGGCGCCCACCCGC‑GTGGCGGGCCCGGGGTGGTCCGGG‑GGGGGGACTT‑CGGGCCACCCggggaaaaccccgctgg                                                                 >  2:1014201/1‑136 (MQ=255)                                                                                      aaCCAGCGGCGCGTCAGGGAGATATGACCCAGTATGGCGGTCCCAATGTCCCGGG‑GTCGGGCGGGGGAAGAAATGGCGGGGCCACCGGC‑GCGGGGTGCCCGGGGTGGTCGGGGGGGGTGGACGC‑TGGGGAACCCGGAAAAAccccgcccgg                                                                 >  2:577366/1‑143 (MQ=255)                                                                                        ccAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGG‑TTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGC‑GTGGCGTGCTCGGTGT‑GTCCGGGCGGGATGACTT‑CGGGCAACCCGGTAAATCCGCTGCTggggg                                                              <  2:109413/151‑1 (MQ=255)                                                                                          aGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGGCCCCTTGTCCAGGG‑GTCGGCGGGGGGAAGAAATTGGGCGCCCCCCCGCGGGGGGGGGCTCGGGGG‑GTCCGGGGGGGGGGACTT‑CGGGCACCCCGGGAAATcccccccggggg                                                               >  2:690425/1‑139 (MQ=255)                                                                                             ggcgCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGG‑TTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGC‑GTGGCGTGCTCGGTGT‑GTCCGGGCGGGATGACTT‑CGGGCAACCCGGTAAATCCGCTGCTGGGGGCGaa                                                          <  2:235471/150‑1 (MQ=255)                                                                                             ggcgCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGG‑TTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGC‑GTGGCGTGCTCGGTGT‑GTCCGGGCGGGATGACTT‑CGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAg                                                         <  2:1066543/151‑1 (MQ=255)                                                                                             ggcgCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGG‑TTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGC‑GTGGCGTGCTCGGTGT‑GTCCGGGCGGGATGACTT‑CGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAg                                                         <  2:295069/151‑1 (MQ=255)                                                                                                  gTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGG‑TTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGC‑GTGGCGTGCTCGGTGT‑GTCCGGGCGGGATGACTT‑CGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGtgct                                                    <  1:160028/151‑1 (MQ=255)                                                                                                      gggAGATATGACTCAGTATGGCGGTCCCATTGGCCAGGG‑TTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGC‑GTGGCGTGCTCGGTGT‑GTCCGGGCGGGATGACTT‑CGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGcc                                                 <  1:1195895/150‑1 (MQ=255)                                                                                                          gaTATGACTCAGTATGGCGGTCCCATTGTCCCAGG‑GTCGGCAGGTGTAAAAATTTGCGCGCCCCCCGGC‑GGGGGGGGGTCGGTGG‑GTCCGGGCGGGGGGAATT‑CCGGCAAACCCGTAAAACCCCCGCGGGGGGCGGAGGGGCTGccccggg                                            >  2:176166/1‑147 (MQ=255)                                                                                                          gaTATGACTCAGTATGGCGGTCCCATTGTCCAGGG‑TTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGC‑GTGGCGTGCTCGGTGT‑GTCCGGGCGGGATGACTT‑CGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGc                                             <  2:982104/150‑1 (MQ=255)                                                                                                                     gTATGGCGGTCCCATTGTCCAGGG‑TTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGC‑GTGGCGTGCTCGGTGT‑GTCCGGGCGGGATGACTT‑CGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCtt                                 <  2:712002/151‑1 (MQ=255)                                                                                                                       aTGGCGGTCCCATTGTCCAGGG‑TTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGC‑GTGGCGTGCTCGGTGT‑GTCCGGGCGGGATGACTG‑CGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTgc                               <  2:956542/151‑1 (MQ=255)                                                                                                                        tGGCGGTCCCATTGTCCAGGG‑TTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGC‑GTGGCGTGCTCGGGGG‑GTCCGGGCGGGGTGACTT‑CGGGCAACCCGGGAAAACCGCCGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGGCCtt                                 >  2:740627/1‑148 (MQ=255)                                                                                                                           cGGTCCCATTGTCCAGGG‑TTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGC‑GCGGCGTGCTCGGCGT‑GTCCGGGCGGGATGACTT‑CGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTg                           >  2:154955/1‑151 (MQ=255)                                                                                                                                ccATTGTCCAGGG‑TTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGC‑GTGGCGTGCTCGGTGT‑GTCCGGGCGGGATGACTT‑CGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTgcccgg                      >  2:123819/1‑151 (MQ=255)                                                                                                                                   ttGTCCAGGG‑TTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGC‑GTGGCGTGCTCGGTGT‑GTCCGGGCGGGATGACTT‑CGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTgcccggccc                   >  1:1152534/1‑151 (MQ=255)                                                                                                                                    tGTCCAGGG‑TTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGC‑GTGGCGTGCTCGGTGT‑GTCCGGGCGGGATGACTT‑CGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTgcccggcccg                  >  1:216536/1‑151 (MQ=255)                                                                                                                                       ccAGGG‑TTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGC‑GTGGCGTGCTCGGTGT‑GTCCGGGCGGGATGACTT‑CGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCt                <  2:1029410/150‑1 (MQ=255)                                                                                                                                              ttCGGCAGGTGTAAGAATTGGCGCGCCCACCGGC‑GTGGCGTGCTCGGTGT‑GTCCGGGCGGGGTGACTT‑CGGGCAACCCGGTAAATCCGCTGCTGGGGGCGGAGGTGCTGCCCGGCGAGGCGGACCTTGCGCTGCCCGGCCCGCTGCCGtt          >  2:123652/1‑150 (MQ=255)                                                                                                                                              ttCGGCAGGTGTAAGAATTGGCGCGCCCACCGGC‑GTGGCGTGCTCGGTGT‑GTCCGGGCGGGATGACTT‑CGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGt           >  2:70627/1‑149 (MQ=255)                                                                                                                                               tCGGCAGGTGTAAGAATTGGCGCGCCCACCGGC‑GTGGCGTGCTCGGTGT‑GTCCGGGCGGGATGACTT‑CGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTc         >  1:721887/1‑150 (MQ=255)                                                                                                                                               tCGGCAGGTGTAAGAATTGGCGCGCCCACCGGC‑GTGGCGTGCTCGGTGT‑GTCCGGGCGGGATGACTT‑CGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTc         <  1:1108405/150‑1 (MQ=255)                                                                                                                                                cGGCAGGTGTAAGAATTGGCGCGCCCACCGGC‑GTGGCGTGCTCGGTGT‑GTCCGGGCGGGATGACTT‑CGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCAt       >  2:835373/1‑151 (MQ=255)                                                                                                                                                  gCAGGTGTAAGAATTGGCGCGCCCACCGGC‑GTGGCGTGCTCGGTGT‑GTCCGGGCGGGATGACTT‑CGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATtc     >  1:1165436/1‑151 (MQ=255)                                                                                                                                                   cAGGTGTAAGAATTGGCGCGCCCACCGGC‑GTGGCGTGCTCGGTGT‑GTCCGGGCGGGATGACTT‑CGGGCAACCCGGTAAATCCGCTGCTGGGGGGGGAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCAtt      >  2:702375/1‑149 (MQ=255)                                                                                                                                                     ggTGTAAGAATTGGCGCGCCCACCGGC‑GTGGCGTGCTCGGTGT‑GTCCGGGCGGGATGACTT‑CGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATtctct  <  2:16597/151‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                                        CTGGTCTTGTGTGTCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGG‑TTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGC‑GTGGCGTGCTCGGTGT‑GTCCGGGCGGGATGACTT‑CGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCT  >  NC_000913/523186‑523483

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