Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,608,695 C→G 12.9% intergenic (+46/+56) 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,6950CG12.9% 1364.6 / 72.2 446intergenic (+46/+56)zntA/tusAzinc, cobalt and lead efflux system/mnm(5)‑s(2)U34‑tRNA 2‑thiolation sulfurtransferase
Reads supporting (aligned to +/- strand):  ref base C (224/158);  new base G (6/56);  total (230/216)
Fisher's exact test for biased strand distribution p-value = 9.59e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 467 total aligned reads displayed.

CCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATA  >  NC_000913/3608546‑3608844
                                                                                                                                                     |                                                                                                                                                     
 cacGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCctct                                                                                                                                                      >  2:2405631/1‑150 (MQ=255)
 cacGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCctct                                                                                                                                                      >  1:2553997/1‑150 (MQ=255)
  acGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTg                                                                                                                                                     <  2:2126707/150‑1 (MQ=255)
   cGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTgg                                                                                                                                                    <  1:615725/150‑1 (MQ=255)
   cGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTgg                                                                                                                                                    <  2:2799486/150‑1 (MQ=255)
   cGCTGTTAGGGATGACCGGGTTGTGGCTGGCAATGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTgg                                                                                                                                                    <  1:7622658/150‑1 (MQ=255)
    gCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTggg                                                                                                                                                   <  2:2717605/150‑1 (MQ=255)
       gtaaGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGgag                                                                                                                                                >  2:8251440/4‑150 (MQ=255)
        ttAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGgaga                                                                                                                                               <  2:3291025/150‑1 (MQ=255)
          tgggATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAgg                                                                                                                                             <  1:3663310/149‑1 (MQ=255)
          aGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCTCCCCAGAGGGGCGAgg                                                                                                                                             <  1:4789583/150‑1 (MQ=255)
            ggATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGt                                                                                                                                           >  2:4506452/1‑150 (MQ=255)
                gACCGGGTGGTGGCGGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAgg                                                                                                                                       <  1:684177/150‑1 (MQ=255)
                  ncGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGGTAgggg                                                                                                                                     >  1:1399728/2‑149 (MQ=255)
                  ccGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGt                                                                                                                                     >  2:3965921/1‑150 (MQ=255)
                    gggTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGa                                                                                                                                   >  1:4543053/1‑150 (MQ=255)
                     ggTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAg                                                                                                                                  <  1:2768516/150‑1 (MQ=255)
                     ggTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAg                                                                                                                                  <  1:1889286/150‑1 (MQ=255)
                     ggTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAg                                                                                                                                  >  1:7665404/1‑150 (MQ=255)
                          tggctggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGTTAAGGCAAACCGATCTCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGTTGAGGGGaa                                                                                                                             <  2:7713033/150‑1 (MQ=255)
                          tggctggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGaa                                                                                                                             <  2:904508/150‑1 (MQ=255)
                           ggctggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAATATTGTTGCGCTGGAGATAAGGCAAACCGATCGCAACTTTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGaaa                                                                                                                            <  1:1583270/150‑1 (MQ=255)
                            gctggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGaaaa                                                                                                                           >  1:1728671/1‑150 (MQ=255)
                            gctggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGaaaa                                                                                                                           <  2:4165613/150‑1 (MQ=255)
                                gcAGTGCTGGCATATACGGGGGCGACGGTTCTTGTGACAGCGAATGCGTTAATATTGTTGCGCAGGAGATAAGGCAAACCGATCTAAACATTTAGCGCTATCGGTCCCCTCGCCCCTCTGGGTAGAGGGTTAGGGTGAGGGGAAAAggcg                                                                                                                       <  2:5590197/150‑1 (MQ=38)
                                    tGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCAt                                                                                                                   >  2:6405273/1‑150 (MQ=255)
                                    tGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCAt                                                                                                                   >  1:5888576/1‑150 (MQ=255)
                                       tGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGa                                                                                                                >  1:425382/1‑150 (MQ=255)
                                       tGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGa                                                                                                                >  1:5469796/1‑150 (MQ=255)
                                         gCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAg                                                                                                              >  1:1995872/1‑150 (MQ=255)
                                         gCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAg                                                                                                              >  1:3215084/1‑150 (MQ=255)
                                          ctgATAGGGGGGGGTGGGTGGTGGTGCCAGCGAATGCGTTAAGATTGTTGCGCCGGTTTTATGGCAAACCGATCGCACCATTGAGGGCGTTCGGTCCCCTCGCCCATCTGGGGTGTGGTTTAGGGTTAGGGGAAAAGGCGGCATCGAAGc                                                                                                             <  2:6032067/148‑1 (MQ=255)
                                           aGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAgg                                                                                                                                 >  1:4401492/1‑129 (MQ=255)
                                           aGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGcc                                                                                                            >  1:4117146/1‑150 (MQ=255)
                                                 gggggCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAg                                                                                                      >  1:7402541/1‑150 (MQ=255)
                                                   gggCGACGGTGCGGGTGACAGCGAATGCGTTATGATTGTTGCGCAGGATATAAGGCACACCGATCGCAACATTGCGCGCTATCGGCCCCCTCGCCCCTCTGGGGAGTGGGTTTGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGcc                                                                                                    <  2:182508/150‑1 (MQ=38)
                                                    ggCGCCGGTGCTGGTGACAGCGAATGCGTTATGATTTTTGCGCTGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGccc                                                                                                   <  2:5164167/150‑1 (MQ=255)
                                                             gcaggTGACAGCGAATGCGTTAAGATTGTTTCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAAcc                                                                                          >  1:7211623/4‑150 (MQ=255)
                                                              cTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCg                                                                                         >  1:3728262/1‑150 (MQ=255)
                                                               tGGTGACTGCTTTTGGGTTATGATTGTGGCGCTGGTGTTCTGGCAACCGGCTGGCAATTTTGCGGGCGCTGGGTCCCCTGGCCCGTGGGGGGTGTGGGTTTGGGTGTGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGc                                                                                        <  2:5050004/150‑1 (MQ=255)
                                                               tGGTGACAGCGAATGCGTTAAGATTGTTGCGCCGGATATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGTGAGGTTTCGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGc                                                                                        <  2:8123410/150‑1 (MQ=39)
                                                                ggTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGcc                                                                                       >  1:6280964/1‑150 (MQ=255)
                                                                 gTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCt                                                                                      >  2:2570169/1‑150 (MQ=255)
                                                                  tGACCGCGATGGCGTTACGTTTGTGGCGCAGGAGATACGGCAACCCGATCGCAACTTTGGGCGCGATGGGGCCCCTCGCCCCTCTGGGGAGTGGGTTTGGGTGCGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCtt                                                                                     <  2:5886723/150‑1 (MQ=255)
                                                                  tGACCGCGAATGCGTTAAGTGTGTTGCGCAGGCTAAACGGCAAACCGACCGCAACATTGCGGGCGTTCGGTCCCCTCGCCCCTCTGGGTATAGGGTTCGGGTGAGGGGACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCtt                                                                                     <  2:7299465/150‑1 (MQ=255)
                                                                    aCAGCGAATGGTTTAATTTTGTTGCGCTGGTTAAATGCCAACCCTTTCGCAAATTTGAGCGCGTTGGGTCCCCTGGCCCCTCTGGGGTTTGGGTTGGGGTGCGGGTAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTa                                                                                   <  2:3919802/150‑1 (MQ=255)
                                                                       gctaaTGCGTTATGATTGTTGCGCAGGATATATGGCAAACCGATCGAAAATTTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGTGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGa                                                                                <  2:1109745/147‑1 (MQ=38)
                                                                        cGACTGCGTTAAGATTGTTGCGCAGGAGATACGGCAACCCGATCGCAACTTTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGCGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGaa                                                                               <  2:4259465/150‑1 (MQ=39)
                                                                         gAATGCTTTAAGATTGTTGCGCAGGATATAAGGCAAACCGATGGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGTGAGGGTTAGGGTGTGGGTACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAAt                                                                              <  2:4821505/150‑1 (MQ=38)
                                                                          tttgggTTAATTTTGTTGCTCAGTATATACGCAACACCTATCCAAAATTTGCGCGCGAGGGGCCCCCCCGCCCGTGGGGGTTGTGGTTTTGGGTGCGGGCAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATc                                                                             <  1:7020223/145‑1 (MQ=255)
                                                                          aaTGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAACCAGCCCCTATCAACCGCCTTTACGAATc                                                                             >  1:8251862/1‑150 (MQ=255)
                                                                               gTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCACCATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATa                                                                        >  1:3945258/1‑150 (MQ=255)
                                                                               gTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATa                                                                        >  1:7121253/1‑150 (MQ=255)
                                                                                ttatgtttgtggCGCAGTAGATACGGAAAACCGATCGCAACTTTGCGCGCGATGGGTCCCCTCGCCCCTGTGGGGTGTGGGTTTGGGTGCGGGGACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATaa                                                                       <  2:5770710/144‑1 (MQ=255)
                                                                                     aTTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATa                                                                  >  2:6661276/1‑150 (MQ=255)
                                                                                      ttgttgCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATaa                                                                 >  1:5687850/1‑150 (MQ=255)
                                                                                      ttgttgCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATaa                                                                 >  1:5405375/1‑150 (MQ=255)
                                                                                       tgttgCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAg                                                                >  2:6902762/1‑150 (MQ=255)
                                                                                        gttgCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAgg                                                               >  1:8162721/1‑150 (MQ=255)
                                                                                                  agaTACGGCAACCCGATCGCAACTTTGTGCGCGATCGGGCCCCTCGCCCCTGGGGGGAGTGGGTTGGGGTGCGGGGACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATCCGGCAGTCCATcc                                                     <  2:4050613/150‑1 (MQ=255)
                                                                                                  agaTAAGGAAAACCGATCGAAACATTGAGCGCGATCGGCCCCCGCGCCCCTCGGGGGTGTGGGTTTGGGTGCGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATcc                                                     <  1:5991975/150‑1 (MQ=255)
                                                                                                      aaGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGtct                                                 >  1:3847338/1‑150 (MQ=255)
                                                                                                      aaGGAAAACCGATCGCAACATTGCGCGCTAGGGGTCCCCGCGCCCCTGGGGGGTGTGGTTTTGGGTGCGGGGACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGtct                                                 <  2:3044704/150‑1 (MQ=255)
                                                                                                        gccAAACCGACCGCAACTTTGAGCGCTCTGGGCCCCCTCGCCCCTGTGGGTTGTGGTTTCGGGTGCGGGGACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTctct                                               <  1:1053719/148‑1 (MQ=255)
                                                                                                           aaccctaTCTCAACTTTGGGGGCGTTGGGTCCCCTCGCCCGTGTGGGGTTTGGTTTTGGGTGCGGGGACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTa                                            <  1:1510368/144‑1 (MQ=255)
                                                                                                           aaaCCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTTCGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTa                                            >  1:1814447/1‑150 (MQ=255)
                                                                                                           aaaCCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTa                                            >  1:1119569/1‑150 (MQ=255)
                                                                                                           aaaCCGATCGCAACATGGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGTGAGGGTTAGGGGGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTa                                            <  1:5027597/150‑1 (MQ=255)
                                                                                                            aaCCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAg                                           >  2:8213115/1‑150 (MQ=255)
                                                                                                            aaCCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGCCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCATTAg                                           >  2:7403331/1‑150 (MQ=255)
                                                                                                              ccGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCa                                         >  1:491879/1‑150 (MQ=255)
                                                                                                                gATCGCACCTTTGCGCGCGATGGGCCCCCGCGCCCCTGTGGGGTGTGGTTTGGGGTGTGGTCACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAAc                                       <  1:2871315/150‑1 (MQ=255)
                                                                                                                gATCGCAACTTTGCGCGCGTGGGGCCCCCGCGCCCCTGTGGGGTGTGGGTTTGGGTGCGGGGACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAAc                                       <  1:6587105/150‑1 (MQ=255)
                                                                                                                gATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTGTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAAc                                       <  2:5491918/150‑1 (MQ=255)
                                                                                                                   cgatataTTGAGCGGGATGGGTCCCCTCGGCCCTCTGTGGTTTGGTTTAGTGTTTGGGTAAAAGGCGGCATCGAAGCATATCAGCCCCTATCAACCGCCTTTACGAATCTAATAACGTTAAGGCAGTCCATCCGTCTCTTTAGCAACCAg                                    <  1:3438880/144‑1 (MQ=255)
                                                                                                                   cgaaaCTTTGAGCGCGCTCGGTCCCCCGGCCCCTGTGGGTAGTGGGTTTGGGTGTGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAg                                    <  2:338548/147‑1 (MQ=255)
                                                                                                                   cGCAACTGTGCGCGCTTTCGGCCCCCTCGCCCCTCGGGGGTGAGGGTTCGGGTGAGGGTAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAg                                    <  1:6480723/150‑1 (MQ=255)
                                                                                                                    gCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGt                                   >  1:2980077/1‑150 (MQ=255)
                                                                                                                     cAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGTGTGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGtt                                  <  2:2023786/150‑1 (MQ=255)
                                                                                                                     cAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGtt                                  >  2:5280509/1‑150 (MQ=255)
                                                                                                                      aaCATTGCGCGCGACGGGTCCCACCGCCCATTTGGGTTTTGGGTTCGGTTGCGGGTAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTAAGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTc                                 <  2:56740/150‑1 (MQ=255)
                                                                                                                         aTTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCgtg                              >  2:7028341/1‑150 (MQ=255)
                                                                                                                         aTTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCgtg                              >  2:3432259/1‑150 (MQ=255)
                                                                                                                          ttGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCgtgt                             >  1:911730/1‑150 (MQ=255)
                                                                                                                          ttGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCgtgt                             >  1:5270575/1‑150 (MQ=255)
                                                                                                                           tGAGCGCGATCGGCCCCCTCGCCCCTCTGGGGAGTGGGTTAGGGTGCGGGCAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGtt                            <  2:1493261/150‑1 (MQ=255)
                                                                                                                               cgcgATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCAt                        >  1:3544459/1‑150 (MQ=255)
                                                                                                                                  gagcGGTCCCCTCGCCCCGCGGGGGAGCGGGTTTGGGTGAGGGGACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATaaa                     <  2:1629583/147‑1 (MQ=255)
                                                                                                                                  gATCGGTCCCCTCGCCCCTCGGGGGTGTGGGTTCGGGTGTGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATaaa                     <  2:7857313/150‑1 (MQ=255)
                                                                                                                                   tgggtcccccTCGCCCCTGTGGGTTGTGGTTTCGGGTGCGGGCACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAg                    <  2:7665404/144‑1 (MQ=255)
                                                                                                                                    tCGGTCCCCTCGCCCCTCTGGGGAGTGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAgg                   <  1:4827815/150‑1 (MQ=255)
                                                                                                                                     cctaacccTCTCCCCAGAGGGGCGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGt                  >  2:3664230/6‑150 (MQ=255)
                                                                                                                                      ctaacccTCTCCCCAGAGGGGCGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTa                 >  1:166628/5‑150 (MQ=255)
                                                                                                                                           ccTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACaaaa            >  2:3562002/1‑150 (MQ=255)
                                                                                                                                           ccTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACaaaa            >  2:6574596/1‑150 (MQ=255)
                                                                                                                                            cTCGCACCTCTGGGGAGAGGGATAGGGTGAGGGGAAAAGTCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGATCGTGTTCCATCAAGGTACAAAAc           >  2:124263/1‑150 (MQ=255)
                                                                                                                                               gCCCCTGTGGGGTGTGGTTTTGGGTGTGGGCACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCa        <  2:6280964/150‑1 (MQ=255)
                                                                                                                                               gCCCCTCTGGGGTGTGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCa        <  2:4676336/150‑1 (MQ=255)
                                                                                                                                               gCCCCTCTGGGGAGTGGGTTTGGGTGTGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCa        <  2:8018862/150‑1 (MQ=255)
                                                                                                                                                     ctGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAata  >  1:4181268/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATA  >  NC_000913/3608546‑3608844

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