Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,378,977 T→A 6.0% intergenic (‑33/+30) ampC ← / ← frdD penicillin‑binding protein; beta‑lactamase, intrinsically weak/fumarate reductase (anaerobic), membrane anchor subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,378,9770TA 6.0% 789.0 / 23.2 242intergenic (‑33/+30)ampC/frdDpenicillin‑binding protein; beta‑lactamase, intrinsically weak/fumarate reductase (anaerobic), membrane anchor subunit
Reads supporting (aligned to +/- strand):  ref base T (95/130);  new base A (11/6);  total (106/136)
Fisher's exact test for biased strand distribution p-value = 8.09e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 242 total aligned reads displayed.

TTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGATCGTGCATCGCGTGGTGCATACGGT  >  NC_000913/4378831‑4379125
                                                                                                                                                  |                                                                                                                                                    
ttgCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAgg                                                                                                                                                   >  1:3624605/1‑150 (MQ=255)
ttgCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAgg                                                                                                                                                   >  1:2318690/1‑150 (MQ=255)
 tgCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGc                                                                                                                                                  >  2:1160181/1‑150 (MQ=255)
  gCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCg                                                                                                                                                 <  1:2690258/150‑1 (MQ=255)
   ctctATATGGGGGTTATTTGTGCGATGCACAATATCTTTGATTTGTTGAGGGGCAGAAAATGTGGAGCAAGAGGCTTTAATTAATACGGCGCAGAGCGTCGTTTTGAACATCGGGTCTGGTTTCAATACAAAACGGCCCGCCATAGGCgg                                                                                                                                                <  1:1175125/150‑1 (MQ=255)
   ctctATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCgg                                                                                                                                                <  2:2674976/150‑1 (MQ=255)
   ctctATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCgg                                                                                                                                                <  2:2094858/150‑1 (MQ=255)
    tctATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCggg                                                                                                                                               >  2:2484875/1‑150 (MQ=255)
      tataAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGcc                                                                                                                                             >  2:3400517/1‑150 (MQ=255)
      tataAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGcc                                                                                                                                             >  2:1593782/1‑150 (MQ=255)
      tataAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTGTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGcc                                                                                                                                             >  2:184364/1‑150 (MQ=255)
       ataAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCg                                                                                                                                            <  2:1701151/150‑1 (MQ=255)
       ataAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCg                                                                                                                                            <  2:1670748/150‑1 (MQ=255)
         aaGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCTATGGCGGGCCgtt                                                                                                                                          <  2:1085387/150‑3 (MQ=255)
         aaGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGa                                                                                                                                          <  2:2339048/150‑1 (MQ=255)
          aGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCAGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGAt                                                                                                                                         >  2:1105194/1‑150 (MQ=255)
           gCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGAtt                                                                                                                                        >  1:2211518/1‑150 (MQ=255)
             ggggTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTa                                                                                                                                      <  2:2406935/150‑1 (MQ=255)
             ggggTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTa                                                                                                                                      <  2:1400717/150‑1 (MQ=255)
               cgataTTGTGCATCGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACa                                                                                                                                    >  2:3726585/5‑150 (MQ=255)
                     tgtgCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGGACTTAATAAGGCGCAGAGCGTCGGTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCg                                                                                                                              >  2:2192879/1‑150 (MQ=255)
                      gtgCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGa                                                                                                                             <  2:1617704/150‑1 (MQ=255)
                      gtgCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGAATTACATTGGCGa                                                                                                                             <  2:2046907/150‑1 (MQ=255)
                       tgCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTGTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGAt                                                                                                                            >  2:1638360/1‑150 (MQ=255)
                         cGATGCAAAATATCGTTGATTTGTTGTGGGGCAGCAATTGTGGAGCAAGTGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATCGGCGGGCCGGATTTACATTGGCGATGc                                                                                                                          <  1:767758/150‑1 (MQ=255)
                               acaATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTCATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGa                                                                                                                    >  2:3577848/1‑150 (MQ=255)
                               acaATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGa                                                                                                                    >  1:2033950/1‑150 (MQ=255)
                               acaATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGa                                                                                                                    >  1:1764401/1‑150 (MQ=255)
                                aaaTATCGTTTTTTTGTTGCGGGGCAGCAAATGTGGAGCAAGTGGCGTTAATTAATACGGCGCTGAGCGTCTTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGAt                                                                                                                   <  1:1394064/149‑1 (MQ=255)
                                 aaTATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGAtt                                                                                                                  >  2:959414/1‑150 (MQ=255)
                                 aaTATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTGTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGAtt                                                                                                                  >  1:2049219/1‑150 (MQ=255)
                                   tatCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGt                                                                                                                >  1:3847903/1‑150 (MQ=255)
                                     tCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTaa                                                                                                              >  2:3367870/1‑150 (MQ=255)
                                     tCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTaa                                                                                                              >  2:375712/1‑150 (MQ=255)
                                      cGTTGATTTGTTGAGGGGCAGCAAATTTGCACCAAGACGCGTTAATTAAAAAGGCGCAGAGCGTAGTTTTCAACATAGGCACTGGTTACAAAACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAAc                                                                                                             <  1:1789245/150‑1 (MQ=255)
                                        ttGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGa                                                                                                           >  1:1674348/1‑150 (MQ=255)
                                        ttGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGa                                                                                                           >  1:1109571/1‑150 (MQ=255)
                                        ttGATTTGGTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGa                                                                                                           >  1:3186987/1‑150 (MQ=255)
                                               gttgAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAAt                                                                                                    <  2:2770122/150‑1 (MQ=255)
                                               gttgAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAAt                                                                                                    >  2:1023062/1‑150 (MQ=255)
                                                ttgAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATc                                                                                                   >  1:2871486/1‑150 (MQ=255)
                                                    ggggCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCg                                                                                               >  2:3477685/1‑150 (MQ=255)
                                                     gggCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGt                                                                                              >  1:3353282/1‑150 (MQ=255)
                                                      ggCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTg                                                                                             >  2:3100396/1‑150 (MQ=255)
                                                       gcagcaAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGa                                                                                            >  2:2745276/1‑150 (MQ=255)
                                                       gcagcaAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGa                                                                                            >  2:3812299/1‑150 (MQ=255)
                                                              aTGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGt                                                                                     <  2:800460/150‑1 (MQ=255)
                                                                gtgGAGCAAGCGGCGGTAATTAATACGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCa                                                                                   <  1:799644/150‑1 (MQ=255)
                                                                gtgGAGCAAGAGGCGTTAATTAATACGGCGCAGAGCGTCGTTTTGAACATTGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCa                                                                                   <  1:1568954/150‑1 (MQ=255)
                                                                 tgGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAg                                                                                  >  2:1880889/1‑150 (MQ=255)
                                                                  ggAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTGTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCCTAGGCGGGCCGGATTTACATTTGCGATGCGTTAGATTGTAACGACACCAAACAGCGTGACAACTGTCAgg                                                                                 >  1:1815694/1‑150 (MQ=255)
                                                                   gAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGa                                                                                >  1:2153205/1‑150 (MQ=255)
                                                                   gAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGa                                                                                <  2:2225651/150‑1 (MQ=255)
                                                                     gCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTCTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATa                                                                              <  2:2304478/150‑1 (MQ=255)
                                                                        agagGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATagca                                                                           <  2:2442319/150‑1 (MQ=255)
                                                                         gagGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACAAAGGGTCTGGTTTCCATACAAAAAGGCCCCCCAAAGGCGGGCCGGATTTAAATTTGCGATTCGTTAGATTGTAACGACACCAAAAAGCGTGACAACAGTCAGGATagcag                                                                          >  1:273388/1‑150 (MQ=255)
                                                                           ggCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGcc                                                                        <  1:2649913/150‑1 (MQ=255)
                                                                            gCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCa                                                                       <  2:3352934/150‑1 (MQ=255)
                                                                            gCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAAAGGCCCCCCAAAGGCGGGCCGGATTTACAATGGCGATGCGTTAGATTGTAAAGACACCAAACAGCGTGACAACAGTAAGGATAGCAGccc                                                                       >  1:2113033/1‑149 (MQ=255)
                                                                              ggTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATCGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGa                                                                     <  2:1805886/150‑1 (MQ=255)
                                                                              ggTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGa                                                                     <  1:2484875/150‑1 (MQ=255)
                                                                               gTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAATCGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGAc                                                                    >  2:562943/1‑150 (MQ=255)
                                                                                 aattaatAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCg                                                                  <  1:1514825/150‑1 (MQ=255)
                                                                                 aattaatAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCg                                                                  <  1:643427/150‑1 (MQ=255)
                                                                                  attaatAAGGCGCAGAGCGTCGTTTTGAATATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGt                                                                 <  2:3662789/150‑1 (MQ=255)
                                                                                   ttaatAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTa                                                                <  1:3119547/150‑1 (MQ=255)
                                                                                     aataaGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGa                                                              <  1:2247651/150‑1 (MQ=255)
                                                                                      ataaGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGACCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGaa                                                             >  2:1422185/1‑150 (MQ=255)
                                                                                       taaGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGaaa                                                            <  2:1264931/150‑1 (MQ=255)
                                                                                        aaGGCGCAGAGCGTCGTTTTGACCATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGaaaa                                                           <  2:2500540/150‑1 (MQ=255)
                                                                                        aaGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGaaaa                                                           <  1:3959925/150‑1 (MQ=255)
                                                                                        aaGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGaaaa                                                           <  1:3726585/150‑1 (MQ=255)
                                                                                          gtcgcAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGTATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAAcc                                                         <  1:767779/148‑1 (MQ=255)
                                                                                           gcgcAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAAccc                                                        <  2:4099908/150‑1 (MQ=255)
                                                                                                     tcgtTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGc                                              >  2:3028296/1‑150 (MQ=255)
                                                                                                     tcgtTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACAATTGCGATGCGTTATATTGTAAAGAAACCAATCAGCGTGACAAATGTAAGGATAGCAGCCAGACCGTAGAAAACCCATTTTccccc                                              >  1:3196491/1‑148 (MQ=255)
                                                                                                      cgtTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCa                                             <  1:2598325/150‑1 (MQ=255)
                                                                                                                 tAGGGTCTGGTTTCCATACAAAACGGCCCGCCTTAGGCGTGCCGGATTTACATTTGAGATGCGTTAGATTTTAACGACACCAATCAGCGTGACAACTGTCAGGTTAGCAGCCAGACCGTAGAATACCCATTTGCCCGCAGGTACGTGGAt                                  <  1:1000209/150‑1 (MQ=255)
                                                                                                                   gggTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGAttt                                >  2:2936015/1‑150 (MQ=255)
                                                                                                                     gTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACCACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCAATTGCCCGCAGGTACGTGGAtttta                              >  1:1026375/1‑149 (MQ=255)
                                                                                                                     gTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTc                              <  1:3744719/150‑1 (MQ=255)
                                                                                                                     gTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTc                              <  1:2990396/150‑1 (MQ=255)
                                                                                                                      tCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCa                             <  1:3166272/150‑1 (MQ=255)
                                                                                                                        tGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGa                           <  1:502618/150‑1 (MQ=255)
                                                                                                                        tGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGa                           <  1:3028296/150‑1 (MQ=255)
                                                                                                                            ttCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGATCGt                       >  2:2448986/1‑150 (MQ=255)
                                                                                                                            ttCCATACAAAACGGCCCCCCATAGGCGGGCCGGATTTAAATTGGCGATGCGTTAGATTGTAACGACACCAAACAGCGTGACAACAGTCAGGATAGCAGCCAGACCGTAGAAAAACCAATTGCCCGCAGGTACGTGGATTTTCAGATCGt                       >  1:876740/1‑150 (MQ=255)
                                                                                                                                    aaaaCGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGATCGTGCATcgcg               >  2:570026/1‑150 (MQ=255)
                                                                                                                                    aaaaCGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGATCGTGCATcgcg               >  2:4019062/1‑150 (MQ=255)
                                                                                                                                     aaaCGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGATCGTGCATCGCgt              >  2:796055/1‑150 (MQ=255)
                                                                                                                                      aaCGGCCCCCCATAGGCGGGCCGGATTTACATTTGCGATGCGTTAGATTTTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCAATTGCCCCCAGGTACGTGGATTTTCAGATCGTGCATCGCgtg             >  1:2823481/1‑150 (MQ=255)
                                                                                                                                         ggCCCGCCTATGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGATCGTGCATCGCgtggtg          >  2:1961156/1‑150 (MQ=255)
                                                                                                                                         ggCCCGCCTATGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGATCGTGCATCGCgtggtg          >  2:2379370/1‑150 (MQ=255)
                                                                                                                                            ccGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCCTTTGCCCGCAGGTCCGTGGCTTTTCAGATCGTGCATCGCGTGGTGcct       >  1:3493983/1‑148 (MQ=255)
                                                                                                                                             cGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGATCGTGCATCGCGTGGTGCATa      <  2:1175227/150‑1 (MQ=255)
                                                                                                                                             cGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGATCGTGCATCGCGTGGTGCATa      <  1:2824886/150‑1 (MQ=255)
                                                                                                                                               ccATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGATCGTGCATCGCGTGGTGCATACg    <  1:1248123/150‑1 (MQ=255)
                                                                                                                                                cATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGATCGTGCATCGCGTGGTGCATACgg   <  2:887564/150‑1 (MQ=255)
                                                                                                                                                cATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGATCGTGCATCGCGTGGTGCATACgg   <  1:1593782/150‑1 (MQ=255)
                                                                                                                                                 aTAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGATCGTGCATCGCGTGGTGCATACGgt  <  2:3186987/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
TTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGATCGTGCATCGCGTGGTGCATACGGT  >  NC_000913/4378831‑4379125

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: