Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,608,695 C→G 6.4% 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,6950CG 6.4% 247.1 / 7.7 91intergenic (+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 (38/47);  new base G (4/2);  total (42/49)
Fisher's exact test for biased strand distribution p-value = 4.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.19e-01

CCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGA  >  NC_000913/3608546‑3608841
                                                                                                                                                     |                                                                                                                                                  
ccacGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCctc                                                                                                                                                    >  2:619933/1‑150 (MQ=255)
   cGCTGTTAGGGCTGACCGGGGGGTGGCTGGCAGTGCCGGCAGATACGGGGGCGACGGCGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGAGCGGTCCCCTCGCCCCTCTgg                                                                                                                                                 <  1:292074/150‑1 (MQ=255)
   cGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTgg                                                                                                                                                 >  1:1278864/1‑150 (MQ=255)
    gCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTggg                                                                                                                                                <  2:956769/150‑1 (MQ=255)
    gCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTggg                                                                                                                                                <  2:798846/150‑1 (MQ=255)
     ctggcaGGGAGGACCGGGCTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGGGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTgggg                                                                                                                                               <  1:45450/145‑1 (MQ=255)
     cTGTTCGGGAGGACCGGGGTGTGGCGGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTgggg                                                                                                                                               <  1:40001/150‑1 (MQ=255)
     cTGTTAGGGATGACCGGGTTGTGGCTGGCAGTTCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTgggg                                                                                                                                               <  1:907470/150‑1 (MQ=255)
      tGTTACGGATCACCGGGTTGTGGCTGGAAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGga                                                                                                                                              <  2:848277/150‑1 (MQ=255)
          aGGGATGACCGGGTTGGGGCGGGCCGGGCTGGCAGAAACGGGGGCGACGGGGCTGGTGACAGCTAAAGCGTTAAGATTGTTGCGCAGGATAAAAGGCAAACCGAACGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAgg                                                                                                                                          <  1:1505745/150‑1 (MQ=255)
          aGGGATGACCGGGGTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACTGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAgg                                                                                                                                          <  1:1162407/150‑1 (MQ=255)
            ggATGCCCGGGGGGTGGCTGGCAGTGCGGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCTCCCCAGAGGGGCGAgggg                                                                                                                                        <  1:549012/150‑2 (MQ=255)
             gATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGtt                                                                                                                                       >  1:490487/1‑150 (MQ=255)
                gACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGAAACGGCAAACCGATCGCAACATTGAGCGCGAGCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAgg                                                                                                                                    <  1:1291739/150‑1 (MQ=255)
                 cccGGGTTGTGGCTGGCCGTGCTGGCAGAGACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTACGCTGGATGCGCAGGAGATACGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAggg                                                                                                                                   <  1:85610/149‑1 (MQ=38)
                  ccGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGTGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGt                                                                                                                                  >  1:1626226/1‑150 (MQ=255)
                  ccGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTgggg                                                                                                                                               >  1:642973/1‑137 (MQ=255)
                  ccGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTgggg                                                                                                                                               <  2:642973/137‑1 (MQ=255)
                         gtgGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGa                                                                                                                           >  2:1451139/1‑150 (MQ=255)
                           ggctggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGaaa                                                                                                                         <  2:679004/150‑1 (MQ=255)
                               ggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACAATGAGCGCGATCGGTCCCCTCGCCCCTCAGGGGAGAGGGATAGGCTGAGGGGAAAAggc                                                                                                                     >  2:1002107/1‑150 (MQ=255)
                                gcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAggcg                                                                                                                    <  1:619933/150‑1 (MQ=255)
                                 cAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGACCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAggcgg                                                                                                                   <  1:63404/150‑1 (MQ=255)
                                         gCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAg                                                                                                           <  1:1218828/150‑1 (MQ=255)
                                          cAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTTCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGACATCGAAGc                                                                                                          >  1:250295/1‑150 (MQ=255)
                                          cAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGc                                                                                                          <  2:490487/150‑1 (MQ=255)
                                           aGATACGGGGGCGCCGGTGCTGGGGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCTACATTGAGCGCGATCGGTCCCCTCGCCCCGCGGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGcc                                                                                                         <  1:931752/150‑1 (MQ=39)
                                           aGATACGGGGGCGACGGTGCTGGTGACTGCGAATCGGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATGGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCTTCGAAGcc                                                                                                         <  2:705573/150‑1 (MQ=38)
                                                  ggggCGACGGTGCAGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCACGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGc                                                                                                  >  1:1641995/1‑150 (MQ=255)
                                                    ggCGCCGGTGCTGGTGACAGCGAATGCGGTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATGGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGTGAGGGTGAGGGTGCGGGGAAAAGGCGGCATGGAGGCGAATCAGccc                                                                                                >  2:672405/1‑150 (MQ=38)
                                                    ggCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGccc                                                                                                >  1:1648683/1‑150 (MQ=255)
                                                    ggCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGccc                                                                                                >  1:1263465/1‑150 (MQ=255)
                                                     gCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCCACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGcccc                                                                                               <  2:275216/150‑1 (MQ=255)
                                                             gcgggggACAGCGAATGCGTTAATATTGGTGCGCCGGATATAAGGAAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAAcc                                                                                       <  2:1278864/144‑1 (MQ=38)
                                                             gCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAAcc                                                                                       <  2:1182026/150‑1 (MQ=255)
                                                              cTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCg                                                                                      <  2:971655/150‑1 (MQ=255)
                                                              cTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCg                                                                                      >  1:1079504/1‑150 (MQ=255)
                                                                ggTGACAGCGAATGCTTTAAGATTGTTGCGCTGTAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGTGTGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGcc                                                                                    <  2:591466/150‑1 (MQ=39)
                                                                  tGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGCGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCtt                                                                                  >  1:1528364/1‑150 (MQ=255)
                                                                   gACAGCGAATGCGTTAAGATGGTGGCGCAGGTGATAAGCCAACCCGATCGCAACATTCAGCGCGATCGGTCCCCTCGCCCCGCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCttt                                                                                 <  1:391648/150‑1 (MQ=38)
                                                                    aCAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGTGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACAGCCTTTa                                                                                >  2:1515393/1‑150 (MQ=255)
                                                                    aCAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTa                                                                                <  2:1626226/150‑1 (MQ=255)
                                                                    aCAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCACGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAAACAGCCCCTATCAACAGCCtttc                                                                                >  2:1521938/1‑149 (MQ=255)
                                                                     cAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTAc                                                                               >  2:628282/1‑150 (MQ=255)
                                                                     cAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGAAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGTTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCGTATAAACCGCATTTAc                                                                               >  2:836696/1‑150 (MQ=39)
                                                                      aGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACg                                                                              <  1:604294/150‑1 (MQ=255)
                                                                       gCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGa                                                                             >  2:1188446/1‑150 (MQ=255)
                                                                        cGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGaa                                                                            >  1:1296701/1‑150 (MQ=255)
                                                                         gAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAAt                                                                           >  1:606264/1‑150 (MQ=255)
                                                                              cGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGAGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGATAGGGTGAGGGGAAAAGGCGGCATCGAAGCAAATCAGCCCCTATCAACCGCCTTTACGAATCaaaa                                                                      >  1:29642/1‑149 (MQ=255)
                                                                                   agagTGTTCCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGa                                                                 <  1:1678689/150‑1 (MQ=255)
                                                                                     aTTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATa                                                               >  2:585047/1‑150 (MQ=255)
                                                                                     aTTGTGGCGCCGGAGTAACGGCAAACCGAGCGCAACATTGCGCGCGAGCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGGGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATa                                                               <  2:1296701/150‑1 (MQ=38)
                                                                                            cgcAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGt                                                        >  2:1157300/1‑150 (MQ=255)
                                                                                              cAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTcc                                                      <  2:1185419/150‑1 (MQ=255)
                                                                                                ggAGATAAGGCAAACCGATCGCAACATTGAGCGGGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTTCGAATAAAATAAAGATAAGGCAGTCAAt                                                    >  1:1630040/1‑150 (MQ=38)
                                                                                                   taaaaGGCAACCCGATCGCACCATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGGGAGTGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTAGCAACCGCCTTTACGAATCAAATAACGATAAGGCAGGCCATCCg                                                 <  2:125824/147‑1 (MQ=38)
                                                                                                        ggCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATAAAATAACGATAAGGCAGTCCATCCggatct                                            >  2:723681/1‑145 (MQ=255)
                                                                                                         gCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCtt                                           >  2:1621718/1‑150 (MQ=255)
                                                                                                           aaaCCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGTGGAGAGGGTTAGGGTGTGGGGAAAAGGCGGCATCGATGCCAATCAGCCCCTATCAACCGCCTTTACTAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTa                                         <  1:1002107/150‑1 (MQ=255)
                                                                                                           aaaCCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTa                                         >  1:1454376/1‑150 (MQ=255)
                                                                                                            aaCCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAg                                        <  2:1641995/150‑1 (MQ=255)
                                                                                                                gATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGACATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAAc                                    <  2:250295/150‑1 (MQ=255)
                                                                                                                 ttCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAAcc                                   <  1:1188446/149‑1 (MQ=255)
                                                                                                                     cAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCACTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGtt                               >  2:460173/1‑150 (MQ=255)
                                                                                                                       aCATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCg                             <  2:903673/150‑1 (MQ=255)
                                                                                                                        aaTTGTGCGCGATCGGTCCCCTCGCCCCTGTGGGGAGAGGGTTAGGGTGAGGGTAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCgt                            <  1:836696/149‑1 (MQ=255)
                                                                                                                         aTTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCgtg                           <  2:1690686/150‑1 (MQ=255)
                                                                                                                           ttaGCGCGATGGGTCCCCTCGCCCCTCTGGGTAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGtt                         <  2:1263465/148‑1 (MQ=255)
                                                                                                                           tGAGCGCGATGTGTCCCCTCGCCCCTCTTGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGtt                         <  2:1648683/150‑1 (MQ=255)
                                                                                                                            gAGCGCGATCGGTCCCCTCGCCACTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCTGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATATCGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTc                        <  2:1008207/150‑1 (MQ=255)
                                                                                                                                 cgATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATTGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATaa                   >  2:993130/1‑150 (MQ=255)
                                                                                                                                 cgATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCGGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATaa                   >  2:242429/1‑150 (MQ=255)
                                                                                                                                 cgATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATaa                   >  1:168050/1‑150 (MQ=255)
                                                                                                                                 cgATCGGTCCCCTCGCCCATCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTAAGAAACAAAAAACGATAAGGCAGTCCAACCGTCTCTTTAGCAACCAGTTCGTGTTCCATaa                   >  2:936443/1‑150 (MQ=39)
                                                                                                                                   atgggTCCCCTCGCCCCGCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAg                 <  1:580727/147‑1 (MQ=255)
                                                                                                                                    tCGGTCCCCTCCCCCCTCTGGGGATAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAgg                <  2:375946/150‑1 (MQ=255)
                                                                                                                                     cctaacccTCTCCCCAGAGGGGCGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCAACCGTCTCTTTAGCAACCAGTTCGAGTTCCATAAAGGt               >  2:1614210/6‑150 (MQ=255)
                                                                                                                                       taacccTCTCCCCAGAGGGGCGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTAc             >  2:1473046/4‑150 (MQ=255)
                                                                                                                                       gTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGCAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTAc             <  2:1079504/150‑1 (MQ=255)
                                                                                                                                       gTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCACCCAGTTCGGGTTCCATAAAGGTAc             >  2:1175205/1‑150 (MQ=255)
                                                                                                                                       gTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTAc             <  1:585047/150‑1 (MQ=255)
                                                                                                                                        tCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACa            <  1:1451139/150‑1 (MQ=255)
                                                                                                                                        aacccTCTCCCCAGAGGGGCGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGACTTTACTAATCAAATAACGATAAGGCAGTCCATCCGACTCTTTAGCAACAAGTTCGTGTTCCATAAAGGTACa            >  2:1291857/3‑150 (MQ=255)
                                                                                                                                         acccTCTCCCCAGAGGGGCGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACaa           >  2:419427/2‑150 (MQ=255)
                                                                                                                                           ccTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCAGCCGTCTCTTTAGCAACCAGTTCGTGTTCCTTAAAGGTACaaaa         >  1:393172/1‑150 (MQ=255)
                                                                                                                                              cGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAAccc      <  2:606264/150‑1 (MQ=255)
                                                                                                                                              cGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAAccc      >  2:474264/1‑150 (MQ=255)
                                                                                                                                                ccccTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAg    >  1:544426/1‑150 (MQ=255)
                                                                                                                                                ccccTCTGCGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGAAAAGGCAGTCCATCCGACTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCac    >  2:633370/1‑149 (MQ=255)
                                                                                                                                                  ccTCTGGGGAGTGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGa  <  2:523198/150‑1 (MQ=255)
                                                                                                                                                  ccTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGa  >  1:474317/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
CCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGA  >  NC_000913/3608546‑3608841

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: