Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,561,719 C→T 6.1% A251T (GCT→ACT) ‡ hisZ ← ATP phosphoribosyltransferase regulatory subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,561,7190CT 6.1% 504.4 / 13.8 163A251T (GCT→ACT) ‡hisZATP phosphoribosyltransferase regulatory subunit
Reads supporting (aligned to +/- strand):  ref base C (85/68);  new base T (7/3);  total (92/71)
Fisher's exact test for biased strand distribution p-value = 5.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 166 total aligned reads displayed.

CGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGCCGGGGCGCGACCGAGGCGCACGCGGGCTTCAGCCAGTACTTCCCGACCACCACACAGCTCGACCAGGGCGCGCAGCATGCTGCCCAGGTCTTTCGGCAGGTCGGCCGTCA  >  NC_002947/5561572‑5561868
                                                                                                                                                   |                                                                                                                                                     
cgtcgtAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGa                                                                                                                                                     <  2:2018272/150‑1 (MQ=255)
cgtcgtAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGa                                                                                                                                                     >  1:2720351/1‑150 (MQ=255)
 gtcgtAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGAt                                                                                                                                                    >  2:5090137/1‑150 (MQ=255)
  gccaatcGATCGCCCAGGGCGGTCGCTTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATc                                                                                                                                                   >  1:3488430/7‑150 (MQ=255)
      aGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCca                                                                                                                                               >  1:3866225/1‑150 (MQ=255)
       gCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCcag                                                                                                                                              >  2:3718938/1‑150 (MQ=255)
         gACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCGCTGATCGcctggc                                                                                                                                            <  1:1633581/150‑5 (MQ=255)
         gACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCcagca                                                                                                                                            <  1:5096036/150‑1 (MQ=255)
          aCCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCcagcag                                                                                                                                           >  2:1983306/1‑150 (MQ=255)
               ccTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCAt                                                                                                                                      >  2:2385904/1‑150 (MQ=255)
                cTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATc                                                                                                                                     <  1:2499414/150‑1 (MQ=255)
                cTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATc                                                                                                                                     <  2:2663180/150‑1 (MQ=255)
                     cgatcgatTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAgcg                                                                                                                                >  1:2313488/1‑150 (MQ=255)
                     cgatcgatTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAgcg                                                                                                                                >  2:4498730/1‑150 (MQ=255)
                       atcgatTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGCCGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGcgcc                                                                                                                              >  1:3262924/1‑150 (MQ=255)
                       atcgatTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGcgcc                                                                                                                              >  1:1274863/1‑150 (MQ=255)
                       atcgatTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGcgcc                                                                                                                              >  1:761359/1‑150 (MQ=255)
                               gCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCAc                                                                                                                      >  2:4979976/1‑150 (MQ=255)
                                 cTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGc                                                                                                                    <  2:2489779/150‑1 (MQ=255)
                                 cTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGc                                                                                                                    >  1:929897/1‑150 (MQ=255)
                                     gCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGc                                                                                                                <  2:5116266/150‑1 (MQ=255)
                                     gCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGc                                                                                                                <  2:2387820/150‑1 (MQ=255)
                                         ggCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGCCggg                                                                                                            <  1:5211738/150‑1 (MQ=255)
                                          gCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCGTGCTGGcggcgc                                                                                                           <  1:559643/150‑6 (MQ=255)
                                          gCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGCCgggg                                                                                                           <  1:5139527/150‑1 (MQ=255)
                                          gCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGCCgggg                                                                                                           <  1:364421/150‑1 (MQ=255)
                                            aCGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCTCGCAGCTCGCCGGGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGCCGGGgcg                                                                                                         >  2:1763997/1‑150 (MQ=255)
                                            aCGAACACGGCGAACACCACGCCCGTGGGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGCCGGGgcg                                                                                                         >  2:3227150/1‑150 (MQ=255)
                                             cGAACACGGCGAACACAACGCCCGTGTGGTAGTTGTAACCACGCAGCTCGACGAGGTCGAAGTACAGCCGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGACATCCAGCGCCGCCAGCACGCTGGCCGGGgcgc                                                                                                        >  1:993222/1‑150 (MQ=255)
                                              gAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGCCGGGGcgcg                                                                                                       >  1:2850618/1‑150 (MQ=255)
                                              gAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGCCGGGGcgcg                                                                                                       >  2:3962231/1‑150 (MQ=255)
                                                acacGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGCCGGGGCGCGAc                                                                                                     >  1:1117082/1‑150 (MQ=255)
                                                acacGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGCCGGGGCGCGAc                                                                                                     <  1:4338815/150‑1 (MQ=255)
                                                acacGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGGCATCCAGCGCCGCCCGCACGCTGGCCGGGGCGCGAc                                                                                                     >  1:3263823/1‑150 (MQ=255)
                                                    ggCGAACACCCCGCCCGTGTGGTAGTGGTACCCCCGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGCTCATCCAGCGCCGCCAGCACGCTGGCCGGGGCGCGACCGAg                                                                                                 <  2:4135565/150‑1 (MQ=255)
                                                     gCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGTGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGCCGGGGCGCGACCGAgg                                                                                                <  1:1502707/150‑1 (MQ=255)
                                                      cGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGCCGGGGCGCGACCGAGgc                                                                                               >  1:2617166/1‑150 (MQ=255)
                                                       gAACACCCCGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGCCGGGGCGCGACCGAGgcg                                                                                              <  2:4030129/150‑1 (MQ=255)
                                                        aaCACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGCCGGGGCGCGACCGAGgcgc                                                                                             >  1:422482/1‑150 (MQ=255)
                                                            ccacGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGCCGGGGCGCGACCGAGGCGCAcgc                                                                                         >  2:3559658/1‑150 (MQ=255)
                                                             cacGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGCCGGGGCGCGACCGAGGCGCAcgcg                                                                                        >  1:5039720/1‑150 (MQ=255)
                                                             cacGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGCCGGGGCGCGACCGAGGCGCAcgcg                                                                                        >  2:4844519/1‑150 (MQ=255)
                                                             cacGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGCCGGGGCGCGACCGAGGCGCAcgcg                                                                                        >  2:2746808/1‑150 (MQ=255)
                                                              acGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGCCGGGGCGCGACCGAGGCGCACGCgg                                                                                       <  2:1728954/150‑1 (MQ=255)
                                                              acGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGCCGGGGCGCGACCGAGGCGCACGCgg                                                                                       >  2:521099/1‑150 (MQ=255)
                                                               cGCCCGTGTGGTAGTGGGAACCACGCAGCTCGCCGAGGTCGAAGGACAGCGGCAGGTCCGGGTAGCGCGACGCCACGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCTCCAGCACGCTGGCCGGGGCGCGACCGAGGCGCACGCtgg                                                                                      <  1:5213026/150‑1 (MQ=255)
                                                                gCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGCCGGGGCGCGACCGAGGCGCACGCGGGc                                                                                     >  1:675657/1‑150 (MQ=255)
                                                                gCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGCCGGGGCGCGACCGAGGCGCACGCGGGc                                                                                     <  2:131246/150‑1 (MQ=255)
                                                                  gggcgtggtagtggtaACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTCGCCGGGGAGCGACCGAGGCGCACGCGGGCtt                                                                                   >  1:3139821/5‑150 (MQ=255)
                                                                    gtgtGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGCCGGGGCGCGACCGAGGCGCACGCGGGCTTCa                                                                                 <  1:1798137/150‑1 (MQ=255)
                                                                    gtgtGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGCCGGGGCGCGACCGAGGCGCACGCGGGCTTCa                                                                                 <  2:4920988/150‑1 (MQ=255)
                                                                      gtggtagtggtaACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGCCGGGGCGCGACCGAGGCGCACGCGGGCTTCAGc                                                                               <  1:2562017/150‑1 (MQ=255)
                                                                      gtggtagtggtaACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGCCGGGGCGCGACCGAGGCGCACGCGGGCTTCAGc                                                                               <  2:4390794/150‑1 (MQ=255)
                                                                      gtggtagtggtaACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGCCGGGGCGCGACCGAGGCGCACGCGGGCTTCAGc                                                                               <  1:4151544/150‑1 (MQ=255)
                                                                         gtagtggtaACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGCCGGGGCGCGACCGAGGCGCACGCGGGCTTCAGCCAg                                                                            >  1:1367321/1‑150 (MQ=255)
                                                                           agtggtaACCCCGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGCCGGGGCGCGACCGAGGCGCACGCGGGCTTCAGCCAGTa                                                                          <  2:1860/150‑1 (MQ=255)
                                                                           agtggtaACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGGAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGCCGGGGCGCGACCGAGGCGCACGCGGGCTTCAGCCAGTa                                                                          <  2:2730969/150‑1 (MQ=255)
                                                                             tggtaACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGTTCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGCCGGGGCGCGACCGAGGCGCACGCGGGCTTCAGCCAGTACt                                                                        >  2:4926230/1‑150 (MQ=255)
                                                                             tggtaACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGc                                                                                                                                                       >  1:3134963/1‑71 (MQ=255)
                                                                             tggtaACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGc                                                                                                                                                       <  2:3134963/71‑1 (MQ=255)
                                                                             tggtaACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTg                                                                                                                  <  2:125603/108‑1 (MQ=255)
                                                                             tggtaACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTg                                                                                                                  >  1:125603/1‑108 (MQ=255)
                                                                               gtaACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGCCGGGGCGCGACCGAGGCGCACGCGGGCTTCAGCCAGTACTTc                                                                      >  1:4105038/1‑150 (MQ=255)
                                                                               gtaACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGCCGGGGCGCGACCGAGGCGCACGCGGGCTTCAGCCAGTACTTc                                                                      >  1:4482502/1‑150 (MQ=255)
                                                                               gtaACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGCCGGGGCGCGACCGAGGCGCACGCGGGCTTCAGCCAGTACTTc                                                                      <  1:1983306/150‑1 (MQ=255)
                                                                                taACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGCCGGGGCGCGACCGAGGCGCACGCGGGCTTCAGCCAGTACTTcc                                                                     >  1:1075581/1‑150 (MQ=255)
                                                                                taACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGCCGGGGCGCGACCGAGGCGCACGCGGGCTTCAGCCAGTACTTcc                                                                     >  2:2862835/1‑150 (MQ=255)
                                                                                taACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGCCGGGGCGCGACCGAGGCGCACGCGGGCTTCAGCCAGTACTTcc                                                                     >  2:5087828/1‑150 (MQ=255)
                                                                                 aaCCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGCCGGGGCGCGACCGAGGCGCACGCGGGCTTCAGCCAGTACTTccc                                                                    <  2:4535235/150‑1 (MQ=255)
                                                                                  aCCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGCCGGGGCGCGACCGAGGCGCACGCGGGCTTCAGCCAGTACTTCCCg                                                                   >  2:1693900/1‑150 (MQ=255)
                                                                                    cACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGCCGGGGCGCGACCGAGGCGCACGCGGGCTTCAGCCAGTACTTCCCGac                                                                 >  1:340727/1‑150 (MQ=255)
                                                                                    cACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGCCGGGGCGCGACCGAGGCGCACGCGGGCTTCAGCCAGTACTTCCCGac                                                                 >  1:712954/1‑150 (MQ=255)
                                                                                    cACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGCCGGGGCGCGACCGAGGCGCACGCGGGCTTCAGCCAGTACTTCCCGac                                                                 >  1:1654946/1‑150 (MQ=255)
                                                                                     ccGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGCCGGGGcgcg                                                                                                       >  2:808584/2‑111 (MQ=255)
                                                                                     aCGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCGGGCCGGGGCGCGACCGAGGCGCACGCGGGCTTCAGCCAGTACTTCCCGacc                                                                <  1:1467575/150‑1 (MQ=255)
                                                                                              gCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGCCGGGGCGCGACCGAGGCGCACGCGGGCTTCAGCCAGTACTTCCCGACCACCACACAg                                                       <  2:214196/150‑1 (MQ=255)
                                                                                              gCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGCCGGGGCGCGACCGAGGCGCACGCGGGCTTCAGCCAGTACTTCCCGACCACCACACAg                                                       >  1:50987/1‑150 (MQ=255)
                                                                                                        aaGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGCCGGGGCGCGACCGAGGCGCACGCGGGCTTCAGCCAGTACTTCCCGACCACCACACAGCTCGACCAgg                                             <  2:670335/150‑1 (MQ=255)
                                                                                                        aaGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGCCGGGGCGCGACCGAGGCGCACGCGGGCTTCAGCCAGTACTTCCCGACCACCACACAGCTCGACCAgg                                             >  2:4798111/1‑150 (MQ=255)
                                                                                                           tACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGCCGGGGCGCGACCGAGGCGCACGCGGGCTTCAGCCAGTACTTCCCGACCACCACACAGCTCGACCAGGgcg                                          <  2:1367321/150‑1 (MQ=255)
                                                                                                               gCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGCCGGGGCGCGACCGAGGCGCACGCGGGCTTCAGCCAGTACTTCCCGACCACCACACAGCTCGACCAGGGCGCgca                                      >  1:3433660/1‑150 (MQ=255)
                                                                                                                    aGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGCCGGGGCGCGACCGAGGCGCACGCGGGCTTCAGCCAGTACTTCCCGACCACCACACAGCTCGACCAGGGCGCGCAGCAtg                                 <  2:3207479/150‑1 (MQ=255)
                                                                                                                    aGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGCCGGGGCGCGACCGAGGCGCACGCGGGCTTCAGCCAGTACTTCCCGACCACCACACAGCTCGACCAGGGCGCGCAGCAtg                                 >  1:3313461/1‑150 (MQ=255)
                                                                                                                    aGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGCCGGGGCGCGACCGAGGCGCACGCGGGCTTCAGCCAGTACTTCCCGACCACCACACAGCTCGACCAGGGCGCGCAGCAtg                                 <  2:2617166/150‑1 (MQ=255)
                                                                                                                           ggTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGCCGGGGCGCGACCGAGGCGCACGCGGGCTTCAGCCAGTACTTCCCGACCACCACACAGCTCGACCAGGGCGCGCAGCATGCTGCCCa                          <  2:918251/150‑1 (MQ=255)
                                                                                                                            gTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGCCGGGGCGCGACCGAGGCGCACGCGGGCTTCAGCCAGTACTTCCCGACCACCACACAGCTCGACCAGGGCGCGCAGCATGCTGCCCAg                         <  2:1075581/150‑1 (MQ=255)
                                                                                                                                 gcgACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGCCGGGGCGCGACCGAGGCGCACGCGGGCTTCAGCCAGTACTTCCCGACCACCACACAGCTCGACCAGGGCGCGCAGCATGCTGCCCAGGTCtt                    >  2:3812345/1‑150 (MQ=255)
                                                                                                                                 gcgACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGCCGGGGCGCGACCGAGGCGCACGCGGGCTTCAGCCAGTACTTCCCGACCACCACACAGCTCGACCAGGGCGCGCAGCATGCTGCCCAGGTCtt                    >  1:1439256/1‑150 (MQ=255)
                                                                                                                                 gcgACGCCAGGAGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGCCGGGGCGCGACCGAGGCGCACGCGGGCTTCAGCCAGTACTTCCCGACCACCACACAGCTCGACCAgg                                             <  2:4678953/125‑1 (MQ=255)
                                                                                                                                 gcgACGCCAGGAGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGCCGGGGCGCGACCGAGGCGCACGCGGGCTTCAGCCAGTACTTCCCGACCACCACACAGCTCGACCAgg                                             >  1:4678953/1‑125 (MQ=255)
                                                                                                                                   gACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGCCGGGGCGCGACCGAGGCGCACGCGGGCTTCAGCCAGTACTTCCCGACCACCACACAGCTCGACCAGGGCGCGCAGCATGCTGCCCAGGTCTTTc                  <  1:3718938/150‑1 (MQ=255)
                                                                                                                                   gACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGCCGGGGCGCGACCGAGGCGCACGCGGGCTTCAGCCAGTACTTCCCGACCACCACACAGCTCGACCAGGGCGCGCAGCATGCTGCCCAGGTCTTTc                  <  1:5006930/150‑1 (MQ=255)
                                                                                                                                   gACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGCCGGGGCGCGACCGAGGCGCACGCGGGCTTCAGCCAGTACTTCCCGACCACCACACAGCTCGACCAGGGCGCGCAGCATGCTGCCCAGGTCTTTc                  <  2:5427859/150‑1 (MQ=255)
                                                                                                                                      tgctggCGATCGCTGATCGCCTGCAGGTCATCCAGCGCCGCCAGCACGCTGGCCGGGGCGCGACCGAGGCGCACGCGGGCTTCAGCCAGTACTTCCCGACCACCACACAGCTCGACCAGGGCGCGCAGCATGCTGCCCAGGTCTTTCGGc               >  2:5212854/5‑150 (MQ=255)
                                                                                                                                         tggCGATCGCTGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGCCGGGGCGCGACCGAGGCGCACGCGGGCTTCAGCCAGTACTTCCCGACCACCACACAGCTCGACCAGGGCGCGCAGCATGCTGCCCAGGTCTTTCGGCAgg            >  2:2400514/2‑150 (MQ=255)
                                                                                                                                         tggCGATCGCTGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGCCGGGGCGCGACCGAGGCGCACGCGGGCTTCAGCCAGTACTTCCCGACCACCACACAGCTCGACCAGGGCGCGCAGCATGCTGCCCAGGTCTTTCGGCAgg            >  2:2135212/2‑150 (MQ=255)
                                                                                                                                         tggCGATCGCTGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGCCGGGGCGCGACCGAGGCGCACGCGGGCTTCAGCCAGTACTTCCCGACCACCACACAGCTCGACCAGGGCGCGCAGCATGCTGCCCAGGTCTTTCGGCAgg            >  2:2147650/2‑150 (MQ=255)
                                                                                                                                           gCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGCCGGGGCGCGACCGAGGCGCACGCGGGCTTCAGCCAGTACTTCCCGACCACCACACAGCTCGACCAGGGCGCGCAGCATGCTGCCCAGGTCTTTCGGCAGGTc          >  2:4213361/1‑150 (MQ=255)
                                                                                                                                                 aGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGCCGGGGCGCGACCGAGGCGCACGCGGGCTTCAGCCAGTACTTCCCGACCACCACACAGCTCGACCAGGGCGCGCAGCATGCTGCCCAGGTCTTTCGGCAGGTCGGCCGt    >  2:1597046/1‑150 (MQ=255)
                                                                                                                                                   cGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGCCGGGGCGCGACCGAGGCGCACGCGGGCTTCAGCCAGTACTTCCCGACCACCACACAGCTCGACCAGGGCGCGCAGCATGCTGCCCAGGTCTTTCGGCAGGTCGGCCGTCa  >  1:5089344/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGCCGGGGCGCGACCGAGGCGCACGCGGGCTTCAGCCAGTACTTCCCGACCACCACACAGCTCGACCAGGGCGCGCAGCATGCTGCCCAGGTCTTTCGGCAGGTCGGCCGTCA  >  NC_002947/5561572‑5561868

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