Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,673,058 T→C 11.3% intergenic (+84/‑83) 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,673,0580TC11.3% 495.9 / 38.9 168intergenic (+84/‑83)prpF/prpDaconitate isomerase/2‑methylcitrate dehydratase
Reads supporting (aligned to +/- strand):  ref base T (89/60);  new base C (7/12);  total (96/72)
Fisher's exact test for biased strand distribution p-value = 8.34e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.45e-01

Only 100 of 174 total aligned reads displayed.

ATCATGAGCCGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTGGCCGA  >  NC_002947/2672909‑2673205
                                                                                                                                                      |                                                                                                                                                   
aTCATGAGCCGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGttt                                                                                                                                                     <  2:5686285/150‑1 (MQ=255)
aTCATGAGCCGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGttt                                                                                                                                                     <  2:7428897/150‑1 (MQ=255)
 tCATGAGCCGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATATCAGGGTTTg                                                                                                                                                    >  1:6709781/1‑150 (MQ=255)
   aTGAGCCGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGaa                                                                                                                                                  >  2:1574752/1‑150 (MQ=255)
   aTGAGCCGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGaa                                                                                                                                                  >  2:5523126/1‑150 (MQ=255)
    tGAGCCGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAAt                                                                                                                                                 <  2:5865699/150‑1 (MQ=255)
       gCCGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCt                                                                                                                                              <  1:3182865/150‑1 (MQ=255)
        ccGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCtg                                                                                                                                             <  1:4742083/150‑1 (MQ=255)
         cGTAGTGCTCGCGTGCTTATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCtgt                                                                                                                                            <  1:6840339/150‑1 (MQ=255)
         cGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCTATAGCTTCTAAAGCACCACTGGCGCTCGGATTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCtgt                                                                                                                                            <  1:3335417/150‑1 (MQ=255)
         cGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCtgt                                                                                                                                            <  1:4324071/150‑1 (MQ=255)
         cGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCtgt                                                                                                                                            <  2:7197840/150‑1 (MQ=255)
          gTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCtgtg                                                                                                                                           >  2:2236564/1‑150 (MQ=255)
                  cgcgTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCt                                                                                                                                   <  2:2181188/150‑1 (MQ=255)
                  cgcgTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCt                                                                                                                                   <  2:273650/150‑1 (MQ=255)
                   gcgTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGGTTGAATGCTGTGATATcccgg                                                                                                                                  >  2:836047/1‑148 (MQ=255)
                     gTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCtgtg                                                                                                                                <  1:932055/150‑1 (MQ=255)
                     gTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCtgtg                                                                                                                                <  1:5495687/150‑1 (MQ=255)
                      tGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTgg                                                                                                                               <  2:6085705/150‑1 (MQ=255)
                      tGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGGATATCACAGCATTCAAACCCTGTGATTACCTGTgg                                                                                                                               <  1:4105743/150‑1 (MQ=18)
                        cTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCTGCTGTGCAGCCCTTCGCGGGCTTGCCCGCTCCCACAGGGATATCACAGCATTCAAACCCTGTGATTACCTGTGGGa                                                                                                                             <  1:3133688/150‑1 (MQ=17)
                        cTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGa                                                                                                                             <  2:3089536/150‑1 (MQ=255)
                        cTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGa                                                                                                                             <  2:624608/150‑1 (MQ=255)
                        cTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGGATATCACAGCATTCAAACCCTGTGATTACCTGTGGGa                                                                                                                             <  1:5223456/150‑1 (MQ=18)
                         tGATGTAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAg                                                                                                                            <  2:3450985/150‑1 (MQ=255)
                         tGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACACGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAg                                                                                                                            >  2:4524010/1‑150 (MQ=255)
                            tGGCGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCgg                                                                                                                         <  1:6920238/150‑1 (MQ=255)
                            tGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCgg                                                                                                                         <  1:6324035/150‑1 (MQ=255)
                            tGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCgg                                                                                                                         <  1:3613429/150‑1 (MQ=255)
                            tGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGGATATCACAGCATTCAAACCCTGTGATTACCTGTGGGAGCgg                                                                                                                         <  1:1491487/150‑1 (MQ=18)
                            tGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGGATATCACAGCATTCAAACCCTGTGATTACCTGTGGGAGCgg                                                                                                                         <  1:5562702/150‑1 (MQ=18)
                             ggAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCggg                                                                                                                        <  2:2199501/150‑1 (MQ=255)
                             ggAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCGGTGGGAGCggg                                                                                                                        >  2:3718444/1‑150 (MQ=255)
                                 ggCTGGGTTCGCGTGCCTGGCGATGGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGGATATCACAGCATTCAAACCCTGTGATTACCTGTGGGa                                                                                                                             <  1:3140677/141‑1 (MQ=17)
                                 ggCTGGGTTCGCGTGCCTGGCGATGGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGGATATCACAGCATTCAAACCCTGTGATTACCTGTGGGa                                                                                                                             >  2:3140677/1‑141 (MQ=17)
                                  gCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGGATATCACAGCATTCAAACCCTGTGATTACCTGTGGGAGCGGGCAAGc                                                                                                                   <  2:337918/150‑1 (MQ=17)
                                       gTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGGATATCACAGCATTCAAACCCTGTGATTACCTGTGGGAGCgg                                                                                                                         >  1:2176746/1‑139 (MQ=17)
                                           gcgTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTTAATTGTGTGATATACCTGTTGGAGCGGGCAAGCCCGCGaatg                                                                                                          >  2:2167496/1‑148 (MQ=38)
                                                   ggCGTTCGCTTCTAACGCACCACGGGCGCGGGGGTGCAGGGGGGGCCGGTTTGCGGCCCTTGGGGGGGTGGCCCGCGCCCCCAGTTAAGCACGGGGTTTGATTGGTGGGTTATCCGGGTGGGCGCGGGCACGCCCGCGAAGGGCTGcaca                                                                                                  <  2:4005136/150‑1 (MQ=255)
                                                       aTAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGGAATCACAGGGTTTGAATGGGGTGATATCCCTGTGGGAGCGGGCCAGCCCGCGGAGGGCTGCAcagcag                                                                                              >  2:6341503/1‑150 (MQ=38)
                                                          gCTTCTAACGCACCACGGGCGCTCGGCTTCACGCGGGGCCGCTGTGCGGCCCTTCGCGGGCTTGCCGGCTCCCCCCGGTAATCCCAGGGTTTGCATGCTGTGATATCCGTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGccc                                                                                           <  2:2399647/150‑1 (MQ=25)
                                                            ttCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCCCAGGGAAACACAGGGTGTGTATGCGGTGATATACCGGTGGGGGCGGGCACGCCCGCGAAGGGCGGCACAGCAGCCCCg                                                                                         >  1:2353933/1‑150 (MQ=37)
                                                             tCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCgg                                                                                        >  2:1226197/1‑150 (MQ=255)
                                                             tCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCCCTCCCAAAGGTAATCAAACGGTTTTAATGCTGTGGTATAACGGTGGGGGGGGGGAAGCCCGCGAAAGGGGGCACACCAGcccccg                                                                                        >  1:1624981/1‑148 (MQ=255)
                                                                 aaGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGGTTGAATGCGGGGATATCCCGGTGGGAGCGGGCCAGCCCGCGCAGGGCTGCACCGCAGCCCCGGccct                                                                                    >  2:5881392/1‑147 (MQ=38)
                                                                  aGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTCATCACAGGGTTTGAATGCTGTGATATCCCGGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCCAtt                                                                                   >  2:2874774/1‑150 (MQ=255)
                                                                   gCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCaa                                                                                                                     >  2:5949939/1‑115 (MQ=255)
                                                                    caccacTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATc                                                                                                                                      <  1:3448998/97‑1 (MQ=255)
                                                                       cacTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCCATTCTGaa                                                                              >  2:1208669/1‑150 (MQ=255)
                                                                               cTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACCGGGGTTGAATGCTGTGATATCCCGGTGGGAGCGGGCAAGCCCGCGTCGGGCTGCACAGCAGCCCCGGGCATTCTGCAGGGAGCTg                                                                      >  2:6579085/1‑150 (MQ=37)
                                                                                tCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCGGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCCAATCGGACGGGAGCTgg                                                                     >  2:3854718/1‑150 (MQ=37)
                                                                                      tcatttggggCCGCTTTGGGGCCCTTCGGGGGCTTGCCCGCTCCCACAGTTAATCACTGGGTTTTATTGCTTTTTTATCCGTGTGGTAGGGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGa                                                               <  2:225121/144‑1 (MQ=255)
                                                                                            ggggCCGCTTGGCGGCCCTTCGCGGGCTGGCCCGCTCCCACAGGTAATCACAGGTTTTGAATGCTGTGATATCCCTGTGGGTGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCa                                                         <  2:3992680/150‑1 (MQ=255)
                                                                                               gCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCGGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCCATTCTGAAGGGAGCTGGATATGAGCGCCAACg                                                      >  2:5419976/1‑150 (MQ=255)
                                                                                                ccGCTTTGCGGCCCTTCGCGGGCTTGCCCCCTCCCACAGGTAAACAAAGGGTTTTAAAGCTGTGATATACCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACACCAGCCCCCGCCATTCTGAAGGGAGCTGGATATGAGCGCCAAcct                                                     >  2:747874/1‑148 (MQ=37)
                                                                                                   cTTTGCGGCCCTTCGCGGGCTTGCCCGCCCCCCCCGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGa                                                  <  2:6885203/150‑1 (MQ=255)
                                                                                                           cccTTCGCGGGCGTGCCCCCTCCCACACGTAATCACAGGGTTTGAATGCTGTGATATACCGGTGGGACCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCCGCAATTGTGAAGGGTCCTGGATATGAGCGCCAACGTAGACCTCAACg                                          >  2:1740055/1‑150 (MQ=37)
                                                                                                             cTTCGCGGGCGTGCCCGCTCCCACAGGTAAACACAGGGTTTTAATTCTGTGATATACCGGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGCAGGGAGCTGGATATGAGCGCCCACGGAGACCGCAACGAt                                        >  2:4681092/1‑150 (MQ=37)
                                                                                                                   gggCTTGCCCGCTCCCCCCGTTACACACCGGGTTTGATTGCTTTGATAACCCTGGGGGAGGGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCg                                  <  1:1226197/150‑1 (MQ=255)
                                                                                                                   gggCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGCATGCTGTGATATCCCGGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACCGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCg                                  >  2:5388682/1‑150 (MQ=255)
                                                                                                                     gCTGGCCCGCTCCCACTGTTAATCACCGGGTTTGAATGCTGTGATAGCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGa                                <  1:4926923/150‑1 (MQ=39)
                                                                                                                       ttGCCCGCTCCCACAGGTAATCACAGTGTTTGAATGCTGTGATATACCTGTGGGAGCGGGCAAGCCAGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAATGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCTCCCGGAtt                              >  1:206569/1‑150 (MQ=38)
                                                                                                                       ttGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGAtt                              >  2:5488442/1‑150 (MQ=255)
                                                                                                                         gCCCGCTCCCACAGGGATATCACAGCATTCAAACCCTGTGATTACCTGTGGGAGCGGGCAAGCCCGCGAAGGGCCGCAAAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTAc                            >  1:7386490/1‑150 (MQ=17)
                                                                                                                          cccGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACg                           >  1:3732418/1‑150 (MQ=255)
                                                                                                                           ccGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGa                          >  2:3497358/1‑150 (MQ=255)
                                                                                                                            cGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGAc                         >  2:2161619/1‑150 (MQ=255)
                                                                                                                             gCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGAcc                        >  1:5096217/1‑150 (MQ=255)
                                                                                                                               tCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAgg                                                                                                          >  1:1731728/1‑66 (MQ=255)
                                                                                                                               tCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAgg                                                                                                          <  2:1731565/66‑1 (MQ=255)
                                                                                                                               tCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAgg                                                                                                          <  2:1731728/66‑1 (MQ=255)
                                                                                                                               tCCCACAGGGATATCACAGCATTCAAACCCTGTGATTACCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCgg                      >  2:7181893/1‑150 (MQ=18)
                                                                                                                                cccACTGTTAATCACAGGTTTTGATTGCTGTGATAGCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCggg                     <  2:1486388/150‑1 (MQ=39)
                                                                                                                                cccACTGTTAAACACCGGGTTTGATTGGTGTTATATCCCTGTGGGCGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCggg                     <  1:5861687/150‑1 (MQ=255)
                                                                                                                                  cccaGTTATTCACTGGTTTTGTATGCTGGGATTGCGCTGTGGGAGCGGGCACGCCCGCGAAGGGCTGCACAGCAGCCCCGGCACTTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTg                   <  2:4359785/148‑1 (MQ=37)
                                                                                                                                  cacaGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGAc                                                 <  2:5629098/120‑1 (MQ=255)
                                                                                                                                  cacaGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTg                   >  1:1164348/1‑150 (MQ=255)
                                                                                                                                  cacaGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTg                   >  2:4626843/1‑150 (MQ=255)
                                                                                                                                  cacaGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTg                   >  2:4626651/1‑150 (MQ=255)
                                                                                                                                  cacaGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTg                   >  1:1164262/1‑150 (MQ=255)
                                                                                                                                   actggTAATAACAGGTTTTTATTGCTGTTATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGc                  <  1:4864967/147‑1 (MQ=39)
                                                                                                                                   acaGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCgg                                                                                        <  2:5166449/80‑1 (MQ=255)
                                                                                                                                   acaGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGc                  >  1:1939523/1‑150 (MQ=255)
                                                                                                                                     aGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTa                >  2:1939867/1‑150 (MQ=255)
                                                                                                                                      ggTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTAc               >  1:3936858/1‑150 (MQ=255)
                                                                                                                                        tAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCggg                     <  2:3436813/142‑1 (MQ=255)
                                                                                                                                           aTCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAAc           >  1:4674471/1‑150 (MQ=255)
                                                                                                                                            tCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAAcc          >  2:5126928/1‑150 (MQ=255)
                                                                                                                                            tCACAGGGTATGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAAATCTGAAGGGAGCTGGATAATAGAGACAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAAcc          >  1:4872657/1‑150 (MQ=39)
                                                                                                                                             cacaGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAAccc         >  1:386643/1‑150 (MQ=255)
                                                                                                                                              acaGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCt        >  1:5440449/1‑150 (MQ=255)
                                                                                                                                              acaGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCt        >  1:6913904/1‑150 (MQ=255)
                                                                                                                                              acaGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCt        >  1:2682341/1‑150 (MQ=255)
                                                                                                                                               caGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTg       >  2:6786712/1‑150 (MQ=255)
                                                                                                                                               caGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTg       >  2:4159846/1‑150 (MQ=255)
                                                                                                                                                aGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCg                                  >  1:4000986/1‑122 (MQ=255)
                                                                                                                                                aGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTgg      >  2:650992/1‑150 (MQ=255)
                                                                                                                                                aGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTgg      >  2:86862/1‑150 (MQ=255)
                                                                                                                                                  ggTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTGGcc    >  2:3606506/1‑150 (MQ=255)
                                                                                                                                                    tttGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTGGCCGa  >  1:6098788/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                   
ATCATGAGCCGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTGGCCGA  >  NC_002947/2672909‑2673205

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: