Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,670,032 G→A 11.4% intergenic (+85/‑82) 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,0320GA11.4% 404.8 / 33.1 151intergenic (+85/‑82)prpF/prpDaconitate isomerase/2‑methylcitrate dehydratase
Reads supporting (aligned to +/- strand):  ref base G (53/78);  new base A (6/11);  total (60/91)
Fisher's exact test for biased strand distribution p-value = 7.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.99e-01

Only 100 of 151 total aligned reads displayed.

CATGAGCCGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTGGCCGA  >  NC_002947/2669884‑2670178
                                                                                                                                                     |                                                                                                                                                  
cATGAGCCGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGa                                                                                                                                                   <  2:5794411/150‑1 (MQ=255)
 aTGAGCCGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCTCGTGCCTGGCGATAGCTTCTAAAGCACCACGGGCGCTCGGCTTCAAGAGGGGCCTCTTTGCGGCCCTTCTCGGGCTTGACCGCTCCCTCAGGTTATCACAGGGTTTGaa                                                                                                                                                  <  2:4560140/150‑1 (MQ=255)
 aTGAGCCGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGaa                                                                                                                                                  <  1:5139203/150‑1 (MQ=255)
 aTGAGCCGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCCCTCCCACACGTAATCACAGGGTTTGaa                                                                                                                                                  >  1:3538806/1‑150 (MQ=255)
    aGCCGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCTCTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGc                                                                                                                                               <  1:6894515/150‑1 (MQ=255)
     gCCGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCt                                                                                                                                              <  1:2636277/150‑1 (MQ=255)
     gCCGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCt                                                                                                                                              <  1:7774743/150‑1 (MQ=255)
     gCCGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCt                                                                                                                                              <  2:7854433/150‑1 (MQ=255)
      ccGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCtg                                                                                                                                             <  1:6748941/150‑1 (MQ=255)
       cGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCtgt                                                                                                                                            <  2:6992314/150‑1 (MQ=255)
       cGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCtgt                                                                                                                                            <  1:4609112/150‑1 (MQ=255)
       cGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCtgt                                                                                                                                            <  1:4160724/150‑1 (MQ=255)
        gTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCtgtg                                                                                                                                           <  1:1150761/150‑1 (MQ=255)
        gTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCtgtg                                                                                                                                           <  1:5623093/150‑1 (MQ=255)
          aGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGat                                                                                                                                         <  1:3690032/150‑1 (MQ=255)
          aGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGat                                                                                                                                         <  1:1575272/150‑1 (MQ=255)
              cTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATcc                                                                                                                                     <  2:2311843/150‑1 (MQ=255)
               tCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACCGGGTGTGAATGCTGTGATATccc                                                                                                                                    >  2:8855358/1‑150 (MQ=255)
                cgcgTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCt                                                                                                                                   <  2:874999/150‑1 (MQ=255)
                cgcgTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCt                                                                                                                                   <  2:1684430/150‑1 (MQ=255)
                cgcgTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCt                                                                                                                                   <  2:5767420/150‑1 (MQ=255)
                 gcgTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCtg                                                                                                                                  <  2:3782190/150‑1 (MQ=255)
                 gcgTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCtg                                                                                                                                  <  2:4752530/150‑1 (MQ=255)
                  cgTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCtgt                                                                                                                                 <  2:8224094/150‑1 (MQ=255)
                  cgTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCtgt                                                                                                                                 <  1:3320515/150‑1 (MQ=255)
                  cgTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCtgt                                                                                                                                 <  2:625651/150‑1 (MQ=255)
                   gTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCtgtg                                                                                                                                <  1:7331163/150‑1 (MQ=255)
                   gTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCtgtg                                                                                                                                <  1:8719157/150‑1 (MQ=255)
                     gCTGATGGAGGGCTGGGGTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGGAATCACAGGGTTTGAATGCTGGGATATACCTGgggg                                                                                                                              >  1:2973270/1‑150 (MQ=39)
                      cTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGGATATCACAGCATTCAAACCCTGTGATTACCTGTGGGa                                                                                                                             <  2:1207065/150‑1 (MQ=18)
                      cTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGGATATCACAGCATTCAAACCCTGTGATTACCTGTGGGa                                                                                                                             <  2:5222915/150‑1 (MQ=18)
                      cTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGGATATCACAGCATTCAAACCCTGTGATTACCTGTGGGa                                                                                                                             <  2:1841683/150‑1 (MQ=18)
                        gATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCACGAGGGGCCGCTTTGCGGCCCATCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGc                                                                                                                           <  2:3839522/150‑1 (MQ=255)
                        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:2913745/150‑1 (MQ=255)
                               ggCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAg                                                                                                                    <  1:4693896/150‑1 (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)
                                                    gctAGCTTCTACAGCACAAGTGGCGCTCGGCTTCCAGAGGGGCCGCTTTGCTCCCATTGGCGGGATTCACCGCTCCGACAGTTAACCCATGTGTTTGATTGCTGTTATATACCTGGGGTAGCGGTCATGCCCGCGAAGGGCTGCAcagca                                                                                               <  2:7406123/148‑1 (MQ=255)
                                                     aTAGCTTCTAACGCACCACTGGCGCTCGGCTTCAAGAGGGGCGGCTTGGCGGCCCTTCGCGGGCTGGCCCGCTCCCCCAGGTAATCACTGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCAcagcag                                                                                              <  2:6557346/150‑1 (MQ=38)
                                                         cTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCCCCGGGAATCACCGGGGTTGTATGCGGGGGGATCCCGGGGGGGGCGGGCAAGCCCGCGCAGGGCGGCGCAGCCGcccc                                                                                          >  1:8903177/1‑150 (MQ=255)
                                                           tCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCgg                                                                                        >  2:4814877/1‑150 (MQ=255)
                                                            caaaaGCACCACTGGCGCTCGGCTACAAGAGGGGCCGCTTGGCGGCCCTTCGGGGGCTTGCCCGCTCCCACAGTTAACCACAGGTTTTGAATGCTGGGATAACCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGc                                                                                       <  2:5500684/148‑1 (MQ=37)
                                                              aaaGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCCTCGCGGGCTTGCCCGCTCCCACCGGTAATCACAGGGTTTGAATGCTGTTATACCCCGGTGGGCGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGcccttcgcg                                                                                     <  2:2747607/150‑7 (MQ=38)
                                                               aaGCACCAGTGGCGCTCGGCTACAGGAGGGGCCGCTTGGCGGCCCTGCGGGGGGTTGCCCGCTCCCACTGGTAATCACTGGGTTGGAAGGCTGTTATATCCGTGTGGGTGCGGGCAAGCCCGCGAAGGGCTGCTGTGCAGcccttcgcgg                                                                                    <  2:7858575/150‑8 (MQ=255)
                                                                aGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACCGGGTTTGAATGCGGTGATATCCCGGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCaaat                                                                                   >  2:1540986/1‑148 (MQ=255)
                                                                    ccacTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTGGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGa                                                                               <  2:8921123/150‑1 (MQ=255)
                                                                     cacTGGCGCTCGGATTCATGCGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCACCCACCGTTAATCACAGGTTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGaa                                                                              <  2:2567575/150‑1 (MQ=37)
                                                                        tGGCTCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAggg                                                                           <  2:5143410/150‑1 (MQ=255)
                                                                              tCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTgg                                                                     >  2:7459876/1‑150 (MQ=255)
                                                                                    tttagagGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCCATTCTGAAGGGAGCTGGATATGa                                                               >  2:3641757/4‑150 (MQ=255)
                                                                                      atgcggggCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCCCCCACAGTAAATAACCGGTTTTGATTGCTGTCATATCCCGGGGGGCGGGGGCACGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAgc                                                             <  1:4466228/146‑1 (MQ=255)
                                                                                       agagGGGCCGCTTTGCGGCCCTTCGCGGGCTTTCCCCCACCCACAGGTAATAACAGGGTTTGGATGCTGTGATAAACCTGTGGGAGCGGGCCAGCCCGCGAAGGGCTGCACCGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGGgcg                                                            >  2:3185917/1‑150 (MQ=37)
                                                                                          ggggCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGTTTTGATTGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCa                                                         <  2:562817/150‑1 (MQ=255)
                                                                                             gCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAAACACAGGGTTTGAATGCTGTGATATCCCGGTGGGAGCGGGCAAGCCCGCGAAGGGCGGCACCGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACg                                                      >  2:4917513/1‑150 (MQ=255)
                                                                                                gCTTTGCGGCCCTTCGCGGGCGTGCCCGCTCCCACAGGTAAACACAGGGTTTGAATGCTGTGATATACCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAAATCTGAAGGGAGCTGGATATGAGCGCCAACGTAg                                                   >  2:6406397/1‑150 (MQ=255)
                                                                                                           cTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGAt                                        >  2:9026722/1‑150 (MQ=255)
                                                                                                            ttCGCGGGCTTGCCCCCACCCCCAGGTAATCACAGGGTTTGCATTGTGTGATATACCTGTGGGGGCGGGCAAGCCCGCCGATGGCTGCAAAGCCGCCCCGGCAAATCTGCAGGGAGCTGGATATGAGCGCCAACGGAGACCTCAACGATc                                       >  2:1946305/1‑150 (MQ=255)
                                                                                                                cgGGCTGGCCCGCTCCCCCCGGTACGCCCAGGGTTTGCATGCTGTGATATCCCTGGGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGccc                                   <  1:4475976/150‑1 (MQ=255)
                                                                                                                 gggCTTGCCCGCTCCCACAGGGATATCACAGCATTCAAACCCTGTGATTACCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCACAACGATCGCCCg                                  >  1:7333969/1‑150 (MQ=17)
                                                                                                                   gCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGa                                >  1:817257/1‑150 (MQ=255)
                                                                                                                   gCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGa                                >  2:705014/1‑150 (MQ=255)
                                                                                                                     ttGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGAtt                              >  1:1326344/1‑150 (MQ=255)
                                                                                                                      tGCCCGCTCCCACAGGTAACCACCGGGTTTGAATGCTGTGATATCCCTGGGGGCGGGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTa                             <  1:2210728/150‑1 (MQ=255)
                                                                                                                        cccGCTCCCACAGGGATATCACAGCATTCAAACCCTGTGATTACCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACg                           >  1:2837963/1‑150 (MQ=18)
                                                                                                                         ccGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGcccc                                                                                          <  2:4525817/86‑1 (MQ=255)
                                                                                                                         ccGCTCCCACAGGGATATCACAGCATTCAAACCCTGTGATTACCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGa                          >  2:6774989/1‑150 (MQ=18)
                                                                                                                             tCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCgg                      >  2:1616541/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)
                                                                                                                                 acaGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGccc                                   <  2:2697831/133‑1 (MQ=255)
                                                                                                                                     gTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACa              >  1:5521444/1‑150 (MQ=255)
                                                                                                                                            acaGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCt        >  1:7060944/1‑150 (MQ=255)
                                                                                                                                            acaGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCt        >  2:4551572/1‑150 (MQ=255)
                                                                                                                                            acaGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCt        >  2:2678249/1‑150 (MQ=255)
                                                                                                                                            acaGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCAGGTGCTACAAACCCt        >  1:1980654/1‑150 (MQ=255)
                                                                                                                                             cggggTTTTAATGCGGTCATCTCCCGGTGGGCGCGGCCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTg       <  1:8855358/148‑1 (MQ=255)
                                                                                                                                             ccgggTTTTAATGCTGTGATATCCCTGGGGGTGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTg       <  2:1206670/148‑1 (MQ=255)
                                                                                                                                             caGGGTTTGAATGCTGTGATATCCCTGTGTGATCGGGCAAGCCCGCGAAGGGCTGCACATCAGACACGGCAAATATGAAGTGAGCTGGATATGAGCGCCAACGTAGACATCAACGAACGCCCGGATTACTAACGGGTGCTACAAACCCTg       >  1:1286945/1‑150 (MQ=255)
                                                                                                                                             caGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGTTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTg       >  2:4344846/1‑150 (MQ=255)
                                                                                                                                             caGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTg       >  2:7914685/1‑150 (MQ=255)
                                                                                                                                              cgggTTCGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTgg      <  2:3538806/149‑1 (MQ=255)
                                                                                                                                              aGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGGGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTgg      >  1:7107797/1‑150 (MQ=255)
                                                                                                                                              aGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTgg      >  1:7520399/1‑150 (MQ=255)
                                                                                                                                                ggTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTGGcc    >  1:4263329/1‑150 (MQ=255)
                                                                                                                                                 gTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTGGCCg   >  2:4225152/1‑150 (MQ=255)
                                                                                                                                                  tttGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTGGCCGa  >  1:9135824/1‑150 (MQ=255)
                                                                                                                                                  tttGAATGCTGGGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTGGCCgt  >  1:9135770/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
CATGAGCCGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTGGCCGA  >  NC_002947/2669884‑2670178

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: