Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,318,512 T→C 10.7% E186E (GAA→GAG ilvC ← ketol‑acid reductoisomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,318,5120TC10.7% 933.6 / 73.8 326E186E (GAA→GAGilvCketol‑acid reductoisomerase
Reads supporting (aligned to +/- strand):  ref base T (157/134);  new base C (21/14);  total (178/148)
Fisher's exact test for biased strand distribution p-value = 5.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.19e-01

Only 100 of 358 total aligned reads displayed.

TCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCGTTGCCGGAAGCGTCCTGGTAGATAGCGATCAGGTCAGGGATGCCGCCGCCTTTGACGAACTCGGAGC  >  NC_002947/5318363‑5318660
                                                                                                                                                     |                                                                                                                                                     
tCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGtcggtttcggtt                                                                                                                                                       >  1:2194611/1‑150 (MQ=255)
tCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGtcggtttcggtt                                                                                                                                                       <  1:3004278/150‑1 (MQ=255)
  tACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCggtttcggtttc                                                                                                                                                     <  2:6002901/150‑1 (MQ=255)
  tACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCggtttcggtttc                                                                                                                                                     <  2:8379682/150‑1 (MQ=255)
  tACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCggtttcggtttc                                                                                                                                                     <  1:4498234/150‑1 (MQ=255)
   aCGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGgtttcggtttcg                                                                                                                                                    <  1:7502745/150‑1 (MQ=255)
   aCGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGgtttcggtttcg                                                                                                                                                    <  2:1066363/150‑1 (MQ=255)
      aTCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTcc                                                                                                                                                 <  1:539795/150‑1 (MQ=255)
           cttcagctCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGaa                                                                                                                                            >  1:6824613/1‑150 (MQ=255)
            ttcatttCGTGCAGGCACTCGAAGTAAGCCATTTCTGGTGCGTAACCAGCTTCGACCAGAACTTCAAAACCTGTGCGAATGAGTTCAACAACGCCACCGCAAAGAACAGCCTGCTCACCGAACAGGTCAGTTTCGCACTCTTCCTTGAAg                                                                                                                                           <  2:4375415/144‑1 (MQ=255)
            ttcagctCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAg                                                                                                                                           >  1:556913/1‑150 (MQ=255)
            ttcagctCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAg                                                                                                                                           >  1:2458588/1‑150 (MQ=255)
            ttcagctCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAg                                                                                                                                           >  1:4831854/1‑150 (MQ=255)
            ttcagctCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAg                                                                                                                                           <  2:271810/150‑1 (MQ=255)
               agctCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAggtg                                                                                                                                        <  2:6594333/150‑1 (MQ=255)
                gctCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAggtgg                                                                                                                                       >  1:4674177/1‑150 (MQ=255)
                  tCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGtt                                                                                                                                     >  1:2350029/1‑150 (MQ=255)
                  tCGTGCAGGCACTCGAAGTAAGCCATTTCTGGTGCGTAACCAGCTTCGACCAGAACTTCAAAACCTGTGCGAATGAGTTCAACAACGCCACCGCAAAGAACAGCCTGCTCACCGAACAGGTCAGTTTCGCACTCTTCCTTGAAGGTGGtt                                                                                                                                     <  2:8095905/150‑1 (MQ=255)
                   cGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGttt                                                                                                                                    <  2:6281520/150‑1 (MQ=255)
                    gTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTc                                                                                                                                   >  2:341935/1‑150 (MQ=255)
                    gTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTc                                                                                                                                   >  2:622947/1‑150 (MQ=255)
                    gTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTc                                                                                                                                   >  2:2047394/1‑150 (MQ=255)
                      gCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCga                                                                                                                                 >  1:8562860/1‑150 (MQ=255)
                      gCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCga                                                                                                                                 <  2:1300963/150‑1 (MQ=255)
                       cAGGCACTCGAAGTAGGCCATTTCCGGAGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCgat                                                                                                                                <  1:8722061/150‑1 (MQ=255)
                         ggCACTCGAAGTAAGCCATTTCTGGTGCGTAACCAGCTTCGACCAGAACTTCAAAACCTGTGCGAATGAGTTCAACAACGCCACCGCAAAGAACAGCCTGCTCACCGAACAGGTCAGTTTCGCACTCTTCCTTGAAGGTGGTTTCaataa                                                                                                                              >  1:334862/1‑148 (MQ=255)
                              tCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCg                                                                                                                         <  2:7032317/150‑1 (MQ=255)
                              tCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCg                                                                                                                         >  1:6529701/1‑150 (MQ=255)
                              tCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCg                                                                                                                         >  2:4226363/1‑150 (MQ=255)
                                   gTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACg                                                                                                                    <  1:3878425/150‑1 (MQ=255)
                                    tAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACgg                                                                                                                   >  2:8421947/1‑150 (MQ=255)
                                    tAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACgg                                                                                                                   >  2:6414052/1‑150 (MQ=255)
                                       gCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGcca                                                                                                                >  2:4431872/1‑150 (MQ=255)
                                        ccATTTCCGGAGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGATTCGTCCTGGAAGGTGGTTTCGATGATACCGGAACGGccac                                                                                                               >  2:3985864/1‑150 (MQ=255)
                                         cATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGccacc                                                                                                              >  1:3820278/1‑150 (MQ=255)
                                          aTTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGccacca                                                                                                             >  1:6195373/1‑150 (MQ=255)
                                             tCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCt                                                                                                          <  1:848207/150‑1 (MQ=255)
                                               cGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTAc                                                                                                        >  2:300555/1‑150 (MQ=255)
                                               cGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTAc                                                                                                        <  1:1223367/150‑1 (MQ=255)
                                                 gcgcGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGc                                                                                                      >  1:2617965/1‑150 (MQ=255)
                                                  cgcgTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGcc                                                                                                     <  2:4261286/150‑1 (MQ=255)
                                                   gcgTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGccc                                                                                                    >  1:5265679/1‑150 (MQ=255)
                                                              ttCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGa                                                                                         >  1:3669048/1‑150 (MQ=255)
                                                                cGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACa                                                                                       <  1:5041111/150‑1 (MQ=255)
                                                                  aCCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGt                                                                                     <  1:1812741/150‑1 (MQ=255)
                                                                   ccAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTg                                                                                    >  1:6668333/1‑150 (MQ=255)
                                                                      gggTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTgtcc                                                                                 >  2:5519785/1‑147 (MQ=255)
                                                                      gggTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGa                                                                                 >  2:7511609/1‑150 (MQ=255)
                                                                            cGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCt                                                                           <  2:489244/150‑1 (MQ=255)
                                                                             gAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCtt                                                                          <  2:8248979/150‑1 (MQ=255)
                                                                             gAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAACCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCtt                                                                          >  2:3760845/1‑150 (MQ=255)
                                                                              aaaCCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTg                                                                         <  1:5635820/150‑1 (MQ=255)
                                                                              aaaCCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTg                                                                         >  2:5831919/1‑150 (MQ=255)
                                                                                  cGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCGt                                                                    >  2:5245604/1‑150 (MQ=255)
                                                                                     cTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCGTTGc                                                                 <  2:4135010/150‑1 (MQ=255)
                                                                                     cTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCGTTGc                                                                 <  1:2019232/150‑1 (MQ=255)
                                                                                     cTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCGTTGc                                                                 <  1:2059179/150‑1 (MQ=255)
                                                                                      tttGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCGTTGcc                                                                <  2:1509390/150‑1 (MQ=255)
                                                                                          aCCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCGTTGCCGGaa                                                            <  1:1853530/150‑1 (MQ=255)
                                                                                          aCCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCGTTGCCGGaa                                                            <  1:1391448/150‑1 (MQ=255)
                                                                                          aCCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCGTTGCCGGaa                                                            <  2:8672939/150‑1 (MQ=255)
                                                                                           ccAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTgg                                                                       <  2:2608840/138‑1 (MQ=255)
                                                                                            cAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCGTTGCCGGAAGc                                                          >  1:1495100/1‑150 (MQ=255)
                                                                                                tCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCGTTGCCGGAAGCGTcc                                                      <  2:5298240/150‑1 (MQ=255)
                                                                                                tCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCGTTGCCGGAAGCGTcc                                                      >  2:862617/1‑150 (MQ=255)
                                                                                                tCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCGTTGCCGGAAGCGTcc                                                      >  1:4120151/1‑150 (MQ=255)
                                                                                                 cGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCGTTGCCGGAAGCGTCCt                                                     <  1:9042970/150‑1 (MQ=255)
                                                                                                 cGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCGTTGCCGGAAGCGTCCt                                                     <  1:803887/150‑1 (MQ=255)
                                                                                                  gACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCGTTGCCGGAAGCGTCCTg                                                    >  2:3576280/1‑150 (MQ=255)
                                                                                                  gACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCGTTGCCGGAAGCGTCCTg                                                    >  1:3977737/1‑150 (MQ=255)
                                                                                                       tACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCGTTGCCGGAAGCGTCCTGGTAGa                                               >  2:1407395/1‑150 (MQ=255)
                                                                                                         ccGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCGTTGCCGGAAGCGTCCTGGTAGATa                                             <  2:8562860/150‑1 (MQ=255)
                                                                                                          cGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCGTTGCCGGAAGCGTCCTGGTAGATAg                                            >  2:7500002/1‑150 (MQ=255)
                                                                                                           gCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCGTTGCCGGAAGCGTCCTGGTAGATAGc                                           >  2:3437067/1‑150 (MQ=255)
                                                                                                           accgcaaaGAACAGCCTGCTCACCGAACAGGTCAGTTTCGCACTCTTCCTTGAAGGTGGTTTCAATAACGCCCGAACGACCGCCACCAACGCCCGAAGCGTAAGACAGAGCGAGATCATGTGCGTTGCCCGAAGCATCCTGATGGATCGc                                           >  2:761613/8‑150 (MQ=255)
                                                                                                           accgcaaaGAACAGCCTGCTCACCGAACAGGTCAGTTTCGCACTCTTCCTTGAAGGTGGTTTCAATAACGCCCGAACGACCGCCACCAACGCCCGAAGCGTAAGACAGAGCGAGATCATGTGCGTTGCCCGAAGCATCCTGATGGATCGc                                           <  1:8959033/143‑1 (MQ=255)
                                                                                                                acaGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCGTTGCCGGAAGCGTCCTGGTAGATAGCGAtca                                      <  2:5426427/150‑1 (MQ=255)
                                                                                                                aaaGAACAGCCTGCTCACCGAACAGGTCAGTTTCGCACTCTTCCTTGAAGGTGGTTTCAATAACGCCCGAACGACCGCCACCAACGCCCGAAGCGTAAGACAGAGTGAGATCATGTGCGTTGCCCGAAGCATCCTGATGGATCGCGAtaa                                      >  2:2908265/3‑148 (MQ=255)
                                                                                                                   gAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCGTTGCCGGAAGCGTCCTGGTAGATAGCGAtcaggt                                   <  1:22659/150‑1 (MQ=255)
                                                                                                                   gAACAGCCTGCTCACCGAACAGGTCAGTTTCGCACTCTTCCTTGAAGGTGGTTTCAATAACGCCCGAACGACCGCCACCAACGCCCGAAGCGTAAGACAGAGCGAGATCATGTGCGTTGCCCGAAGCATCCTGATGGATCGCGAtaaggc                                   <  2:1930986/150‑2 (MQ=255)
                                                                                                                   gAACAACCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCGTTGCCGGAAGCGTCCTGGTAGATAGCGAtcaggt                                   <  1:3760845/150‑1 (MQ=255)
                                                                                                                     aCAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCGTTGCCGGAAGCGTCCTGGTAGATAGCGAtcaggtca                                 <  1:7086217/150‑1 (MQ=255)
                                                                                                                       aGCCTGCTCACCGAACAGGTCAGTTTCGCACTCTTCCTTGAAGGTGGTTTCAATAACGCCCGAACGACCGCCACCAACGCCCGAAGCGTAAGACAGAGCGAGATCATGTGCGTTGCCCGAAGCATCCTGATGGATCGCGAtaaggcaagg                               >  2:2235839/1‑145 (MQ=255)
                                                                                                                       aGCCTGCTCACCGAACAGGTCAGTTTCGCACTCTTCCTTGAAGGTGGTTTCAATAACGCCCGAACGACCGCCACCAACGCCCGAAGCGTAAGACAGAGCGAGATCATGTGCGTTGCCCGAAGCATCCTGATGGATCGCGAtaaggcaagg                               >  1:757323/1‑145 (MQ=255)
                                                                                                                        gCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCGTTGCCGGAAGCGTCCTGGTAGATAGCGATCAGGTCAggg                              >  1:106546/1‑150 (MQ=255)
                                                                                                                         ccTGCTCACCGAACAGGTCAGTTTCGCACTCTTCCTTGAAGGTGGTTTCAATAACGCCCGAACGACCGCCACCAACGCCCGAAGCGTAAGACAGAGCGAGATCATGTGCGTTGCCCGAAGCATCCTGATGGATCGCGAtaaggcaaggaa                             >  2:2337775/1‑143 (MQ=255)
                                                                                                                          cTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCGTTGCCGGAAGCGTCCTGGTAGATAGCGATCAGGTCAGGGAt                            <  2:8951236/150‑1 (MQ=255)
                                                                                                                             cTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCGTTGCCGGAAGCGTCCTGGTAGATAGCGATCAGGTCAGGGATgcc                         <  2:7452572/150‑1 (MQ=255)
                                                                                                                               cACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCGTTGCCGGAAGCGTCCTGGTAGATAGCGATCAGGTCAGGGATgccgc                       >  1:7259570/1‑150 (MQ=255)
                                                                                                                               cACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCGTTGCCGGAAGCGTCCTGGTAGATAGCGATCAGGTCAGGGATgccgc                       >  1:6410544/1‑150 (MQ=255)
                                                                                                                                aCCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCGTTGCCGGAAGCGTCCTGGTAGATAGCGATCAGGTCAGGGATgccgcc                      >  2:4608307/1‑150 (MQ=255)
                                                                                                                                aCCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCGTTGCCGGAAGCGTCCTGGTAGATAGCGATCAGGACAGGGATgccgcc                      >  2:1193534/1‑150 (MQ=255)
                                                                                                                                 ccGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCGTTGCCGGAAGCGTCCTGGTAGATAGCGATCAGGTCAGGGATgccgccg                     >  2:1789614/1‑150 (MQ=255)
                                                                                                                                    aaCAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCGTTGCCGGAAGCGTCCTGGTAGATAGCGATCAGGTCAGGGATGCCGCCGCCt                  <  1:831428/150‑1 (MQ=255)
                                                                                                                                         gTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCGTTGCCGGAAGCGTCCTGGTAGATAGCGATCAGGTCAGGGATGCCGCCGCCTTTGAc             >  2:1680350/1‑150 (MQ=255)
                                                                                                                                         gTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCGTTGCCGGAAGCGTCCTGGTAGATAGCGATCAGGTCAGGGATGCCGCCGCCTTTGAc             <  2:2458588/150‑1 (MQ=255)
                                                                                                                                                tcggtttcgTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCGTTGCCGGAAGCGTCCTGGTAGATAGCGATCAGGTCAGGGATGCCGCCGCCTTTGACGAACTCg      <  2:1574246/150‑1 (MQ=255)
                                                                                                                                                  ggtttcgTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCGTTGCCGGAAGCGTCCTGGTAGATAGCGATCAGGTCAGGGATGCCGCCGCCTTTGACGAACTCGGa    <  2:7248383/150‑1 (MQ=255)
                                                                                                                                                  ggtttcgTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCGTTGCCGGAAGCGTCCTGGTAGATAGCGATCAGGTCAGGGATGCCGCCGCCTTTGACGAACTCGGa    <  1:4370708/150‑1 (MQ=255)
                                                                                                                                                    tttcgTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCGTTGCCGGAAGCGTCCTGGTAGATAGCGATCAGGTCAGGGATGCCGCCGCCTTTGACGAACTCGGAGc  >  1:5956399/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCGTTGCCGGAAGCGTCCTGGTAGATAGCGATCAGGTCAGGGATGCCGCCGCCTTTGACGAACTCGGAGC  >  NC_002947/5318363‑5318660

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: