Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,670,044 Δ1 bp 8.2% intergenic (+97/‑70) 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,0440A. 8.2% 262.2 / 18.9 85intergenic (+97/‑70)prpF/prpDaconitate isomerase/2‑methylcitrate dehydratase
Reads supporting (aligned to +/- strand):  ref base A (41/36);  new base . (4/3);  total (46/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.70e-01

GTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTGGCCGATTACGCCCTCGGCTA  >  NC_002947/2669895‑2670193
                                                                                                                                                      |                                                                                                                                                     
gTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGata                                                                                                                                                       <  2:4763085/150‑1 (MQ=255)
gTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGata                                                                                                                                                       <  2:4763263/150‑1 (MQ=255)
     cgcgTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCt                                                                                                                                                  <  2:518837/150‑1 (MQ=255)
       cgTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCtgt                                                                                                                                                <  2:397025/150‑1 (MQ=255)
       cgTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCtgt                                                                                                                                                <  1:4187656/150‑1 (MQ=255)
        gTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCtgtg                                                                                                                                               <  1:4838700/150‑1 (MQ=255)
         tGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTgg                                                                                                                                              <  1:1182575/150‑1 (MQ=255)
               tGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCgg                                                                                                                                        <  2:2350790/150‑1 (MQ=255)
               tGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGGATATCACAGCATTCAAACCCTGTGATTACCTGTGGGAGCgg                                                                                                                                        <  2:242530/150‑1 (MQ=18)
               tGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGGATATCACAGCATTCAAACCCTGTGATTACCTGTGGGAGCgg                                                                                                                                        <  1:2389049/150‑1 (MQ=18)
                 gAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGACCCCGACCACACGGAATAAAACGGGTTGAAATCTTTGATAAACCAGGGGGGGCgggg                                                                                                                                      >  1:4143739/1‑149 (MQ=255)
                             cgcgTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCTCTCCCACAGGGATATCACAGCATTCAAACCCTGTGATTACCTGTGGGAGCGGGCAAGCCCGCGAAg                                                                                                                          <  1:3137705/150‑1 (MQ=17)
                                    cTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCGTGCCCGCTCCCCCAGGGAATCACAGGGGTTGAAAGGTGTGATATCCCGGTGGGGGCGGGCAAGCCCGCGAAGGGCTGcc                                                                                                                   >  2:1984693/1‑149 (MQ=255)
                                     tGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGGAATCACAGGGTTTGAATGCTGTGATAACCCTGGGGGAGCGGGCAAGCCCGCGAAGGGCTGcac                                                                                                                  >  2:5361063/1‑150 (MQ=255)
                                     tGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCGTGGCCGCGCCCCCACGTAAAAACAGGGGTTTTATTGGGTGTTATACCGGTGGGGGCGGGGCAGCCCGCGACGCGCGGcac                                                                                                                  >  1:3213934/1‑150 (MQ=255)
                                          aTAGCTTCTAAAGCACCACTGGCGGTCGGCTTCAAGTGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGTTAATCACAGTGTTTGAATGCTTTTATATCCCTGTGGGAGCGGGAAAGCCCGCGAAGGGCTGCAcagcag                                                                                                             <  2:4314938/150‑1 (MQ=38)
                                              cTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACCGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGcccc                                                                                                         >  2:3333994/1‑150 (MQ=255)
                                              cTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTgg                                                                                                                                              >  1:3713480/1‑113 (MQ=255)
                                              cTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTgg                                                                                                                                              <  2:3713480/113‑1 (MQ=255)
                                               ttCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTGCCACAGGTAATAAAAGGGTTTGAATTCTTTGATATGCCTGTGGGTGCGGGCAAGCCCGCGTAGGGCTGGACAGCAGCCCCg                                                                                                        >  1:868967/1‑150 (MQ=38)
                                               ttCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCGTGCCCGCTCCCCCACGGAAACACAGGGGTTGAATGCGGGGCTATCCCGCTGGGAGCCGGCCAGCCCCCGCACGGCAGCACACCCGccccc                                                                                                        >  2:2806904/1‑149 (MQ=255)
                                                   aaaGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCACCCACCGGTAATCACAGGGTTTGAATGCTGTGGTATCCCGGTGGGAGCGGGCCAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCac                                                                                                    >  2:1033546/1‑149 (MQ=38)
                                                          cccgggCGCTCGGCTTCATGAGGGGCCGGTTGGCGGCCCTTCGGGGGCTTGCCCGCTCCCACCGTTAATCACTGGTTTTTATTGCTGTGATATCCCTGGGGGAGGGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGaa                                                                                             <  2:2742314/146‑1 (MQ=255)
                                                                  cTCGGCTCCATGTGGGGCCGCTTGGCGGCCCTTCGCGGGCTTGCCCGACCCCACTGTTAAACACGGGTTTTTATTGGTGTTTTAGCCCTTTGGGTGGGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTg                                                                                     <  2:1814405/150‑1 (MQ=255)
                                                                      gCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACCGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGata                                                                                 <  2:2668365/150‑1 (MQ=255)
                                                                      gCTTCAAGAGGGGCCGCTGTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACCGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGata                                                                                 <  2:5437876/150‑1 (MQ=255)
                                                                         tCACGAGGGGCCGCTGTGGGGCCCTTCGCGGGCTTGCCCGCTCCCACCGGTAATCACCGGGTTTGACGGCTGTGAGATCCCGGTGGGAGCGGGCACGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGa                                                                              <  2:4775415/150‑1 (MQ=37)
                                                                                gggCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACCGGTAATCACCGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCACGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCaa                                                                       <  2:4771821/150‑1 (MQ=255)
                                                                                  gCCGCTTGGCGGCCCTTCGCGGGCTGGCCCGCCCCCACCGGTAACCACCGGGTTTGAAGGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACg                                                                     <  2:137012/150‑1 (MQ=38)
                                                                                   ccGCTTTGCGGCCCTGCGCGGGCTTGCCCGCTCCCACCGGTAATCACCGGGTTTGACTGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGt                                                                    <  1:3991677/150‑1 (MQ=255)
                                                                                      cTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGa                                                                 <  2:3608901/150‑1 (MQ=255)
                                                                                           cGGCCCTTCGCGGGCTTGCCCCCTCCCACCGGTAAACACAGGGGTTGAATGCGGGGATATACCGGTGGGACCGGGCAAGCCCCCGAACGGGGGCACACCCCCCCCCGCCAATATGCAGGGAGCTGGGAATGAGCGCCAACGGACAccgcc                                                            >  2:5050103/1‑147 (MQ=255)
                                                                                             gCCCTTGGCGGGGTTGCCCGCTCCCACCGTTAATCACCGGTTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAAc                                                          <  2:1911619/150‑1 (MQ=38)
                                                                                               ccTTCGGGGGCTTGCCCGCTCCCACAGGTAATCACTGGGTTTGAATGCTGTGATATCCCTGTTGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGa                                                        <  1:2721749/150‑1 (MQ=255)
                                                                                               ccTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTAGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGa                                                        <  1:5444995/150‑1 (MQ=255)
                                                                                                 ttCGCGGGCTTCCCCGCTCCCACAGGTAAGCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATc                                                      <  2:3248367/150‑1 (MQ=255)
                                                                                                        gattGCCCGCTCCCACAGTTAATCACAGGGTTTGAATGCTGTTATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGa                                               <  1:3600241/148‑1 (MQ=255)
                                                                                                          ttGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGTAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGAtt                                             >  2:2472293/1‑150 (MQ=255)
                                                                                                          ttGCCCGCTCCCACAGGGATATCACAGCATTCAAACCCTGTGATTACCTGTGGGAGCGGGCAAGCCCGCGAAGGGCCGCAAAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGAtt                                             >  2:1616719/1‑150 (MQ=14)
                                                                                                            gCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGAt                                              >  1:1291288/1‑147 (MQ=255)
                                                                                                            gCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGAt                                              <  2:1291288/147‑1 (MQ=255)
                                                                                                            gCCCGCTCCCACAGGGATATCACAGCATTCAAACCCTGTGATTACCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGAGTAc                                           >  2:5066603/1‑150 (MQ=17)
                                                                                                             cccGCTCCCACAGGGATATCACAGCATTCAAACCCTGTGAGTACCTGTGGGAGCGGGCAAGCCCGCGAAGGGCCGCACAGCAGCCCCGGCAAGTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACg                                          >  2:2542293/1‑150 (MQ=17)
                                                                                                              ccGCTCCCACAGGGATATCACAGCATTCAAACCCTGTGATTACCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGa                                         >  1:4361214/1‑150 (MQ=18)
                                                                                                                gCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGAcc                                       >  1:2402884/1‑150 (MQ=255)
                                                                                                                gCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGAcc                                       >  2:4456918/1‑150 (MQ=255)
                                                                                                                  tCCCACCGTTAATCACTGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCgg                                     <  1:5344581/150‑1 (MQ=255)
                                                                                                                  tCCCACAGGTAATCCCAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCgg                                     <  2:2570793/150‑1 (MQ=255)
                                                                                                                  tCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCgg                                     >  2:2053313/1‑150 (MQ=255)
                                                                                                                  tCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCgg                                     >  2:4099242/1‑150 (MQ=255)
                                                                                                                  tCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCgg                                     >  2:2053182/1‑150 (MQ=255)
                                                                                                                  tCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCgg                                     >  2:2052969/1‑150 (MQ=255)
                                                                                                                  tCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCgg                                     >  1:2363627/1‑150 (MQ=255)
                                                                                                                  tCCCACAGGGATATCACAGCATTCAAACCCTGTGATTACCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCgg                                     >  1:941755/1‑150 (MQ=18)
                                                                                                                    ccACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGt                                   >  1:4392596/1‑150 (MQ=255)
                                                                                                                      acaGTTAATCCCCGGGTTTGAATGCTGTTATATCCCTGTGGGCGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGc                                 <  1:1785633/150‑1 (MQ=255)
                                                                                                                      acaGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTg                                                                                                                     <  2:473681/66‑1 (MQ=255)
                                                                                                                      acaGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTg                                                                                                                     >  1:473681/1‑66 (MQ=255)
                                                                                                                         ggTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGaa                                                                                             >  1:2741031/1‑87 (MQ=255)
                                                                                                                         ggTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGaa                                                                                             <  2:2741031/87‑1 (MQ=255)
                                                                                                                          gTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACa                             >  1:2212451/1‑150 (MQ=255)
                                                                                                                               tCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACg                                                         >  1:3874954/1‑118 (MQ=255)
                                                                                                                               tCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACg                                                         <  2:3874954/118‑1 (MQ=255)
                                                                                                                               tCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTATGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTAAACGATCGCCCGGAATACGACCGGGTGCTACAAAcc                         >  1:4549255/1‑150 (MQ=255)
                                                                                                                                cactGGTTTTGATTGCTGTTATATCCGTGGGGGCGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAAccc                        <  1:335561/150‑1 (MQ=255)
                                                                                                                                cacaGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAAccc                        >  1:4114615/1‑150 (MQ=255)
                                                                                                                                cacaGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAAccc                        >  2:682786/1‑150 (MQ=255)
                                                                                                                                 ccaGGTTTTTAATGCTGTTAAATCCCTGTGGGAGGGGGAATGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCt                       <  2:97019/149‑1 (MQ=255)
                                                                                                                                 actgggTTTGATGGGTGTGATATCCGTGGGGGCGGGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCt                       <  2:1070969/147‑1 (MQ=255)
                                                                                                                                 acaGGGTTTGAATGCTGTGATATCCCTGTGGGATCGGGCAAGCCCGAGAAGGGCTGCACAGCAGCCCCGGCACTTCTGAAGGGAGCTGGATATGAGAGCCAACGTAGACCTCAACGATCGCCCGGATTACGAACGGGTGCTACAAACCCt                       >  2:4226584/1‑150 (MQ=39)
                                                                                                                                 acaGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCt                       >  2:2876283/1‑150 (MQ=255)
                                                                                                                                 acaGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCt                       >  2:948721/1‑150 (MQ=255)
                                                                                                                                 acaGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCt                       >  2:1660160/1‑150 (MQ=255)
                                                                                                                                  caGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTg                      >  1:5233777/1‑150 (MQ=255)
                                                                                                                                  caGGGTGTGTAGGCTGTGATATCCCTGTGGGTGGGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTg                      <  2:4253357/150‑1 (MQ=255)
                                                                                                                                     ggTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCgg                                     <  1:3008725/132‑1 (MQ=255)
                                                                                                                                     ggTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCgg                                     >  2:3008725/1‑132 (MQ=255)
                                                                                                                                       tttGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAAc                                                                      >  1:2072887/1‑97 (MQ=255)
                                                                                                                                       tttGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAAc                                                                      <  2:2072887/97‑1 (MQ=255)
                                                                                                                                             tGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTGGCCGATTACGc           >  2:4346981/1‑150 (MQ=255)
                                                                                                                                                tgtgATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTGGCCGATTACGCCCt        >  1:2648452/1‑150 (MQ=255)
                                                                                                                                                tgtgATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTGGCCGATTACGCCCt        >  1:2987422/1‑150 (MQ=255)
                                                                                                                                                tgtgATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTGGCCGATTACGCCCt        >  1:763638/1‑150 (MQ=255)
                                                                                                                                                tgtgATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTGGCCGATTACGCCCt        >  1:1505437/1‑150 (MQ=255)
                                                                                                                                                tgtgATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTGGCCGATTACGCCCt        >  1:318441/1‑150 (MQ=255)
                                                                                                                                                tgtgATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGAACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTGGCCGATTACGCCCt        >  2:2565745/1‑150 (MQ=255)
                                                                                                                                                 gtgATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTGGCCGATTACGCCCTc       <  2:4114615/150‑1 (MQ=255)
                                                                                                                                                  ttatatCCCTGGGGGTGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTGGCCGATTACGCCCTCg      <  1:5361063/148‑1 (MQ=255)
                                                                                                                                                    atatCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTGGCCGATTACGCCCTCGGc    >  2:2362086/1‑150 (MQ=255)
                                                                                                                                                     tatCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTGGCCGATTACGCCCTCGGCt   >  1:4327396/1‑150 (MQ=255)
                                                                                                                                                      atCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTa                                                                   >  2:2189371/1‑85 (MQ=255)
                                                                                                                                                      atCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTa                                                                   <  1:2189371/85‑1 (MQ=255)
                                                                                                                                                      atCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGAc                                                                <  2:440674/88‑1 (MQ=255)
                                                                                                                                                      atCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGAc                                                                >  1:440674/1‑88 (MQ=255)
                                                                                                                                                      atCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTGGCCGATTACGCCCTCGGCTa  >  2:1680385/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
GTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTGGCCGATTACGCCCTCGGCTA  >  NC_002947/2669895‑2670193

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: