Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,318,413 C→T 10.8% A219A (GCG→GCA ilvC ← ketol‑acid reductoisomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,318,4130CT10.8% 464.1 / 36.0 167A219A (GCG→GCAilvCketol‑acid reductoisomerase
Reads supporting (aligned to +/- strand):  ref base C (77/72);  new base T (9/9);  total (86/81)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01

Only 100 of 168 total aligned reads displayed.

TTCGTTGATGACTTCCGGGCCAGTGACGTACTCACCGTATTCGGCGTTGTTGGAGATCGAGTAGTTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACG  >  NC_002947/5318266‑5318560
                                                                                                                                                   |                                                                                                                                                   
      gatgaCTTCCGGGCCAGTGACGTACTCACCGTATTCGGCGTTGTTGGAGATCGAGTAGTTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGcc                                                                                                                                             >  1:3367093/1‑150 (MQ=255)
       atgaCTTCCGGGCCAGTGACGTACTCACCGTATTCGGCGTTGTTGGAGATCGAGTAGTTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCa                                                                                                                                            <  1:3680177/150‑1 (MQ=255)
        tgaCTTCCGGGCCAGTGACGTACTCACCGTATTCGGCGTTGTTGGAGATCGAGTAGTTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAg                                                                                                                                           >  2:2613583/1‑150 (MQ=255)
        tctgcttcGGGCCAGTGACGTACTCACCGTATTCGGCGTTGTTGGACATCGAGTAGTTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAg                                                                                                                                           >  2:1578309/8‑150 (MQ=255)
         gatgcgcGGACCAGTGACGTATTCGCCCCACTCAGCAGTGTTGGAGATCGAGTAGTTCATGTTCGCAATGCCGCCTTCATAGATCAGGTCTACGATGAGCTTCATTTCGTGCAGGCACTCGAAGTAAGCCATTTCTGGTGCGTAACCAGc                                                                                                                                          >  2:3724958/7‑150 (MQ=255)
         gaCTTCCGGGCCAGTGACGTACTCACCGTATTCGGCGTTGTTGGAGATCGAGTAGTTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGc                                                                                                                                          >  2:1188839/1‑150 (MQ=255)
            ttCCGGGCCAGTGACGTACTCACCGTATTCGGCGTTGTTGGAGATCGAGTAGTTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTc                                                                                                                                       >  1:3349717/1‑150 (MQ=255)
              gcGGACCAGTGACGTATTCGCCCCACTCAGCAGTGTTGGAGATCGAGTAGTTCATGTTCGCAATGCCGCCTTCATAGATCAGGTCTACGATGAGCTTCATTTCGTGCAGGCACTCGAAGTAAGCCATTTCTGGTGCGTAACCAGCTTCGa                                                                                                                                     >  2:2596658/2‑150 (MQ=255)
                 gaccAGTGACGTATTCGCCCCACTCAGCAGTGTTGGAGATCGAGTAGTTCATGTTCGCAATGCCGCCTTCATAGATCAGGTCTACGATGAGCTTCATTTCGTGCAGGCACTCGAAGTAAGCCATTTCTGGTGCGTAACCAGCTTCAACCa                                                                                                                                  <  2:1902722/148‑1 (MQ=255)
                  gCCAGTGACGTACTCACCGTATTCGGCGTTGTTGGAGATCGAGTAGTTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAg                                                                                                                                 <  1:1829714/150‑1 (MQ=255)
                   ccAGTGACGTACTCACCGTATTCGGCGTTGTTGGAGATCGAGTAGTTCTTGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAgg                                                                                                                                <  2:3367059/150‑1 (MQ=255)
                   ccAGTGACGTACTCACCGTATTCGGCGTTGTTGGAGATCGAGTAGTTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAgg                                                                                                                                <  2:3367093/150‑1 (MQ=255)
                        gACGTACTCACCGTATTCGGCGTTGTTGGAGATCGAGTAGTTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTc                                                                                                                           <  2:35737/150‑1 (MQ=255)
                            tACTCACCGTATTCGGCGTTGTTGGAGATCGAGTAGTTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGaaa                                                                                                                       <  1:984666/150‑1 (MQ=255)
                              cTCACCGTATTCGGCGTTGTTGGAGATCGAGTAGTTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAAcc                                                                                                                     >  2:1672160/1‑150 (MQ=255)
                               tCACCGTATTCGGCGTTGTTGGAGATCGAGTAGTTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCg                                                                                                                    >  2:4239190/1‑150 (MQ=255)
                                cACCGTATTCGGCGTTGTTGGAGATCGAGTAGTTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCgg                                                                                                                   >  2:3521265/1‑150 (MQ=255)
                                 aCCGTATTCGGCGTTGTTGGAGATCGAGTAGTTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGc                                                                                                                  >  1:1040194/1‑150 (MQ=255)
                                  ccGTATTCGGCGTTGTTGGAGATCGAGTAGTTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCt                                                                                                                 <  2:2865559/150‑1 (MQ=255)
                                  ccGTATTCGGCGTTGTTGGAGATCGAGTAGTTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCt                                                                                                                 <  2:623788/150‑1 (MQ=255)
                                   cGTATTCGGCGTTGTTGGAGATCGAGTAGTTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCtt                                                                                                                <  2:2333246/150‑1 (MQ=255)
                                       ttCGGCGTTGTTGGAGATCGAGTAGTTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGaaa                                                                                                                       >  1:2875184/1‑139 (MQ=255)
                                          ggCGTTGTTGGAGATCGAGTAGTTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAg                                                                                                         >  2:1503691/1‑150 (MQ=255)
                                           gCGTTGTTGGAGATCGAGTAGTTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGc                                                                                                        <  2:4173233/150‑1 (MQ=255)
                                             gttgttGGAGATCGAGTAGTTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTc                                                                                                      >  1:1128528/1‑150 (MQ=255)
                                               tgttgGAGATCGAGTAGTTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGa                                                                                                    <  2:787068/150‑1 (MQ=255)
                                                gttgGAGATCGAGTAGTTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCACGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGAc                                                                                                   <  2:279866/150‑1 (MQ=255)
                                                 ttgGAGATCGAGTAGTTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACg                                                                                                  <  1:1197878/150‑1 (MQ=255)
                                                 ttgGAGATCGAGTAGTTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACg                                                                                                  >  2:3648815/1‑150 (MQ=255)
                                                  tgGAGATCGAGTAGTTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACgg                                                                                                 >  2:5082179/1‑150 (MQ=255)
                                                    gagaTCGAGTAGTTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTa                                                                                               >  1:4742690/1‑150 (MQ=255)
                                                            gtagtTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCcaca                                                                                       >  2:4179322/1‑150 (MQ=255)
                                                              agtTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCcacaca                                                                                     <  1:3366371/150‑1 (MQ=255)
                                                                   aTGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACa                                                                                >  1:2078989/1‑150 (MQ=255)
                                                                   aTGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACa                                                                                <  2:3895742/150‑1 (MQ=255)
                                                                   aTGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTATGACCAGCTCGACGGTACCGCCACACAGAACa                                                                                >  2:4054298/1‑150 (MQ=255)
                                                                     gTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGc                                                                              <  1:4706471/150‑1 (MQ=255)
                                                                      ttGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCGATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGcc                                                                             >  1:4487612/1‑150 (MQ=255)
                                                                      ttGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGcc                                                                             >  1:1914189/1‑150 (MQ=255)
                                                                        ggCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTg                                                                           <  2:309895/150‑1 (MQ=255)
                                                                        ggCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTg                                                                           >  2:4892864/1‑150 (MQ=255)
                                                                         gCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGc                                                                          >  1:5023367/1‑150 (MQ=255)
                                                                              gccgccTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCAcc                                                                     <  2:3456570/150‑1 (MQ=255)
                                                                              gccgccTTCGTACATGAGGACTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCAcc                                                                     >  1:1442211/1‑150 (MQ=255)
                                                                                cgccTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGa                                                                   >  2:3513871/1‑150 (MQ=255)
                                                                                cgccTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGa                                                                   >  1:3866563/1‑150 (MQ=255)
                                                                                 gccTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGaa                                                                  <  2:751011/150‑1 (MQ=255)
                                                                                 gccTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGaa                                                                  >  1:3006970/1‑150 (MQ=255)
                                                                                  ccTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAAc                                                                 <  2:2078989/150‑1 (MQ=255)
                                                                                    ttCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAg                                                               >  1:135631/1‑150 (MQ=255)
                                                                                       gTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGtc                                                            <  2:3929772/150‑1 (MQ=255)
                                                                                         aCATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGtcgg                                                          >  1:3452868/1‑150 (MQ=255)
                                                                                          cATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGaa                                                                  >  1:4816081/1‑141 (MQ=255)
                                                                                          cATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGaa                                                                  <  2:4816081/141‑1 (MQ=255)
                                                                                          cATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGtcggt                                                         <  2:4596881/150‑1 (MQ=255)
                                                                                           atcaGGTCTACGATGAGCTTCATTTCGTGCAGGCACTCGAAGTAAGCCATTTCTGGTGCGTAACCAGCTTCGACCAGAACTTCAAAACCTGTGCGAATGAGTTCAACAACGCCACCGCAAAGAACAGCCTGCTCACCGAACAGGtcagtt                                                        <  2:3191563/147‑1 (MQ=255)
                                                                                            tcaGGTCTACGATGAGCTTCATTTCGTGCAGGCACTCGAAGTAAGCCATTTCTGGTGCGTAACCAGCTTCGACCAGAACTTCAAAACCTGTGCGAATGAGTTCAACAACGCCACCGCAAAGAACAGCCTGCTCACCGAACAGGtcagttt                                                       >  2:3157125/3‑150 (MQ=255)
                                                                                              aGGTCTACGATGAGCTTCATTTCGTGCAGGCACTCGAAGTAAGCCATTTCTGGTGCGTAACCAGCTTCGACCAGAACTTCAAAACCTGTGCGAATGAGTTCAACAACGCCACCGCAAAGAACAGCCTGCTCACCGAACAGGtcagtttcg                                                     <  1:4270451/150‑1 (MQ=255)
                                                                                               ggTCTACGATCGGCTTCAGCTCGGGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACCGAACAGCCTGCTCACCGAACAGGtcggtttcgg                                                    <  1:3360417/150‑1 (MQ=255)
                                                                                                gTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGtcggtttcggt                                                   <  2:4298343/150‑1 (MQ=255)
                                                                                                  cTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCATAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTcggtttcggttt                                                 >  1:3188964/1‑150 (MQ=255)
                                                                                                   tACGATGAGCTTCATTTCGTGCAGGCACTCGAAGTAAGCCATTTCTGGTGCGTAACCAGCTTCGACCAGAACTTCAAAACCTGTGCGAATGAGTTCAACAACGCCACCGCAAAGAACAGCCTGCTCACCGAACAGGtcagtttcgcactc                                                >  1:1583464/1‑145 (MQ=255)
                                                                                                   tACGATGAGCTTCATTTCGTGCAGGCACTCGAAGTAAGCCATTTCTGGTGCGTAACCAGCTTCGACCAGAACTTCAAAACCTGTGCGAATGAGTTCAACAACGCCACCGCAAAGAACAGCCTGCTCACCGAACAGGtcagtgtcgcactc                                                >  2:1568382/1‑145 (MQ=255)
                                                                                                   tACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCggtttcggtttc                                                >  2:2652573/1‑150 (MQ=255)
                                                                                                   tACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCggtttcggtttc                                                <  2:3048505/150‑1 (MQ=255)
                                                                                                     cGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGt                                              <  2:319590/150‑1 (MQ=255)
                                                                                                        tcagcttcagctCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCt                                           >  2:653967/1‑150 (MQ=255)
                                                                                                           gcttcatttCGTGCAGGCACTCGAAGTAAGCCATTTCTGGTGCGTAACCAGCTTCGACCAGAACTTCAAAACCTGTGCGAATGAGTTCAACAACGCCACCGCAAAGAACAGCCTGCTCACCGAACAGGTCAGTTTCGCACTCTTCCTTGa                                        >  1:4475178/1‑150 (MQ=255)
                                                                                                           gcttcagctCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGa                                        <  2:1128528/150‑1 (MQ=255)
                                                                                                            cttcagctCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTcc                                            >  2:1783472/1‑145 (MQ=255)
                                                                                                            cttcagctCGTGCAGGCACTCGAAGTAGGCCATGTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGaa                                       <  1:1020700/150‑1 (MQ=255)
                                                                                                               cagctCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAggt                                    <  1:1421472/150‑1 (MQ=255)
                                                                                                                aggtCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAggtg                                   >  1:2176490/4‑150 (MQ=255)
                                                                                                                agctCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTGGAAggtg                                   >  1:2176525/1‑150 (MQ=255)
                                                                                                                 gctCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACa                                                                <  1:4271682/120‑1 (MQ=255)
                                                                                                                 gctCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAggtgg                                  <  1:653967/150‑1 (MQ=255)
                                                                                                                 gctCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAggtgg                                  >  2:3591611/1‑150 (MQ=255)
                                                                                                                   tCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGtt                                >  1:1946566/1‑150 (MQ=255)
                                                                                                                       gCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCga                            <  2:2822347/150‑1 (MQ=255)
                                                                                                                       gCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCga                            <  2:4453179/150‑1 (MQ=255)
                                                                                                                       gCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCga                            >  2:2576131/1‑150 (MQ=255)
                                                                                                                       gCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCga                            <  1:3648563/150‑1 (MQ=255)
                                                                                                                        cAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCgat                           <  1:4597963/150‑1 (MQ=255)
                                                                                                                        cAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCgat                           <  1:4597938/150‑1 (MQ=255)
                                                                                                                        cAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCgat                           <  1:3521265/150‑1 (MQ=255)
                                                                                                                         aGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGtcg                                                           >  1:617685/1‑117 (MQ=255)
                                                                                                                            cACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATa                       <  1:1182891/150‑1 (MQ=255)
                                                                                                                            cACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATa                       <  1:4534619/150‑1 (MQ=255)
                                                                                                                            cACTCGAAGTAAGCCATTTCTGGTGCGTAACCAGCTTCGACCAGAACTTCAAAACCTGTGCGAATGAGTTCAACAACGCCACCGCAAAGAACAGCCTGCTCACCGAACAGGTCAGTTTCGCACTCTTCCTTGAAGGTGGTTTCaataacg                       <  1:3788304/150‑5 (MQ=255)
                                                                                                                              cTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATAcc                     >  2:2899673/1‑150 (MQ=255)
                                                                                                                               tCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCg                    >  2:3186035/1‑150 (MQ=255)
                                                                                                                               tCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCg                    >  2:1997077/1‑150 (MQ=255)
                                                                                                                                 gAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGt                  >  1:2947086/1‑150 (MQ=255)
                                                                                                                                         ccATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGccac          >  1:958368/1‑150 (MQ=255)
                                                                                                                                         ccATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGccac          <  1:4374386/150‑1 (MQ=255)
                                                                                                                                         ccATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGccac          <  1:140227/150‑1 (MQ=255)
                                                                                                                                           aTTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGccacca        <  1:2613583/150‑1 (MQ=255)
                                                                                                                                           aTTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGccacca        <  1:2380967/150‑1 (MQ=255)
                                                                                                                                           aTTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGccacca        >  1:1601019/1‑150 (MQ=255)
                                                                                                                                                cGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTAc   >  1:2310440/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TTCGTTGATGACTTCCGGGCCAGTGACGTACTCACCGTATTCGGCGTTGTTGGAGATCGAGTAGTTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACG  >  NC_002947/5318266‑5318560

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: