Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,318,452 G→A 8.6% V206V (GTC→GTT) ‡ ilvC ← ketol‑acid reductoisomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,318,4520GA 8.6% 545.0 / 30.3 186V206V (GTC→GTT) ‡ilvCketol‑acid reductoisomerase
Reads supporting (aligned to +/- strand):  ref base G (85/85);  new base A (12/4);  total (97/89)
Fisher's exact test for biased strand distribution p-value = 6.86e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01

Only 100 of 186 total aligned reads displayed.

ATTCGGCGTTGTTGGAGATCGAGTAGTTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCGTTGCC  >  NC_002947/5318304‑5318598
                                                                                                                                                    |                                                                                                                                                  
 ttCGGCGTTGTTGGAGATCGAGTAGTTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGAc                                                                                                                                                  >  2:2074943/1‑150 (MQ=255)
 ttCGGCGTTGTTGGAGATCGAGTAGTTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGAc                                                                                                                                                  <  1:2844013/150‑1 (MQ=255)
 ttCGGCGTTGTTGGAGATCGAGTAGTTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGAc                                                                                                                                                  >  2:4390429/1‑150 (MQ=255)
 ttCGGCGTTGTTGGAGATCGAGTAGTTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGAc                                                                                                                                                  >  2:1247356/1‑150 (MQ=255)
     gCGTTGTTGGAGATCGAGTAGTTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTACGACCAGGGTTTCGAAACCGGCTTTGACCAGc                                                                                                                                              >  2:618358/1‑150 (MQ=255)
       gttgttGGAGATCGAGTAGTTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCTGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTc                                                                                                                                            <  2:4487396/150‑1 (MQ=255)
          gttgGAGATCGAGTAGTTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGCTCAGCTTCAGCACGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGAc                                                                                                                                         <  2:1185388/150‑1 (MQ=255)
             ggAGATCGAGTAGTTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGt                                                                                                                                      >  1:2276987/1‑150 (MQ=255)
                gaTCGAGTAGTTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTAcc                                                                                                                                   <  1:982093/150‑1 (MQ=255)
                    gAGTAGTTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCca                                                                                                                               >  2:5173927/1‑150 (MQ=255)
                         gtTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAg                                                                                                                          >  2:4997096/1‑150 (MQ=255)
                           tCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGaa                                                                                                                        <  1:1139063/150‑1 (MQ=255)
                            cATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAAc                                                                                                                       >  2:3027973/1‑150 (MQ=255)
                             aTGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACa                                                                                                                      <  1:1546770/150‑1 (MQ=255)
                             aTGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACa                                                                                                                      <  2:116588/150‑1 (MQ=255)
                               gTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGc                                                                                                                    <  1:1612041/150‑1 (MQ=255)
                               gTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGc                                                                                                                    <  1:3111522/150‑1 (MQ=255)
                               gTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGc                                                                                                                    <  1:3011980/150‑1 (MQ=255)
                                ttGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCca                                                                                                                               <  1:1821262/138‑1 (MQ=255)
                                ttGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCGCACAGAACAGcc                                                                                                                   >  2:4637366/1‑150 (MQ=255)
                                ttGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGcc                                                                                                                   <  1:3487039/150‑1 (MQ=255)
                                ttGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGcc                                                                                                                   >  1:1071164/1‑150 (MQ=255)
                                  ggCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTg                                                                                                                 <  2:2112033/150‑1 (MQ=255)
                                           gccTTCATAGATCAGGTCTACGATGAGCTTCATTTCGTGCAGGCACTGGAAGTAAGCCATTTCTGGTGCGTAACCAGCTTCGACCAGAACTTCAAAACCTGTGCGAATGAGTTCAACAACGCCACCGCAAAGAACAGCCTGCTCACCGaa                                                                                                        <  1:4378368/150‑1 (MQ=255)
                                              ttCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAg                                                                                                     <  1:3378314/150‑1 (MQ=255)
                                              ttCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAg                                                                                                     >  2:533989/1‑150 (MQ=255)
                                                cGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGt                                                                                                   <  1:3823723/150‑1 (MQ=255)
                                                 gTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGtc                                                                                                  <  1:3085549/150‑1 (MQ=255)
                                                  tagatcaGGTCTACGATGAGCTTCATTTCGTGCAGGCACTCGAAGTAAGCCATTTCTGGTGCGTAACCAGCTTCGACCAGAACTTCAAAACCTGTGCGAATGAGTTCAACAACGCCACCGCAAAGAACAGCCTGCTCACCGAACAGGtca                                                                                                 >  2:2376485/7‑149 (MQ=255)
                                                     atcaGGTCTACGATGAGCTTCATTTCGTGCAGGCACTCGAAGTAAGCCATTTCTGGTGCGTAACCAGCTTCGACCAGAACTTCAAAACCTGTGCGAATGAGTTCAACAACGCCACCGCAAAGAACAGCCTGCTCACCGAACAgg                                                                                                    <  1:3545579/141‑1 (MQ=255)
                                                     atcaGGTCTACGATGAGCTTCATTTCGTGCAGGCACCCGAAGTAAGCCATTTCTGGTGCGTAACCAGCTTCGACCAGAACTTCAAAACCTGTGCGAATGAGTTCAACAACGCCACCGCAAAGAACAGCCTGCTCACCGAACAGGtcagtt                                                                                              >  2:3392950/4‑150 (MQ=255)
                                                        aGGTCTACGATGAGCTTCATTTCGTGCAGGCACTCGAAGTAAGCCATTTCTGGTGCGTAACCAGCTTCGACCAGAACTTCAAAACCTGTGCGAATGAGTTCAACAACGCCACCGCAAAGAACAGCCTGCTCACCGAACAGGtcagtttcg                                                                                           >  2:1985383/1‑150 (MQ=255)
                                                        aGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGtcggtttcg                                                                                           <  2:1888686/150‑1 (MQ=255)
                                                        aGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGtcggtttcg                                                                                           <  1:2010160/150‑1 (MQ=255)
                                                        aGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGtcggtttcg                                                                                           <  2:5365406/150‑1 (MQ=255)
                                                        aGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTACGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGtcggtttcg                                                                                           <  2:5090922/150‑1 (MQ=255)
                                                          gTCTACGATGAGCTTCATTTCGTGCAGGCACTCGAAGTAAGCCATTTCTGGTGCGTAACCAGCTTCGACCAGAACTTCAAAACCTGTGCGAATGAGTTCAACAACGCCACCGCAAAGAACAGCCTGCTCACCGAACAGGtcagtttcgca                                                                                         >  2:1449409/1‑148 (MQ=255)
                                                          gTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGtcggtttcggt                                                                                         <  1:1731064/150‑1 (MQ=255)
                                                           tCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGtcggtttcggtt                                                                                        >  2:723549/1‑150 (MQ=255)
                                                           tCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGtcggtttcggtt                                                                                        >  2:2969845/1‑150 (MQ=255)
                                                             tACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGACCAGCCTGCTCACCGAACAGGTCggtttcggtttc                                                                                      <  2:3279630/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)
                                                                     gcttcatttCGTGCAGGCACTCGAAGTAAGCCATTTCTGGTGCGTAACCAGCTTCGACCAGAACTTCAAAACCTGTGCGAATGAGTTCAACAACGCCACCGCAAAGAACAGCCTGCTCACCGAACAGGTCAGTTTCGCACTCTTCCTTGa                                                                              >  2:3812120/1‑150 (MQ=255)
                                                                      cttcatttCGTGCAGGCACTCGAAGTAAGCCATTTCTGGTGCGTAACCAGCTTCGACCAGAACTTCAAAACCTGTGCGAATGAGTTCAACAACGCCACCGCAAAGAACAGCCTGCTCACCGAACAGGTCAGTTTCGCACTCTTCCTTGaa                                                                             <  1:1673985/150‑1 (MQ=255)
                                                                      cttcagctCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGaa                                                                             >  2:2526920/1‑150 (MQ=255)
                                                                       ttcatttCGTGCAGGCACTCGAAGTAAGCCATTTCTGGTGCGTAACCAGCTTCGACCAGAACTTCAAAACCTGTGCGAATGAGTTCAACAACGCCACCGCAAAGAACAGCCTGCTCACCGAACAGGTCAGTTTCGCACTCTTGCTTGAAg                                                                            >  1:3377839/7‑150 (MQ=255)
                                                                         cagctCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGCACAGGTCGGTTTCGGTTTCGTCCTTGAAggt                                                                          <  2:2030001/150‑1 (MQ=255)
                                                                           gctCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAggtgg                                                                        <  2:3091615/150‑1 (MQ=255)
                                                                            ctCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAggtggt                                                                       >  2:5445438/1‑150 (MQ=255)
                                                                             tCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGtt                                                                      >  2:2601683/1‑150 (MQ=255)
                                                                                 gCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCga                                                                  <  1:2767457/150‑1 (MQ=255)
                                                                                  cAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCgat                                                                 <  1:4332755/150‑1 (MQ=255)
                                                                                  cAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCgat                                                                 <  1:5187505/150‑1 (MQ=255)
                                                                                  cAGGCACTCGAAGTAAGCCATTTCTGGTGCGTAACCAGCTTCGACCAGAACTTCAAAACCTGTGCGAATGAGTTCAACAACGCCACCGCAAAGAACAGCCTGCTCACCGAACAGGTCAGTTTCGCACTCTTCCTTGAAGGTGGTTTCaat                                                                 >  2:824567/1‑150 (MQ=255)
                                                                                     gCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCgatgat                                                              >  2:2420655/1‑150 (MQ=255)
                                                                                       aCTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATAc                                                            >  2:520900/1‑150 (MQ=255)
                                                                                         tCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCg                                                          >  2:147376/1‑150 (MQ=255)
                                                                                         tCGAAGTAAGCCATTTCTGGTGCGTAACCAGCTTCGACCAGAACTTCAAAACCTGTGCGAATGAGTTCAACAACGCCACCGCAAAGAACAGCCTGCTCACCGAACAGGTCAGTTTCGCACTCTTCCTTGAAGGTGGTTTCaataacgccc                                                          >  2:2851012/1‑143 (MQ=255)
                                                                                             aGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTAc                                                      <  2:2692975/150‑1 (MQ=255)
                                                                                                  gCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGcca                                                 >  2:3753865/1‑150 (MQ=255)
                                                                                                   ccATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGccac                                                <  2:3091523/150‑1 (MQ=255)
                                                                                                    cATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGccacc                                               >  2:4330766/1‑150 (MQ=255)
                                                                                                    cATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGccacc                                               >  1:1326425/1‑150 (MQ=255)
                                                                                                     aTTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGccacca                                              >  1:4300867/1‑150 (MQ=255)
                                                                                                      tttCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCcaccac                                             >  2:3701208/1‑150 (MQ=255)
                                                                                                              gcgTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGccc                                     <  1:2045342/150‑1 (MQ=255)
                                                                                                                 tAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGaa                                  <  2:2713580/150‑1 (MQ=255)
                                                                                                                 tAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGaa                                  >  2:1584760/1‑150 (MQ=255)
                                                                                                                 tAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGaa                                  <  2:178253/150‑1 (MQ=255)
                                                                                                                    ccAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCa                               >  1:3219769/1‑150 (MQ=255)
                                                                                                                     cAGCTTCGACCCGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCAt                              <  2:4830273/150‑1 (MQ=255)
                                                                                                                     cAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCAt                              >  1:2508607/1‑150 (MQ=255)
                                                                                                                     cAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCAt                              <  1:1206640/150‑1 (MQ=255)
                                                                                                                       gCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAggtg                                                                         >  1:3369438/1‑105 (MQ=255)
                                                                                                                       gCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAg                            <  1:2500454/150‑1 (MQ=255)
                                                                                                                       gCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAg                            <  1:1637259/150‑1 (MQ=255)
                                                                                                                       gCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAg                            <  1:2757706/150‑1 (MQ=255)
                                                                                                                        cTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAgg                           >  1:4451517/1‑150 (MQ=255)
                                                                                                                         ttCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGa                          <  1:1835761/150‑1 (MQ=255)
                                                                                                                         ttCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCAGAGGa                          >  2:3430821/1‑150 (MQ=255)
                                                                                                                            gACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCAAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAg                       <  2:2276987/150‑1 (MQ=255)
                                                                                                                             aCCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGt                      <  1:1247356/150‑1 (MQ=255)
                                                                                                                              ccAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTg                                                                                                                 <  2:857086/58‑1 (MQ=255)
                                                                                                                              ccAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTg                                                                                                                 >  1:857086/1‑58 (MQ=255)
                                                                                                                               cAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGc                    >  2:4995417/1‑150 (MQ=255)
                                                                                                                                 gggTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGa                  >  1:2423158/1‑150 (MQ=255)
                                                                                                                                    tttCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGt               >  1:1053168/1‑150 (MQ=255)
                                                                                                                                        gAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCtt           >  1:4781587/1‑150 (MQ=255)
                                                                                                                                        gAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCtt           <  2:4945226/150‑1 (MQ=255)
                                                                                                                                         aaaCCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTg          <  1:976972/150‑1 (MQ=255)
                                                                                                                                         aaaCCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTg          >  2:5155450/1‑150 (MQ=255)
                                                                                                                                            ccGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCg       >  1:5209976/1‑150 (MQ=255)
                                                                                                                                             cGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCGt      <  2:929090/150‑1 (MQ=255)
                                                                                                                                              ggCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCGtt     <  2:3505792/150‑1 (MQ=255)
                                                                                                                                               gCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCGTTg    >  1:5040472/1‑150 (MQ=255)
                                                                                                                                                cTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCGTTGc   >  1:4529643/1‑150 (MQ=255)
                                                                                                                                                 tttGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCGTTGcc  <  2:2579641/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
ATTCGGCGTTGTTGGAGATCGAGTAGTTCATGTTGGCGATGCCGCCTTCGTACATGAGGTCTACGATCAGCTTCAGCTCGTGCAGGCACTCGAAGTAGGCCATTTCCGGCGCGTAGCCAGCTTCGACCAGGGTTTCGAAACCGGCTTTGACCAGCTCGACGGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCGTTGCC  >  NC_002947/5318304‑5318598

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: