Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,670,027 G→C 11.8% intergenic (+80/‑87) 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,0270GC11.8% 437.7 / 40.4 144intergenic (+80/‑87)prpF/prpDaconitate isomerase/2‑methylcitrate dehydratase
Reads supporting (aligned to +/- strand):  ref base G (50/76);  new base C (6/11);  total (56/88)
Fisher's exact test for biased strand distribution p-value = 7.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.22e-01

Only 100 of 150 total aligned reads displayed.

GGCAATCATGAGCCGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTGGCC  >  NC_002947/2669878‑2670176
                                                                                                                                                      |                                                                                                                                                      
    aTCATGAGCCGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGttt                                                                                                                                                    <  2:4208774/150‑1 (MQ=255)
    aTCATGAGCCGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGttt                                                                                                                                                    <  1:7988753/150‑1 (MQ=255)
      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)
           gCCGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCt                                                                                                                                            <  1:2636277/150‑1 (MQ=255)
            ccGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCtg                                                                                                                                           <  1:6748941/150‑1 (MQ=255)
            ccGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCtg                                                                                                                                           <  1:6613259/150‑1 (MQ=255)
            ccGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCtg                                                                                                                                           <  1:8812577/150‑1 (MQ=255)
             cGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACTACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCtgt                                                                                                                                          <  1:8475594/150‑1 (MQ=255)
             cGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCtgt                                                                                                                                          <  1:4160724/150‑1 (MQ=255)
             cGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCtgt                                                                                                                                          <  1:4609112/150‑1 (MQ=255)
             cGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCtgt                                                                                                                                          <  2:8490370/150‑1 (MQ=255)
              gTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCtgtg                                                                                                                                         <  1:1150761/150‑1 (MQ=255)
              gTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCtgtg                                                                                                                                         <  1:5623093/150‑1 (MQ=255)
              gTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCtgtg                                                                                                                                         <  1:735895/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:1684430/150‑1 (MQ=255)
                      cgcgTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCt                                                                                                                                 <  2:5767420/150‑1 (MQ=255)
                       gcgTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCtg                                                                                                                                <  2:4752530/150‑1 (MQ=255)
                       gcgTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCtg                                                                                                                                <  2:3782190/150‑1 (MQ=255)
                        cgTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCtgt                                                                                                                               <  1:3320515/150‑1 (MQ=255)
                        cgTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCtgt                                                                                                                               <  2:8224094/150‑1 (MQ=255)
                         gTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCtgtg                                                                                                                              <  1:7331163/150‑1 (MQ=255)
                           gCTGATGGAGGGCTGGGGTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGGAATCACAGGGTTTGAATGCTGGGATATACCTGgggg                                                                                                                            >  1:2973270/1‑150 (MQ=39)
                            cTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGGATATCACAGCATTCAAACCCTGTGATTACCTGTGGGa                                                                                                                           <  2:1841683/150‑1 (MQ=18)
                            cTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGGATATCACAGCATTCAAACCCTGTGATTACCTGTGGGa                                                                                                                           <  2:5222915/150‑1 (MQ=18)
                            cTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGGATATCACAGCATTCAAACCCTGTGATTACCTGTGGGa                                                                                                                           <  2:1207065/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)
                              gATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGGATATCACAGCATTCAAACCCTGTGATTACCTGTGGGAGc                                                                                                                         <  2:5900703/150‑1 (MQ=18)
                                tGGAGGGCTGGGTTCGCGTGCCTGGCGCTAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGGATATCACAGCATTCAAACCCTGTGATTACCTGTGGGAGCgg                                                                                                                       <  2:1310739/150‑1 (MQ=18)
                                tGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCgg                                                                                                                       <  1:8142900/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)
                                     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)
                                                           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)
                                                                     aaGCACCACGGGCGCTCGGCTTCAAGTGGGGCCGCTTGGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACCGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCTAGCCCGCGAAGGGCTGCTGTGCAGcccttcgcgg                                                                                  <  2:7858636/150‑8 (MQ=37)
                                                                      aGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACCGGGTTTGAATGCGGTGATATCCCGGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCaaat                                                                                 >  2:1540986/1‑148 (MQ=255)
                                                                           cccgggggcTGGGCTTCACGCGGGGCCGGTTTGGGGCCCTGGGGGGGCTTGCCCGCGCCCACGGTTAAGCCCTGGGTTTGTTTGCTGTTTTTTCCCTGTGGGCGGGGGAACGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGaa                                                                            <  2:2566853/143‑1 (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:7340876/150‑1 (MQ=255)
                                                                                    tCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTgg                                                                   >  2:7459876/1‑150 (MQ=255)
                                                                                    tCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTgg                                                                   <  2:7588144/150‑1 (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)
                                                                                                      gCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACCGGTAATCACAGGGTTTGAATGCGGGGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCGGCCCAGCCGCCCCGGCACTTCTGAAGGGAGCTGGATATGTGCGCCAACGTAg                                                 >  2:6716790/1‑150 (MQ=37)
                                                                                                      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)
                                                                                                                         gCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGa                              >  1:817257/1‑150 (MQ=255)
                                                                                                                         gCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGa                              >  2:6834373/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)
                                                                                                                                   tCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCgg                    >  2:1616541/1‑150 (MQ=255)
                                                                                                                                   acccACCGTTAATCACTGGTTTTGATTGGTTTTTTACCCCTGTGGTTGGGGGAAACCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCgg                    <  1:8773172/149‑1 (MQ=255)
                                                                                                                                       acaGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGccc                                 >  1:2697831/1‑133 (MQ=255)
                                                                                                                                       acaGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGccc                                 <  2:2697831/133‑1 (MQ=255)
                                                                                                                                       acaGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGc                >  2:6427687/1‑150 (MQ=255)
                                                                                                                                         cggTACTCACAGGGTTTGAAGGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTa              <  2:8903177/149‑1 (MQ=255)
                                                                                                                                           gTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACa            >  1:5521444/1‑150 (MQ=255)
                                                                                                                                                 cactGGGTTTGAATGCTGTGATTTCCCTGTGGGCGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAAccc       <  2:831733/150‑1 (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)
                                                                                                                                                   ccgggTTTTAATGCTGTGATATCCCTGGGGGTGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTg     <  2:1206670/148‑1 (MQ=255)
                                                                                                                                                   caGGGTTTGAATGCTGTGATATCCCTGTGTGATCGGGCAAGCCCGCGAAGGGCTGCACATCAGACACGGCAAATATGAAGTGAGCTGGATATGAGCGCCAACGTAGACATCAACGAACGCCCGGATTACTAACGGGTGCTACAAACCCTg     >  1:1286945/1‑150 (MQ=255)
                                                                                                                                                   caGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGTTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTg     >  2:4344846/1‑150 (MQ=255)
                                                                                                                                                    cgggTTCGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTgg    <  2:3538806/149‑1 (MQ=255)
                                                                                                                                                    aGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTgg    >  1:7174036/1‑150 (MQ=255)
                                                                                                                                                    aGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTgg    >  1:7520399/1‑150 (MQ=255)
                                                                                                                                                      ggTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTGGcc  >  1:4263329/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                      
GGCAATCATGAGCCGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTGGCC  >  NC_002947/2669878‑2670176

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: