Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,608,696 T→G 5.6% intergenic (+47/+55) 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,6960TG 5.6% 720.6 / 4.5 244intergenic (+47/+55)zntA/tusAzinc, cobalt and lead efflux system/mnm(5)‑s(2)U34‑tRNA 2‑thiolation sulfurtransferase
Reads supporting (aligned to +/- strand):  ref base T (124/97);  new base G (3/15);  total (130/114)
Fisher's exact test for biased strand distribution p-value = 2.24e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 258 total aligned reads displayed.

CCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGAT‑CGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATAT  >  NC_000913/3608546‑3608845
                                                                                                                                                       |                                                                                                                                                     
 cacGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGAT‑CGCAACATTGAGCGCGATCGGTCCCCTCGCCCctct                                                                                                                                                       >  2:685940/1‑150 (MQ=255)
 cacGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGAT‑CGCAACATTGAGCGCGATCGGTCCCCTCGCCCctct                                                                                                                                                       >  1:473435/1‑150 (MQ=255)
 cacGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGAT‑CGCAACATTGAGCGCGATCGGTCCCCTCGCCCctct                                                                                                                                                       >  2:2539257/1‑150 (MQ=255)
 cacGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGAT‑CGCAACATTGAGCGCGATCGGTCCCCTCGCCCctct                                                                                                                                                       >  2:2539301/1‑150 (MQ=255)
 cacGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGAT‑CGCAACATTGAGCGCGATCGGTCCCCTCGCCCctct                                                                                                                                                       >  2:3267528/1‑150 (MQ=255)
  acGCTGTTAGGGATGCCCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGAT‑CGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTg                                                                                                                                                      <  1:3525853/150‑1 (MQ=255)
  acGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGAT‑CGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTg                                                                                                                                                      >  2:2181948/1‑150 (MQ=255)
    gCTGTTAGGGATGACCGGGGTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGTTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGAG‑CGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTggg                                                                                                                                                    <  1:36850/150‑1 (MQ=255)
      tGTTAGGGATGCCCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGAT‑CGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGga                                                                                                                                                  <  1:2460153/150‑1 (MQ=255)
      tGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGAT‑CGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGga                                                                                                                                                  <  1:3099584/150‑1 (MQ=255)
      tGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGAT‑CGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGga                                                                                                                                                  >  1:951297/1‑150 (MQ=255)
          aGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGAT‑CGCAACATTGAGCGCGATCGGTCCCCTCTCCCCAGAGGGGCGAgg                                                                                                                                              <  2:159116/150‑1 (MQ=255)
          aGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGAT‑CGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGCGAgg                                                                                                                                              <  2:1292648/150‑1 (MQ=255)
          aGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGAT‑CGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAgg                                                                                                                                              <  1:2224792/150‑1 (MQ=255)
            ggATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGAT‑CGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGt                                                                                                                                            <  1:1620079/150‑1 (MQ=255)
            ggATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGAT‑CGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGt                                                                                                                                            <  1:2939914/150‑1 (MQ=255)
              aTGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGAT‑CGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTa                                                                                                                                          <  1:3821092/150‑1 (MQ=255)
                  ccGGGATGAGGCTTGCAGTGCTGGCAGATACGGGGCCGACGATGCTGGTGACAGCGAATGCGGTAAGATTGTGGCGCGGGAGATAAGGCAAACCGAT‑CCCAACATAGAGCTCGATCGGACCCCACGCCCCGCGGGGGAGCGGGTGAGGGt                                                                                                                                      >  2:3788024/1‑150 (MQ=255)
                   cGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGAT‑CGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTg                                                                                                                                     >  1:3164666/1‑150 (MQ=255)
                     ggTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGAT‑CGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAg                                                                                                                                   <  2:2817763/150‑1 (MQ=255)
                     ggTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGAT‑CGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAg                                                                                                                                   <  2:3917581/150‑1 (MQ=255)
                      gTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGAT‑CGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAgg                                                                                                                                  <  2:2881990/150‑1 (MQ=255)
                         gtgGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGAT‑CGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGa                                                                                                                               >  2:3143179/1‑150 (MQ=255)
                          tggctggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGAT‑CGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGaa                                                                                                                              >  1:2588458/1‑150 (MQ=255)
                            gctggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGAT‑CGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGaaaa                                                                                                                            >  1:553774/1‑150 (MQ=255)
                            gctggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGAT‑CGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGaaaa                                                                                                                            >  1:1785801/1‑150 (MQ=255)
                              tggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGAT‑CGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAgg                                                                                                                          >  1:1499239/1‑150 (MQ=255)
                                gcAGTGCTGGCAGCTACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTTAGATTGTTGCGCAGGAGATAAGGCAAACCGAT‑CGCCACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGTGGAGAGGGTTAGGGTGAGGGGAAAAggcg                                                                                                                        <  2:1732639/150‑1 (MQ=255)
                                  cgTGCTGGCAGATACGGGGGCGACGGTGCGGGTGACCGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAACCCGAT‑CGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTTGGTTGAGGGGAAAAggcggc                                                                                                                      <  2:513736/149‑1 (MQ=39)
                                  aGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGAT‑CGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAggcggc                                                                                                                      >  2:1305493/1‑150 (MQ=255)
                                  aGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGAT‑CGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTGAGGGTGAGGGGAACAggcggc                                                                                                                      >  1:3266289/1‑150 (MQ=255)
                                       tGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGAT‑CGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGa                                                                                                                 >  1:1175971/1‑150 (MQ=255)
                                         gCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTATTTTTGTTGCGATGGAGATAAGGCAAACCGAT‑CGCAACATTGCGCGCTATCGGTCCCCTCGCCCCTCTGGGGAGTGGTTTAGGGTGAGGGGAAAAGGCGGCATCGAAg                                                                                                               <  2:2745669/150‑1 (MQ=38)
                                          cAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGAT‑CGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGc                                                                                                              >  1:17408/1‑150 (MQ=255)
                                            gATCCGGGGGCGACGGTGCGGGTGACCGCGAATGCGTTACGATGGTTGCGCAGGTGATAAGGCAACCCGAT‑CGCAACATTGCGCGCGATCGGTCCCCACGCCCCTCGGGGGAGAGGGTTCGGGTGAGGGGAAAAGGCGGCATCGAAGCCa                                                                                                            <  2:4063244/150‑1 (MQ=255)
                                            gATAGGGGGGCGACGGTGCTGGTGACTGCTTATGCGTTAAGATTGTTGCGCTGTAGATACGGAAAACCGAT‑CGCAACATTGAGCGCTTTCGGTCCCCCCGCCCCTGTGGGGTGAGGGTTAGGGTGTGGGGAAAAGGCGGCATCGAAGCCa                                                                                                            <  2:76040/150‑1 (MQ=255)
                                              tACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGAT‑CGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGG‑AAAGGCGGCATCGAAGCCAATc                                                                                                         >  2:177090/1‑150 (MQ=255)
                                               aCGGGGGCGCCGGTGCGGGTGACCGCGCATGCGTTACGATGGTGGCGCAGGAGATACGGCAAACCGAT‑CGCAACATGGCGCGCGATCGGTCCCCGCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATc                                                                                                         <  2:1995019/150‑1 (MQ=37)
                                                  ggggggAGGGTGCTTGTGACAGCTATTGCTTTAAGTTTTTTGGGCAGTTTATAAGGAAAACCGTT‑CGCAACATTGAGCGCTATCGGTCCCCGCGCCCCTCGGGGTTGTGGTTTTGGGTGTGGGGAAAAGGCGGCATCGAAGCCAATCAGc                                                                                                      <  2:1909315/150‑1 (MQ=255)
                                                    ggCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGAT‑CGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGG‑AAAGGCGGCATCGAAGCCAATCAGcccc                                                                                                   >  2:1000073/1‑150 (MQ=255)
                                                                   gACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGAT‑CGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCttt                                                                                     >  1:2605197/1‑150 (MQ=255)
                                                                    aCAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGAT‑CGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTCTCAACCGCCtttc                                                                                    >  1:3541810/1‑149 (MQ=255)
                                                                     cAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGAT‑CGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTAc                                                                                   >  1:2138073/1‑150 (MQ=255)
                                                                      aGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGAT‑CGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACg                                                                                  >  1:117379/1‑150 (MQ=255)
                                                                        cgtattggTTATGTTTTTTGCGCTGCAGATAAGGC‑AACCGCTGCTCAACATTGAGCGCGTTCGGTCCCATCGCCACTCTGGGGTTTGGGTTAGGGTTAGGGTAAAAGGCGGCATCTAAGCCAATCAGCCCCTATCAACCGCCTTTACGaa                                                                                <  2:2401132/143‑1 (MQ=255)
                                                                        cGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGAT‑CGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGaa                                                                                >  1:1764575/1‑150 (MQ=255)
                                                                               ggtAAGATTGTTGCGCAGGAGATAAGGCAAACCGAT‑CGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATa                                                                         >  1:240331/3‑150 (MQ=255)
                                                                               gTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGAT‑CGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATa                                                                         >  1:1315088/1‑150 (MQ=255)
                                                                               gTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGAT‑CGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATa                                                                         >  1:1962108/1‑150 (MQ=255)
                                                                               gTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGAT‑CGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATa                                                                         >  1:2559324/1‑150 (MQ=255)
                                                                                ttAAGATTGTTGCGCAGGAGATAAGGCAAACCGAT‑CGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATaa                                                                        >  1:1358248/1‑150 (MQ=255)
                                                                                  aaGATTGTTGCGCAGGAGATAAGGCAAACCGAT‑CGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACg                                                                      >  2:3807157/1‑150 (MQ=255)
                                                                                       tgttgCGCAGGAGATAAGGCAAACCGAT‑CGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAg                                                                 >  1:22358/1‑150 (MQ=255)
                                                                                        gttgCGCTGGATATAAGGCAAACCTTT‑CGAAACTTTGTGCGCTATGGGGCCCCTCGCCCCTGTGGGTTGTGGTTTTGGGTGTGGGGACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAgg                                                                <  1:2970625/150‑1 (MQ=255)
                                                                                        gttgCGCAGGAGATAAGGCAAACCGAT‑CGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAgg                                                                >  2:704397/1‑150 (MQ=255)
                                                                                          tgCGCTGGCTATAAGGCAACCGGTT‑CGAAACTTTGTGCGCGTTCGGTCCCCTCGCCCCTCTGGGGAGTGGGTTCGGGGGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCa                                                              <  2:2004218/150‑1 (MQ=37)
                                                                                               aGGAGATAAGGCAAACCGAT‑CGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCa                                                         >  2:64609/1‑150 (MQ=255)
                                                                                                  tgaTAAGGCAAACCGAT‑CGCAACTTTGAGCGCGATGGGTCCCCTCGCCCCTGTGGGTAGTGGTTTTGGGTGCGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATcc                                                      <  1:3062478/149‑1 (MQ=255)
                                                                                                  agaTACGCCAAACCGAG‑CGCAACATTGCGCGCGATGGGTCCCCCCCCCCGTCGGGGGTGCGGGTTTGGGTGGGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATcc                                                      <  1:227247/150‑1 (MQ=255)
                                                                                                        ggCAAACCGAT‑CGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTATACGAATCAAATAACGATAAGGCAGTCCATCCGTctct                                                >  2:3320721/1‑150 (MQ=255)
                                                                                                        ggCAAACCGAT‑CGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGGGAGGGGAAAAGGCGGCATCGAAGCCACTCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTctct                                                >  2:3442936/1‑150 (MQ=255)
                                                                                                         gCAAACCGAT‑CGCAACATGGCGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGCGGGGACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCtt                                               <  2:3686778/150‑1 (MQ=255)
                                                                                                          cAAACCGAT‑CGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGTGTGGGTTTGGGTGAGGGGACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCttt                                              <  1:2730988/150‑1 (MQ=255)
                                                                                                           aacccGAT‑CCCACCATTGAGGGCGACCGGCCCCCCCGCCCCTCGGGGGAGTGGGTTTGGGTGCGGGTACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATACCGATAAGGCAGTCCATCCGTCCCTTTa                                             <  1:869115/147‑1 (MQ=255)
                                                                                                           aaaCCGAT‑CGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGTGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTa                                             <  2:2638516/150‑1 (MQ=255)
                                                                                                            acccgac‑cgaAACTTGGCGCGCTCTCGGCCCCCTCGCCCGTGTGGGGTGTGGGTTTGGGGCCGGGTACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAg                                            <  1:1435968/148‑1 (MQ=255)
                                                                                                                gAT‑CGCAACTTTGCGCGGGAGGGGTCCCCTCGCCCGTGTGGGGCGTGGGTTTGGGTGTGGGGACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAAc                                        <  1:4102443/150‑1 (MQ=255)
                                                                                                                gAT‑CGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAAc                                        >  2:2766683/1‑150 (MQ=255)
                                                                                                                gAT‑CGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAAc                                        >  2:475188/1‑150 (MQ=255)
                                                                                                                gAT‑CGCAACATTGAGCGCGATCGGTCCCCGCGCCCCTCTGGGGAGTGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAAc                                        <  1:3048929/150‑1 (MQ=255)
                                                                                                                 gcccGCAACTTGGCGCGCGATCGGCCCCCGCGCCCCTGTGGGGTGTGGGTTGGGGGGTGGGCACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAAc                                        <  1:716982/147‑1 (MQ=255)
                                                                                                                  t‑CGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCa                                      >  2:1719395/1‑150 (MQ=255)
                                                                                                                     gCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAAcccgg                                    >  2:31113/1‑147 (MQ=255)
                                                                                                                      cAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGtt                                   >  2:2050393/1‑150 (MQ=255)
                                                                                                                      cAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGtt                                   >  2:4084926/1‑150 (MQ=255)
                                                                                                                       aaCATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTc                                  <  2:2378134/150‑1 (MQ=255)
                                                                                                                       aaCATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTc                                  >  2:609909/1‑150 (MQ=255)
                                                                                                                       aaCATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTc                                  >  2:590072/1‑150 (MQ=255)
                                                                                                                       aaCATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACAAGTTc                                  >  1:1514136/1‑150 (MQ=255)
                                                                                                                         cATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCgt                                >  1:3442948/1‑150 (MQ=255)
                                                                                                                             gAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGG‑AAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTcc                           <  1:1000073/150‑1 (MQ=255)
                                                                                                                              cgcgctttCGGTCCCCCCGCCCCGGGGGGTAGTGTTTTTGGGGGCGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTcc                           <  2:1901496/143‑1 (MQ=255)
                                                                                                                              aGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTcc                           >  2:442822/1‑150 (MQ=255)
                                                                                                                                cgctATCGGTCCCCTCGCCCCTCTGGGTAGTGGGTTTGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCAt                         <  2:514956/150‑1 (MQ=255)
                                                                                                                                 gggATCGGTCCCCGCGCCCGTGGGGGGTGTGGGTTTGGGTGCGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATa                        <  2:1917328/148‑1 (MQ=255)
                                                                                                                                 gcgATCGGTCCCCTCGCCCCTCGGGGGTGCGGTTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATa                        <  2:575900/150‑1 (MQ=255)
                                                                                                                                    aacgcaccccccccccCTCTGGGCACAGGGTAAGGGTCAGCCCAAAACGCGGCATCGAAGCCAAACAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAACGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAg                     <  2:11840/144‑1 (MQ=255)
                                                                                                                                         tCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACa                >  1:964310/1‑150 (MQ=255)
                                                                                                                                          acccTCTCCCCAGAGGGGCGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACaa               >  1:3887920/2‑150 (MQ=255)
                                                                                                                                           cccTCTCCCCAGAGGGGCGAGGGGTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACaaa              >  2:208778/1‑150 (MQ=255)
                                                                                                                                            ccTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACaaaa             >  2:2939949/1‑150 (MQ=255)
                                                                                                                                            ccTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACaaaa             >  2:664963/1‑150 (MQ=255)
                                                                                                                                              tCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAAcc           >  2:2427488/1‑150 (MQ=255)
                                                                                                                                                gCCCCTCTGGGTTGAGGGTTTGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCa         <  1:2594926/150‑1 (MQ=255)
                                                                                                                                                gCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCa         >  1:1643487/1‑150 (MQ=255)
                                                                                                                                                 ccccTCTGGGTAGTGGGTTTGGGTGTGGGTAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAg        <  1:621119/150‑1 (MQ=255)
                                                                                                                                                 ccccTCTGGGGTGTGGGTTTGGGTGAGGGTAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAg        <  2:3547337/150‑1 (MQ=255)
                                                                                                                                                  cccTCGGGGGAGTGGTTTCGGTGGCGGTGACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAgg       <  1:2929557/150‑1 (MQ=255)
                                                                                                                                                     tctGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAAACAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAACGTACAAAACCCAGGaaa    >  2:2659389/1‑149 (MQ=255)
                                                                                                                                                       tGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAatat  >  1:1043977/1‑150 (MQ=255)
                                                                                                                                                       |                                                                                                                                                     
CCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGAT‑CGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATAT  >  NC_000913/3608546‑3608845

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: ‑