Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,318,371 C→G 8.2% L233L (CTG→CTC ilvC ← ketol‑acid reductoisomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,318,3710CG 8.2% 515.7 / 28.0 159L233L (CTG→CTCilvCketol‑acid reductoisomerase
Reads supporting (aligned to +/- strand):  ref base C (73/73);  new base G (8/5);  total (81/78)
Fisher's exact test for biased strand distribution p-value = 5.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.60e-01

Only 100 of 178 total aligned reads displayed.

TCCTGGATGCGCTTCAGAGCATTGCGCATGGCCTTGCGGGATTCTTCGTTGATGACTTCCGGGCCAGTGACGTACTCACCGTATTCGGCGTTGTTGGAGATCGAGTAGTTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCT  >  NC_002947/5318222‑5318519
                                                                                                                                                     |                                                                                                                                                     
tCCTGGATGCGCGTCAGAGCATTGCGCATGGCCTTGCGGGATTCTTCGTTGATGACTTCCGGGCCAGTGACGTACTCACCGTATTCGGCGTTGTTGGAGATCGAGTAGTTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGAtc                                                                                                                                                       <  2:73433/150‑1 (MQ=255)
  cTGGATGCGCTTCAGAGCATTGCGCATGGCCTTGCGGGATTCTTCGTTGATGACTTCCGGGCCAGTGACGTACTCACCGTATTCGGCGTTGTTGGAGATCGAGTAGTTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGAtcag                                                                                                                                                     >  1:4830273/1‑150 (MQ=255)
  cTGGATGCGCTTCAGAGCATTGCGCATGGCCTTGCGGGATTCTTCGTTGATGACTTCCGGGCCAGTGACGTACTCACCGTATTCGGCGTTGTTGGAGATCGAGTAGTTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGAtcag                                                                                                                                                     >  1:2693104/1‑150 (MQ=255)
   tGGATGCGCTTCAGAGCATTGCGCATGGCCTTGCGGGATTCTTCGTTGATGACTTCCGGGCCAGTGACGTACTCACCGTATTCGGCGTTGTTGGAGATCGAGTAGTTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGAtcagc                                                                                                                                                    <  1:5409846/150‑1 (MQ=255)
     gATGCGCTTCAGAGCATTGCGCATGGCCTTGCGGGATTCTTCGTTGATGACTTCCGGGCCAGTGACGTACTCACCGTATTCGGCGTTGTTGGAGATCGAGTAGTTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGAtcagctt                                                                                                                                                  <  2:3709870/150‑1 (MQ=255)
     gATGCGCTTCAGAGCATTGCGCATGGCCTTGCGGGATTCTTCGTTGATGACTTCCGGGCCAGTGACGTACTCACCGTATTCGGCGTTGTTGGAGATCGAGTAGTTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGAtcagctt                                                                                                                                                  >  1:178253/1‑150 (MQ=255)
      aTGCGCTTCAGAGCATTGCGCATGGCCTTGCGGGATTCTTCGTTGATGACTTCCGGGCCAGTGACGTACTCACCGTATTCGGCGTTGTTGGAGATCGAGTAGTTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGAtcagcttc                                                                                                                                                 <  2:1880980/150‑1 (MQ=255)
      aTGCGCTTCAGAGCATTGCGCATGGCCTTGCGGGATTCTTCGTTGATGACTTCCGGGCCAGTGACGTACTCACCGTATTCGGCGTGGTTGGAGATCGAGTAGTTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGTtcagcttc                                                                                                                                                 >  2:1637259/1‑150 (MQ=255)
        gcgcTTCAGAGCATTGCGCATGGCCGTGCGGGATTCTTCGTTGATGACTTCCGGGCCAGTGACGTACTCACCGTATTCGGCGTTGTTGGAGATCGAGTAGTTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGAtcagcttcag                                                                                                                                               <  1:4635575/150‑1 (MQ=255)
           cTTCAGAGCATTGCGCATGGCCTTGCGGGATTCTTCGTTGATGACTTCCGGGCCAGTGACGTACTCACCGTATTCGGCGTTGTTGGAGATCGAGTAGTTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTc                                                                                                                                            <  2:2413936/150‑1 (MQ=255)
            ttCAGAGCATTGCGCATGGCCTTGCGGGATTCTTCGTTGATGACTTCCGGGCCAGTGACGTACTCACCGTATTCGGCGTTGTTGGAGATCGAGTAGTTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCg                                                                                                                                           >  2:3487039/1‑150 (MQ=255)
                gagCATTGCGCATGGCCTTGCGGGATTCTTCGTTGATGACTTCCGGGCCAGTGACGTACTCACCGTATTCGGCGTTGTTGGAGATCGAGTAGTTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCa                                                                                                                                       <  2:3070531/150‑1 (MQ=255)
                  gCATTGCGCATGGCCTTGCGGGATTCTTCGTTGATGACTTCCGGGCCAGTGACGTACTCACCGTATTCGGCGTTGTTGGAGATCGAGTAGTTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAgg                                                                                                                                     <  1:1766783/150‑1 (MQ=255)
                   cATTGCGCATGGCCTTGCGGGATTCTTCGTTGATGACTTCCGGGCCAGTGACGTACTCACCGTATTCGGCGTTGTTGGAGATCGAGTAGTTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGc                                                                                                                                    >  2:2045342/1‑150 (MQ=255)
                     ttGCGCATGGCCTTGCGGGATTCTTCGTTGATGACTTCCGGGCCAGTGACGTACTCACCGTATTCGGCGTTGTTGGAGATCGAGTAGTTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCAc                                                                                                                                  >  2:2757706/1‑150 (MQ=255)
                     ttGCGCATGGCCTTGCGGGATTCTTCGTTGATGACTTCCGGGCCAGTGACGTACTCACCGTATTCGGCGTTGTTGGAGATCGAGTAGTTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCAc                                                                                                                                  <  2:509644/150‑1 (MQ=255)
                     ttGCGCATGGCCTTGCGGGATTCTTCGTTGATGACTTCCGGGCCAGTGACGTACTCACCGTATTCGGCGTTGTTGGAGATCGAGTAGTTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCAc                                                                                                                                  >  2:2500454/1‑150 (MQ=255)
                           aTGGCCTTGCGGGATTCTTCGTTGATGACTTCCGGGCCAGTGACGTACTCACCGTATTCGGCGTTGTTGGAGATCGAGTAGTTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAg                                                                                                                            >  1:2114364/1‑150 (MQ=255)
                             ggCCTTGCGGGATTCTTCGTTGATGACTTCCGGGCCAGTGACGTACTCACCGTATTCGGCGTTGTTGGAGATCGAGTAGTTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTa                                                                                                                          <  1:2026343/150‑1 (MQ=255)
                              gCCTTGCGGGATTCTTCGTTGATGACTTCCGGGCCAGTGACGTACTCACCGTATTCGGCGTTGTTGGAGATCGAGTAGTTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAg                                                                                                                         <  1:2880803/150‑1 (MQ=255)
                               ccTTGCGGGATTCTTCGTTGATGACTTCCGGGCCAGTGACGTACTCACCGTATTCGGCGTTGTTGGAGATCGAGTAGTTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAgg                                                                                                                        <  2:2025421/150‑1 (MQ=255)
                                cTTGCGGGATTCTTCGTTGATGACTTCCGGGCCAGTGACGTACTCACCGTATTCGGCGTTGTTGGAGATCGAGTAGTTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGc                                                                                                                       <  2:4839731/150‑1 (MQ=255)
                                   gCGGGATTCTTCGTTGATGACTTCCGGGCCAGTGACGTACTCACCGTATTCGGCGTTGTTGGAGATCGAGTAGTTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCAt                                                                                                                    >  2:91738/1‑150 (MQ=255)
                                    cGGGATTCTTCGTTGATGACTTCCGGGCCAGTGACGTACTCACCGTATTCGGCGTTGTTGGAGATCGAGTAGTTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCAtt                                                                                                                   >  2:4269499/1‑150 (MQ=255)
                                     gggATTCTTCGTTGATGACTTCCGGGCCAGTGACGTACTCACCGTATTCGGCGTTGTTGGAGATCGAGTAGTTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCAttt                                                                                                                  >  2:1206640/1‑150 (MQ=255)
                                         ttcttcGTTGATGACTTCCGGGCCAGTGACGTACTCACCGTATTCGGCGTTGTTGGAGATCGAGTAGTTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCgg                                                                                                              <  1:1469641/150‑1 (MQ=255)
                                            ttcGTTGATGACTTCCGGGCCAGTGACGTACTCACCGTATTCGGCGTTGTTGGAGATCGAGTAGTTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGgcgc                                                                                                           >  1:3279630/1‑150 (MQ=255)
                                             tcGTTGATGACTTCCGGGCCAGTGACGTACTCACCGTATTCGGCGTTGTTGGAGATCGAGTAGTTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGcgcg                                                                                                          >  2:2767457/1‑150 (MQ=255)
                                               gTTGATGACTTCCGGGCCAGTGACGTACTCACCGTATTCGGCGTTGTTGGAGATCGAGTAGTTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATATCCGGcgcgaa                                                                                                        >  1:116588/1‑148 (MQ=255)
                                                 tgatgaCTTCCGGGCCAGTGACGTACTCACCGTATTCGGCGTTGTTGGAGATCGAGTAGTTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGc                                                                                                      >  2:4634147/1‑150 (MQ=255)
                                                  gatgaCTTCCGGGCCAGTGACGTACTCACCGTATTCGGCGTTGTTGGAGATCGAGTAGTTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGcc                                                                                                     >  2:183896/1‑150 (MQ=255)
                                                  gatgaCTTCCGGGCCAGTGACGTACTCACCGTATTCGGCGTTGTTGGAGATCGAGTAGTTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGcc                                                                                                     <  2:4790820/150‑1 (MQ=255)
                                                      atgcgcGGACCAGTGACGTATTCGTCCCACTCAGCAGTGTTGGAGATCGAGTAGTTCATGTTCGCAATGCCGCCTTCATAGATCAGGTCTACGATGAGCTTCATTTCGTGCAGGCACTCGAAGTAAGCCATTTCTGGTGCGTAACCAGCt                                                                                                 >  2:1970114/6‑150 (MQ=255)
                                                        ttCCGGGCCAGTGACGTACTCACCGTATTCGGCGTTGTTGGAGATCGAGTAGTTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTc                                                                                               >  2:4154795/1‑150 (MQ=255)
                                                        gcgcGGACCAGTGACGTATTCGCCCCACTCAGCAGTGTTGGAGATCGAGTAGTTCATGTTCGCAATGCCGCCTTCATAGATCAGGTCTACGATGAGCTTCATTTCGTGCAGGCACTCGAAGTAAGCCATTTCTGGTGCGTAACCAGCTTc                                                                                               <  1:3806661/147‑1 (MQ=255)
                                                           cGGGCCAGTGACGTACTCACCGTATTCGGCGTTGTTGGAGATCGAGTAGTTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGAc                                                                                            <  1:3569508/150‑1 (MQ=255)
                                                           cGGACCAGTGACGTATTCGCCCCACTCAGCAGTGTTGGAGATCGAGTAGTTCATGTTCGCAATGCCGCCTTCATAGATCAGGTCTACGATGAGCTTCATTTCGTGCAGGCACTCGAAGTAAGCCATTTCTGGTGCGTAACCAGCTTCGAc                                                                                            <  1:4736552/150‑1 (MQ=255)
                                                             ggCCAGTGACGTACTCACCGTATTCGGCGTTGTTGGAGATCGAGTAGTTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCa                                                                                          <  2:387634/150‑1 (MQ=255)
                                                             ggCCAGTGACGTACTCACCGTATTCGGCGTTGTTGGAGATCGAGTAGTTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCa                                                                                          >  1:3091523/1‑150 (MQ=255)
                                                               ccAGTGACGTACTCACCGTATTCGGCGTTGTTGGAGATCGAGTAGTTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAgg                                                                                        <  2:3022230/150‑1 (MQ=255)
                                                                 aGTGACGTACTCACCGTATTCGGCGTTGTTGGAGATCGAGTAGTTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGt                                                                                      <  2:2425867/150‑1 (MQ=255)
                                                                    gACGTACTCACCGTATTCGGCGTTGTTGGAGATCGAGTAGTTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTc                                                                                   <  1:5040462/150‑1 (MQ=255)
                                                                             cccGTATTCGGCGTTGTTGGAGATCGAGTAGTTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGc                                                                          >  2:976972/2‑150 (MQ=255)
                                                                             aCCGTATTCGGCGTTGTTGGAGATCGAGTAGTTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGc                                                                          >  2:1679591/1‑150 (MQ=255)
                                                                             aCCGTATTCGGCGTTGTTGGACATCGAGTAGTTCATGTAGGCGAAGCCGCCTTCGTACATGAGGTCTACCATCAGCATCAGCACGTGCAGGCACACGAAGTAGGCCAAATCCGGCCCGTAGCCAGCTTCGACCAGCGTTTCGAAACCGGc                                                                          >  1:928939/1‑150 (MQ=255)
                                                                              ccGTATTCGGCGTTGTTGGAGATCGAGTAGTTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTc                                                                                                                                            >  1:707455/1‑83 (MQ=255)
                                                                              ccGTATTCGGCGTTGTTGGAGATCGAGTAGTTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTc                                                                                                                                            <  2:707455/83‑1 (MQ=255)
                                                                              ccGTATTCGGCGTTGTTGGAGATCGAGTAGTTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCt                                                                         <  2:2869922/150‑1 (MQ=255)
                                                                               cGTATTCGGCGTTGTTGGAGATCGAGTAGTTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCtt                                                                        >  1:848359/1‑150 (MQ=255)
                                                                               cGTATTCGGCGTTGTTGGAGATCGAGTAGTTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCtt                                                                        >  1:2713580/1‑150 (MQ=255)
                                                                                  aTTCGGCGTTGTTGGAGATCGAGTAGTTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGATCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGa                                                                     <  2:3675935/150‑1 (MQ=255)
                                                                                   ttCGGCGTTGTTGGAGATCGAGTAGTTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGAc                                                                    >  2:1247356/1‑150 (MQ=255)
                                                                                   ttCGGCGTTGTTGGAGATCGAGTAGTTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGAc                                                                    >  2:4390429/1‑150 (MQ=255)
                                                                                       gCGTTGTTGGAGATCGAGTAGTTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTACGACCAGGGTTTCGAAACCGGCTTTGACCAGc                                                               >  2:618358/1‑150 (MQ=255)
                                                                                         gttgttGGAGATCGAGTAGTTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCTGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTc                                                             <  2:4487396/150‑1 (MQ=255)
                                                                                            gttgGAGATCGAGTAGTTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGAc                                                          <  2:512099/150‑1 (MQ=255)
                                                                                               ggAGATCGAGTAGTTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCa                                                                                                    >  2:1185384/1‑106 (MQ=255)
                                                                                               ggAGATCGAGTAGTTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCa                                                                                                    <  1:1185384/106‑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:5174328/1‑150 (MQ=255)
                                                                                                       agtagtTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCcac                                               >  1:3505792/1‑150 (MQ=255)
                                                                                                        gtagtTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCcaca                                              >  1:112899/1‑150 (MQ=255)
                                                                                                           gtTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAg                                           >  2:4997096/1‑150 (MQ=255)
                                                                                                             tCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGaa                                         <  1:4269499/150‑1 (MQ=255)
                                                                                                               aTGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACa                                       <  2:116588/150‑1 (MQ=255)
                                                                                                                 gTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGc                                     <  1:3011980/150‑1 (MQ=255)
                                                                                                                 gTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGc                                     <  1:3111522/150‑1 (MQ=255)
                                                                                                                  ttGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCca                                                >  2:1821262/1‑138 (MQ=255)
                                                                                                                  ttGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGcc                                    >  1:1071164/1‑150 (MQ=255)
                                                                                                                    ggCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTg                                  <  2:2112033/150‑1 (MQ=255)
                                                                                                                      cGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCt                                <  1:4598804/150‑1 (MQ=255)
                                                                                                                             gccTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGaa                         >  1:176049/1‑150 (MQ=255)
                                                                                                                             gccTTCATAGATCAGGTCTACGATGAGCTTCATTTCGTGCAGGCACTGGAAGTAAGCCATTTCTGGTGCGTAACCAGCTTCGACCAGAACTTCAAAACCTGTGCGAATGAGTTCAACAACGCCACCGCAAAGAACAGCCTGCTCACCGaa                         <  1:4378368/150‑1 (MQ=255)
                                                                                                                                ttCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAg                      <  1:1316464/150‑1 (MQ=255)
                                                                                                                                ttCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAg                      <  1:3378314/150‑1 (MQ=255)
                                                                                                                                 tCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAgg                     >  1:1568816/1‑150 (MQ=255)
                                                                                                                                   gTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGtc                   <  1:1217649/150‑1 (MQ=255)
                                                                                                                                   gTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGtc                   <  1:3085549/150‑1 (MQ=255)
                                                                                                                                    tagatcaGGTCTACGATGAGCTTCATTTCGTGCAGGCACTCGAAGTAAGCCATTTCTGGTGCGTAACCAGCTTCGACCAGAACTTCAAAACCTGTGCGACTGAGTTCAACAACGCCACCGCAAAGAACAGCCTGCTCACCGAACAGGtca                  >  2:2376439/7‑149 (MQ=255)
                                                                                                                                       atcaGGTCTACGATGAGCTTCATTTCGTGCAGGCACTCGAAGTAAGCCATTTCTGGTGCGTAACCAGCTTCGACCAGAACTTCAAAACCTGTGCGAATGAGTTCAACAACGCCACCGCAAAGAACAGCCTGCTCACCGAACAgg                     <  1:3545579/141‑1 (MQ=255)
                                                                                                                                       atcaGGTCTACGATGAGCTTCATTTCGTGCAGGCACTCGAAGTAAGCCATTTCTGGTGCGTAACCAGCTTCGACCAGAACTTCAAAACCTGTGCGAATGAGTTCAACAACGCCACCGCAAAGAACAGCCTGCTCACCGAACAgg                     >  2:3545579/4‑144 (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)
                                                                                                                                            gTCTACGATGAGCTTCATTTCGTGCAGGCACTCGAAGTAAGCCATTTCTGGTGCGTAACCAGCTTCGACCAGAACTTCAAAACCTGTGCGAATGAGTTCAACAACGCCACCGCAAAGAACAGCCTGCTCACCGAACAGGtcagtttcgca          >  2:1449409/1‑148 (MQ=255)
                                                                                                                                            gTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGtcggtttcggt          <  1:91738/150‑1 (MQ=255)
                                                                                                                                            gTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGtcggtttcggt          <  1:1731064/150‑1 (MQ=255)
                                                                                                                                             tCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGtcggtttcggtt         >  2:2969845/1‑150 (MQ=255)
                                                                                                                                               tACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCCGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAAc                                        <  2:5272556/117‑1 (MQ=255)
                                                                                                                                               tACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGACCAGCCTGCTCACCGAACAGGTCggtttcggtttc       <  2:3279630/150‑1 (MQ=255)
                                                                                                                                               tACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAAc                                        >  1:5272556/1‑117 (MQ=255)
                                                                                                                                               tACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCggtttcggtttc       <  1:250249/150‑1 (MQ=255)
                                                                                                                                               gaCGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCggtttcggtttc       <  2:3608723/149‑1 (MQ=255)
                                                                                                                                                  gATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTc    >  1:4247407/1‑150 (MQ=255)
                                                                                                                                                    tcagcttcagctCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCt  >  1:2589816/1‑150 (MQ=255)
                                                                                                                                                    tcagcttcagctCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCCGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCt  >  1:1967893/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TCCTGGATGCGCTTCAGAGCATTGCGCATGGCCTTGCGGGATTCTTCGTTGATGACTTCCGGGCCAGTGACGTACTCACCGTATTCGGCGTTGTTGGAGATCGAGTAGTTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCT  >  NC_002947/5318222‑5318519

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: