Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913523,3350TG38.1% 156.0 / 15.0 84G25G (GGT→GGG) rhsDRhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor
Reads supporting (aligned to +/- strand):  ref base T (9/43);  new base G (32/0);  total (41/43)
Fisher's exact test for biased strand distribution p-value = 2.24e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCTTGTGTGTCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGC‑GTGCTCGGTGT‑GT‑CCGGGCGGGAT‑GA‑CTTC‑GGGC‑AACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCT  >  NC_000913/523190‑523483
                                                                                                                                                 |                                                                                                                                                           
tCTTGTGTGTCATTTTA‑TTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTTGGCAGgggt                                                                                                                                                           >  2:870012/1‑147 (MQ=255)
      gtgtCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTg                                                                                                                                                   <  2:578248/150‑1 (MQ=255)
       tgtCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTg                                                                                                                                                   <  2:1264562/149‑1 (MQ=255)
         tCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGGCCCCTTGTCCCGGGGTCCGCGGGGGGAAGAATTGgcg                                                                                                                                                >  2:403975/1‑150 (MQ=255)
         tCATTTTATTTTTTTTTGTTGCGAACACAGGGATGTGAACACTAACTAAAAGGACTCTTTATATGAGCGATAAACCAGGGGCGCGTCAGGGTGATATGACTCAGTATGGCGGGCCCACTGTAGCGGGTTCGGCTGGGGTAGAAATTgtg                                                                                                                                                 >  1:608961/1‑147 (MQ=255)
           attttattttTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGCAGGTGTAAGAATTGgcgcg                                                                                                                                              >  2:805708/1‑150 (MQ=255)
                   tttttttGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGGGTAAGAATTGGCGCGCCCAcccg                                                                                                                                      >  2:139052/1‑148 (MQ=255)
                    ttttttGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCCTTGTCCAGGGGTCGGCAGGTGTAAGAATTTGGGCGCCCACCggc                                                                                                                                     >  2:1387727/1‑150 (MQ=255)
                     tttttGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCCTTGTGCAGGGGTCGGGAGGGGTAAGAAATTGCGGGGCCACCcgcgg                                                                                                                                   >  1:712721/1‑150 (MQ=255)
                       tttGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCggcgtg                                                                                                                                  <  2:914595/150‑1 (MQ=255)
                       tttGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGCACCCTCGTCCCGGGGTGGGGAGGGGTTAAAAAGGGCGCGGCCACCggcggcg                                                                                                                                 >  1:1076852/1‑148 (MQ=255)
                        ttgttgCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCggcgtgg                                                                                                                                 <  2:320969/150‑1 (MQ=255)
                             gCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCCTTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGgcgtggc‑gtg                                                                                                                            >  2:166613/1‑149 (MQ=255)
                             gCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCCGGGGTCGGCAGGTGTAAGAAATGGGGCGCCCACCcgcgtgg                                                                                                                                 >  1:299982/1‑145 (MQ=255)
                              cTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCCTTGTCCAGGGGTCGGCAGGGGGAAGAATTTGGGCGCCCACCggcgtggc‑g                                                                                                                              >  2:950672/1‑146 (MQ=255)
                                aaCACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGGCCCCTTGTGCCGGGGTCGGGCGGGGTTTGAAATGGGGCGCCCACCCGGGTCGC‑GGGcc                                                                                                                          >  2:889672/1‑147 (MQ=255)
                                 acacAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGC‑GTGCTCGgt                                                                                                                      <  1:1094415/151‑1 (MQ=255)
                                 acacAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGCAGGGGTAAGAATTGGCGCGCCCACCGGCGTGGC‑GGGCTCgg                                                                                                                       >  2:1232544/1‑150 (MQ=255)
                                             gAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTTCCGGGGTCGGCAGGGGGAAGAATTGGCGCGCCCCGCGGCGGGGG‑GGGCTCCGGGT‑GT‑Cccggggg                                                                                                        >  2:782526/1‑151 (MQ=255)
                                              aaCAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGC‑GTGCTCGGTGT‑GT‑Ccgggcggg                                                                                                       <  2:115959/151‑1 (MQ=255)
                                              aaCAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGGGTAAGAATTGGCGCGCCCACCGGCGTGGC‑GTGGTCGGGGT‑GT‑Ccgggcggg                                                                                                       >  2:497605/1‑151 (MQ=255)
                                               aCAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGC‑GTGCTCGGTGT‑GT‑CCGGGCGGGa                                                                                                      <  2:1195522/151‑1 (MQ=255)
                                               aCAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGCAGGGGGAAGAATTGGCGCGCCCACCCGGGGGGG‑GGGCTCCGTGTGGT‑Ccggggggg                                                                                                       >  1:730764/1‑151 (MQ=255)
                                                cAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATGTGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGC‑GTGCTCGGTGT‑GT‑CCGGGCGGGAt                                                                                                     <  1:665076/151‑1 (MQ=255)
                                                cAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCCGGGGTAAGAATTGGGGCCCCCCCCCGCCGGGC‑GGGGGCGGGGT‑GG‑Ccggggggg                                                                                                       >  1:265937/1‑149 (MQ=255)
                                                cAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGC‑GTGCTCGGTGT‑GT‑CCGGGCGGGAt                                                                                                     <  2:507084/151‑1 (MQ=255)
                                                cAATAACTAAAAGGGCACTTTATATGAGCAGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGC‑GTGCTCGGTGT‑GT‑CCGGGCGGGAt                                                                                                     <  2:40087/151‑1 (MQ=255)
                                                  aTAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGC‑GTGCTCGGTGT‑GT‑CCGGGCGGGAT‑Ga                                                                                                  <  1:5138/151‑1 (MQ=255)
                                                  aTAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGCAGGGGTAAGAAATGGCGCGCCCACCGGCGGGGC‑GTGCCCGGGGG‑GT‑Cggggggggggg                                                                                                    >  2:56306/1‑147 (MQ=255)
                                                     aCTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTGCCACAGGTCCGGGGTCGGCAGGGGGAAAAAATTGCGCTCCCCCCCGCGCGGC‑GGGCTCGgcgt‑gggccgggggggtg                                                                                                     >  1:1410648/1‑134 (MQ=255)
                                                       tAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGCAGGGGTAAGAATTGGCGCGCCCACCGGCGTGGG‑GGTCTCGGGGT‑GT‑TCCGGCGGGGT‑g                                                                                                   >  2:439972/1‑145 (MQ=255)
                                                       tAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGCAGGGGTAAGAATTGGCGCGCCCACCGGCGTGGG‑GGGCTCGGGGT‑GT‑CCGGGGGGGGT‑GA‑Att                                                                                              >  2:1303630/1‑149 (MQ=255)
                                                        aaaaGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGC‑GTGCTCGGTGT‑GT‑CCGGGCGGGAT‑GA‑CTTC‑gg                                                                                          <  1:468039/151‑1 (MQ=255)
                                                         aaaGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGC‑GTGCTCGGTGT‑GT‑CCGGGCGGGAT‑GA‑CTTC‑ggg                                                                                         <  2:583265/151‑1 (MQ=255)
                                                         aaaGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGC‑GTGCTCGGTGT‑GT‑CCGGGCGGGAT‑GA‑CTTC‑ggg                                                                                         <  2:111523/151‑1 (MQ=255)
                                                           aGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGC‑GTGCTCGGTGT‑GT‑CCGGGCGGGAT‑GA‑CTTC‑GGGC‑a                                                                                      <  1:982877/151‑1 (MQ=255)
                                                            gggCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGC‑GTGCTCGGTGT‑GT‑CCGGGCGGGAT‑GA‑CTTC‑GGGC‑a                                                                                      <  2:470396/150‑1 (MQ=255)
                                                              gCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGC‑GTGCTCGGTGT‑GT‑CCGGGCGGGAT‑GA‑CTTC‑GGGC‑AAcc                                                                                   <  2:374048/151‑1 (MQ=255)
                                                               cACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGC‑GTGCTCGGTGT‑GT‑CCGGGCGGGAT‑GA‑CTTC‑GGGC‑AAccc                                                                                  <  1:443205/151‑1 (MQ=255)
                                                                aCTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGGGTAAGAATTGGCGCGCCCACCCGCGGGGG‑GTGCTCGGGGT‑GT‑CCGGGGGGGGG‑GA‑CTTC‑CGGG‑AAacc                                                                                  >  2:166711/1‑150 (MQ=255)
                                                                 cTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCCGGGGTCGGCCGGGGTAAGAATTGGCGCGCCCACCCGCGGGGG‑GTGGCCGGGGG‑GG‑CCGGGGGGGGTGGC‑CTTC‑GGGC‑AACCCg                                                                                 >  2:378537/1‑151 (MQ=255)
                                                                  tttATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTTTCCAGGGGTCGGCAGGGGTAAGAATTGGCGCGCCCACCGGCGGGGC‑GGGGTCGGGGT‑GT‑CCCGGGGGGGT‑GACCTTC‑GGGA‑AACCggg                                                                                >  2:699027/1‑151 (MQ=255)
                                                                   ttATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGC‑GTGCTCGGTGT‑GT‑CCGGGCGGGAT‑GA‑CTTC‑GGGC‑AACCCGGt                                                                               <  2:1386780/150‑1 (MQ=255)
                                                                   ttATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGC‑GTGCTCGGTGT‑GT‑CCGGGCGGGAT‑GA‑CTGC‑GGGC‑AACCCgg                                                                                <  1:966383/149‑1 (MQ=255)
                                                                     atatGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGGCCCCTTGTTCAGGGGTCGGCCGGGGTAAAAAATGGCGGGCCCCCCCGGGGGGG‑GGGGCCGGGGG‑GTGCCGGGGGGGGT‑GA‑CCTC‑GGGG‑AACCCggg                                                                               >  1:778825/1‑148 (MQ=255)
                                                                       atGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGC‑GTGCTCGGTGT‑GT‑CCGGGCGGGAT‑GA‑CTTC‑AGGC‑AACCCGGTAAAt                                                                           <  1:800481/150‑1 (MQ=255)
                                                                       atGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGGGTAAGAATTGGCGCGCCCACCGGCGTGGC‑GTGGTCGGGGT‑GT‑CCGGGCGGGGT‑GA‑CTTC‑GGGC‑AACCCGGGaaa                                                                            >  2:298992/1‑149 (MQ=255)
                                                                       atGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGCCGGGGTAAGAAATGGCGCGCCCACCGGCGTGGCGGGGCCGGGGGG‑GT‑CCGGG‑GGGGG‑GG‑CTTC‑CGGC‑AACCCCGGaaa                                                                            >  1:657050/1‑149 (MQ=255)
                                                                        tGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCCTTGTCCAGGGGTCGGCAGGGGTAAGAATTGGCGCGCCCACCCGCGTGGG‑GGGCTCGGGGT‑GT‑CCGGGGGGGGT‑GA‑CTTCGGGGC‑AACCCGGGaaacc                                                                          >  2:685117/1‑149 (MQ=255)
                                                                         gAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCGCACCGGCGTGGC‑GTGCTCGGTGT‑GT‑CCGGGCGGGAT‑GA‑CTTC‑GGGC‑AACCCGTTAAATCCg                                                                        <  1:668682/151‑1 (MQ=255)
                                                                          aGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCCTTGTCCCGGGGTCGGCCGGGGTGAGAATTGGGGCGCCCCCCCGCGGGGG‑GTGGTCGGGGT‑GG‑CCGGGGGGGGT‑GG‑CTTC‑CGGC‑AACCCCGGAAAACCgc                                                                       >  1:223333/1‑151 (MQ=255)
                                                                          aGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGC‑GTGCTCGGTGT‑GT‑CCGGGCGGGAT‑GA‑CTTC‑GGGC‑AACCCGGTAAATCCgc                                                                       <  1:1076522/151‑1 (MQ=255)
                                                                           gCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCCTTGTCCAGGGGTCGGCAGGGGTAAGAATTGGGGCGCCCACCGGCGTGGC‑GTGCTCGGGGT‑GG‑CCGGGGGGGGT‑GG‑CCTC‑GGGC‑AACCCGGGAAATcc                                                                         >  2:1122501/1‑148 (MQ=255)
                                                                           gCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGC‑GTGCTCGGTGT‑GT‑CCGGGCGGGAT‑GA‑CTTC‑GGGC‑AACCCGGTAAATCCgct                                                                      <  1:765382/151‑1 (MQ=255)
                                                                              gAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGC‑GTGCTCGGTGT‑GT‑CCGGGCGGGAT‑GA‑CTTC‑GGGC‑AACCCGGTAAATCCgctgct                                                                   <  2:902139/151‑1 (MQ=255)
                                                                                aaaCCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGCAGGGGTAAGAATTGGCGCGCCCACCGGCGTGGC‑GTGGTCGGGGT‑GT‑CCCGGGGGGGG‑GA‑CTTCGGGGA‑AACCCGGAAAAcccccgccgg                                                                  >  2:36319/1‑144 (MQ=255)
                                                                                  aCCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGC‑GTGCTCGGTGT‑GT‑CCGGGCGGGAT‑GA‑CTTC‑GGGC‑AACCCGGTAAATCCGCTGCTgggg                                                               <  1:231972/151‑1 (MQ=255)
                                                                                   ccAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGC‑GTGCTCGGTGT‑GT‑CCGGGCGGGAT‑GA‑CTTC‑GGGC‑AACCCGGTAAATCCGCTGCTggggg                                                              <  1:1102595/151‑1 (MQ=255)
                                                                                    cAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGCAGGGGTAAGAATTGGCGCGCCCACCGGCGTGGC‑GTGCTCGGGGT‑GT‑CCGGGCGGGGT‑GA‑CTTC‑GGGC‑AACCCGGGAAATCCCCTGCTggggg                                                              >  2:1250385/1‑150 (MQ=255)
                                                                                    cAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGCAGGGGGAAGAATTGGCGCGCCCACCGGCGGGGC‑GTGCTCGGGGG‑GT‑CCGGGGGGGAT‑GA‑CTTC‑GGGC‑AACCCGGGAAATCCGCTGCTggggg                                                              >  2:405582/1‑150 (MQ=255)
                                                                                        ggcgCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGC‑GTGCTCGGTGT‑GT‑CCGGGCGGGAT‑GA‑CTTC‑GGGC‑AACCCGGTAAATCCGCTGCTGGGGGCg                                                            <  1:1386033/148‑1 (MQ=255)
                                                                                        ggcgCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGC‑GTGCTCGGTGT‑GT‑CCGGGCGGGAT‑GA‑CTTC‑GGGC‑AACCCGGTAAATCCGCTGCTGGGGGCGAAg                                                         <  2:847156/151‑1 (MQ=255)
                                                                                          cgcgTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGC‑GTGCTCGGTGT‑GT‑CCGGGCGGGAT‑GA‑CTTC‑GGGC‑AACCCGGTAAATCCGCTGCTGGGGGCGAAGGt                                                       <  2:721601/151‑1 (MQ=255)
                                                                                          cgcgTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGC‑GTGCTCGGTGT‑GT‑CCGGGCGGGAT‑GA‑CTTC‑GGGC‑AACCCGGTAAATCCGCTGCTGGGGGCGAAGGt                                                       <  1:368728/151‑1 (MQ=255)
                                                                                           gcgTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGC‑GTGCTCGGTGT‑GT‑CCGGGCGGGAT‑GA‑CTTC‑GGGC‑AACCCGGTAAATCCGCTGCTGGGGGCGAAGGtg                                                      <  2:899683/151‑1 (MQ=255)
                                                                                            cgTCAGGGAGATATGACTCAGTATGGCGGTCCCCTTGTCCCGGGGTCGGGAGGGGGCAGAATTTGCGCGCCCACCTGCGGGGC‑GGGCTCCGGGTGGT‑CCGGGGGGGGG‑GG‑TTTC‑GGGC‑CACCCGGAAAAACGGCT‑CTGGGGGCggaag                                                        >  1:746671/1‑144 (MQ=255)
                                                                                             gTCAGGGAGATATGACTCAGTATGGCGGTCCAATTTTCCAGGGGTCGGGAGGGGTTAGAAATGTCGGCCCCCCCCCCGGGGG‑GCGATTGGTTT‑GT‑CCGGGCGGGGG‑AA‑CTTC‑GGGCAAACCCGGAAAACCCCCGGCGCGGGGGGAAgg                                                        >  1:1166272/1‑148 (MQ=255)
                                                                                               cAGGGAGATATGACTCAGTATGGCGGTCCCACTGTCCAGGGGTCGGCAGGGGTAAGAATTGGCGCGCCCACCGGCGGGGC‑GTGCGCGGGGT‑GT‑CCGGGGGGGGT‑GA‑CTTCGGGGC‑ACCCCGGGAAA‑CCCCCGCTGGGGGGGGAgggg                                                      >  1:421659/1‑145 (MQ=255)
                                                                                                aGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGGGTAAGAATTGGCGCGCCCACCGGCGTGGC‑GTGCTCGGTGT‑GT‑CCGGGGGGGGT‑GA‑CATCGGGGA‑AACCCGGAAAAT‑CGCCGCTGGGGGGGGAGGGGGTGcc                                                 >  2:1351597/1‑151 (MQ=255)
                                                                                                      atatGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGC‑GTGCTCGGTGT‑GT‑CCGGGCGGGAT‑GA‑CTTC‑GGGC‑AACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCga                                           <  2:1328625/151‑1 (MQ=255)
                                                                                                           aCTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGC‑GTGCTCGGTGT‑GT‑CCGGGCGGGAT‑GA‑CTTC‑GGGC‑AACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACgg                                      <  2:457212/151‑1 (MQ=255)
                                                                                                           aCTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGC‑GTGCTCGGTGT‑GT‑CCGGGCGGGAT‑GA‑CTTC‑GGGC‑AACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACgg                                      <  1:1113313/151‑1 (MQ=255)
                                                                                                              cAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGC‑GTGCTCGGTGT‑GT‑CCGGGCGGGAT‑GA‑CTTC‑GGGC‑AACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGAc                                    <  1:1252269/150‑1 (MQ=255)
                                                                                                                  aTGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGC‑GTGCTCGGTGT‑GT‑CCGGGCGGGAT‑GA‑CTTC‑GGGC‑AACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTgc                               <  2:1069885/151‑1 (MQ=255)
                                                                                                                    ggCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGC‑GTGCTCGGTGT‑GT‑CCGGGCGGGAT‑GA‑CTTC‑GGGC‑AACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTgc                               >  2:481819/1‑149 (MQ=255)
                                                                                                                     gCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGC‑GTGCTCGGTGT‑GT‑CCGGGCGGGAT‑GA‑CTTC‑GGGC‑AACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCt                            <  2:390365/151‑1 (MQ=255)
                                                                                                                         tCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGC‑GTGCTCGGTGT‑GT‑CCGGGCGGGAT‑GA‑CTTC‑GGGC‑AACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTgccc                        >  1:200209/1‑151 (MQ=255)
                                                                                                                          cccATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGC‑GTGCTCGGTGT‑GT‑CCGGGCGGGAT‑GA‑CTTC‑GGGC‑AACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTgcccg                       <  2:821779/151‑1 (MQ=255)
                                                                                                                            cATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGC‑GTGCTCGGTGT‑GT‑CCGGGCGGGAT‑GA‑CTTC‑GGGC‑AACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTgcccgg                      <  2:85821/150‑1 (MQ=255)
                                                                                                                                    aGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGC‑GTGCTCGGTGT‑GT‑CCGGGCGGGAT‑GA‑CTTC‑GGGC‑AACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGcc             >  1:174722/1‑151 (MQ=255)
                                                                                                                                    aGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGC‑GTGCTCGGTGT‑GT‑CCGGGCGGGAT‑GA‑CTTC‑GGGC‑AACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGcc             <  1:533617/151‑1 (MQ=255)
                                                                                                                                     gggTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGC‑GTGCTCGGTGT‑GT‑CCGGGCGGGAT‑GA‑CTTC‑GGGC‑AACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCg            <  2:561671/151‑1 (MQ=255)
                                                                                                                                           ggCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGC‑GTGCTCGGTGT‑GT‑CCGGGCGGGAT‑GA‑CTTC‑GGGC‑AACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCAtt      <  1:1033475/151‑1 (MQ=255)
                                                                                                                                            gCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGC‑GTGCTCGGTGT‑GT‑CCGGGCGGGAT‑GA‑CTTC‑GGGC‑AACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCAtt      >  2:744770/1‑150 (MQ=255)
                                                                                                                                              aGGTGTAAGAATTGGCGCGCCCACCGGCGTGGC‑GTGCTCGGTGT‑GT‑CCGGGCGGGAT‑GA‑CTTC‑GGGC‑AACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATtctc   >  2:1017174/1‑151 (MQ=255)
                                                                                                                                                gtgtAAGAATTGGCGCGCCCACCGGCGTGGC‑GTGCTCGGTGT‑GT‑CCGGGCGGGAT‑GA‑CTTC‑GGGC‑AACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATtctct  <  1:1203938/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                           
TCTTGTGTGTCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGC‑GTGCTCGGTGT‑GT‑CCGGGCGGGAT‑GA‑CTTC‑GGGC‑AACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCT  >  NC_000913/523190‑523483

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