Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	523,335	0	T	G	55.9%	37.3 / 67.9	60	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 (9/17);  new base G (33/0);  total (43/17)
Fisher's exact test for biased strand distribution p-value = 1.13e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.83e-07

GGTCTTGTGTGTCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGC‑TCGGTGTGT‑CCGGGCGGGATGA‑CTTCGGGCAACCCGGT‑AAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCT  >  NC_000913/523188‑523483                                                                                                                                                    |                                                                                                                                                         ggTCTTGTGTGTCATTTTATTTTTTTTTGTTGCTAACACAGGGATTTGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAAATATGACTCAGTATGGCGGTCCCCTTGTCACGGGTTCGGGAGGGGTa                                                                                                                                                       >  2:268472/1‑151 (MQ=255) ggTCTTGTGTGTCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCCTTGTCCAGGGTTCGGCAGGGGTa                                                                                                                                                       >  1:611757/1‑151 (MQ=255)     ttGTGTGTCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGaa                                                                                                                                                   <  1:542185/151‑1 (MQ=255)     ttGTGTGTCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGaa                                                                                                                                                   <  1:430928/151‑1 (MQ=255)     ttGTGTGTCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGaa                                                                                                                                                   <  2:70228/151‑1 (MQ=255)     ttGTGTGTCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGATCGGCAGGTGTAAGaa                                                                                                                                                   <  2:459918/151‑1 (MQ=255)         gtgtCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTTCCATTTTCCCGGGTTTGGGCGGTGTAAGAAATgg                                                                                                                                               >  2:33453/1‑151 (MQ=255)         gtgtCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCCTTGTTCAGGGTTCGGCAGGTGTAAGAATTgg                                                                                                                                               >  2:1105533/1‑151 (MQ=255)         gtgtCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGCAGGGGGAAGAATTgg                                                                                                                                               >  1:472364/1‑151 (MQ=255)            tCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGgcg                                                                                                                                             <  2:48353/150‑1 (MQ=255)                          tttGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCggcgtgg                                                                                                                              <  1:719008/151‑1 (MQ=255)                            tgttgCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCggcgtggcg                                                                                                                            <  1:1094107/151‑1 (MQ=255)                               tgCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTTGGCAGGGGGAAGAATTGGCGCGCCCCCCggcgtggcggg                                                                                                                          >  1:722121/1‑148 (MQ=255)                               tgCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATCGTCCCAGGTTCTGCACGTGTAATAATTTGCGTGCCCCCCggccgcgcctg                                                                                                                          >  2:636580/1‑142 (MQ=255)                                gCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTTTCCAGGGGTCGGCAGGGGGAAGAATTGGCGCGCCCACCggcggggcggg                                                                                                                          >  1:404758/1‑147 (MQ=255)                                gCTAACAAAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCCATGTCCAGGGGTCGGCACGGGTAAATATTTGCCCCCCCCCCCgccttgcgtg                                                                                                                          >  1:696375/1‑149 (MQ=255)                                     cacaGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCCTTTTCCCGGGGTCGTCAGGGGTTAGAATTTGCGCGCCCACCCGCGTGGGGGGC‑GCGgtg                                                                                                                  >  2:634769/1‑151 (MQ=255)                                     cacaGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCCCTTTCCAGGGTTCGGGGGGGGTAAGAATTGGGGCCCCCACCcgggtgggggtttccgcg                                                                                                                   >  1:963980/1‑140 (MQ=255)                                         gggATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGCAGGGGTAAGAATTGGCGCGCCCACCGGCGGGGGGGGC‑TCGgggtgt                                                                                                               >  1:602893/1‑150 (MQ=255)                                           gATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGCCGGGGTAAGAAATTGGGCGCCCACCCGCGGGGCGGGG‑TCgggtggttcg                                                                                                            >  2:294804/1‑143 (MQ=255)                                           gATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGCAGGGGGAAGAATTGGCGCGCCCACCGGCGGGGCGGGG‑TCgggggggccg                                                                                                            >  1:400652/1‑143 (MQ=255)                                           gATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCACTGTCCACGGTTCTGCAGGGGGTAGAATTGGGGCGCCCACCCGCGGGGCGGGC‑GCTGTGTGT‑Gcg                                                                                                           >  2:363927/1‑151 (MQ=255)                                            atatGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGC‑TCGGTGTGT‑Ccg                                                                                                           <  1:1135643/150‑1 (MQ=255)                                                 aaCAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGGGTAAGAATTGGGGGGCCCACCGGCGGGGCGGGC‑TCCGGGTGT‑Cccgggggg                                                                                                     >  1:520172/1‑151 (MQ=255)                                                   cAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGGGTAAGAATTGGCGCGCCCACCGGCGGGGGGGGC‑TCCGGGGGT‑Ccggggggggg                                                                                                   >  1:753831/1‑149 (MQ=255)                                                       aaCTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGGGTAAGAAATTGGGCGCCCACCGGCGTGGCGGGC‑TCGGTGTGT‑CCGGGGGGGGTgg                                                                                                 >  2:270189/1‑148 (MQ=255)                                                        aCTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGC‑TCGGTGTGT‑CCGGGCGGGATGA‑Ctt                                                                                             <  1:268556/151‑1 (MQ=255)                                                         cTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTGCCATTTTCCAGGGGTCGGCAGGGGTAAAAAATGGCGGGCCCACCGGCGTGGCGTGCTTCGGGGGGT‑Cggggggggtgaact                                                                                               >  2:757446/1‑144 (MQ=255)                                                                ggCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGC‑TCGGTGTGT‑CCGGGCGGGATGA‑CTTCGGGCa                                                                                       <  1:320144/149‑1 (MQ=255)                                                                    cTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCCGGGTTCGGCAGGGGTAAGAATTGGCGCGCCCACCCGCGTGGCGGGC‑CCGGGGTGT‑CCGGGCGGGGTGA‑CTTCGGGCAAACCgg                                                                                 >  1:946132/1‑151 (MQ=255)                                                                    cTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGGAGGGGTAAGAATTGGCGCGCCCACCGGGGGGGCGGGG‑CCGGTGTGT‑CCGGGCGGGGGGA‑GTTCGGGGAAACCgg                                                                                 >  1:1055647/1‑151 (MQ=255)                                                                    cTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGC‑TCGgtg                                                                                                                  <  2:622748/120‑1 (MQ=255)                                                                    cTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGGAGGGGTAAGAAATGGCGCGCCCACCGGCGTGGCGGGC‑CCGGGGGGT‑CCCGGCCGGGTGA‑CTTTGGGGAAacccg                                                                                 >  1:1021368/1‑149 (MQ=255)                                                                    cTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGCAGGGGTAAGAATTTGGGCGCCCACCGGCGTGGCGGGC‑GCGGGGTGT‑CCGGGGGGGGTGG‑CTTCGGGCAAccccg                                                                                 >  1:217395/1‑149 (MQ=255)                                                                    cTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGCAGGGGTAAGAATTGGGGCGCCCCCCGGCGTGGCGGGC‑TCGGGGTTT‑GCGGGCGGGGTGA‑GTTCGGGCaaacagg                                                                                 >  2:974815/1‑146 (MQ=255)                                                                    cTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGCAGGGGTAAGAATTGGCGCGCCCACCGGCGGGGCGGGC‑TCGGGGTGT‑CCCGGGGGGGTGA‑CTTCGGGGAACCCgg                                                                                 >  1:195643/1‑151 (MQ=255)                                                                    cTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGG‑TCGGCAGGGGTAAGAATTGGCGCGCCCCCCGGCGGGGCGGGC‑CCGGGGGGT‑CCGGGGGGGGTGA‑CTTCGGGCAACCCggg                                                                                >  1:898273/1‑150 (MQ=255)                                                                    cTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCAATGTCCCAGGGTCGGCAGGGGTAAAAATTGGCGCGCCCACCGGCGTGGGCTGG‑CCGGGTTGT‑CCTGGCGGGGAGACCTTCGGGAAAccgg                                                                                  >  2:854398/1‑149 (MQ=255)                                                                    cTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCAATGTCCAGGGGTCGGGAGGGGTAAGAATTTGCGCGCCCACCGGCGTGGCGGGC‑TCGGTGTGG‑CCGGGGGGGGGGG‑CGTCCGGGAACCCgg                                                                                 >  2:216124/1‑151 (MQ=255)                                                                    cTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCAATGTCCAGGGGTCGGCAGGTGTAAGGATTGGCGCGCCCACCCGCGTGGCCTGC‑TCGGTGTGT‑GCCGGGGGGGTTC‑CATTGGGCAAccccg                                                                                 >  2:526544/1‑149 (MQ=255)                                                                      ttATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGGACCTTCTTTCAGGGGTGTGTAGGGGTGAGAATGTGGGCCCCCACCGGTGCTGCGTGC‑CCGGGGCGT‑CGCGGGGGGGTGA‑GGTTGGGGAACCCggga                                                                               >  2:455016/1‑149 (MQ=255)                                                                       tataTGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGC‑TCGGTGTGT‑CCGGGCGGGATGA‑CTTCGGGCAACCCGGT‑aa                                                                             <  2:800009/151‑1 (MQ=255)                                                                       tataTGAGCGGAAAACCAGCGGCGCGTCAGAGAGATAAGACTATGTATGGCATTCCCATGGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTTC‑TCGGTGTGT‑CCGGGCGGGATGA‑CTTCGGGCAACCCGGT‑aa                                                                             <  2:517269/151‑1 (MQ=255)                                                                          atGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGC‑TCGGTGTGT‑CCGGGCGGGATGA‑CTTCGGGCAACCCGGT‑AAATc                                                                          <  2:312813/151‑1 (MQ=255)                                                                                 gAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCCTTGTCCCGGGATCGGCAGGGGTAAGAAATTGCGCGCCCCCCGGTGTGGCGTGG‑TCGGGGAGG‑CCGGGCGGGGTGA‑CCTTGGGGAACACGGTGAAACCCgcgc                                                                     >  2:1144441/1‑148 (MQ=255)                                                                                 gAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCCATGTCCAGGGGTCCGCAGGGGTAAGAATTGGGGCGCCCACCGGCGTGGCGTGC‑TCGGTGTGT‑CCGGGCGGGGTGA‑CTTCGGGCAACCCGGT‑GAATCCgcttc                                                                    >  2:105951/1‑148 (MQ=255)                                                                                                       agaTATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGCAGGGGTAAGAATTGGCGCGCCCACCGGCGTGGCGGGC‑TCGGGGGGT‑CCGGGCGGGGTGA‑CTTCGGGGAAACCGGG‑AAAACCGCTTCTGGGGGGGGAGGGGGCGCCCgg                                              >  1:627824/1‑150 (MQ=255)                                                                                                        gaTATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGGGTAAGAATTGGCGCGCCCACCGGCGTGGCGGGC‑TCGGGGGGT‑CCGGGCGGGGTGA‑CTTCGGGGAACCCGGG‑AAATCCGCCGCTGGGGGGGGAGGGGCTGCCCGGCg                                            >  1:940126/1‑151 (MQ=255)                                                                                                         atatGACTCAGTATGGCGGTCCCCTTGTCCCGGGGTCGGCCGGTGTAAGAATTTGCGCGCCCACCCGCCGGGCCTGC‑GCCGGGGGG‑CCCGGCGGGGTGA‑CTTCGGGGAACCCCGT‑AAAACCCCTGCTGGGGGGGGAGGGGTGGGCCGGCgg                                           >  2:791824/1‑150 (MQ=255)                                                                                                         atatGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGGGTAAGAATTGGCGCGCCCACCGGCGGGGCGTGC‑TCGGGGTGT‑CCGGGCGGGGTGA‑CTTCGGGCAACCCGGT‑AAATCCCCTGCTGGGGGGGGAGGGGCTGCCCGGCgg                                           >  1:488359/1‑150 (MQ=255)                                                                                                              aCTCAGTATGGCGGTCCCATTTTCCAGGGTTCGGCAGGGGTAAGAATTGGCGCGCCCCCCGGCGTGGGGGGC‑TCGGTGTGTGCGGGGGGGGGGGA‑TTCCGGGCACCCCGG‑‑AAATTCGCTGCTGGGGGGGGAGGGGGTGCCCGGGGGGACgg                                      >  2:140720/1‑151 (MQ=255)                                                                                                              aCTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGC‑TCGGTGTGT‑CCGGGCGGGGTGA‑CTTCGGGCAACCCGGT‑AAATCCGCTGCTGGGGGCGAAGGGGCTGCCCGGCGGGACgg                                      >  1:946786/1‑151 (MQ=255)                                                                                                                tgaGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACGGGCGTGGCGTGC‑GCGGTGTGT‑CCGGGGGGGATGA‑CTTCGGGCAACCCGGT‑AAATCCGCTGCTGGGGGCGAGGGTGCTGCCCGGCGAGACGGAc                                    <  2:852622/149‑1 (MQ=255)                                                                                                                  aGTATGGCGGTCCCATTGTCCAGGGGTCGGCAGGGGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGC‑TCGGGGTGT‑CCGGGGGGGGGGA‑ATTCGGGGAACCCGGG‑AAATCCCCTGCTGGGGGGGGAGGGGGTGCCCgggggggcgggc                                    >  1:689092/1‑139 (MQ=255)                                                                                                                      tGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGC‑TCGGTGTGT‑CCGGGCGGGATGA‑CTTCGGGCAACCCGGT‑AAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTgcg                              <  1:33459/151‑1 (MQ=255)                                                                                                                           gTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCCCCGGCGTGGCGTGC‑TCGGTGTGT‑TCGGGCGGGATGA‑CCTCGGGCAACCCGGC‑AAATCCGCTGTTGGGGGCGAAGGCTCCGCTCGGCGAGCCCGACCTTGCGCTgc                          >  2:1142248/1‑150 (MQ=255)                                                                                                                             cccATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGC‑TCGGTGTGT‑CCGGGCGGGATGA‑CTTCGGGCAACCCGGT‑AAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTgcccg                       >  1:174827/1‑151 (MQ=255)                                                                                                                             cccATTGTCCAGGGTTCGGCAGGCGTAAGAATTGGCGCCCCCACCGGCGCGGCGCGC‑TCGGTGTGT‑CCGGGCGGGAGGA‑CCTCGGGCAACCCGAT‑AAATCCGCGGCTGGGGGCGAAGGTGATTCCCGGCGAGCCGGACCTTGCGCTccccg                       >  2:490293/1‑151 (MQ=255)                                                                                                                                   gTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGC‑TCGGTGTGT‑CCGGGCGGGATGA‑CTTCGGGCAACCCGGT‑AAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTgcccggcccg                  >  1:952370/1‑150 (MQ=255)                                                                                                                                           ttCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGC‑TCGGTGTGT‑CCGGGCGGGATGA‑CTTCGGGCAACCCGGT‑AAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTc         <  1:757783/151‑1 (MQ=255)                                                                                                                                                   gtgtAAGAATTGGCGCGCCCACCGGCGTGGCGTGC‑TCGGTGTGT‑CCGGGCGGGATGA‑CTTCGGGCAACCCGGT‑AAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATtctct  >  1:34279/1‑150 (MQ=255)                                                                                                                                                    |                                                                                                                                                         GGTCTTGTGTGTCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGC‑TCGGTGTGT‑CCGGGCGGGATGA‑CTTCGGGCAACCCGGT‑AAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCT  >  NC_000913/523188‑523483

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