Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,608,721 A→C 6.2% intergenic (+72/+30) zntA → / ← tusA zinc, cobalt and lead efflux system/mnm(5)‑s(2)U34‑tRNA 2‑thiolation sulfurtransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,608,7210AC 6.2% 721.2 / 9.9 241intergenic (+72/+30)zntA/tusAzinc, cobalt and lead efflux system/mnm(5)‑s(2)U34‑tRNA 2‑thiolation sulfurtransferase
Reads supporting (aligned to +/- strand):  ref base A (151/70);  new base C (2/18);  total (153/88)
Fisher's exact test for biased strand distribution p-value = 3.70e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 258 total aligned reads displayed.

TGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCGGATCGTCGG  >  NC_000913/3608572‑3608869
                                                                                                                                                     |                                                                                                                                                    
tggctggcAGTGCTGGCAGATACGGGGGCGACGGTGCTTGTGACAGCGAATGCGTTAAGATTGTTGAGCAGGAGATAATGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCATCTGGGGAGAGGGTTGGGGTGAGGGGaa                                                                                                                                                      >  2:1741844/1‑150 (MQ=255)
tggctggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGaa                                                                                                                                                      >  1:2988081/1‑150 (MQ=255)
tggctggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGaa                                                                                                                                                      >  1:3548345/1‑150 (MQ=255)
tggctggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGaa                                                                                                                                                      >  1:3548370/1‑150 (MQ=255)
  gctggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGaaaa                                                                                                                                                    >  2:389725/1‑150 (MQ=255)
  gctggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGaaaa                                                                                                                                                    >  1:594962/1‑150 (MQ=255)
  gctggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGaaaa                                                                                                                                                    <  2:3685258/150‑1 (MQ=255)
   ctggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAg                                                                                                                                                   >  1:1728945/1‑150 (MQ=255)
    tggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCACGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAgg                                                                                                                                                  >  2:2634462/1‑150 (MQ=255)
     ggcAGGGCTGGCAGATAGGGGGGCGACGGTGCTGGTTACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTTAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGTTTAGGGTGAGGGGAAAAggc                                                                                                                                                 <  2:3978307/150‑1 (MQ=39)
            cTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCg                                                                                                                                          >  1:2695799/1‑150 (MQ=255)
             tGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTGAGGGTGAGGGGAAAAGGCGGCATCGa                                                                                                                                         >  2:397061/1‑150 (MQ=255)
              ggCTGATACGGGGGCGACGGTGCTGGTGACAGCGATTGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGaa                                                                                                                                        <  2:3126272/150‑1 (MQ=255)
                   aTACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCaa                                                                                                                                   >  2:2948772/1‑150 (MQ=255)
                      cGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCa                                                                                                                                >  2:840484/1‑150 (MQ=255)
                         gggCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGcc                                                                                                                             >  1:2568073/1‑150 (MQ=255)
                          ggCGACGGTGCTGGTGCCAGCGAATGCGTTATGATTGTTGCGCTGGAGATAAGCCAACCCGATCGAAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGccc                                                                                                                            <  2:471327/150‑1 (MQ=38)
                          ggCGACGGTGCTGGTGACAGCTATTGCGTTAATTTTGTTGCGCAGTATATAAGGCAACCCGATCGAAAATTTGAGCGCGATGGGTCCCCTCGCCCCTCTGGGGATTGGTTTAGGGTGCGGGGAAAAGGCGGCATCGAAGCCAATCAGccc                                                                                                                            <  2:4235263/150‑1 (MQ=255)
                          ggCGACGGTGCTGGTCACACCGAATGCTATAACATTGTTCCGCAGCACAACCGCCAAACCCAACGCAACATTCAGCGACAACGGCCCCCACGCCCCACTGGGGAGAGGTTACGGGCGAGGGGAAAAGGCGGCATCGAAGCCAATCAGccc                                                                                                                            <  2:2008560/150‑1 (MQ=255)
                             gACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTa                                                                                                                         >  1:1146947/1‑150 (MQ=255)
                                ggTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCa                                                                                                                      >  1:706801/1‑150 (MQ=255)
                                   gCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGCGAGGGGTAGGGGGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCaaac                                                                                                                   >  2:4612201/1‑148 (MQ=255)
                                        tGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCtt                                                                                                              >  2:4614718/1‑150 (MQ=255)
                                         gACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCttt                                                                                                             >  1:2379137/1‑150 (MQ=255)
                                         gACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCGGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCttt                                                                                                             >  1:929041/1‑150 (MQ=255)
                                         gACAGCGAATGCGTTAAGATTGTTGCACAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCttt                                                                                                             >  1:3108576/1‑150 (MQ=255)
                                           cAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTAc                                                                                                           >  2:1072727/1‑150 (MQ=255)
                                           cAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTAc                                                                                                           >  2:3486388/1‑150 (MQ=255)
                                            aGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACg                                                                                                          >  1:3163809/1‑150 (MQ=255)
                                            aGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACg                                                                                                          >  1:2422504/1‑150 (MQ=255)
                                              cGAATGGTTTAAGTTTGTTGCTCTGTATAAACGCCAAACCGTCCCAAACTGTGCGCGCGATCGGCCCCCCCCCCCCGCGGGGGAGGGGGTTCGGGTGGGGGCACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGaa                                                                                                        <  1:3880370/150‑1 (MQ=255)
                                                   gCGTTAAGTTTGTTGGGCTGGATATAAGGCAACCGGATCGAAACATTGCGCGCGCTGGGCCCCCTCGCCCCGGTGGGTTGTGGGTTTGGGTGTGGGGACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCaaa                                                                                                   <  2:485171/150‑1 (MQ=255)
                                                     gTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATa                                                                                                 >  1:4414706/1‑150 (MQ=255)
                                                     gTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATa                                                                                                 >  1:4587796/1‑150 (MQ=255)
                                                     gTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATa                                                                                                 >  1:696990/1‑150 (MQ=255)
                                                     gTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATa                                                                                                 >  2:2321043/1‑150 (MQ=255)
                                                     gTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACAGTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGCGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATa                                                                                                 >  2:4321703/1‑150 (MQ=255)
                                                      ttAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATaa                                                                                                >  1:2644963/1‑150 (MQ=255)
                                                      ttAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATaa                                                                                                >  1:2782724/1‑150 (MQ=255)
                                                         aGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCACCATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGTGGTTTAGGGTGCGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGa                                                                                             <  2:1111702/150‑1 (MQ=255)
                                                           tttgttgCGCTGGATATAAGGCAAACCGATCGCAACATTGCGCGCGATCGGTCCCCTCGCCCATCTGGGGAGTGGGTTCGGGTGCGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATa                                                                                           <  2:592117/149‑1 (MQ=255)
                                                            ttgttgCGCTGGTGTTAAGGAAAACCGATCGCAACATTTCGCGCGATCGGTCCCCTCGCCCCTGTGGGGTGTGGGTTCGGGGGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATaa                                                                                          <  2:4687340/150‑1 (MQ=255)
                                                            ttgttgCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGGGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATaa                                                                                          >  2:3331057/1‑150 (MQ=255)
                                                                tgCGCTGGAGATAAGGAAAACCGATCGAAACATTGAGCGCTATCGGTCCCCTCGCCCCTCTGGGGAGTGGTTTTGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCa                                                                                      <  2:289944/150‑1 (MQ=39)
                                                                 gcgcAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAg                                                                                     >  2:839700/1‑150 (MQ=255)
                                                                   gcAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTc                                                                                   >  1:1443421/1‑150 (MQ=255)
                                                                     aGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCa                                                                                 >  1:2909840/1‑150 (MQ=255)
                                                                           tAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGtc                                                                           >  1:1908862/1‑150 (MQ=255)
                                                                              ggCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTctct                                                                        >  2:1421027/1‑150 (MQ=255)
                                                                              ggCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTctct                                                                        >  1:4209954/1‑150 (MQ=255)
                                                                               gCAAACCTATGGCAACTTTGCGCGCGATCGGTCCCCGCGCCCCTCTGGGGAGTGGGTTCGGGTGAGGCGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCtt                                                                       <  2:4188782/150‑1 (MQ=255)
                                                                               gCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCtt                                                                       >  2:369675/1‑150 (MQ=255)
                                                                               gCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCAATCCGTCTCtt                                                                       >  1:220763/1‑150 (MQ=255)
                                                                               ccAACCCGACCGCAACATGGCGCGCGAGCGGCCCCCTCGCCCCCGGGGGTAGCGGGTTCGGGGGAGGGCACACGGCGGCATCGCCGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCtt                                                                       <  1:228343/149‑1 (MQ=255)
                                                                                cAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCttt                                                                      <  2:2422504/150‑1 (MQ=255)
                                                                                cAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCGGGGGAGAGGGTGAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTAGCAACCGCCTTTACGACTCAAATAACGATAAGGCAGTCCATCCGTCTCgtt                                                                      >  2:3464537/1‑150 (MQ=255)
                                                                                 aacccGATCGCAACATTGAGCGCGTGCGGGCCCCTCGCCCCTGGGGGGAGTGGGTTCGGGTGGGGGGAACAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTa                                                                     <  2:1519976/147‑1 (MQ=255)
                                                                                  aaCCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAg                                                                    >  2:1298977/1‑150 (MQ=255)
                                                                                          gCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGt                                                            >  2:1087957/1‑150 (MQ=255)
                                                                                           cAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGtt                                                           >  2:3641204/1‑150 (MQ=255)
                                                                                            aataTTGCGCGCGAGGGGTCCCCTCGCCCCTCTGGGGTGTGGGTTAGGGTGTGGGTACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTc                                                          <  2:4532517/147‑1 (MQ=255)
                                                                                            aaCATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTc                                                          >  2:3690229/1‑150 (MQ=255)
                                                                                              ctttGAGCGCGATGGGCCCCCTCGCCCCTGGGGGTTTTGGTTTTGGGTGCGGGTAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCgt                                                        <  2:1046710/148‑1 (MQ=255)
                                                                                               aTTGCGCGCGATCGGTCCCCTCGCCCCTGTGGGGAGTGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCgtg                                                       <  1:726238/150‑1 (MQ=255)
                                                                                               aTTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCgtg                                                       >  1:4029369/1‑150 (MQ=255)
                                                                                                  gcgcgcGATCGGTCCCCGCGCCCCTGGGGGGTGTGGTTTGGGGTGCGGGCACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTc                                                    <  2:1713212/148‑1 (MQ=255)
                                                                                                        gATGGGTCCCGCCGCCCCTCGGGGTAGAGGGTTCGGGTGCGGGGACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATaaa                                              <  1:1520450/150‑1 (MQ=255)
                                                                                                        gATCGGCCCCCTCGCCCCTCTGGGGTGTGGGTTTGGGTGTGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATaaa                                              <  1:4229647/150‑1 (MQ=255)
                                                                                                        gATCGCACCCCCCCCCCCTCTGGGGAGTGGTTTAGGCCGAGGCCAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATaaa                                              <  1:1361790/150‑1 (MQ=255)
                                                                                                         aTCGGTCCCCCCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAg                                             <  1:2634462/150‑1 (MQ=255)
                                                                                                             taacccTCTCCCCAGAGGGGCGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTAc                                         >  2:1145252/4‑150 (MQ=255)
                                                                                                                  cTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAAc                                    >  1:3345294/1‑150 (MQ=255)
                                                                                                                   tCGCCCCTGTGGGGTGTGGTTTTGGGTGTGGGTAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAAcc                                   <  2:2695799/150‑1 (MQ=255)
                                                                                                                   tCGCCCCTCTGGGGAGTGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAAcc                                   <  2:2993793/150‑1 (MQ=255)
                                                                                                                   tCGCCCCTCTGGGGAGAGGGTTCGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTCTCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAAcc                                   <  2:3517830/150‑1 (MQ=255)
                                                                                                                   accccccACTGGGGAGAGGGTAAGGCTGAGCCCAAAAGGCGCCATCGAACCCAATCAGCCCCTATCAACCGCCATTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAAcc                                   <  2:3306514/147‑1 (MQ=255)
                                                                                                                     gCCCCTCTGGGTTGAGTTTTTGGGTGAGGGCAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCa                                 <  1:1965117/150‑1 (MQ=255)
                                                                                                                     gCCCCTCAGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACTATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTACGTGTTCCATAAAGGTACAAAACCCa                                 >  1:2837575/1‑150 (MQ=255)
                                                                                                                      cccgtggggggTGTGGGTTTGGTGGCGGGCACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAg                                <  1:2099615/143‑1 (MQ=255)
                                                                                                                      ccccTCTGGGGTGTGGGTTTGGGTGTGGGTAAAAGGCGGCATCGAAGCCATTCAGCCCCTATCAACGGCCTTTACGAATCAAATAACTATATGGATTTCAATCCTTCTCTTTAGCAACATTTTGGTTTTCAATAAAGGTAAAAAACCCAg                                <  2:706801/150‑1 (MQ=255)
                                                                                                                      ccccTCTGGGGAGAGGGTTAGGGTGCGGGTAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAg                                <  1:289418/150‑1 (MQ=255)
                                                                                                                       cccTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAgg                               >  2:1907142/1‑150 (MQ=255)
                                                                                                                          tctGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAat                            >  2:947854/1‑150 (MQ=255)
                                                                                                                          tctGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAat                            >  2:3189098/1‑150 (MQ=255)
                                                                                                                           ctGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAata                           >  1:2548528/1‑150 (MQ=255)
                                                                                                                             ggggAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATc                         >  2:1332096/1‑150 (MQ=255)
                                                                                                                              gggtttgggTTCGGGTGCGGGCACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATcg                        <  1:840484/144‑1 (MQ=255)
                                                                                                                                gagaGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATcgcg                      >  2:2105686/1‑150 (MQ=255)
                                                                                                                                gagaGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATcgcg                      >  2:1956772/1‑150 (MQ=255)
                                                                                                                                    gggTTTGGGGGCGGGCACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTa                  <  1:1674274/150‑1 (MQ=255)
                                                                                                                                     ggTTAGGGTGAGGGGAAAAGGCGGTATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAg                 >  2:1143851/1‑150 (MQ=255)
                                                                                                                                     ggTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAg                 >  1:324289/1‑150 (MQ=255)
                                                                                                                                      gTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGAGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGt                >  2:1548901/1‑150 (MQ=255)
                                                                                                                                       ttAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTg               >  1:3674775/1‑150 (MQ=255)
                                                                                                                                           ggTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCg           <  1:471068/150‑1 (MQ=255)
                                                                                                                                               aGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCGGAtc       >  2:222516/1‑150 (MQ=255)
                                                                                                                                                   gAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCGGAtcgtcg   >  1:4092673/1‑150 (MQ=255)
                                                                                                                                                   gAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCGGAtcgtcg   <  2:4587796/150‑1 (MQ=255)
                                                                                                                                                    aaaaGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCGGATCGTCgg  <  1:3735333/150‑1 (MQ=255)
                                                                                                                                                    aaaaGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCGGATCGTCgg  <  1:101124/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCGGATCGTCGG  >  NC_000913/3608572‑3608869

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