Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,359,703 T→A 5.0% intergenic (+89/+24) PP_2949 → / ← PP_2950 permease/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,359,7030TA 5.0% 740.4 / 17.2 213intergenic (+89/+24)PP_2949/PP_2950permease/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (100/102);  new base A (4/7);  total (104/109)
Fisher's exact test for biased strand distribution p-value = 5.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 220 total aligned reads displayed.

GATGCTGCTCGTGACGCTGGGTGCTGTCATCTGTGTTTACGGGGCCAGGAAATATGCCTAAGCACAGCGACGTGCCCAGTCAAACAGCTGGGTATAGCCCCCGCCGTCTGCCCCCGGGCTTAGTTGCACAAAACACCTGATCCTTTATGTAAAAGGATCAGGTTTAACCTGACTTACAGGCCCAGTTCAGTCAGCCCTGGGTGATCATCCGGACGGCGACCCAGTGGCCAGTGGAACTTGCGCTCGCTTTCCTTGATGGGCATGTCGTTGATGCAGGCGAAACGGCGGTGCATCAAGCC  >  NC_002947/3359554‑3359852
                                                                                                                                                     |                                                                                                                                                     
gATGCTGCTCGTGACGCTGGGTGCTGTCATCTGTGTTTACGGGGCCAGGAAATATGCCTAAGCACAGCGACGTGCCCAGTCAATCAGCTGGGTATAGCCCCCGCCGTCTGCCCCCGGGCTTAGTTGCACTAATCACCTGATCCTTTATGt                                                                                                                                                       >  2:2644653/1‑150 (MQ=255)
gATGCTGCTCGTGACGCTGGGTGCTGTCATCTGTGTTTACGGGGCCAGGAAATATGCCTAAGCACAGCGACGTGCCCAGTCAAACAGCTGGGTATAGCCCCCGCCGTCTGCCCCCGGGCTTAGTTGCACAAAACACCTGATCCTTTATGt                                                                                                                                                       >  1:4666703/1‑150 (MQ=255)
gATGCTGCTCGTGACGCTGGGTGCTGTCATCTGTGTTTACGGGGCCAGGAAATATGCCTAAGCACAGCGACGTGCCCAGTCAAACAGCTGGGTATAGCCCCCGCCGTCTGCCCCCGGGCTTAGTTGCACAAAACACCTGATCCTTTATGt                                                                                                                                                       >  1:3547400/1‑150 (MQ=255)
    ctgctCGTGACGCTGGGTGCTGTCATCTGTGTTTACGGGGCCAGGAAATATGCCTAAGCACAGCGACGTGCCCAGTCAAACAGCTGGGTATAGCCCCCGCCGTCTGCCCCCGGGCTTAGTTGCACAAAACACCTGATCCTTTATGTaaaa                                                                                                                                                   <  2:3624493/150‑1 (MQ=255)
     tgctCGTGACGCTGGGTGCTGTCATCTGTGTTTACGGGGCCAGGAAATATGCCTAAGCACAGCGACGTGCCCAGTCAAACAGCTGGGTATAGCCCCCGCCGTCTGCCCCCGGGCTTAGTTGCACAAAACACCTGATCCTTTATGTAAAAg                                                                                                                                                  >  1:5369939/1‑150 (MQ=255)
      gctCGTGACGCTGGGTGCTGTCATCTGTGTTTACGGGGCCAGGAAATATGCCTAAGCACAGCGACGTGCCCAGTCAAACAGCTGGGTATAGCCCCCGCCGTCTGCCCCCGGGCTTAGTTGCACAAAACACCTGATCCTTTATGTAAAAgg                                                                                                                                                 >  1:791712/1‑150 (MQ=255)
       ctCGTGACGCTGGGTGCTGTCATCTGTGTTTACGGGGCCAGGAAATATGCCTAAGCACAGCGACGTGCCCAGTCAAACAGCTGGGTATAGCCCCCGCCGTCTGCCCCCGGGCTTAGTTGCACAAAACACCTGATCCTTTATGTAAAAGGa                                                                                                                                                >  1:1855957/1‑150 (MQ=255)
        tCGTGACGCTGGGTGCTGTCATCTGTGTTTACGGGGCCAGGAAATATGCCTAAGCACAGCGACGTGCCCAGTCAAACAGCTGGGTATAGCCCCCGCCGTCTGCCCCCGGGCTTAGTTGCACAAAACACCTGATCCTTTATGTAAAAGGAt                                                                                                                                               <  1:2837832/150‑1 (MQ=255)
        tCGTGACGCTGGGTGCTGTCATCTGTGTTTACGGGGCCAGGAAATATGCCTAAGCACAGCGACGTGCCCAGTCAAACAGCTGGGTATAGCCCCCGCCGTCTGCCCCCGGGCTTAGTTGCACAAAACACCTGATCCTTTATGTAAAAGGAt                                                                                                                                               >  1:1013454/1‑150 (MQ=255)
        tCGTGACGCTGGGTGCTGTCATCTGTGTTTACGGGGCCAGGAAATATGCCTAAGCACAGCGACGTGCCCAGTCAAACAGCTGGGTATAGCCCCCGCCGTCTGCCCCCGGGCTTAGTTGCACAAAACACCTGATCCTTTATGTAAAAGGAt                                                                                                                                               <  2:5474551/150‑1 (MQ=255)
          gTGACGCTGGGTGCTGTCATCTGTGTTTACGGGGCCAGGAAATATGCCTAAGCACAGCGACGTGCCCAGTCAAACAGCTGGGTATTGGGCCCGCCGTCTGCCCCCGGGCTTAGTTGCACAAAACACCTGATCCTTTATGTAAAAGGATCa                                                                                                                                             <  2:1974173/150‑1 (MQ=255)
          gTGACGCTGGGTGCTGTCATCTGTGTTTACGGGGCCAGGAAATATGCCTAAGCACAGCGACGTGCCCAGTCAAACAGCTGGGTATAGCCCCCGCCGTCTGCCCCCGGGCTTAGTTGCACAAAACACCTGATCCTTTATGTAAAAGGATCa                                                                                                                                             <  1:2314955/150‑1 (MQ=255)
          gTGACGCTGGGTGCTGTCATCTGTGTTTACGGGGCCAGGAAATATGCCTAAGCACAGCGACGTGCCCAGTCAAACAGCTGGGTATAGCCCCCGCCGTCTGCCCCCGGGCTTAGTTGCACAAAACACCTGATCCTTTATGTAAAAGGATCa                                                                                                                                             >  2:2848127/1‑150 (MQ=255)
          gTGACGCTGGGTGCTGTCATCTGTGTTTACGGGGCCAGGAAATATGCCTAAGCACAGCGACGTGCCCAGTCAAACAGCTGGGTATAGCCCCCGCCGTCTGCCCCCGGGCTTAGTTGCACAAAACACCTGATCCTTTATGTAAAAGGATCa                                                                                                                                             <  1:531311/150‑1 (MQ=255)
           tGACGCTGGGTGCTGTCATCTGTGTTTACGGGGCCAGGAAATATGCCTAAGCACAGCGACGTGCCCAGTCAAACAGCTGGGTATAGCCCCCGCCGTCTGCCCCCGGGCTTAGTTGCACAAAACACCTGATCCTTTTACATAAAGGATCAg                                                                                                                                            <  1:3691133/150‑1 (MQ=255)
           tGACGCTGGGTGCTGTCATCTGTGTTTACGGGGCCAGGAAATATGCCTAAGCACAGCGACGTGCCCAGTCAAACAGCTGGGTATAGCCCCCGCCGTCTGCCCCCGGGCTTAGTTGCACAAAACACCTGATCCTTTATGTAAAAGGATCAg                                                                                                                                            >  1:1681009/1‑150 (MQ=255)
           tGACGCTGGGTGCTGTCATCTGTGTTTACGGGGCCAGGAAATATGCCTAAGCACAGCGACGTGCCCAGTCAAACAGCTGGGTATAGCCCCCGCCGTCTGCCCCCGGGCTTAGTTGCACAAAACACCTGATCCTTTATGTAAAAGGATCAg                                                                                                                                            <  2:3103747/150‑1 (MQ=255)
            gACGCTGGGTGCTGTCATCTGTGTTTACGGGGCCAGGAAATATGCCTAAGCACAGCGACGTGCCCAGTCAAACAGCTGGGTATAGCCCCCGCCGTCTGCCCCCGGGCTTAGTTGCACAAAACACCTGATCCTTTTACATAAAGGATCAgg                                                                                                                                           <  1:1543277/150‑1 (MQ=255)
                  gggTGCTGTCATCTGTGTTTACGGGGCCAGGAAATATGCCTAAGCACAGCGACGTGCCCAGTCAAACAGCTGGGTATAGCCCCCGCCGTCTGCCCCCGGGCTTAGTTGCACAAAACACCTGATCCTTTATGTAAAAGGATCAGGTTTAAc                                                                                                                                     <  1:2295686/150‑1 (MQ=255)
                      gCTGTCATCTGTGTTTACGGGGCCAGGAAATATGCCTAAGCACAGCGACGTGCCCAGTCAAACAGCTGGGTATAGCCCCCGCCGTCTGCCCCCGGGCTTAGTTGCACAAAACACCTGATCCTTTATGTAAAAGGATCAGGTTTAACCTgn                                                                                                                                 <  1:3321748/150‑2 (MQ=255)
                      gCTGTCATCTGTGTTTACGGGGCCAGGAAATATGCCTAAGCACAGCGACGTGCCCAGTCAAACAGCTGGGTATAGCCCCCGCCGTCTGCCCCCGGGCTTAGTTGCACAAAACACCTGATCCTTTATGTAAAAGGATCAGGTTTAACCTGa                                                                                                                                 <  1:5610104/150‑1 (MQ=255)
                       cTGTCATCTGTGTTTACGGGGCCAGGAAATATGCCTAAGCACAGCGACGTGCCCAGTCAAACAGCTGGGTATAGCCCCCGCCGTCTGCCCCCGGGCTTAGTTGCACAAAACACCTGATCCTTTATGTAAAAGGATCAGGTTTAACCTGAc                                                                                                                                <  1:25803/150‑1 (MQ=255)
                       cTGTCATCTGTGTTTACGGGGCCAGGAAATATGCCTAAGCACAGCGACGTGCCCAGTCAAACAGCTGGGTATAGCCCCCGCCGTCTGCCCCCGGGCTTAGTTGCACAAAACACCTGATCCTTTATGTAAAAGGATCAGGTTTAACCTGAc                                                                                                                                >  2:5620277/1‑150 (MQ=255)
                       cTGTCATCTGTGTTTACGGGGCCAGGAAATATGCCTAAGCACAGCGACGTGCCCAGTCAAACAGCTGGGTATAGCCCCCGCCGTCTGCCCCCGGGCTTAGTTGCACAAAACACCTGATCCTTTATGTAAAAGGATCAGGTTTAACCTGAc                                                                                                                                <  1:2848127/150‑1 (MQ=255)
                         gTCATCTGTGTTTACGGGGCCAGGAAATATGCCTAAGCACAGCGACGTGCCCAGTCAAACAGCTGGGTATAGCCCCCGCCGTCTGCCCCCGGGCTTAGTTGCACAAAACACCTGATCCTTTATGTAAAAGGATCAGGTTTAACCTGACtt                                                                                                                              >  1:5068696/1‑150 (MQ=255)
                            aTCTGTGTTTACGGGGCCAGGAAATATGCCTAAGCACAGCGACGTGCCCAGTCAAACAGCTGGGTATAGCCCCCGCCGTCTGCCCCCGGGCTTAGTTGCACAAAACACCTGATCCTTTATGTAAAAGGATCAGGTTAAACCTGatccttt                                                                                                                           <  2:537490/150‑7 (MQ=255)
                             tCTGTGTTTACGGGGCCAGGAAATATGCCTAAGCACAGCGACGTGCCCAGTCAAACAGCTGGGTATAGCCCCCGCCGTCTGCCCCCGGGCTTAGTTGCACAAAACACCTGATCCTTTATGTAAAAGGATCAGGTTTAACCTGACTTACAg                                                                                                                          >  1:207847/1‑150 (MQ=255)
                                gtgtTTACGGGGCCAGGAAATATGCCTAAGCACAGCGACGTGCCCAGTCAAACAGCTGGGTATAGCCCCCGCCGTCTGCCCCCGGGCTTAGTTGCACAAAACACCTGATCCTTTATGTAAAAGGATCAGGTTTAACCTGACTTACAGGcc                                                                                                                       <  2:4195955/150‑1 (MQ=255)
                                      aCGGGGCCAGGAAATATGCCTAAGCACAGCGACGTGCCCAGTCAAACAGCTGGGTATAGCCCCCGCCGTCTGCCCCCGGGCTTAGTTGCACAAAACACCTGATCCTTTATGTAAAAGGATCAGGTTTAACCTGACTTACAGGCCCAGTtc                                                                                                                 <  2:2228108/150‑1 (MQ=255)
                                       cGGGGCCAGGAAATATGCCTAAGCACAGCGACGTGCCCAGTCAAACAGCTGGGTATAGCCCCCGCCGTCTGCCCCCGGGCTTAGTTGCACAAAACACCTGATCCTTTATGTAAAAGGATCAGGTTTAACCTGACTTACAGGCCCAGTtca                                                                                                                <  2:1801006/150‑1 (MQ=255)
                                       cGGGGCCAGGAAATATGCCTAAGCACAGCGACGTGCCCAGTCAAACAGCTGGGTATAGCCCCCGCCGTCTGCCCCCGGGCTTAGTTGCACAAAACACCTGATCCTTTATGTAAAAGGATCAGGTTTAACCTGACTTACAGGCCCAGTtca                                                                                                                <  1:5038311/150‑1 (MQ=255)
                                        ggggCCAGGAAATATGCCTAAGCACAGCGACGTGCCCAGTCAAACAGCTGGGTATAGCCCCCGCCGTCTGCCCCCGGGCTTAGTTGCACAAAACACCTGATCCTTTATGTAAAAGGATCAGGTTTAACCTGACTTACAGGCCCAGTtcag                                                                                                               >  2:4621050/1‑150 (MQ=255)
                                         gggCCAGGAAATATGCCTAAGCACAGCGACGTGCCCAGTCAAACAGCTGGGTATAGCCCCCGCCGTCTGCCCCCGGGCTTAGTTGCACAAAACACCTGATCCTTTATGTAAAAGGATCAGGTTTAACCTGACTTACAGGCCCAGTtcagt                                                                                                              <  1:212482/150‑1 (MQ=255)
                                          ggCCAGGAAATATGCCTAAGCACAGCGACGTGCCCAGTCAAACAGCTGGGTATAGCCCCCGCCGTCTGCCCCCGGGCTTAGTTGCACAAAACACCTGATCCTTTATGTAAAAGGATCAGGTTTAACCTGACTTACAGGCCCAGTtcagtc                                                                                                             <  2:4679879/150‑1 (MQ=255)
                                           gCCAGGAAATATGCCTAAGCACAGCGACGTGCCCAGTCAAACAGCTGGGTATAGCCCCCGCCGTCTGCCCCCGGGCTTAGTTGCACAAAACACCTGATCCTTTATGTAAAAGGATCAGGTTTAACCTGACTTACAGGCCCAGTtcagtca                                                                                                            >  1:1997370/1‑150 (MQ=255)
                                             cAGGAAATATGCCTAAGCACAGCGACGTGCCCAGTCAAACAGCTGGGTATAGCCCCCGCCGTCTGCCCCCGGGCTTAGTTGCACAAAACACCTGATCCTTTATGTAAAAGGATCAGGTTTAACCTGACTTACAGGCCCAGTTCAGTCAGc                                                                                                          <  2:2011729/150‑1 (MQ=255)
                                             cAGGAAATATGCCTAAGCACAGCGACGTGCCCAGTCAAACAGCTGGGTATAGCCCCCGCCGTCTGCCCCCGGGCTTAGTTGCACAAAACACCTGATCCTTTATGTAAAAGGATCAGGTTTAACCTGACTTACAGGCCCAGTTCAGTCAGc                                                                                                          <  1:4160729/150‑1 (MQ=255)
                                             cAGGAAATATGCCTAAGCACAGCGACGTGCCCAGTCAAACAGCTGGGTATAGCCCCCGCCGTCTGCCCCCGGGCTTAGTTGCACAAAACACCTGATCCTTTATGTAAAAGGATCAGGTTTAACCTGACTTACAGGCCCAGTTCAGTCAGc                                                                                                          <  2:3888937/150‑1 (MQ=255)
                                             cAGGAAATATGCCTAAGCACAGCGACGTGCCCAGTCAAACAGCTGGGTATAGCCCCCGCCGTCTGCCCCCGGGCTTAGTTGCACAAAACACCTGATCCTTTATGTAAAAGGATCAGGTTTAACCTGACTTACAGGCCCAGTTCAGTCAGc                                                                                                          >  2:1274562/1‑150 (MQ=255)
                                                  aaTATGCCTAAGCACAGCGACGTGCCCAGTCAAACAGCTGGGTATAGCCCCCGCCGTCTGCCCCCGGGCTTAGTTGCACAAAACACCTGATCCATTATGAAAAAGGATCAGGATTAACCTGACTTACAGGCCCAGTTCAGTCAGCACTgg                                                                                                     <  2:5101830/150‑1 (MQ=255)
                                                      tGCCTAAGCACAGCGACGTGCCCAGTCAAACAGCTGGGTATAGCCCCCGCCGTCTGCCCCCGGGCTTAGTTGCACAAAACACCTGATCCTTTATGTAAAAGGATCAGGTTTAACCTGACTTACAGGCCCAGTTCAGTCAGCCCTGGGTGa                                                                                                 <  1:987589/150‑1 (MQ=255)
                                                          tAAGCACAGCGACGTGCCCAGTCAAACAGCTGGGTATAGCCCCCGCCGTCTGCCCCCGGGCTTAGTTGCACAAAACACCTGATCCTTTATGTAAAAGGATCAGGTTTAACCTGACTTACAGGCCCAGTTCAGTCAGCCCTGGGTGatcat                                                                                             >  2:2585206/1‑150 (MQ=255)
                                                            aGCACAGCGACGTGCCCAGTCAAACAGCTGGGTATAGCCCCCGCCGTCTGCCCCCGGGCTTAGTTGCACAAAACACCTGATCCTTTATGTAAAAGGATCAGGTTTAACCTGACTTACAGGCCCAGTTCAGTCAGCCCTGGGTGATCATcc                                                                                           <  1:4829587/150‑1 (MQ=255)
                                                            aGCACAGCGACGTGCCCAGTCAAACAGCTGGGTATAGCCCCCGCCGTCTGCCCCCGGGCTTAGTTGCACAAAACACCTGATCCTTTATGTAAAAGGATCAGGTTTAACCTGACTTACAGGCCCAGTTCAGTCAGCCCTGGGTGATCATcc                                                                                           <  2:3944489/150‑1 (MQ=255)
                                                            aGCACAGCGACGTGCCCAGTCAAACAGCTGGGTATAGCCCCCGCCGTCTGCCCCCGGGCTTAGTTGCACAAAACACCTGATCCTTTATGTAAAAGGATCAGGTTTAACCTGACTTACAGGCCCAGTTCAGTCAGCCCTGGGTGATCATcc                                                                                           <  2:1454545/150‑1 (MQ=255)
                                                               acaGCGACGTGCCCAGTCAAACAGCTGGGTATAGCCCCCGCCGTCTGCCCCCGGGCTTAGTTGCACAAAACACCTGATCCTTTATGTAAAAGGATCAGGTTTAACCTGACTTACAGGCCCAGTTCAGTCAGCCCTGGGTGATCATCCGGa                                                                                        >  1:2399033/1‑150 (MQ=255)
                                                               acaGCGACGTGCCCAGTCAAACAGCTGGGTATAGCCCCCGCCGTCTGCCCCCGGGCTTAGTTGCACAAAACACCTGATCCTTTATGTAAAAGGATCAGGTTTAACCTGACTTACAGGCCCAGTTCAGTCAGCCCTGGGTGATCATCCGGa                                                                                        <  1:2420322/150‑1 (MQ=255)
                                                                  gCGACGTGCCCAGTCAAACAGCTGGGTATAGCCCCCGCCGTCTGCCCCCGGGCTTAGTTGCACAAAACACCTGATCCTTTATGTAAAAGGATCAGGTTTAACCTGACTTACAGGCCCAGTTCAGTCAGCCCTGGGTGATCATCCGGACgg                                                                                     <  1:252636/150‑1 (MQ=255)
                                                                     aCGTGCCCAGTCAAACAGCTGGGTATAGCCCCCGCCGTCTGCCCCCGGGCTTAGTTGCACAAAACACCTGATCCTTTATGTAAAAGGATCAGGTTTAACCTGACTTACAGGCCCAGTTCAGTCAGCCCTGGGTGATCATCCGGACGGCGa                                                                                  <  2:791712/150‑1 (MQ=255)
                                                                      cGTGCCCAGTCAAACAGCTGGGTATAGCCCCCGCCGTCTGCCCCCGGGCTTAGTTGCACAAAACACCTGATCCTTTATGTAAAAGGATCAGGTTTAACCTGACTTACAGGCCCAGTTCAGTCAGCCCTGGGTGATCATCCGGACGGCGAc                                                                                 <  2:1544421/150‑1 (MQ=255)
                                                                       gTGCCCAGTCAAACAGCTGGGTATAGCCCCCGCCGTCTGCCCCCGGGCTTAGTTGCACAAAACACCTGATCCTTTATGTAAAAGGATCAGGTTTAACCTGACTTACAGGCCCAGTTCAGTCAGCCCTGGGTGATCATCCGGACGGCGAcc                                                                                <  2:5415399/150‑1 (MQ=255)
                                                                       gTGCCCAGTCAAACAGCTGGGTATAGCCCCCGCCGTCTGCCCCCGGGCTTAGTTGCACAAAACACCTGATCCTTTATGTAAAAGGATCAGGTTTAACCTGACTTACAGGCCCAGTTCAGTCAGCCCTGGGTGATCATCCGGACGGCGAcc                                                                                <  1:1558237/150‑1 (MQ=255)
                                                                       gTGCCCAGTCAAACAGCTGGGTATAGCCCCCGCCGTCTGCCCCCGGGCTTAGTTGCACAAAACACCTGATCCTTTATGTAAAAGGATCAGGTTTAACCTGACTTACAGGCCCAGTTCAGTCAGCCCTGGGTGATCATCCGGACGGCGAcc                                                                                >  1:5136875/1‑150 (MQ=255)
                                                                       gTGCCCAGTCAAACAGCTGGGTATAGCCCCCGCCGTCTGCCCCCGGGCTTAGTTGCACAAAACACCTGATCCTTTATGTAAAAGGATCAGGTTTAACCTGACTTACAGGCCCAGTTCAGTCAGCCCTGGGTGATCATCCGGACGGCGAcc                                                                                >  2:682088/1‑150 (MQ=255)
                                                                         gCCCAGTCAAACAGCTGGGTATAGCCCCCGCCGTCTGCCCCCGGGCTTAGGTGCACAAAACACCTGATCCTTTAGGTAAAAGGATCAGGTTTAACCTGACTTACAGGCCCAGTTCAGTCAGCCCTGGGTGATCATCCGGACGGCGACcca                                                                              <  2:3442463/150‑1 (MQ=255)
                                                                              gTCAAACAGCTGGGTATAGCCCCCGCCGTCTGCCCCCGGGCTTAGTTGCACAAAACACCTGATCCTTTATGTAAAAGGATCAGGTTTAACCTGACTTACAGGCCCAGTTCAGTCAGCCCTGGGTGATCATCCGGACGGCGACccagtggc                                                                         >  2:3381539/1‑150 (MQ=255)
                                                                               tCAAACAGCTGGGTATAGCCCCCGCCGTCTGCCCCCGGGCTTAGTTGCACAAAACACCTGATCCTTTATGTAAAAGGATCAGGTTTAACCTGACTTACAGGCCCAGTTCAGTCAGCCCTGGGTGATCATCCGGACGGCGACccagtggcc                                                                        >  2:494384/1‑150 (MQ=255)
                                                                               tCAAACAGCTGGGTATAGCCCCCGCCGTCTGCCCCCGGGCTTAGTTGCACAAAACACCTGATCCTTTATGTAAAAGGATCAGGTTAACCTGACTTACAGGCCCAGTTCAGTCAGCCCTGGGTGATCATCCGGACGGCGACccagtggcca                                                                       >  1:1501959/1‑150 (MQ=255)
                                                                                 aaaCAGCTGGGTATAGCCCCCGCCGTCTGCCCCCGGGCTTAGTTGCACAAAACACCTGATCCTTTATGTAAAAGGATCAGGTTTAACCTGACTTACAGGCCCAGTTCAGTCAGCCCTGGGTGATCAGCCGGACGGCGACccagtggccag                                                                      >  2:4531526/1‑150 (MQ=255)
                                                                                  aaCAGCTGGGTATAGCCCCCGCCGTCTGCCCCCGGGCTTAGTTGCACAAAACACCTGATCCTTTATGTAAAAGGATCAGGTTTAACCTGACTTACAGGCCCAGTTCAGTCAGCCCTGGGTGATCATCCGGACGGCGACccagtggccagt                                                                     <  1:1370102/150‑1 (MQ=255)
                                                                                  aaCAGCTGGGTATAGCCCCCGCCGTCTGCCCCCGGGCTTAGTTGCACAAAACACCTGATCCTTTATGTAAAAGGATCAGGTTTAACCTGACTTACAGGCCCAGTTCAGTCAGCCCTGGGTGATCATCCGGACGGCGACccagtggccagt                                                                     <  1:2589892/150‑1 (MQ=255)
                                                                                   aCAGCTGGGTATAGCCCCCGCCGTCTGCCCCCGGGCTTAGTTGCACAAAACACCTGATCCTTTATGTAAAAGGATCAGGTTTAACCTGACTTACAGGCCCAGTTCAGTCAGCCCTGGGTGATCATCCGGACGGCGACccagtggccagtg                                                                    <  1:4234797/150‑1 (MQ=255)
                                                                                   aCAGCTGGGTATAGCCCCCGCCGTCTGCCCCCGGGCTTAGTTGCACAAAACACCTGATCCTTTATGTAAAAGGATCAGGTTTAACCTGACTTACAGGCCCAGTTCAGTCAGCCCTGGGTGATCATCCGGACGGCGACccagtggccagtg                                                                    <  1:3436209/150‑1 (MQ=255)
                                                                                      gCTGGGTATAGCCCCCGCCGTCTGCCCCCGGGCTTAGTTGCACAAAACACCTGATCCTTTATGTAAAAGGATCAGGTTTAACCTGACTTACAGGCCCAGTTCAGTCAGCCCTGGGTGATCATCCGGACGGCGACCCAGTGGCCAGTGGaa                                                                 <  2:2401378/150‑1 (MQ=255)
                                                                                         gggTATAGCCCCCGCCGTCTGCCCCCGGGCTTAGTTGCACAAAACACCTGATCCTTTATGTAAAAGGATCAGGTTTAACCTGACTTACAGGCCCAGTTCAGTCAGCCCTGGGTGATCATCCGGACGGCGACCCAGTGGCCAGTGGAACtt                                                              <  2:1959918/150‑1 (MQ=255)
                                                                                           gTATAGCCCCCGCCGTCTGCCCCCGGGCTTAGTTGCACAAAACACCTGATCCTTTATGTAAAAGGATCAGGTTTAACCTGACTTACAGGCCCAGTTCAGTCAGCCCTGGGTGATCATCCGGACGGCGACCCAGTGGCCAGTGGAACTTgc                                                            <  2:5595470/150‑1 (MQ=255)
                                                                                           gTATAGCCCCCGCCGTCTGCCCCCGGGCTTAGTTGCACAAAACACCTGATCCTTTATGTAAAAGGATCAGGTTTAACCTGACTTACAGGCCCAGTTCAGTCAGCCCTGGGTGATCATCCGGACGGCGACCCAGTGGCCAGTGGAACTTgc                                                            >  2:2135392/1‑150 (MQ=255)
                                                                                                 cccccGCCGTCTGCCCCCGGGCTTAGTTGCACAAAACACCTGATCCTTTATGTAAAAGGATCAGGTTTAACCTGACTTACAGGCCCAGTTCAGTCAGCCCTGGGTGATCATCCGGACGGCGACCCAGTGGCCAGTGGAACTTGcgctcgc                                                      <  1:577923/150‑1 (MQ=255)
                                                                                                   cccGCCGTCTGCCCCCGGGCTTAGTTGCACAAAACACCTGATCCTTTATGTAAAAGGATCAGGTTTAACCTGACTTACAGGCCCAGTTCAGTCAGCCCTGGGTGATCATCCGGACGGCGACCCAGTGGCCAGTGGAACTTGCGCTCGCtt                                                    >  1:3443635/1‑150 (MQ=255)
                                                                                                      gccgTCTGCCCCCGGGCTTAGTTGCACAAAACACCTGATCCTTTATGTAAAAGGATCAGGTTTAACCTGACTTACAGGCCCAGTTCAGTCAGCCCTGGGTGATCATCCGGACGGCGACCCAGTGGCCAGTGGAACTTGCGCTCGCTTTcc                                                 >  2:150604/1‑150 (MQ=255)
                                                                                                            tGCCCCCGGGCTTAGTTGCACAAAACACCTGATCCTTTATGTAAAAGGATCAGGTTTAACCTGACTTACAGGCCCAGTTCAGTCAGCCCTGGGTGATCATCCGGACGGCGACCCAGTGGCCAGTGGAACTTGCGCTCGCTTTCCTTGATg                                           >  2:1709788/1‑150 (MQ=255)
                                                                                                             gCCCCCGGGCTTAGTTGCACAAAACACCTGATCCTTTATGTAAAAGGATCAGGTTTAACCTGACTTACAGGCCCAGTTCAGTCAGCCCTGGGTGATCATCCGGACGGCGACCCAGTGGCCAGTGGAACTTGCGCTCGCTTTCCTTGATgg                                          >  2:1861263/1‑150 (MQ=255)
                                                                                                               ccccGGGCTTAGTTGCACAAAACACCTGATCCTTTATGTAAAAGGATCAGGTTTAACCTGACTTACAGGCCCAGTTCAGTCAGCCCTGGGTGATCATCCGGACGGCGACCCAGTGGCCAGTGGAACTTGCGCTCGCTTTCCTTGATGGGc                                        >  1:90006/1‑150 (MQ=255)
                                                                                                                cccGGGCTTAGTTGCACAAAACACCTGATCCTTTATGTAAAAGGATCAGGTTTAACCTGACTTACAGGCCCAGTTCAGTCAGCCCTGGGTGATCATCCGGACGGCGACCCAGTGGCCAGTGGAACTTGCGCTCGCTTTCCTTGATGGGCa                                       >  2:5007671/1‑150 (MQ=255)
                                                                                                                cccGGGCTTAGTTGCACAAAACACCTGATCCTTTATGTAAAAGGATCAGGTTTAACCTGACTTACAGGCCCAGTTCAGTCAGCCCTGGGTGATCATCCGGACGGCGACCCAGTGGCCAGTGGAACTTGCGCTCGCTTTCCTTGATGGGCa                                       >  1:325019/1‑150 (MQ=255)
                                                                                                                 ccGGGCTTAGTTGCACAAAACACCTGATCCTTTATGTAAAAGGATCAGGTTTAACCTGACTTACAGGCCCAGTTCAGTCAGCCCTGGGTGATCATCCGGACGGCGACCCAGTGGCCAGTGGAACTTGCGCTCGCTTTCCTTGATGGGCAt                                      >  2:4749837/1‑150 (MQ=255)
                                                                                                                  cGGGCTTAGTTGCACAAAACACCTGATCCTTTATGTAAAAGGATCAGGTTTAACCTGACTTACAGGCCCAGTTCAGTCAGCCCTGGGTGATCATCCGGACGGCGACCCAGTGGCCAGTGGAACTTGCGCTCGCTTTCCTTGATGGGCATg                                     <  2:4123413/150‑1 (MQ=255)
                                                                                                                  cGGGCTTAGTTGCACAAAACACCTGATCCTTTATGTAAAAGGATCAGGTTTAACCTGACTTACAGGCCCAGTTCAGTCAGCCCTGGGTGATCATCCGGACGGCGACCCAGTGGCCAGTGGAACTTGCGCTCGCTTTCCTTGATGGGCATg                                     >  1:4459179/1‑150 (MQ=255)
                                                                                                                       ttAGTTGCACAAAACACCTGATCCTTTATGTAAAAGGATCAGGTTTAACCTGACTTACAGGCCCAGTTCAGTCAGCCCTGGGTGATCATCCGGACGGCGACCCAGTGGCCAGTGGAACTTGCGCTCGCTTTCCTTGATGGGCATGTCGtt                                >  1:2242913/1‑150 (MQ=255)
                                                                                                                       ttAGTTGCACAAAACACCTGATCCTTTATGTAAAAGGATCAGGTTTAACCTGACTTACAGGCCCAGTTCAGTCAGCCCTGGGTGATCATCCGGACGGCGACCCAGTGGCCAGTGGAACTTGCGCTCGCTTTCCTTGATGGGCATGTCGtt                                <  2:3338512/150‑1 (MQ=255)
                                                                                                                          gTTGCACAAAACACCTGATCCTTTATGTAAAAGGATCAGGTTTAACCTGACTTACAGGCCCAGTTCAGTCAGCCCTGGGTGATCATCCGGACGGCGACCCAGTGGCCAGTGGAACTTGCGCTCGATTTCCTTGATGGGCATGTCGTTGAt                             >  1:3142980/1‑150 (MQ=255)
                                                                                                                          gTTGCACAAAACACCTGATCCTTTATATAAAAGGATCAGGTTTAACCTGACTTACAGGCCCAGTTCAGTCAGCCCTGGGTGATCATCCGGACGGCGACCCAGTGGCCAGTGGAACTTACGCTCGCTTTCCTTGATGGGCATGTCGTTGAt                             <  1:2632950/150‑1 (MQ=255)
                                                                                                                           ttGCACAAAACACCTGATCCTTTATGTAAAAGGATCAGGTTTAACCTGACTTACAGGCCCAGTTCAGTCAGCCCTGGGTGATCATCCGGACGGCGACCCAGTGGCCAGTGGAACTTGCGCTCGCTTTc                                                  >  2:767350/1‑128 (MQ=255)
                                                                                                                             gCACAAAACACCTGATCCTTTATGTAAAAGGATCAGGTTTAACCTGACTTACAGGCCCAGTTCAGTCAGCCCTGGGTGATCATCCGGACGGCGACCCAGTGGCCAGTGGAACTTGCGCTCGCTTTCCTTGATGGGCATGTCGTTGATGCa                          <  1:2310176/150‑1 (MQ=255)
                                                                                                                              cacaAAACACCTGATCCTTTATGTAAAAGGATCAGGTTTAACCTGACTTACAGGCCCAGTTCAGTCAGCCCTGGGTGATCATCCGGACGGCGACCCAGTGGCCAGTGGAACTTGCGCTCGCTTTCCTTGATGGGCATGTCGTTGATGCAg                         >  2:3088826/1‑150 (MQ=255)
                                                                                                                               acaAAACACCTGATCCTTTATGTAAAAGGATCAGGTTTAACCTGACTTACAGGCCCAGTTCAGTCAGCCCTGGGTGATCATCCGGACGGCGACCCAGTGGCCAGTGGAACTTGCGCTCGCTTTCCTTGATGGGCATGTCGTTGATGCAgg                        >  1:4030282/1‑150 (MQ=255)
                                                                                                                                 aaaaCACCTGATCCTTTATGTAAAAGGATCAGGTTTAACCTGACTTACAGGCCCAGTTCAGTCAGCCCTGGGTGATCATCCGGACGGCGACCCAGTGGCCAGTGGAACTTGCGCTCGCTTTCCTTGATGGGCATGTCGTTGATGCAGGCg                      >  1:860130/1‑150 (MQ=255)
                                                                                                                                   aaCACCTGATCCTTTATGTAAAAGGATCAGGTTTAACCTGACTTACAGGCCCAGTTCAGTCAGCCCTGGGTGATCATCCGGACGGCGACCCAGTGGCCAGTGGAACTTGCGCTCGCTTTCCTTGATGGGCATGTCGTTGATGCAGGCGaa                    >  1:1728930/1‑150 (MQ=255)
                                                                                                                                      acCTGATCCTTTATGTAAAAGGATCAGGTTTAACCTGACTTACAGGCCCAGTTCAGTCAGCCCTGGGTGATCATCCGGACGGCGACCCAGTGGCCAGTGGAACTTGCGCTCGCTTTCCTTGATGGGCATGTCGTTGATGCAGGCGAAAcg                 >  2:2927483/1‑150 (MQ=255)
                                                                                                                                        cTGATCCTTTATGTAAAAGGATCAGGTTTAACCTGACTTACAGGCCCAGTTCAGTCAGCCCTGGGTGATCATCCGGACGGCGACCCAGTGGCCAGTGGAACTTGCGCTCGCTTTCCTTGATGGGCATGTCGTTGATGCAGGCGAAAcggc               <  1:1936121/150‑1 (MQ=255)
                                                                                                                                           aTCCTTTATGTAAAAGGATCAGGTTTAACCTGACTTACAGGCCCAGTTCAGTCAGCCCTGGGTGATCATCCGGACGGCGACCCAGTGGCCAGTGGAACTTGCGCTCGCTTTCCTTGATGGGCATGTCGTTGATGCAGGCGAAACGGCGGt            >  1:1408136/1‑150 (MQ=255)
                                                                                                                                             ccTTTATGTAAAAGGATCAGGTTTAACCTGACTTACAGGCCCAGTTCAGTCAGCCCTGGGTGATCATCCGGACGGCGACCCAGTGGCCAGTGGAACTTGCGCTCGCTTTCCTTGATGGGCATGTCGTTGATGCAGGCGAAACGGCGGTGc          <  1:4723168/150‑1 (MQ=255)
                                                                                                                                              cTTTATGTAAAAGGATCAGGTTTAACCTGACTTACAGGCCCAGTTCAGTCAGCCCTGGGTGATCATCCGGACGGCGACCCAGTGGCCAGTGGAACTTGCGCTCGCTTTCCTTGATGGGCATGTCGTTGATGCAGGCGAAACGGCGGTGCa         >  1:5425161/1‑150 (MQ=255)
                                                                                                                                               tttATGTAAAAGGATCAGGTTTAACCTGACTTACAGGCCCAGTTCAGTCAGCCCTGGGTGATCATCCGGACGGCGACCCAGTGGCCAGTGGAACTTGCGCTCGCTTTCCTTGATGGGCATGTCGTTGATGCAGGCGAAACGGCGGTGCAt        >  1:744933/1‑150 (MQ=255)
                                                                                                                                               tttATGTAAAAGGATCAGGTTTAACCTGACTTACAGGCCCAGTTCAGTCAGCCCTGGGTGATCATCCGGACGGCGACCCAGTGGCCAGTGGAACTTGCGCTCGCTTTCCTTGATGGGCATGTCGTTGATGCAGGCGAAACGGCGGTGCAt        >  1:2195176/1‑150 (MQ=255)
                                                                                                                                                   tGTAAAAGGATCAGGTTTAACCTGACTTACAGGCCCAGTTCAGTCAGCCCTGGGTGATCATCCGGACGGCGACCCAGTGGCCAGTGGAACTTGCGCTCGCTTTCCTTGATGGGCATGTCGTTGATGCAGGCGAAACGGCGGTGCATCAAg    >  1:2899331/1‑150 (MQ=255)
                                                                                                                                                   tGTAAAAGGATCAGGTTTAACCTGACTTACAGGCCCAGTTCAGTCAGCCCTGGGTGATCATCCGGACGGCGACCCAGTGGCCAGTGGAACTTGCGCTCGCTTTCCTTGATGGGCATGTCGTTGATGCAGGCGAAACGGCGGTGCATCAAg    <  1:1274562/150‑1 (MQ=255)
                                                                                                                                                    gTAAAATGATCAGGTTTAACCTGTCGTACAGGCCCCGTTCAGTCCGCCCTGGGTGTTCATCCGGTCGGGGACCCAGGGGCCAGTGGCACTTGGGGTCGCTTTCCTTGATGGGCATGTCGTTGATGCAGGCGAAACGGCGGTGCATCAAGc   <  2:5068762/150‑1 (MQ=255)
                                                                                                                                                    gTAAAAGGATCAGGTTTAACCTGACTTACAGGCCCAGTTCAGTCAGCCCTGGGTGATCATCCGGACGGCGACCCAGTGGCCAGTGGAACTTGCGCTCGCTTTCCTTGATGGGCATGTCGTTGATGCAGGCGAAACGGCGGTGCATCAAGc   >  2:291550/1‑150 (MQ=255)
                                                                                                                                                     tAAAAGGATCAGGTTTAACCTGACTTACAGGCCCAGTTCAGTCAGCCCTGGGTGATCATCCg                                                                                          >  2:493310/1‑62 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GATGCTGCTCGTGACGCTGGGTGCTGTCATCTGTGTTTACGGGGCCAGGAAATATGCCTAAGCACAGCGACGTGCCCAGTCAAACAGCTGGGTATAGCCCCCGCCGTCTGCCCCCGGGCTTAGTTGCACAAAACACCTGATCCTTTATGTAAAAGGATCAGGTTTAACCTGACTTACAGGCCCAGTTCAGTCAGCCCTGGGTGATCATCCGGACGGCGACCCAGTGGCCAGTGGAACTTGCGCTCGCTTTCCTTGATGGGCATGTCGTTGATGCAGGCGAAACGGCGGTGCATCAAGCC  >  NC_002947/3359554‑3359852

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