Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	523,326	0	T	G	44.1%	81.3 / 10.6	60	G22G (GGT→GGG)	rhsD	Rhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor
Reads supporting (aligned to +/- strand):  ref base T (10/23);  new base G (26/0);  total (37/23)
Fisher's exact test for biased strand distribution p-value = 8.86e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACATTAGCTGGTCTTGTGTGTCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGG‑TTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGC‑GTGGCGTGCTCGGTGT‑GTCCGGGCGGGATGACTT‑CGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTT  >  NC_000913/523179‑523475                                                                                                                                                     |                                                                                                                                                         aCATTAGCTGGTCTTGTGTGTCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGG‑TTc                                                                                                                                                        <  1:1231886/150‑1 (MQ=255) aCATTAGCTGGTCTTGTGTGTCATTTTA‑TTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGG‑TTCg                                                                                                                                                       <  2:619713/150‑1 (MQ=255)        cTGGTCTTGTGTGTCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTTCAGGG‑GTTGGCAGgtg                                                                                                                                                >  1:186677/1‑151 (MQ=255)          gggcTTGTGTGTCATTTTATTTTTTTTTGTTGCTAACACAGAGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGTCGGTCCCATTGTCCAGGG‑TTCGGCAGGTGTa                                                                                                                                              <  1:799659/148‑1 (MQ=255)               tgtgtgTCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCCTTGTCCCGGG‑GTCGGCAgggggtagaa                                                                                                                                          >  2:314883/1‑142 (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)                                                                                                                                                     |                                                                                                                                                         ACATTAGCTGGTCTTGTGTGTCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGG‑TTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGC‑GTGGCGTGCTCGGTGT‑GTCCGGGCGGGATGACTT‑CGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTT  >  NC_000913/523179‑523475

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