Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,670,019:1 +T 12.9% intergenic (+72/‑95) prpF → / → prpD aconitate isomerase/2‑methylcitrate dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,670,0191.T12.9% 534.9 / 79.5 139intergenic (+72/‑95)prpF/prpDaconitate isomerase/2‑methylcitrate dehydratase
Reads supporting (aligned to +/- strand):  ref base . (41/80);  new base T (7/11);  total (48/91)
Fisher's exact test for biased strand distribution p-value = 7.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.00e-01

Only 100 of 139 total aligned reads displayed.

CCAAGGCAATCATGAGCCGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACA  >  NC_002947/2669874‑2670166
                                                                                                                                                  |                                                                                                                                                    
ccAAGGCAATCATGAGCCGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATca                                                                                                                                                  >  1:2122133/1‑150 (MQ=255)
 cAAGGCAATCATGAGCCGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCGTGCCCGCTCCCCCAGGTAATcac                                                                                                                                                 >  1:3558987/1‑150 (MQ=255)
  aaGGCAATCATGAGCCGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATcaca                                                                                                                                                <  2:2235717/150‑1 (MQ=255)
  aaGGCAATCATGAGCCGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATcaca                                                                                                                                                <  2:6857755/150‑1 (MQ=255)
  aaGGCAATCATGAGCCGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATcaca                                                                                                                                                >  1:107173/1‑150 (MQ=255)
    ggCAATCATGAGCCGTATTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTTAAGCACCACTGGCGCTCGGCTTCAATAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGATCCCACAGGTAATCACAgg                                                                                                                                              <  2:2412211/150‑1 (MQ=255)
    ggCAATCATGAGCCGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGGATATCACAg                                                                                                                                               >  1:6494570/1‑150 (MQ=255)
        aTCATGAGCCGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGttt                                                                                                                                          <  1:7988753/150‑1 (MQ=255)
        aTCATGAGCCGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGttt                                                                                                                                          <  2:4208774/150‑1 (MQ=255)
           aTGAGCCGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCTCGTGCCTGGCGATAGCTTCTAAAGCACCACGGGCGCTCGGCTTCAAGAGGGGCCTCTTTGCGGCCCTTCTCGGGCTTGACCGCTCCCTCAGGTTATCACAGGGTTTGaa                                                                                                                                      <  2:4560140/150‑1 (MQ=255)
           aTGAGCCGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGaa                                                                                                                                      <  1:5139203/150‑1 (MQ=255)
               gCCGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCt                                                                                                                                  <  2:7854433/150‑1 (MQ=255)
               gCCGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCt                                                                                                                                  <  1:7774743/150‑1 (MQ=255)
                ccGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCtg                                                                                                                                 <  1:6748941/150‑1 (MQ=255)
                ccGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCtg                                                                                                                                 <  1:8812577/150‑1 (MQ=255)
                ccGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCtg                                                                                                                                 <  1:6613259/150‑1 (MQ=255)
                 cGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCtgt                                                                                                                                <  2:5369033/150‑1 (MQ=255)
                 cGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCtgt                                                                                                                                <  2:6992314/150‑1 (MQ=255)
                 cGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCtgt                                                                                                                                <  1:4609112/150‑1 (MQ=255)
                  gTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCtgtg                                                                                                                               <  1:5623093/150‑1 (MQ=255)
                  gTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCtgtg                                                                                                                               <  1:5308489/150‑1 (MQ=255)
                  gTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCtgtg                                                                                                                               <  1:1150761/150‑1 (MQ=255)
                  gTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCtgtg                                                                                                                               <  1:735895/150‑1 (MQ=255)
                    aGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGat                                                                                                                             <  1:4365499/150‑1 (MQ=255)
                    aGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGat                                                                                                                             <  1:3690032/150‑1 (MQ=255)
                    aGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGat                                                                                                                             <  1:1575272/150‑1 (MQ=255)
                        cTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATcc                                                                                                                         <  2:107173/150‑1 (MQ=255)
                        cTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATcc                                                                                                                         <  2:2311843/150‑1 (MQ=255)
                         tCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACCGGGTGTGAATGCTGTGATATccc                                                                                                                        >  2:8855358/1‑150 (MQ=255)
                          cgcgTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCt                                                                                                                       <  2:1684286/150‑1 (MQ=255)
                          cgcgTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCt                                                                                                                       <  2:5767420/150‑1 (MQ=255)
                          cgcgTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCt                                                                                                                       <  2:874999/150‑1 (MQ=255)
                           gcgTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCtg                                                                                                                      <  2:3782190/150‑1 (MQ=255)
                           gcgTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCtg                                                                                                                      <  2:4752530/150‑1 (MQ=255)
                            cgTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCtgt                                                                                                                     <  1:3320515/150‑1 (MQ=255)
                            cgTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCtgt                                                                                                                     <  1:3445542/150‑1 (MQ=255)
                            cgTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCtgt                                                                                                                     <  2:5755069/150‑1 (MQ=255)
                            cgTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCtgt                                                                                                                     <  2:4901230/150‑1 (MQ=255)
                            cgTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCtgt                                                                                                                     <  1:3842530/150‑1 (MQ=255)
                            cgTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCtgt                                                                                                                     <  2:625651/150‑1 (MQ=255)
                             gTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCtgtg                                                                                                                    <  1:8719157/150‑1 (MQ=255)
                             gTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCtgtg                                                                                                                    <  1:7331163/150‑1 (MQ=255)
                               gCTGATGGAGGGCTGGGGTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGGAATCACAGGGTTTGAATGCTGGGATATACCTGgggg                                                                                                                  >  1:2973270/1‑150 (MQ=39)
                                cTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGGATATCACAGCATTCAAACCCTGTGATTACCTGTGGGa                                                                                                                 <  2:1207065/150‑1 (MQ=18)
                                cTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGGATATCACAGCATTCAAACCCTGTGATTACCTGTGGGa                                                                                                                 <  2:1841683/150‑1 (MQ=18)
                                  gATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCCCCGGGAAACACACGGGTTGAAAGCAGAGATATACCAGTGGGAGc                                                                                                               >  1:5138150/1‑150 (MQ=255)
                                  gATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCtg                                                                                                                                 <  2:3145280/132‑1 (MQ=255)
                                  gATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCtg                                                                                                                                 >  1:3145280/1‑132 (MQ=255)
                                  gATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGc                                                                                                               <  2:8342836/150‑1 (MQ=255)
                                  gATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGGATATCACAGCATTCAAACCCTGTGATTACCTGTGGGAGc                                                                                                               <  2:5900703/150‑1 (MQ=18)
                                    tGGAGGGCTGGGTTCGCGTGCCTGGCGCTAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGGATATCACAGCATTCAAACCCTGTGATTACCTGTGGGAGCgg                                                                                                             <  2:1310739/150‑1 (MQ=18)
                                    tGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCgg                                                                                                             <  1:8142914/150‑1 (MQ=255)
                                    tGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGGATATCACAGCATTCAAACCCTGTGATTACCTGTGGGAGCgg                                                                                                             <  2:5298316/150‑1 (MQ=18)
                                     ggAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCggg                                                                                                            <  2:907383/150‑1 (MQ=255)
                                      gAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCCCAGGGATATCACAGCATTCAAACCCTGTGATTACCTGTGGGAGCGGGc                                                                                                           <  2:6814019/150‑1 (MQ=18)
                                        gggCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCaa                                                                                                         <  2:2578738/150‑1 (MQ=255)
                                        gggCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCaa                                                                                                         <  2:2913745/150‑1 (MQ=255)
                                         ggCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAg                                                                                                        <  1:4693896/150‑1 (MQ=255)
                                         ggCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCGCCCCCACGTAAACACACGGTTTGGAATGTGGGTTATCACGGGGGGAGCGGGCacg                                                                                                        >  1:7034878/1‑148 (MQ=255)
                                            tGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGGTATCCCTGTGGGAGCGGGCAAGccc                                                                                                     >  2:1273294/1‑150 (MQ=255)
                                            tGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGccc                                                                                                     <  2:4287624/150‑1 (MQ=255)
                                                   gcgTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGGGGGCGCGGGCAAGCCCGCGAAgg                                                                                              >  2:8451293/1‑150 (MQ=255)
                                                          tGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGcac                                                                                       >  2:5919184/1‑150 (MQ=255)
                                                           ggCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACCGGGTTTGAATGCTGTTATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGcaca                                                                                      <  2:3749592/150‑1 (MQ=255)
                                                              gctAGCTTCTACAGCACAAGTGGCGCTCGGCTTCCAGAGGGGCCGCTTTGCTCCCATTGGCGGGATTCACCGCTCCGACAGTTAACCCATGTGTTTGATTGCTGTTATATACCTGGGGTAGCGGTCATGCCCGCGAAGGGCTGCAcagca                                                                                   <  2:7406123/148‑1 (MQ=255)
                                                               aTAGCTTCTAACGCACCACTGGCGCTCGGCTTCAAGAGGGGCGGCTTGGCGGCCCTTCGCGGGCTGGCCCGCTCCCCCAGGTAATCACTGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCAcagcag                                                                                  <  2:6557346/150‑1 (MQ=38)
                                                                     tCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCgg                                                                            >  2:4814877/1‑150 (MQ=255)
                                                                      gccgaGCGCCAGTGGTGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCACCCACAGGTAATAACAGGGTTTGAATTCTGTGATATCCCAGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGc                                                                           >  2:190250/5‑150 (MQ=37)
                                                                      caaaaGCACCACTGGCGCTCGGCTACAAGAGGGGCCGCTTGGCGGCCCTTCGGGGGCTTGCCCGCTCCCACAGTTAACCACAGGTTTTGAATGCTGGGATAACCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGc                                                                           <  2:5500684/148‑1 (MQ=37)
                                                                        aaaGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCCTCGCGGGCTTGCCCGCTCCCACCGGTAATCACAGGGTTTGAATGCTGTTATACCCCGGTGGGCGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGcccttcgcg                                                                         <  2:2747607/150‑7 (MQ=38)
                                                                         aaGCACCAGTGGCGCTCGGCTACAGGAGGGGCCGCTTGGCGGCCCTGCGGGGGGTTGCCCGCTCCCACTGGTAATCACTGGGTTGGAAGGCTGTTATATCCGTGTGGGTGCGGGCAAGCCCGCGAAGGGCTGCTGTGCAGcccttcgcgg                                                                        <  2:7858575/150‑8 (MQ=255)
                                                                         aaGCACCACGGGCGCTCGGCTTCAAGTGGGGCCGCTTGGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACCGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCTAGCCCGCGAAGGGCTGCTGTGCAGcccttcgcgg                                                                        <  2:7858636/150‑8 (MQ=37)
                                                                          aGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACCGGGTTTGAATGCGGTGATATCCCGGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCaaat                                                                       >  2:1540986/1‑148 (MQ=255)
                                                                               cacTGGCGCTCGGATTCATGCGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCACCCACCGTTAATCACAGGTTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGaa                                                                  <  2:2567575/150‑1 (MQ=37)
                                                                                  tGGCTCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAggg                                                               <  2:5143410/150‑1 (MQ=255)
                                                                                        tgggCTTCATGTGGGGCGGCTTTGCGGCCCTTGGCGGGGTTGCCGGCTCCCACTGGTAATCACTGGTTTTGAATGCTGTGATTGCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTgg                                                         <  1:9024983/148‑1 (MQ=25)
                                                                                        tCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTgg                                                         >  2:7459876/1‑150 (MQ=255)
                                                                                              tttagagGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCCATTCTGAAGGGAGCTGGATATGa                                                   >  2:3641757/4‑150 (MQ=255)
                                                                                                atgcggggCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCCCCCACAGTAAATAACCGGTTTTGATTGCTGTCATATCCCGGGGGGCGGGGGCACGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAgc                                                 <  1:4466228/146‑1 (MQ=255)
                                                                                                aatagGGGCAGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAAAAAAAGGGTTAGAATGCGGTGATATCCCTGTGGGACCCGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAgc                                                 >  2:5175912/4‑150 (MQ=37)
                                                                                                 agagGGGCCGCTTTGCGGCCCTTCGCGGGCTTTCCCCCACCCACAGGTAATAACAGGGTTTGGATGCTGTGATAAACCTGTGGGAGCGGGCCAGCCCGCGAAGGGCTGCACCGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGGgcg                                                >  2:3185917/1‑150 (MQ=37)
                                                                                                         cGCTTTGCGGCCCTTCGCGGGCTTTCCCGCACCCAAAGGGAATAACCGGGTTTTCATGCGGTGATAAACCTGGGGGTGCCGGCAAACCCGCGAACGGCCGCAAAGAAGCCACCGCCAATGTGTAGGGAGCTGGATATGCTCGcccaactg                                        >  2:4058435/1‑144 (MQ=255)
                                                                                                                     cTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGAt                            >  2:9026722/1‑150 (MQ=255)
                                                                                                                           gggCTTGCCCGCTCCCACAGGGATATCACAGCATTCAAACCCTGTGATTACCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCg                      >  1:2800268/1‑150 (MQ=17)
                                                                                                                             gCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGa                    >  1:817257/1‑150 (MQ=255)
                                                                                                                             gCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGa                    >  2:705014/1‑150 (MQ=255)
                                                                                                                                tGCCCGCTCCCACAGGTAACCACCGGGTTTGAATGCTGTGATATCCCTGGGGGCGGGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTa                 <  1:2210728/150‑1 (MQ=255)
                                                                                                                                   ccGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGcccc                                                                              <  2:4525817/86‑1 (MQ=255)
                                                                                                                                   ccGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGcccc                                                                              >  1:4525817/1‑86 (MQ=255)
                                                                                                                                   ccGCTCCCACAGGGATATCACAGCATTCAAACCCTGTGATTACCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGa              >  2:6774989/1‑150 (MQ=18)
                                                                                                                                   cagCTCCCACAGGGATATCACAGCATTCAAACCCTGTGATTAACTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCACCGGCAATTCTGAAGGGAGCTGGATATGAGCGACAACGAAGACCTCAACGATCGCACGGATTACGa              >  2:2740586/3‑150 (MQ=17)
                                                                                                                                     gCTCCCACAGGGATATCACAGCATTCAAACCCTGTGATTACCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGAcc            >  1:1294990/1‑150 (MQ=18)
                                                                                                                                       tCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCgg          >  2:1616541/1‑150 (MQ=255)
                                                                                                                                       tCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTATGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCgg          >  1:3083138/1‑150 (MQ=255)
                                                                                                                                       acccACCGTTAATCACTGGTTTTGATTGGTTTTTTACCCCTGTGGTTGGGGGAAACCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCgg          <  1:8773172/149‑1 (MQ=255)
                                                                                                                                         ccACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGt        >  2:6670781/1‑150 (MQ=255)
                                                                                                                                           acaGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGccc                       >  1:2697831/1‑133 (MQ=255)
                                                                                                                                           acaGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGc      >  2:6427687/1‑150 (MQ=255)
                                                                                                                                             cggTACTCACAGGGTTTGAAGGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTa    <  2:8903177/149‑1 (MQ=255)
                                                                                                                                               gTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACa  >  1:5521444/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
CCAAGGCAATCATGAGCCGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACA  >  NC_002947/2669874‑2670166

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: