Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 349,717 T→C 5.8% intergenic (+145/‑295) prpB → / → prpC 2‑methylisocitrate lyase/2‑methylcitrate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913349,7170TC 5.8% 2024.6 / 69.7 697intergenic (+145/‑295)prpB/prpC2‑methylisocitrate lyase/2‑methylcitrate synthase
Reads supporting (aligned to +/- strand):  ref base T (172/440);  new base C (28/10);  total (242/455)
Fisher's exact test for biased strand distribution p-value = 2.40e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 716 total aligned reads displayed.

AATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACT‑GAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGT  >  NC_000913/349568‑349866
                                                                                                                                                     |                                                                                                                                                      
aaTAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCAGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACAGTAGGTCGGATAAGACGGATGGCGTCGCATCCCACAACCGATGCCTGATGCGCCGCTGACGt                                                                                                                                                        >  1:10502360/1‑150 (MQ=255)
   aaaaaaCGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGAc                                                                                                                                                     <  1:4104834/150‑1 (MQ=255)
        aCGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATc                                                                                                                                                <  2:4934481/150‑1 (MQ=255)
         cGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCa                                                                                                                                               <  2:7151621/150‑1 (MQ=255)
          gCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAg                                                                                                                                              <  1:3157214/150‑1 (MQ=255)
           cccGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAgg                                                                                                                                             <  2:4518931/150‑1 (MQ=255)
                 gATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACg                                                                                                                                       <  2:6150897/150‑1 (MQ=255)
                   ttGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACggg                                                                                                                                     <  1:10917444/150‑1 (MQ=255)
                   ttGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACggg                                                                                                                                     <  2:9574665/150‑1 (MQ=255)
                    tGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGCATGGCGTCGCATCCGACAACCGAAGCCTCAAGCGCCGCTGACGCCACTTATCAGGCCAACtggg                                                                                                                                    >  2:8002512/1‑150 (MQ=14)
                     gTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGt                                                                                                                                   <  2:9807824/150‑1 (MQ=255)
                     gTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGt                                                                                                                                   <  2:4095074/150‑1 (MQ=255)
                      tATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTg                                                                                                                                  <  2:4622342/150‑1 (MQ=255)
                       aTTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGc                                                                                                                                 <  1:3779144/150‑1 (MQ=255)
                        ttCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCa                                                                                                                                <  2:354645/150‑1 (MQ=255)
                               gCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACT‑GAACt                                                                                                                        <  1:9039603/150‑1 (MQ=255)
                                ccGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACT‑GAACTg                                                                                                                       >  2:9006732/1‑150 (MQ=255)
                                      ccTGATGCGTCGCTGACGCGAGTTATCTGGCCTACGATGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACT‑GAACTGTAGGTc                                                                                                                 <  2:7439727/150‑1 (MQ=255)
                                       cTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACT‑GAACTGTAGGTCg                                                                                                                <  2:2398462/150‑1 (MQ=255)
                                       cTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGCATGGCGTCGCATCCGACAAACGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACT‑GTACTGTAGGTCg                                                                                                                >  2:10444303/1‑150 (MQ=255)
                                         gATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACT‑GAACTGTAGGTCGGa                                                                                                              <  1:4642892/150‑1 (MQ=255)
                                            gCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACT‑GAACTGTAGGTCGGATaa                                                                                                           <  1:6038200/150‑1 (MQ=255)
                                            gCGTCGCTGACGCGACTTATCAGGCCAACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACT‑GAACTGTAGGTCGGATaa                                                                                                           >  1:5791133/1‑150 (MQ=255)
                                             cGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACT‑GAACTGTAGGTCGGATAAg                                                                                                          <  2:523817/150‑1 (MQ=255)
                                             cGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACT‑GAACTGTAGGTCGGATAAg                                                                                                          <  1:5735609/150‑1 (MQ=255)
                                              gTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACT‑GAACTGTAGGTCGGATAAGa                                                                                                         <  1:296500/150‑1 (MQ=255)
                                               tCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACT‑GAACTGTAGGTCGGATAAGAc                                                                                                        >  2:2586902/1‑150 (MQ=255)
                                                      cgcgACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACT‑GAACTGTAGGTCGGATAAGACGCATAGc                                                                                                 <  1:1754786/150‑1 (MQ=255)
                                                        cgACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACT‑GAACTGTAGGTCg                                                                                                                >  1:5148128/1‑133 (MQ=40)
                                                         gACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACT‑GAACTGTAGGTCGGATAAGACGCATAGCGTc                                                                                              <  1:5975865/150‑1 (MQ=39)
                                                         gACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACT‑GAACTGTAGGTCGGATAAGACGCATAGCGTc                                                                                              >  1:6466448/1‑150 (MQ=39)
                                                         gACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACT‑GAACTGTAGGTCGGATAAGACGCATAGCGTc                                                                                              <  2:4803471/150‑1 (MQ=39)
                                                         gACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACT‑GAACTGTAGGTCGGATAAGACGCATAGCGTc                                                                                              <  1:7420179/150‑1 (MQ=39)
                                                          aCTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACT‑GAACTGTAGGTCGGATAAGACGCATAGCGTCg                                                                                             <  2:11071791/150‑1 (MQ=39)
                                                            ttATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGAGACTTATCAGGCCTACGGGGTGAACT‑Gaa                                                                                                                          >  1:8683901/1‑119 (MQ=255)
                                                              aTCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACT‑GAACTGTAGGTCGGATAAGACGCATAGCGTCGCATc                                                                                         <  2:3722508/150‑1 (MQ=39)
                                                               tCAGGCCTACGAGGTGCACTGAACTCTAGGTCGGATAAGACGGATGGCGACGCATCCGACAACCGATGCCTGATGCGCCGCAGACGTGACTTATCAGGCCTACGGGGTGCACT‑GAACTGTAGGTCGGATAAGACGCATAGCGTCGCATcc                                                                                        >  1:10680031/1‑150 (MQ=255)
                                                                 aGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACT‑GAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGa                                                                                      <  2:1470584/150‑1 (MQ=39)
                                                                    ccTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACT‑GAACTGTAg                                                                                                                    >  2:5502438/1‑117 (MQ=25)
                                                                     cTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGAACT‑GAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGa                                                                                      >  1:5419259/1‑146 (MQ=37)
                                                                     cTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACT‑GAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAAc                                                                                  >  2:7061987/1‑150 (MQ=21)
                                                                       aCGTGGTGCACTGAACTGTTGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACT‑GAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCg                                                                                <  1:7845207/150‑1 (MQ=255)
                                                                       aCGAGGTGCACTTAAGTGTTGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGGTGACGTGACTTATCCGGCCTACGGGGTGAACT‑GAACTGTAGGTCGGATGCGACGCTATGCGTCTTATCCGAcctaca                                                                                <  2:8616332/150‑6 (MQ=11)
                                                                        cGAGGTGCACTGAAGTGTAGGTCGGATAAGACGGTTGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACT‑GAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGa                                                                               <  1:4934544/150‑1 (MQ=255)
                                                                          aGGTGCACACAACACTAGCACGCATAACACCCATGGCCTCGCATACGACAACCGATGCCTGAAGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACT‑GAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATg                                                                             <  2:3941897/150‑1 (MQ=11)
                                                                            gTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACT‑GAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGcc                                                                           >  1:1032007/1‑150 (MQ=37)
                                                                                aCTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACT‑GAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGAt                                                                       >  1:10026125/1‑150 (MQ=37)
                                                                                aCTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACT‑GAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGAt                                                                       <  2:1755964/150‑1 (MQ=37)
                                                                                 cTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACT‑GAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATg                                                                      <  2:1641668/150‑1 (MQ=37)
                                                                                 cTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACAGAATACAGATGCACAGCGGACGTGACTTAAAGAGCCTACGGGGCGCACTAGAGCT‑AAGGTCGCGGAGGACGAATAGCGTCGCATCCGGCCTCCGGTGCCTgggg                                                                      >  2:6553469/1‑147 (MQ=2)
                                                                                  tGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGAGCAGCTGAAGTGACTTATCAGGCCTCCGGGGTGAACG‑GACCTGGATGACGGATAAGACGCATAGCGTCGCAGCCGACAACCGATGCCTGATgc                                                                     >  1:7058720/1‑150 (MQ=11)
                                                                                    aaCTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACT‑GAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATgcgc                                                                   <  2:8202173/150‑1 (MQ=37)
                                                                                      cTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACT‑GAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGcgccg                                                                 <  1:1139371/150‑1 (MQ=37)
                                                                                       tGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACT‑GAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGcgccgc                                                                <  1:8264066/150‑1 (MQ=37)
                                                                                          aGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACT‑GAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGa                                                             <  2:5748431/150‑1 (MQ=37)
                                                                                          aGGTCGGAAAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACT‑GAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGa                                                             <  1:6889348/150‑1 (MQ=37)
                                                                                           ggTCGGATAAGACGGATGGCGTCGCATCCGACACCCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACT‑GAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGAc                                                            <  2:210931/150‑1 (MQ=37)
                                                                                           ggTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACT‑GAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGAc                                                            >  2:533492/1‑150 (MQ=37)
                                                                                             tCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACT‑GAACTGTAGGTCGGATAAGACGCATAGCGTc                                                                                              >  2:3345286/1‑114 (MQ=39)
                                                                                              cGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACT‑GAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGAcgcg                                                         >  2:9450394/1‑150 (MQ=37)
                                                                                               tgATAATACGGATGTCGTCGCATCCGACATCCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGTGGTGCACT‑GAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGa                                                        <  2:2918019/149‑1 (MQ=255)
                                                                                                  tAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACT‑GAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACtt                                                     >  1:2217505/1‑150 (MQ=37)
                                                                                                    aGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACT‑GAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTAt                                                   <  1:4520912/150‑1 (MQ=17)
                                                                                                     gACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACT‑GAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATc                                                  <  2:10808653/150‑1 (MQ=37)
                                                                                                      aCGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACT‑GAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCa                                                 <  2:2027089/150‑1 (MQ=37)
                                                                                                      aCGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACT‑GAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCa                                                 >  2:2808172/1‑150 (MQ=37)
                                                                                                      aCGGAGGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACT‑GAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCa                                                 <  2:3175547/150‑1 (MQ=25)
                                                                                                         gATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACT‑GAACTGTAGGTCGGATAAGACGCATAGCGTCGCCTCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGc                                              <  2:1209070/150‑1 (MQ=25)
                                                                                                               gTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACT‑GAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACgg                                        <  1:3962627/150‑1 (MQ=32)
                                                                                                                tCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACT‑GAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACggg                                       <  1:1610895/150‑1 (MQ=32)
                                                                                                                 cGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACT‑GAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACgggg                                      <  1:1905303/150‑1 (MQ=32)
                                                                                                                  gCATCCGACACCCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACT‑GAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGt                                     <  2:10086003/150‑1 (MQ=17)
                                                                                                                  gCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACT‑GAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGt                                     <  1:2622135/150‑1 (MQ=32)
                                                                                                                    aTCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACT‑GAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTCTCGGGCCTACGGGGTGc                                   >  1:8026270/1‑150 (MQ=17)
                                                                                                                     tCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACT‑GAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGcc                                  >  1:7535269/1‑149 (MQ=17)
                                                                                                                      ccGCCAACCGATGCCTGATGCGCCGCTGCCGTGACTTATCAGGCCTACGGGGTGCACT‑GAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCAc                                 <  1:1422051/150‑1 (MQ=21)
                                                                                                                      ccGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACT‑GAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCAc                                 <  2:5963745/150‑1 (MQ=34)
                                                                                                                       cGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACT‑GAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACt                                <  1:2843663/150‑1 (MQ=34)
                                                                                                                       cGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACT‑GAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACt                                <  1:9609878/150‑1 (MQ=34)
                                                                                                                          cAACCGATGCCTGATGCGCCGCTGACTTGTCTGATCCTGACTAGGGGGTGCACT‑GCACTGTAGGTATTATAAGACGCAGTGGGTGGCATCCGACACCCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGaa                             <  2:9321009/150‑1 (MQ=11)
                                                                                                                          cAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACT‑GAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGaa                             <  1:6334483/150‑1 (MQ=34)
                                                                                                                           aaCCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACT‑GAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAAc                            >  1:9415029/1‑150 (MQ=34)
                                                                                                                            aCCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACT‑GAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACt                           <  2:6457048/150‑1 (MQ=34)
                                                                                                                             ccGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACT‑GAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTg                          <  1:6558192/150‑1 (MQ=34)
                                                                                                                             ccGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACT‑GAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTg                          >  2:6895052/1‑150 (MQ=34)
                                                                                                                              cgtggCCTGCTGCGCGGCGGCCGGGACTTATCAGGCCTACGGGGTGCACT‑GAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGt                         <  2:5862392/146‑1 (MQ=11)
                                                                                                                              cGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACT‑GAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTga                         >  2:4397598/1‑149 (MQ=21)
                                                                                                                              cGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACT‑GAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGt                         <  2:6341692/150‑1 (MQ=34)
                                                                                                                               gATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACT‑GAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTa                        <  1:2421272/150‑1 (MQ=34)
                                                                                                                                aTGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACT‑GAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAg                       <  1:8540709/150‑1 (MQ=34)
                                                                                                                                aTGCCCGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACT‑GAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAg                       <  2:10614884/150‑1 (MQ=21)
                                                                                                                                 tGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACT‑GAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACACGACTTATCAGGCCTACGGGGTGCACTGAACTGTAgg                      >  1:9167805/1‑150 (MQ=21)
                                                                                                                                      gATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACT‑GAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGa                 <  1:4419861/150‑1 (MQ=34)
                                                                                                                                       aTGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACT‑GAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGAt                <  2:1089681/150‑1 (MQ=34)
                                                                                                                                        tGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACT‑GAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATa               <  1:3418011/150‑1 (MQ=34)
                                                                                                                                          cgccgcTGACGTGACTTATCAGGCCTACGGGGTGCACT‑GAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAg             <  2:8911099/150‑1 (MQ=34)
                                                                                                                                           gccgcTGACGTGACTTATCAGGCCTAAGGGGTGCACT‑GACCTGTAGGACGGATAAGACGCATAGCGACGCATCCGACACCCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTCAACTGTAGGTCGGATAAGa            <  2:8452675/150‑1 (MQ=11)
                                                                                                                                               cTGACGTGACTTATCAGGCCTACGGGGTGCACT‑GAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCa        <  1:9941019/150‑1 (MQ=34)
                                                                                                                                                     tGACTTATCAGGCCTACGGGGTGCACT‑GAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGt  >  2:7728211/1‑150 (MQ=32)
                                                                                                                                                     tGACTTATCAGGCCTACGGGGTGCACT‑GAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGt  >  1:621638/1‑150 (MQ=32)
                                                                                                                                                     |                                                                                                                                                      
AATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACT‑GAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGT  >  NC_000913/349568‑349866

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: ‑