Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,318,446 G→T 8.8% A208A (GCC→GCA) ‡ ilvC ← ketol‑acid reductoisomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,318,4460GT 8.8% 621.5 / 36.2 192A208A (GCC→GCA) ‡ilvCketol‑acid reductoisomerase
Reads supporting (aligned to +/- strand):  ref base G (90/85);  new base T (13/4);  total (103/89)
Fisher's exact test for biased strand distribution p-value = 7.25e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 195 total aligned reads displayed.

CACCGTATTCGGCGTTGTTGGAGATCGAGTAGTTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCGTT  >  NC_002947/5318298‑5318595
                                                                                                                                                    |                                                                                                                                                     
cACCGTATTCGGCGTTGTTGGAGATCGAGTAGTTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCgg                                                                                                                                                      >  2:5187505/1‑150 (MQ=255)
cACCGTATTCGGCGTTGTTGGAGATCGAGTAGTTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCgg                                                                                                                                                      <  2:2766884/150‑1 (MQ=255)
 aCCGTATTCGGCGTTGTTGGAGATCGAGTAGTTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGc                                                                                                                                                     >  2:1679591/1‑150 (MQ=255)
 aCCGTATTCGGCGTTGTTGGACATCGAGTAGTTCATGTAGGCGAAGCCGCCTTCGTACATGAGGTCTACCATCAGCATCAGCACGTGCAGGCACACGAAGTAGGCCAAATCCGGCCCGTAGCCAGCTTCGACCAGCGTTTCGAAACCGGc                                                                                                                                                     >  1:928939/1‑150 (MQ=255)
   cGTATTCGGCGTTGTTGGAGATCGAGTAGTTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCtt                                                                                                                                                   >  1:2990578/1‑150 (MQ=255)
   cGTATTCGGCGTTGTTGGAGATCGAGTAGTTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCtt                                                                                                                                                   >  1:2713580/1‑150 (MQ=255)
    gcataccgCGTCGTTGGAGAGCGAGCACCGCCTGTTGGCGAAGCCGCCCACGCCCATGAGGTCGACGATCAGCGTCAGCTCGTGCAGGCACGCCAAGAAGGCCAGTTCCGTCGCGTAGCCAGCTTCGACCAGGGTATCGAAACCGGCttt                                                                                                                                                  <  2:3577818/143‑1 (MQ=255)
    gTATTCGGCGTTGTTGGAGATCGAGTAGTTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCttt                                                                                                                                                  >  1:2579641/1‑150 (MQ=255)
      aTTCGGCGTTGTTGGAGATCGAGTAGTTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGATCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGa                                                                                                                                                <  2:3675935/150‑1 (MQ=255)
       ttCGGCGTTGTTGGAGATCGAGTAGTTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGAc                                                                                                                                               <  1:2844013/150‑1 (MQ=255)
       ttCGGCGTTGTTGGAGATCGAGTAGTTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGAc                                                                                                                                               >  2:1247356/1‑150 (MQ=255)
       ttCGGCGTTGTTGGAGATCGAGTAGTTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGAc                                                                                                                                               >  2:2074943/1‑150 (MQ=255)
        tCGGCGTTGTTGGAGATCGAGTAGTTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGAACAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGAcc                                                                                                                                              <  2:930441/150‑1 (MQ=255)
          ggCGTTGTTGGAGATCGAGTAGTTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCCGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAg                                                                                                                                            >  1:3123344/1‑150 (MQ=255)
           gCGTTGTTGGAGATCGAGTAGTTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGc                                                                                                                                           <  1:4390429/150‑1 (MQ=255)
           gCGTTGTTGGAGATCGAGTAGTTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTACGACCAGGGTTTCGAAACCGGCTTTGACCAGc                                                                                                                                           >  2:618358/1‑150 (MQ=255)
                gttgGAGATCGAGTAGTTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGCTCAGCTTCAGCACGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGAc                                                                                                                                      <  2:1185388/150‑1 (MQ=255)
                   ggAGATCGAGTAGTTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGt                                                                                                                                   >  1:2276987/1‑150 (MQ=255)
                      gaTCGAGTAGTTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTAcc                                                                                                                                >  2:3058779/1‑150 (MQ=255)
                      gaTCGAGTAGTTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTAcc                                                                                                                                <  1:982093/150‑1 (MQ=255)
                          gAGTAGTTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCca                                                                                                                            >  2:2873147/1‑150 (MQ=255)
                          gAGTAGTTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCca                                                                                                                            >  2:5173927/1‑150 (MQ=255)
                               gtTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAg                                                                                                                       >  2:4997096/1‑150 (MQ=255)
                                 tCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGaa                                                                                                                     <  1:1139063/150‑1 (MQ=255)
                                   aTGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACa                                                                                                                   <  1:1546770/150‑1 (MQ=255)
                                   aTGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACa                                                                                                                   <  2:116588/150‑1 (MQ=255)
                                     gTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGc                                                                                                                 <  1:1612041/150‑1 (MQ=255)
                                      ttGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGcc                                                                                                                >  1:1071164/1‑150 (MQ=255)
                                      ttGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGcc                                                                                                                <  1:3487039/150‑1 (MQ=255)
                                          cGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCt                                                                                                            >  2:484813/1‑150 (MQ=255)
                                           gATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTc                                                                                                           >  2:8404/1‑150 (MQ=255)
                                                 gccTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGaa                                                                                                     >  1:176049/1‑150 (MQ=255)
                                                    ttCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAg                                                                                                  <  1:3378314/150‑1 (MQ=255)
                                                    ttCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAg                                                                                                  >  2:4332755/1‑150 (MQ=255)
                                                      cGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGt                                                                                                <  1:3823723/150‑1 (MQ=255)
                                                      cGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGt                                                                                                <  1:3095211/150‑1 (MQ=255)
                                                        tagatcaGGTCTACGATGAGCTTCATTTCGTGCAGGCACTCGAAGTAAGCCATTTCTGGTGCGTAACCAGCTTCGACCAGAACTTCAAAACCTGTGCGAATGAGTTCAACAACGCCACCGCAAAGAACAGCCTGCTCACCGAACAGGtca                                                                                              >  2:2376485/7‑149 (MQ=255)
                                                           atcaGGTCTACGATGAGCTTCATTTCGTGCAGGCACTCGAAGTAAGCCATTTCTGGTGCGTAACCAGCTTCGACCAGAACTTCAAAACCTGTGCGAATGAGTTCAACAACGCCACCGCAAAGAACAGCCTGCTCACCGAACAgg                                                                                                 <  1:3545579/141‑1 (MQ=255)
                                                           atcaGGTCTACGATGAGCTTCATTTCGTGCAGGCACCCGAAGTAAGCCATTTCTGGTGCGTAACCAGCTTCGACCAGAACTTCAAAACCTGTGCGAATGAGTTCAACAACGCCACCGCAAAGAACAGCCTGCTCACCGAACAGGtcagtt                                                                                           >  2:3392950/4‑150 (MQ=255)
                                                             caGGTCTACGATGAGCTTCATTTCGTGCAGGCACTCGAAGTAAGCCATTTCTGGTGCGTAACCAGCTTCGACCAGAACTTCAAAACCTGTGCGAATGAGTTCAACAACGCCACCGCAAAGAACAGCCTGCTCACCGAACAGGtcagtttc                                                                                         <  2:2398123/149‑1 (MQ=255)
                                                              aGGTCTACGATGAGCTTCATTTCGTGCAGGCACTCGAAGTAAGCCATTTCTGGTGCGTAACCAGCTTCGACCAGAACTTCAAAACCTGTGCGAATGAGTTCAACAACGCCACCGCAAAGAACAGCCTGCTCACCGAACAGGtcagtttcg                                                                                        >  2:4724473/1‑150 (MQ=255)
                                                              aGGTCTACGATGAGCTTCATTTCGTGCAGGCACTCGAAGTAAGCCATTTCTGGTGCGTAACCAGCTTCGACCAGAACTTCAAAACCTGTGCGAATGAGTTCAACAACGCCACCGCAAAGAACAGCCTGCTCACCGAACAGGtcagtttcg                                                                                        >  2:1985383/1‑150 (MQ=255)
                                                              aGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGtcggtttcg                                                                                        <  2:5365406/150‑1 (MQ=255)
                                                              aGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGtcggtttcg                                                                                        <  2:1888686/150‑1 (MQ=255)
                                                              aGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTACGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGtcggtttcg                                                                                        <  2:5090922/150‑1 (MQ=255)
                                                                gTCTACGATGAGCTTCATTTCGTGCAGGCACTCGAAGTAAGCCATTTCTGGTGCGTAACCAGCTTCGACCAGAACTTCAAAACCTGTGCGAATGAGTTCAACAACGCCACCGCAAAGAACAGCCTGCTCACCGAACAGGtcagtttcgca                                                                                      >  2:1449409/1‑148 (MQ=255)
                                                                   tACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCCGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAAc                                                                                                                    <  2:5272556/117‑1 (MQ=255)
                                                                   tACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGACCAGCCTGCTCACCGAACAGGTCggtttcggtttc                                                                                   <  2:3279630/150‑1 (MQ=255)
                                                                   tACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAAc                                                                                                                    >  1:5272556/1‑117 (MQ=255)
                                                                   tACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCggtttcggtttc                                                                                   >  2:4448685/1‑150 (MQ=255)
                                                                      gATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTc                                                                                >  1:4247407/1‑150 (MQ=255)
                                                                           gcttcatttCGTGCAGGCACTCGAAGTAAGCCATTTCTGGTGCGTAACCAGCTTCGACCAGAACTTCAAAACCTGTGCGAATGAGTTCAACAACGCCACCGCAAAGAACAGCCTGCTCACCGAACAGGTCAGTTTCGCACTCTTCCTTGa                                                                           >  2:3812120/1‑150 (MQ=255)
                                                                           gcttcagctCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGa                                                                           <  1:3058779/150‑1 (MQ=255)
                                                                            cttcagctCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGaa                                                                          >  2:2526920/1‑150 (MQ=255)
                                                                             ttcatttCGTGCAGGCACTCGAAGTAAGCCATTTCTGGTGCGTAACCAGCTTCGACCAGAACTTCAAAACCTGTGCGAATGAGTTCAACAACGCCACCGCAAAGAACAGCCTGCTCACCGAACAGGTCAGTTTCGCACTCTTGCTTGAAg                                                                         >  1:3377839/7‑150 (MQ=255)
                                                                               cagctCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGCACAGGTCGGTTTCGGTTTCGTCCTTGAAggt                                                                       <  2:2030001/150‑1 (MQ=255)
                                                                               cagctCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAggt                                                                       <  1:2705454/150‑1 (MQ=255)
                                                                                 gctCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGATACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACCGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAggtgg                                                                     <  2:3091606/150‑1 (MQ=255)
                                                                                 gctCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAggtgg                                                                     <  2:156972/150‑1 (MQ=255)
                                                                                  ctCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAggtggt                                                                    >  2:5445438/1‑150 (MQ=255)
                                                                                   tCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGtt                                                                   >  2:2601683/1‑150 (MQ=255)
                                                                                       gCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCga                                                               <  1:2767457/150‑1 (MQ=255)
                                                                                        cAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCgat                                                              >  1:1825654/1‑150 (MQ=255)
                                                                                        cAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCgat                                                              <  1:5187505/150‑1 (MQ=255)
                                                                                           gCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCgatgat                                                           >  2:2420655/1‑150 (MQ=255)
                                                                                               tCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTTGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCg                                                       >  2:2976929/1‑150 (MQ=255)
                                                                                               tCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGtcggtttcg                                                                                        >  1:5171685/1‑117 (MQ=255)
                                                                                               tCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCg                                                       >  2:745997/1‑150 (MQ=255)
                                                                                               tCGAAGTAAGCCATTTCTGGTGCGTAACCAGCTTCGACCAGAACTTCAAAACCTGTGCGAATGAGTTCAACAACGCCACCGCAAAGAACAGCCTGCTCACCGAACAGGTCAGTTTCGCACTCTTCCTTGAAGGTGGTTTCaataacgccc                                                       >  2:2851012/1‑143 (MQ=255)
                                                                                                cGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCgg                                                      <  2:2114364/150‑1 (MQ=255)
                                                                                                   aGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTAc                                                   <  1:1976677/150‑1 (MQ=255)
                                                                                                   aGTAGGCCAGGTCCGGCGCGCACCCAGCTTCGACCAGGGTATCGAAACCGGCTTTGACCAGCTCGACGGTCCCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGCGTCGGTTTCCTCCTTGAAGGTGGTTTCGATGATACCGGTAc                                                   <  2:2693104/150‑1 (MQ=255)
                                                                                                        gCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGcca                                              >  1:1914530/1‑150 (MQ=255)
                                                                                                          cATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGccacc                                            >  2:4330766/1‑150 (MQ=255)
                                                                                                          cATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGccacc                                            >  1:1326425/1‑150 (MQ=255)
                                                                                                           aTTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGccacca                                           >  1:4300867/1‑150 (MQ=255)
                                                                                                            tttCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGtcggtttcgg                                                                                       >  1:4528768/1‑105 (MQ=255)
                                                                                                            tttCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGtcggtttcgg                                                                                       <  2:4528768/105‑1 (MQ=255)
                                                                                                                    gcgTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGccc                                  <  1:2045342/150‑1 (MQ=255)
                                                                                                                       tAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGaa                               <  2:178253/150‑1 (MQ=255)
                                                                                                                          ccAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCa                            >  1:3334651/1‑150 (MQ=255)
                                                                                                                           cAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCAt                           >  1:2508607/1‑150 (MQ=255)
                                                                                                                           cAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCAt                           <  1:1206640/150‑1 (MQ=255)
                                                                                                                            aGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATa                          >  1:4359984/1‑150 (MQ=255)
                                                                                                                             gCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAg                         <  1:1637485/150‑1 (MQ=255)
                                                                                                                             gCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAg                         <  1:2500454/150‑1 (MQ=255)
                                                                                                                               ttCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCAGAGGa                       >  2:3430821/1‑150 (MQ=255)
                                                                                                                                  gACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCAAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAg                    <  2:2276987/150‑1 (MQ=255)
                                                                                                                                   aCCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGt                   <  1:1247356/150‑1 (MQ=255)
                                                                                                                                    ccAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTg                                                                                                              >  1:857086/1‑58 (MQ=255)
                                                                                                                                     cAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGc                 >  2:4995417/1‑150 (MQ=255)
                                                                                                                                       gggTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGa               >  1:2423158/1‑150 (MQ=255)
                                                                                                                                          tttCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGt            >  2:3729820/1‑150 (MQ=255)
                                                                                                                                          tttCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGt            >  1:1053168/1‑150 (MQ=255)
                                                                                                                                              gAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCtt        >  1:4781587/1‑150 (MQ=255)
                                                                                                                                               aaaCCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTg       >  2:5155450/1‑150 (MQ=255)
                                                                                                                                               aaaCCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTg       <  1:976972/150‑1 (MQ=255)
                                                                                                                                                 aCCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGc     <  1:1679591/150‑1 (MQ=255)
                                                                                                                                                   cGGCTTTGACCAGCTCGACGGTACCGCCACACTGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCTAAGCATAGGACAGTGCGACGTTCttttcgt   <  2:928939/150‑6 (MQ=255)
                                                                                                                                                    ggCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCGtt  <  2:3505792/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CACCGTATTCGGCGTTGTTGGAGATCGAGTAGTTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCGTT  >  NC_002947/5318298‑5318595

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: