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→GGGrhsDRhs 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

TCTTGTGTGTCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCT  >  NC_000913/523190‑523483
                                                                                                                                                 |                                                                                                                                                           
tCTTGTGTGTCATTTTATTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTTGGCAGgggt                                                                                                                                                           >  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)
                             gCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCCTTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGgcgtggcgtg                                                                                                                            >  2:166613/1‑149 (MQ=255)
                             gCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCCGGGGTCGGCAGGTGTAAGAAATGGGGCGCCCACCcgcgtgg                                                                                                                                 >  1:299982/1‑145 (MQ=255)
                              cTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCCTTGTCCAGGGGTCGGCAGGGGGAAGAATTTGGGCGCCCACCggcgtggcg                                                                                                                              >  2:950672/1‑146 (MQ=255)
                                aaCACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGGCCCCTTGTGCCGGGGTCGGGCGGGGTTTGAAATGGGGCGCCCACCCGGGTCGCGGGcc                                                                                                                          >  2:889672/1‑147 (MQ=255)
                                 acacAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGgt                                                                                                                      <  1:1094415/151‑1 (MQ=255)
                                 acacAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGCAGGGGTAAGAATTGGCGCGCCCACCGGCGTGGCGGGCTCgg                                                                                                                       >  2:1232544/1‑150 (MQ=255)
                                             gAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTTCCGGGGTCGGCAGGGGGAAGAATTGGCGCGCCCCGCGGCGGGGGGGGCTCCGGGTGTCccggggg                                                                                                        >  2:782526/1‑151 (MQ=255)
                                              aaCAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCcgggcggg                                                                                                       <  2:115959/151‑1 (MQ=255)
                                              aaCAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGGGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGGTCGGGGTGTCcgggcggg                                                                                                       >  2:497605/1‑151 (MQ=255)
                                               aCAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGa                                                                                                      <  2:1195522/151‑1 (MQ=255)
                                               aCAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGCAGGGGGAAGAATTGGCGCGCCCACCCGGGGGGGGGGCTCCGTGTGGTCcggggggg                                                                                                       >  1:730764/1‑151 (MQ=255)
                                                cAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATGTGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGAt                                                                                                     <  1:665076/151‑1 (MQ=255)
                                                cAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCCGGGGTAAGAATTGGGGCCCCCCCCCGCCGGGCGGGGGCGGGGTGGCcggggggg                                                                                                       >  1:265937/1‑149 (MQ=255)
                                                cAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGAt                                                                                                     <  2:507084/151‑1 (MQ=255)
                                                cAATAACTAAAAGGGCACTTTATATGAGCAGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGAt                                                                                                     <  2:40087/151‑1 (MQ=255)
                                                  aTAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGa                                                                                                  <  1:5138/151‑1 (MQ=255)
                                                  aTAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGCAGGGGTAAGAAATGGCGCGCCCACCGGCGGGGCGTGCCCGGGGGGTCggggggggggg                                                                                                    >  2:56306/1‑147 (MQ=255)
                                                     aCTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTGCCACAGGTCCGGGGTCGGCAGGGGGAAAAAATTGCGCTCCCCCCCGCGCGGCGGGCTCGgcgtgggccgggggggtg                                                                                                     >  1:1410648/1‑134 (MQ=255)
                                                       tAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGCAGGGGTAAGAATTGGCGCGCCCACCGGCGTGGGGGTCTCGGGGTGTTCCGGCGGGGTg                                                                                                   >  2:439972/1‑145 (MQ=255)
                                                       tAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGCAGGGGTAAGAATTGGCGCGCCCACCGGCGTGGGGGGCTCGGGGTGTCCGGGGGGGGTGAAtt                                                                                              >  2:1303630/1‑149 (MQ=255)
                                                        aaaaGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCgg                                                                                          <  1:468039/151‑1 (MQ=255)
                                                         aaaGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCggg                                                                                         <  2:583265/151‑1 (MQ=255)
                                                         aaaGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCggg                                                                                         <  2:111523/151‑1 (MQ=255)
                                                           aGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCa                                                                                      <  1:982877/151‑1 (MQ=255)
                                                            gggCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCa                                                                                      <  2:470396/150‑1 (MQ=255)
                                                              gCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAAcc                                                                                   <  2:374048/151‑1 (MQ=255)
                                                               cACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAAccc                                                                                  <  1:443205/151‑1 (MQ=255)
                                                                aCTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGGGTAAGAATTGGCGCGCCCACCCGCGGGGGGTGCTCGGGGTGTCCGGGGGGGGGGACTTCCGGGAAacc                                                                                  >  2:166711/1‑150 (MQ=255)
                                                                 cTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCCGGGGTCGGCCGGGGTAAGAATTGGCGCGCCCACCCGCGGGGGGTGGCCGGGGGGGCCGGGGGGGGTGGCCTTCGGGCAACCCg                                                                                 >  2:378537/1‑151 (MQ=255)
                                                                  tttATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTTTCCAGGGGTCGGCAGGGGTAAGAATTGGCGCGCCCACCGGCGGGGCGGGGTCGGGGTGTCCCGGGGGGGTGACCTTCGGGAAACCggg                                                                                >  2:699027/1‑151 (MQ=255)
                                                                   ttATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGt                                                                               <  2:1386780/150‑1 (MQ=255)
                                                                   ttATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTGCGGGCAACCCgg                                                                                <  1:966383/149‑1 (MQ=255)
                                                                     atatGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGGCCCCTTGTTCAGGGGTCGGCCGGGGTAAAAAATGGCGGGCCCCCCCGGGGGGGGGGGCCGGGGGGTGCCGGGGGGGGTGACCTCGGGGAACCCggg                                                                               >  1:778825/1‑148 (MQ=255)
                                                                       atGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCAGGCAACCCGGTAAAt                                                                           <  1:800481/150‑1 (MQ=255)
                                                                       atGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGGGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGGTCGGGGTGTCCGGGCGGGGTGACTTCGGGCAACCCGGGaaa                                                                            >  2:298992/1‑149 (MQ=255)
                                                                       atGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGCCGGGGTAAGAAATGGCGCGCCCACCGGCGTGGCGGGGCCGGGGGGGTCCGGGGGGGGGGCTTCCGGCAACCCCGGaaa                                                                            >  1:657050/1‑149 (MQ=255)
                                                                        tGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCCTTGTCCAGGGGTCGGCAGGGGTAAGAATTGGCGCGCCCACCCGCGTGGGGGGCTCGGGGTGTCCGGGGGGGGTGACTTCGGGGCAACCCGGGaaacc                                                                          >  2:685117/1‑149 (MQ=255)
                                                                         gAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCGCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGTTAAATCCg                                                                        <  1:668682/151‑1 (MQ=255)
                                                                          aGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCCTTGTCCCGGGGTCGGCCGGGGTGAGAATTGGGGCGCCCCCCCGCGGGGGGTGGTCGGGGTGGCCGGGGGGGGTGGCTTCCGGCAACCCCGGAAAACCgc                                                                       >  1:223333/1‑151 (MQ=255)
                                                                          aGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCgc                                                                       <  1:1076522/151‑1 (MQ=255)
                                                                           gCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCCTTGTCCAGGGGTCGGCAGGGGTAAGAATTGGGGCGCCCACCGGCGTGGCGTGCTCGGGGTGGCCGGGGGGGGTGGCCTCGGGCAACCCGGGAAATcc                                                                         >  2:1122501/1‑148 (MQ=255)
                                                                           gCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCgct                                                                      <  1:765382/151‑1 (MQ=255)
                                                                              gAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCgctgct                                                                   <  2:902139/151‑1 (MQ=255)
                                                                                aaaCCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGCAGGGGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGGTCGGGGTGTCCCGGGGGGGGGACTTCGGGGAAACCCGGAAAAcccccgccgg                                                                  >  2:36319/1‑144 (MQ=255)
                                                                                  aCCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTgggg                                                               <  1:231972/151‑1 (MQ=255)
                                                                                   ccAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTggggg                                                              <  1:1102595/151‑1 (MQ=255)
                                                                                    cAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGCAGGGGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGGGTGTCCGGGCGGGGTGACTTCGGGCAACCCGGGAAATCCCCTGCTggggg                                                              >  2:1250385/1‑150 (MQ=255)
                                                                                    cAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGCAGGGGGAAGAATTGGCGCGCCCACCGGCGGGGCGTGCTCGGGGGGTCCGGGGGGGATGACTTCGGGCAACCCGGGAAATCCGCTGCTggggg                                                              >  2:405582/1‑150 (MQ=255)
                                                                                        ggcgCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCg                                                            <  1:1386033/148‑1 (MQ=255)
                                                                                        ggcgCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAg                                                         <  2:847156/151‑1 (MQ=255)
                                                                                          cgcgTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGt                                                       <  2:721601/151‑1 (MQ=255)
                                                                                          cgcgTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGt                                                       <  1:368728/151‑1 (MQ=255)
                                                                                           gcgTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGtg                                                      <  2:899683/151‑1 (MQ=255)
                                                                                            cgTCAGGGAGATATGACTCAGTATGGCGGTCCCCTTGTCCCGGGGTCGGGAGGGGGCAGAATTTGCGCGCCCACCTGCGGGGCGGGCTCCGGGTGGTCCGGGGGGGGGGGTTTCGGGCCACCCGGAAAAACGGCTCTGGGGGCggaag                                                        >  1:746671/1‑144 (MQ=255)
                                                                                             gTCAGGGAGATATGACTCAGTATGGCGGTCCAATTTTCCAGGGGTCGGGAGGGGTTAGAAATGTCGGCCCCCCCCCCGGGGGGCGATTGGTTTGTCCGGGCGGGGGAACTTCGGGCAAACCCGGAAAACCCCCGGCGCGGGGGGAAgg                                                        >  1:1166272/1‑148 (MQ=255)
                                                                                               cAGGGAGATATGACTCAGTATGGCGGTCCCACTGTCCAGGGGTCGGCAGGGGTAAGAATTGGCGCGCCCACCGGCGGGGCGTGCGCGGGGTGTCCGGGGGGGGTGACTTCGGGGCACCCCGGGAAACCCCCGCTGGGGGGGGAgggg                                                      >  1:421659/1‑145 (MQ=255)
                                                                                                aGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGGGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGGGGGGTGACATCGGGGAAACCCGGAAAATCGCCGCTGGGGGGGGAGGGGGTGcc                                                 >  2:1351597/1‑151 (MQ=255)
                                                                                                      atatGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCga                                           <  2:1328625/151‑1 (MQ=255)
                                                                                                           aCTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACgg                                      <  2:457212/151‑1 (MQ=255)
                                                                                                           aCTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACgg                                      <  1:1113313/151‑1 (MQ=255)
                                                                                                              cAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGAc                                    <  1:1252269/150‑1 (MQ=255)
                                                                                                                  aTGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTgc                               <  2:1069885/151‑1 (MQ=255)
                                                                                                                    ggCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTgc                               >  2:481819/1‑149 (MQ=255)
                                                                                                                     gCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCt                            <  2:390365/151‑1 (MQ=255)
                                                                                                                         tCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTgccc                        >  1:200209/1‑151 (MQ=255)
                                                                                                                          cccATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTgcccg                       <  2:821779/151‑1 (MQ=255)
                                                                                                                            cATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTgcccgg                      <  2:85821/150‑1 (MQ=255)
                                                                                                                                    aGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGcc             >  1:174722/1‑151 (MQ=255)
                                                                                                                                    aGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGcc             <  1:533617/151‑1 (MQ=255)
                                                                                                                                     gggTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCg            <  2:561671/151‑1 (MQ=255)
                                                                                                                                           ggCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCAtt      <  1:1033475/151‑1 (MQ=255)
                                                                                                                                            gCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCAtt      >  2:744770/1‑150 (MQ=255)
                                                                                                                                              aGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATtctc   >  2:1017174/1‑151 (MQ=255)
                                                                                                                                                gtgtAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATtctct  <  1:1203938/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                           
TCTTGTGTGTCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCT  >  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: