Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 587,079 C→G 7.5% P43R (CCG→CGG)  PP_0497 → Cro/CI family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947587,0790CG 7.5% 567.6 / 26.4 174P43R (CCG→CGG) PP_0497Cro/CI family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (93/68);  new base G (9/4);  total (102/72)
Fisher's exact test for biased strand distribution p-value = 5.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

Only 100 of 180 total aligned reads displayed.

TCGCGTCACCTGGATTCCGCATATGAGCGATCACTTCGCTACAAACCTCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCGGTGCGCGGGTATCGCCCCAGGCCCAGCAACCCAGTGACCCGATCAGCGAGTTGTTCC  >  NC_002947/586930‑587228
                                                                                                                                                     |                                                                                                                                                     
  gcgTCACCTGGATTCCGCATATGAGCGATCACTTCGCTACAAACCTCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGa                                                                                                                                                     >  2:648711/1‑150 (MQ=255)
        ccTGGATTCCGCATATGAGCGATCACTTCGCTACAAACCTCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCtt                                                                                                                                               <  1:4245/150‑1 (MQ=255)
         cTGGATTCCGCATATGAGCGATCACTTCGCTACAAACCTCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCttt                                                                                                                                              <  2:4546643/150‑1 (MQ=255)
          tGGATTCCGCATATGAGCGATCACTTCGCTACAAACCTCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTc                                                                                                                                             >  2:1938951/1‑150 (MQ=255)
            gATTCCGCATATGAGCGATCACTTCGCTACAAACCTCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCaa                                                                                                                                           <  1:3908644/150‑1 (MQ=255)
            gATTCCGCATATGAGCGATCACTTCGCTACAAACCTCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCaa                                                                                                                                           <  2:4557318/150‑1 (MQ=255)
             aTTCCGCATATGAGCGATCACTTCGCTACAAACCTCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAAc                                                                                                                                          <  1:4459746/150‑1 (MQ=255)
             aTTCCGCATATGAGCGATCACTTCGCTACAAACCTCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAAc                                                                                                                                          <  2:1951365/150‑1 (MQ=255)
              ttCCGCATATGAGCGATCACTTCGCTACAAACCTCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAAcc                                                                                                                                         >  1:5344034/1‑150 (MQ=255)
                ccGCATATGAGCGATCACTTCGCTACAAACCTCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTg                                                                                                                                       >  1:2401682/1‑150 (MQ=255)
                 cGCATATGAGCGATCACTTCGCTACAAACCTCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGa                                                                                                                                      >  1:3990194/1‑150 (MQ=255)
                  gCATATGAGCGATCACTTCGCTACAAACCTCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGaa                                                                                                                                     <  1:4140617/150‑1 (MQ=255)
                    atatGAGCGATCACTTCGCTACAAACCTCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGc                                                                                                                                   >  2:4379298/1‑150 (MQ=255)
                        gAGCGATCACTTCGCTACAAACCTCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTAt                                                                                                                               >  1:1010280/1‑150 (MQ=255)
                          gCGATCACTTCGCTACAAACCTCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCg                                                                                                                             >  2:1983295/1‑150 (MQ=255)
                              tCACTTCGCTACAAACCTCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGa                                                                                                                         >  1:38373/1‑150 (MQ=255)
                              tCACTTCGCTACAAACCTCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGa                                                                                                                         <  1:408637/150‑1 (MQ=255)
                                      cTACAAACCTCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCg                                                                                                                 <  1:400315/150‑1 (MQ=255)
                                      cTACAAACCTCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCg                                                                                                                 <  2:110552/150‑1 (MQ=255)
                                       tACAAACCTCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCgg                                                                                                                >  1:292968/1‑150 (MQ=255)
                                       tACAAACCTCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCgg                                                                                                                >  2:1555397/1‑150 (MQ=255)
                                       tACAAACCTCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCgg                                                                                                                >  1:2004906/1‑150 (MQ=255)
                                        aCAAACCTCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGgt                                                                                                               >  2:3933183/1‑150 (MQ=255)
                                         cAAACCTCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGgtg                                                                                                              >  1:2623368/1‑150 (MQ=255)
                                          aaaCCTCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGgtgt                                                                                                             >  1:4174889/1‑150 (MQ=255)
                                             ccTCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGa                                                                                                          <  1:4202432/150‑1 (MQ=255)
                                              cTCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAg                                                                                                         >  2:3186047/1‑150 (MQ=255)
                                              cTCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAg                                                                                                         <  1:1714848/150‑1 (MQ=255)
                                               tCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGa                                                                                                          <  2:3899795/148‑1 (MQ=255)
                                               tCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAgg                                                                                                        <  1:278387/150‑1 (MQ=255)
                                               tCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAgg                                                                                                        >  2:4634037/1‑150 (MQ=255)
                                                cAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGa                                                                                                       <  2:2258787/150‑1 (MQ=255)
                                                  aaCTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGAtt                                                                                                     <  2:1396765/150‑1 (MQ=255)
                                                     ttGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACg                                                                                                  >  2:327743/1‑150 (MQ=255)
                                                     ttGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACg                                                                                                  <  2:1081963/150‑1 (MQ=255)
                                                     ttGCCTGCAGTCACTACCGCTCTATTTACGAAGTTTGCATCCAGCTGTATATCAACCGCTCTCAGTTCAACAATTACCTGAGAGTGCAAAGCCGTCCGACGGATTTCAACATGAAGCGTATCGGCTACTTCTTCGGTGTCgatgataaag                                                                                                  >  2:3437259/1‑146 (MQ=255)
                                                       gCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGaa                                                                                                <  1:4297660/150‑1 (MQ=255)
                                                       gCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGaa                                                                                                <  1:3241956/150‑1 (MQ=255)
                                                        ccTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAAc                                                                                               >  2:947336/1‑150 (MQ=255)
                                                        ccTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGAAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAAc                                                                                               <  2:2641049/150‑1 (MQ=255)
                                                         cTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACt                                                                                              >  1:2600442/1‑150 (MQ=255)
                                                               tCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCt                                                                                        >  2:1022541/1‑150 (MQ=255)
                                                                 aCTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAAc                                                                                                                                          <  1:2162446/98‑1 (MQ=255)
                                                                   tACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGccc                                                                                    >  2:4045703/1‑150 (MQ=255)
                                                                   tACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGccc                                                                                    >  1:2102358/1‑150 (MQ=255)
                                                                    gccGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGTAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCccg                                                                                   >  1:829170/2‑150 (MQ=255)
                                                                     ccGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCccgc                                                                                  >  1:1927980/1‑150 (MQ=255)
                                                                      cGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCt                                                                                        >  2:2950553/1‑143 (MQ=255)
                                                                      cGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCccgcc                                                                                 >  2:4842408/1‑150 (MQ=255)
                                                                           tATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCa                                                                            >  2:5223035/1‑150 (MQ=255)
                                                                               tCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTc                                                                        >  1:1501258/1‑150 (MQ=255)
                                                                               tCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTc                                                                        >  2:4413714/1‑150 (MQ=255)
                                                                                 cGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGc                                                                      >  1:3757653/1‑150 (MQ=255)
                                                                                 cGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACCGAACCTGCCCGCCGAGCAGTTCGc                                                                      >  1:640997/1‑150 (MQ=255)
                                                                                  gAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACgg                                                                                                                                                  <  2:2361493/73‑1 (MQ=255)
                                                                                     gTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGGTCGCCCGt                                                                  >  1:4273658/1‑150 (MQ=255)
                                                                                      tttGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGGACCTGCCCGCCGAGCAGTTCGCCCGTc                                                                 <  1:2584503/150‑1 (MQ=255)
                                                                                      tttGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTc                                                                 <  2:3186721/150‑1 (MQ=255)
                                                                                          ccgccAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGAt                                                             <  2:3601630/150‑1 (MQ=255)
                                                                                            gccAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCg                                                           <  2:2093067/150‑1 (MQ=255)
                                                                                             ccAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCgg                                                          >  1:1073038/1‑150 (MQ=255)
                                                                                                 cTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCGGTgcg                                                      >  1:3029839/1‑150 (MQ=255)
                                                                                                  tGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGc                                                                             >  1:5267556/1‑126 (MQ=255)
                                                                                                   gTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCGGTGcgcg                                                    >  2:3484981/1‑150 (MQ=255)
                                                                                                    tCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGTCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCGGTGCTCgg                                                   <  2:3232921/150‑1 (MQ=255)
                                                                                                    tCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCGGTGCGCgg                                                   >  1:2694999/1‑150 (MQ=255)
                                                                                                    tCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCGGTGCGCgg                                                   <  1:4634037/150‑1 (MQ=255)
                                                                                                      gATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCGGTGCGCGGGt                                                 >  1:2213371/1‑150 (MQ=255)
                                                                                                        tCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCa                                                                            >  2:1645446/1‑121 (MQ=255)
                                                                                                        tCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCa                                                                            <  1:1645446/121‑1 (MQ=255)
                                                                                                          aaCCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCGGTGCGCGGGTATCg                                             >  1:4006618/1‑150 (MQ=255)
                                                                                                           aCCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCGGTGCGCGGGTATCGc                                            >  1:3470922/1‑150 (MQ=255)
                                                                                                              gcgcgcAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCGGTGCGCGGGTATCGcccc                                         >  1:3691665/1‑150 (MQ=255)
                                                                                                                gcgcAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCGGTGCGCGGGTATCGCCCCAg                                       <  2:2600442/150‑1 (MQ=255)
                                                                                                                  gcAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCGGTGCGCGGGTATCGCCCCAGGc                                     >  2:586893/1‑150 (MQ=255)
                                                                                                                     gTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCGGTGCGCGGGTATCGCCCCAGGCCCa                                  <  1:1938951/150‑1 (MQ=255)
                                                                                                                     gTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCGGTGCGCGGGTATCGCCCCAGGCCCa                                  <  2:1920573/150‑1 (MQ=255)
                                                                                                                      ttCAATAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCGGTGCGCGGGTATCGCCCCAGGCCCAg                                 >  1:2876649/1‑150 (MQ=255)
                                                                                                                         aacaaGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCGGTGCGCGGGTATCGCCCCAGGCCCAGCaa                              >  1:1147820/1‑150 (MQ=255)
                                                                                                                          acaaGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCGGTGCGCGGGTATCGCCCCAGGCCCAGCAAc                             >  2:3737776/1‑150 (MQ=255)
                                                                                                                             aGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCGGTGCGCGGGTATCGCCCCAGGCCCAGCAACCCa                          <  2:5344034/150‑1 (MQ=255)
                                                                                                                              gTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCGGTGCGCGGGTATCGCCCCAGGCCCAGCAACCCAg                         >  1:5191404/1‑150 (MQ=255)
                                                                                                                               tACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCGGTGCGCGGGTATCGCCCCAGGCCCAGCAACCCAGt                        <  1:3186047/150‑1 (MQ=255)
                                                                                                                                 ccTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCGGTGCGCGGGTATCGCCCCAGGCCCAGCAACCCAGTGa                      <  1:32845/150‑1 (MQ=255)
                                                                                                                                  cTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCGGTGCGCGGGTATCGCCCCAGGCCCAGCAACCCAGTGAc                     <  1:4068955/150‑1 (MQ=255)
                                                                                                                                   tGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCGGTGCGCGGGTATCGCCCCAGGCCCAGCAACCCAGTGAcc                    <  1:4379298/150‑1 (MQ=255)
                                                                                                                                   tGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCGGTGCGCGGGTATCGCCCCAGGCCCAGCAACCCAGTGAcc                    >  2:5279418/1‑150 (MQ=255)
                                                                                                                                         ttgaaaGCCGTCGGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCGGTGCGCGGGTATCGCCCCAGGCCCAGCAACCCAGt                        >  1:102712/4‑140 (MQ=255)
                                                                                                                                         ttgaaaGCCGTCGGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCGGTGCGCGGGTATCGCCCCAGGCCCAGCAACCCAGTGACCCGATCa              >  1:4456811/4‑150 (MQ=255)
                                                                                                                                          gCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCGGTGCGCGGGTATCGCCCCAGGCCCAGCAACCCAGTGACCCGATCAg             >  1:4362974/1‑150 (MQ=255)
                                                                                                                                          gCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCGGTGCGCGGGTATCGCCCCAGGCCCAGCAACCCAGTGACCCGATCAg             <  1:3372805/150‑1 (MQ=255)
                                                                                                                                           cAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCGGTGCGCGGGTATCGCCCCAGGCCCAGCAACCCAGTGACCCGATCAGc            <  2:38373/150‑1 (MQ=255)
                                                                                                                                            aaaGCCGTCGGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCGGTGCGCGGGTATCGCCCCAGGCCCAGCAACCCAGTGACCCGATCAGCg           >  1:4950369/1‑150 (MQ=255)
                                                                                                                                            aaaGCCGTCGGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCGGTGCGCGGGTATCGCCCCAGGCCCAGCAACCCAGTGACCCGATCAGCg           >  1:4950310/1‑150 (MQ=255)
                                                                                                                                            aaaGCCGTCGGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCGGTGCGCGGGTATCGCCCCAGGCCCAGCAACCCAGTGACCCGATCAGCg           >  1:4950183/1‑150 (MQ=255)
                                                                                                                                             aaGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCGGTGCGCGGGTATCGCCCCAGGCCCAGCAACCCAGTGACCCGATCAGCGa          <  2:2488029/150‑1 (MQ=255)
                                                                                                                                              aGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCGGTGCGCGGGTATCGCCCCAGGCCCAGCAACCCAGTGACCCGATCAGCGAg         >  1:2345208/1‑150 (MQ=255)
                                                                                                                                                    ccGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCGGTGCGCGGGTATCGCCCCAGGCCCAGCAACCCAGTGACCCGATCAGCGAGTTGTTc   <  2:3706345/150‑1 (MQ=255)
                                                                                                                                                    ccGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCGGTGCGCGGGTATCGCCCCAGGCCCAGCAACCCAGTGACCCGATCAGCGAGTTGTTc   <  1:648711/150‑1 (MQ=255)
                                                                                                                                                     cGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCGGTGCGCGGGTATCGCCCCAGGCCCAGCAACCCAGTGACCCGATCAGCGAGTTGTTcc  <  1:3837714/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TCGCGTCACCTGGATTCCGCATATGAGCGATCACTTCGCTACAAACCTCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCGGTGCGCGGGTATCGCCCCAGGCCCAGCAACCCAGTGACCCGATCAGCGAGTTGTTCC  >  NC_002947/586930‑587228

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: