Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,608,721 A→C 8.9% 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 8.9% 640.0 / 18.3 218intergenic (+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 (124/71);  new base C (2/21);  total (126/92)
Fisher's exact test for biased strand distribution p-value = 3.36e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 223 total aligned reads displayed.

TGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCGGATCGTCGG  >  NC_000913/3608572‑3608869
                                                                                                                                                     |                                                                                                                                                    
  gctggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGaaaa                                                                                                                                                    <  2:1210466/150‑1 (MQ=255)
   ctggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAg                                                                                                                                                   >  2:3779924/1‑150 (MQ=255)
    tggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAgg                                                                                                                                                  >  1:1098958/1‑150 (MQ=255)
    tggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAgg                                                                                                                                                  >  2:1786016/1‑150 (MQ=255)
    tggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAgg                                                                                                                                                  >  2:3239428/1‑150 (MQ=255)
    tggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAgg                                                                                                                                                  <  2:3310778/150‑1 (MQ=255)
        aGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGTGGAAAAggcggc                                                                                                                                              <  2:546549/150‑1 (MQ=255)
        aGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAAGGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGTGGAAAAggcggc                                                                                                                                              <  2:2486062/150‑1 (MQ=255)
           gCTGGCAGATACGGGGGCGAGGGTGCTGGTGACTGCGAATGCGTTATGATTGTTGCGCAGTAGATAAGGCAAACCGATCGCAACATTGAGCGCTATCGGTCCCCTCGCCCCTCTGGGGTGTGGTTTAGGGTGTGGGGAAAAGGCGGCATc                                                                                                                                           <  2:3121652/150‑1 (MQ=255)
            cTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGCGGGTTCGGGTGAGGGGAAAAGGCGGCATCg                                                                                                                                          >  2:1600461/1‑150 (MQ=255)
             tGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGa                                                                                                                                         >  1:1425342/1‑150 (MQ=255)
              ggCAGATCGGGGGGCGCCGTTGGTGGTGACAGCTCTTGCTTTATTTTTGTTGCGCCGGTTTTAAGGCAAACCGTTCGCAACTTTGGGGGGGTTGGGGCCCCTGGCCCTTGTGGGGTGTGGTTTTGGGTGAGGGGAAAAGGCGGCATCGaa                                                                                                                                        <  2:1928852/150‑1 (MQ=255)
               gCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAg                                                                                                                                       >  2:3509336/1‑150 (MQ=255)
                cagagcgggggggGACGGGGCTGGTGACAGCGACTTCTTTTTGTTTGTTGCGCTGGATATACGGCAACCCGATCGAAACTTTGCGGGCTCTGGGTCCCCTCGCCCCTGGGGGGTGTGGTTTTGTGTGTGGGGAAAAGGCGGCATCGAAGc                                                                                                                                      <  2:1080173/143‑1 (MQ=255)
                cAGATACGGGGGCGAGGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCCGGAGATCAGGCAAACCGATCGCAACATTGAGCGCTATCGGTCCCCTCGCCCCTCTGGGGAGTGGGTTCGGGTGAGGGGAAAAGGCGGCATCGAAGc                                                                                                                                      <  2:2691833/150‑1 (MQ=39)
                 aGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGGTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGcc                                                                                                                                     >  1:4003552/1‑150 (MQ=255)
                     aCGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATc                                                                                                                                 >  2:1848221/1‑150 (MQ=255)
                        ggggCGACGGTGCTGGTGACAGCGAATGCGTTATGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGTTGAGGTTTTGGGTGCGGGGAAAAGGCGGCATCGAAGCCAATCAGc                                                                                                                              <  2:1332918/150‑1 (MQ=255)
                        ggggCGACGGTGCGGGTGACAGCTAATGGTTTAAGATTGTTGCGCTGGATATATGGAAAACCGATCGCAACATTGAGCGCGATGGGTCCCCTCGCCCCTCTGGGGAGTGGGTTTGGGTGAGGGTAAAAGGCGGCATCGAAGCCAATCAGc                                                                                                                              <  2:1646462/150‑1 (MQ=37)
                             gtgggTGCTGGTGACTGGGAATGCGTTATGTTTGTTGCGCTGGTGATACGGCAAACCGATCGCAACTTTGTGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGCGGGTTTGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTa                                                                                                                         <  2:4090059/147‑1 (MQ=37)
                             gagggTGCTGGTGACAGCGAATGCTTTATGATTGTTGCGCTGGATATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGTGTGGGTTTGGGTGCGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTa                                                                                                                         <  2:3660498/147‑1 (MQ=255)
                             gACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTa                                                                                                                         >  2:1353109/1‑150 (MQ=255)
                                ggTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGGTAGGGGGAGGGGACAAGGCGGCATCGAAGCCCATCAGCCCCGATCa                                                                                                                      >  2:3374286/1‑150 (MQ=39)
                                  tGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAAc                                                                                                                    >  1:2130665/1‑150 (MQ=255)
                                     tGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAAccac                                                                                                                 >  1:3421584/1‑148 (MQ=255)
                                     tGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAAccac                                                                                                                 >  1:1642702/1‑148 (MQ=255)
                                        tGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCtt                                                                                                              >  2:2964702/1‑150 (MQ=255)
                                        tGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCtt                                                                                                              <  2:4142350/150‑1 (MQ=255)
                                         gACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCttt                                                                                                             >  1:2743086/1‑150 (MQ=255)
                                         gACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCttt                                                                                                             >  2:2114083/1‑150 (MQ=255)
                                           cAGCTAAGGCGTTAAGATTGTTGCGCAGGAGATACGGCAACCCGAGCGCAACATTGCGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGCGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTAc                                                                                                           <  2:492333/150‑1 (MQ=38)
                                             gCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGa                                                                                                         >  1:2910311/1‑150 (MQ=255)
                                             gCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGa                                                                                                         >  1:2532342/1‑150 (MQ=255)
                                              cGAATGCGTTATGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGaa                                                                                                        <  2:2410207/150‑1 (MQ=255)
                                              cGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGaa                                                                                                        >  1:2652243/1‑150 (MQ=255)
                                                     gTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATa                                                                                                 >  2:86862/1‑150 (MQ=255)
                                                     gTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATa                                                                                                 >  1:581068/1‑150 (MQ=255)
                                                      ttAAGATTGTTGCGCAGGAGATACGGCAAACCGATCGCAACTTTGCGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGTGGGTTTGGGTGAGGGTAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATaa                                                                                                <  2:3382285/150‑1 (MQ=39)
                                                       tatttttgttgCGCTGTTTATAAGGCAAACCGATCGAAACTTTGCGCGCTATGGGTCCCCTCGCCCCTGTGGGGAGTGGGTTCGGGTGTGGGTAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAAc                                                                                               <  1:3074710/145‑1 (MQ=255)
                                                                   gcAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTc                                                                                   >  1:2018599/1‑150 (MQ=255)
                                                                   gcAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTc                                                                                   >  1:2427386/1‑150 (MQ=255)
                                                                   gcAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTc                                                                                   >  1:2989407/1‑150 (MQ=255)
                                                                      ggAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCAt                                                                                >  1:767170/1‑150 (MQ=255)
                                                                        ttaTAAGCCAACCGGATCGAAACTTTGCGCGCGTTGGGCCCCCCCGCCCCTCGGGGGTGTGGGTTGGGGTGGGGGTACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATcc                                                                              <  1:1549090/148‑1 (MQ=255)
                                                                            aaGGAAAACGTATCGAAACATTGCGCGGTATCGGTCCCCTCGCCCCTCGGGGGCGCGGGTTCGGGGGCGGGCACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGtct                                                                          <  1:2504353/150‑1 (MQ=255)
                                                                              ggCAAACCGTTCGCAACTTTGCGCGCGATCGGTCCCCTCGCCCCTCTGGGTTGTGGTTTTGGGTGTGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTctct                                                                        <  2:2001798/150‑1 (MQ=255)
                                                                              ggCAAACCGATCGCAACATTGCGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTTGGGGGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTctct                                                                        <  2:312073/150‑1 (MQ=255)
                                                                              ggCAAACCGATCGCAACATTGAGCGCTATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTctct                                                                        <  2:275701/150‑1 (MQ=255)
                                                                               cccacccGATCGCAACATGGTGCGCGATCGGACCCCGCGCCCCTCTGGGTTTTGGTTTTGGTTGGGGGTAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCtt                                                                       <  2:1440639/145‑1 (MQ=255)
                                                                                 aaaCCGTTCCCAACCTTGCGCGCGCTCGGCCCCCCCGCCCCTGTGGGGTGTGGGTTTGGGGGCGGGCACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTa                                                                     <  1:3881792/150‑1 (MQ=255)
                                                                                 aaaCCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTa                                                                     >  1:3655190/1‑150 (MQ=255)
                                                                                 aaaCCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTa                                                                     >  2:210853/1‑150 (MQ=255)
                                                                                 aaaCCGAGCGAACCTTTGCGCGCGATCGGTCCCCCCGCCCCGGGGGGTAGTGGTTTCGGGTGTGGGCACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTa                                                                     <  1:510724/150‑1 (MQ=255)
                                                                                 aaaCCGACCGCAACTTTGCGGGCGATCGGCCCCCTCGCCCCTCTGGGTTGTGGTTTTGGGTGCGGGGACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTa                                                                     <  1:103173/150‑1 (MQ=255)
                                                                                    ccGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCa                                                                  >  1:3361280/1‑150 (MQ=255)
                                                                                      gATCGCAAGATGGAGCGCGTGCGGACCCCCCGCCCCTCGGGGGCGTGGTTTAGGGTGCGGGTACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAAc                                                                <  2:764813/150‑1 (MQ=255)
                                                                                      gATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAAc                                                                >  1:399331/1‑150 (MQ=255)
                                                                                       aTCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAAcc                                                               >  1:3085941/1‑150 (MQ=255)
                                                                                        tCGCAACATTGAGCGCTATCGGACCCCTCGCCCCTCTGGGGTGTGGGTTTGGGTGTGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCa                                                              <  1:2067906/150‑1 (MQ=255)
                                                                                        tCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGTGGGTTTGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCa                                                              <  1:3853994/150‑1 (MQ=255)
                                                                                         cgaaaCTTTGCGCGCGATGGGCCCCCGCGCCCCTGGGGGGTGTGGGTTCGGGGGCGGGCACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAg                                                             <  2:1098958/147‑1 (MQ=255)
                                                                                              cATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCgt                                                        >  2:4036054/1‑150 (MQ=255)
                                                                                               aTTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCgtg                                                       >  2:357983/1‑150 (MQ=255)
                                                                                                   aGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTcc                                                   <  2:1384419/150‑1 (MQ=255)
                                                                                                     cgcgATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCAt                                                 >  1:1018726/1‑150 (MQ=255)
                                                                                                        gATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATaaa                                              >  2:2231862/1‑150 (MQ=255)
                                                                                                         aTCGGGCCCCTCGCCCCTCTGGGGAGTGGGTTAGGGTGAGGGGACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAg                                             <  2:3859522/150‑1 (MQ=255)
                                                                                                           gggTCCCCGCGCCCGTGTGGGGCGTGGTTTGGGGTGGGGGGACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGt                                           <  2:1867167/149‑1 (MQ=255)
                                                                                                           cctaacccTCTCCCCAGAGGGGCGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGt                                           >  1:2335135/6‑150 (MQ=255)
                                                                                                           cctaacccTCTCCCCAGAGGGGCGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGACCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGt                                           >  1:2335358/6‑150 (MQ=255)
                                                                                                            ggTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTa                                          >  1:1163707/1‑150 (MQ=255)
                                                                                                             gTCCCCGCGCCCCTCTGGGGTGTGGGTTTGGGTGCGGGGACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTAc                                         <  2:1369278/150‑1 (MQ=255)
                                                                                                              tCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACa                                        >  1:99059/1‑150 (MQ=255)
                                                                                                               ccccTCGCCCCTCTGGGGAGTGGGTTAGGGTGCGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACaa                                       <  2:2298068/150‑1 (MQ=255)
                                                                                                                cccTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACaaa                                      >  1:906696/1‑150 (MQ=255)
                                                                                                                cccTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACaaa                                      >  2:2535181/1‑150 (MQ=255)
                                                                                                                cccTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACaaa                                      >  1:3412312/1‑150 (MQ=255)
                                                                                                                cccTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACaaa                                      >  1:345910/1‑150 (MQ=255)
                                                                                                                  cTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAAc                                    >  2:1420860/1‑150 (MQ=255)
                                                                                                                  cTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAAc                                    >  2:1527380/1‑150 (MQ=255)
                                                                                                                  cTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAAc                                    >  2:450648/1‑150 (MQ=255)
                                                                                                                  cTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAAc                                    >  2:73754/1‑150 (MQ=255)
                                                                                                                   tCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAAcc                                   >  1:115518/1‑150 (MQ=255)
                                                                                                                    ccccccACTGGGCAGAGGGTTAGCGTCACCCCAAAACCCCCCATCGAAGCCAATCACCCCCTATCAACCGCCTTTACGAAACAAATAACGATAACGCAGTCCATCCGTCACTTTAGCAACCAGTTCGTGTTCCATAAAGCTACAAAAccc                                  <  2:596358/148‑1 (MQ=255)
                                                                                                                       cccTCTGGGGAGTGGGTTTGGGTGTGGGGACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAgg                               <  2:2089889/150‑1 (MQ=255)
                                                                                                                         ctctGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGaa                             >  2:1561432/1‑150 (MQ=255)
                                                                                                                           ctGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAata                           >  1:80492/1‑150 (MQ=255)
                                                                                                                                 tgtgggTTTGGGTGTGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCgg                     <  2:2156214/147‑1 (MQ=255)
                                                                                                                                 agagGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCgg                     <  1:3128512/150‑1 (MQ=255)
                                                                                                                                  gagGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCggg                    >  1:4117441/1‑150 (MQ=255)
                                                                                                                                    ggtttggggTGCGGGGACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTa                  <  2:3085895/144‑1 (MQ=255)
                                                                                                                                           ggTGCGGGGCCAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCGTTCCGAATCAAACAACGATACGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATACAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCg           <  2:1510739/150‑1 (MQ=255)
                                                                                                                                           ggTGCGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCg           <  2:2754197/150‑1 (MQ=255)
                                                                                                                                            gTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCgg          >  1:738951/1‑150 (MQ=255)
                                                                                                                                               cggggAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCGGAtc       <  2:2922757/149‑1 (MQ=255)
                                                                                                                                               aGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCGGAtc       >  2:2666928/1‑150 (MQ=255)
                                                                                                                                                 gggAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCGGAtcgt     <  2:1630188/150‑1 (MQ=255)
                                                                                                                                                  ggAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCGGAtcgtc    >  2:3479938/1‑150 (MQ=255)
                                                                                                                                                  ggAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCGGAtcgtc    <  1:3827848/150‑1 (MQ=255)
                                                                                                                                                    aaaaGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCGGATCGTCgg  <  2:3968412/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: