Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,349,055 G→T 7.5% P13T (CCT→ACT)  PP_1168 ← TRAP dicarboxylate transporter subunit DctQ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,349,0550GT 7.5% 954.1 / 46.1 286P13T (CCT→ACT) PP_1168TRAP dicarboxylate transporter subunit DctQ
Reads supporting (aligned to +/- strand):  ref base G (134/130);  new base T (10/12);  total (144/142)
Fisher's exact test for biased strand distribution p-value = 6.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 309 total aligned reads displayed.

GGTGCCCAGCACGTAGCGGGCGAAGATGCCCCACGGGATGATCAGCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTAGAACGGCTTCTTGTCGACGGTAA  >  NC_002947/1348906‑1349204
                                                                                                                                                     |                                                                                                                                                     
ggTGCCCAGCACGTAGCGGGCGAAGATGCCCCACGGGATGATCAGCGTCATGGCCCGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAgg                                                                                                                                                       >  2:2893180/1‑150 (MQ=255)
ggTGCCCAGCACGTAGCGGGCGAAGATGCCCCACGGGATGATCAGCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAgg                                                                                                                                                       >  2:1926601/1‑150 (MQ=255)
ggTGCCCAGCACGTAGCGGGCGAAGATGCCCCACGGGATGATCAGCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAgg                                                                                                                                                       <  2:1789873/150‑1 (MQ=255)
     ccAGCACGTAGCGGGCGAAGATGCCCCACGGGATGATCAGCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGggcggc                                                                                                                                                  <  1:3849467/150‑1 (MQ=255)
      cAGCACGTAGCGGGCGAAGATGCCCCACGGGATGATCAGCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGgcggcg                                                                                                                                                 >  1:1946041/1‑150 (MQ=255)
      cAGCACGTAGCGGGCGAAGATGCCCCACGGGATGATCAGCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGgcggcg                                                                                                                                                 >  1:2917260/1‑150 (MQ=255)
        gCACGTAGCGGGCGAAGATGCCCCACGGGATGATCAGCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCggcggg                                                                                                                                               >  2:7281114/1‑150 (MQ=255)
        gCACGTAGCGGGCGAAGATGCCCCACGGGATGATCAGCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCggcggg                                                                                                                                               >  1:4497602/1‑150 (MQ=255)
        gCACGTAGCGGGCGAAGATGCCCCACGGGATGATCAGCGTCATGGCCAGGATCGACAAGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCggcggg                                                                                                                                               <  2:2760907/150‑1 (MQ=255)
          aCGTAGCGGGCGAAGATGCCCCACGGGATGATCAGCGTCATGGCCAGGATCGACGGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCggcgggcg                                                                                                                                             <  1:2521874/150‑1 (MQ=255)
          aCGTAGCGGGCGAAGATGCCCCACGGGATGATCAGCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCggcgggcg                                                                                                                                             <  2:5125406/150‑1 (MQ=255)
            gTAGCGGGCGAAGATGCCCCACGGGATGATCAGCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGcgcc                                                                                                                                           <  2:4132948/150‑1 (MQ=255)
               gCGGGCGAAGATGCCCCACGGGATGATCAGCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCa                                                                                                                                        <  1:5635048/150‑1 (MQ=255)
                 gggCGAAGATGCCCCACGGGATGATCAGCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGTGCGGCGCCCGCCGCCCt                                                                                                                                      <  2:3053875/150‑1 (MQ=255)
                 gggCGAAGATGCCCCACGGGATGATCAGCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGTGCGGCGCCCGCCGCCCt                                                                                                                                      <  2:6853415/150‑1 (MQ=255)
                   gCGAAGATGCCCCACGGGATGATCAGCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGTGCGGCGCCCGCCGCCCTGc                                                                                                                                    <  2:6408901/150‑1 (MQ=255)
                      aaGATGCCCCACTGTATGATCAGCGTCATGGCCAGGATCGACAGGCCCGCGATCCTGATGCAACTGCGGTACAGCGTGTCGTTCTCGCTCAGGAAAAGCGAGTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGa                                                                                                                                 <  1:6085915/150‑1 (MQ=255)
                         aTGCCCCACGGGATGATCAGCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGt                                                                                                                              >  1:1102698/1‑150 (MQ=255)
                            ccccACGGGATGATCAGCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGa                                                                                                                           >  2:6456086/1‑150 (MQ=255)
                             cccACGGGATGATCAGCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGTGCGGCGCCCGCCGCCCTGc                                                                                                                                    <  2:2583538/140‑1 (MQ=255)
                               cACGGGATGATCAGCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAACGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTg                                                                                                                        <  2:3530483/150‑1 (MQ=255)
                               cACGGGATGATCAGCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCCTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTg                                                                                                                        >  2:1729249/1‑150 (MQ=255)
                                aCGGTATGATCAGCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTTTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTgg                                                                                                                       <  1:7195519/150‑1 (MQ=255)
                                aCGGGATGATCAGCGTCATGGCCAGGATTACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAATGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTggg                                                                                                                      <  1:2966993/150‑1 (MQ=255)
                                 cGGGATGATCAGCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTggg                                                                                                                      >  1:4714301/1‑150 (MQ=255)
                                  gggATGATCAGCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGc                                                                                                                     >  1:4594510/1‑150 (MQ=255)
                                    ggtgatCAGCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCaa                                                                                                                   >  2:6630986/3‑150 (MQ=255)
                                    gatgatCAGCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCTATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCaa                                                                                                                   <  2:1871866/150‑1 (MQ=255)
                                         tCAGCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGggcg                                                                                                              >  1:2863402/1‑150 (MQ=255)
                                           aGCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGAGGCAACTGCGGTACAGCGTGTCGTTCCCGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGggcggc                                                                                                            <  2:1978914/150‑1 (MQ=255)
                                            gcttcttGGCCAGGATCGACAGGCCCGCTATCCAGAGGCAAGTGCGGTACCGCGTGTCGTTCACGCTCAGGAAAAGCGATTCCATAGGCATCACCAACGCGGCCGGGCGGGGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCt                                                                                                           <  1:3986973/144‑1 (MQ=255)
                                            gcgcCATGGCCAGGATCGACAGGCCCGCGCTCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCt                                                                                                           <  1:4128262/150‑1 (MQ=255)
                                            gCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCt                                                                                                           >  1:116148/1‑150 (MQ=255)
                                              gtactGGCCAGGATCGACCGGCCCGCGATCCAGATGCAACGGCGTTACCGCGGGTCTTTCACGCCCCGGAAAAGCGATTACATCGGCAACACCAACGCGGCCGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGt                                                                                                         <  1:956770/146‑1 (MQ=255)
                                               tCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGtt                                                                                                        >  1:5558276/1‑150 (MQ=255)
                                                cATGGCCAGGATCGACAGGCCCGCGATCCTGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTa                                                                                                       <  1:6696935/150‑1 (MQ=255)
                                                   ggCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGACAAGCGATGTCAGAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGAGTGGGCACGGGCGGCGGGTAcgg                                                                                                    >  2:2068879/1‑148 (MQ=255)
                                                    gCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAATGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTgg                                                                                                   >  1:528363/1‑150 (MQ=255)
                                                         gATCGACAGGCCCGCGATCAAGATGCACCTGCGGTACAGCGTGTCGTTCACGCTCCGGAAAAGCGTTTTCATAGGCAACACAAACGCGGCTGGGGGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGc                                                                                              <  1:3351069/150‑1 (MQ=255)
                                                          aTCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCg                                                                                             <  2:7213131/150‑1 (MQ=255)
                                                            cGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTc                                                                                           <  2:7268368/150‑1 (MQ=255)
                                                                   cccGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGc                                                                                    >  1:2700994/1‑150 (MQ=255)
                                                                   cccGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGc                                                                                    >  1:6449594/1‑150 (MQ=255)
                                                                     cgcgATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCtt                                                                                  <  1:1696042/150‑1 (MQ=255)
                                                                           ccAGATGCAACGGCGGTACCGCGTGTCGTTCCCGCTCAGCAAAAGCGATTTCATAGGCATCACCAACGGGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAg                                                                            <  1:5764252/150‑1 (MQ=255)
                                                                            cAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAgg                                                                           >  2:4997143/1‑150 (MQ=255)
                                                                            cAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCAGAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAgg                                                                           >  2:3541011/1‑150 (MQ=255)
                                                                             aGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGt                                                                          >  2:4573871/1‑150 (MQ=255)
                                                                             aGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACAGCGTCGATACGCTTCATCAGGt                                                                          >  1:4971235/1‑150 (MQ=255)
                                                                              gATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTc                                                                         >  2:3199237/1‑150 (MQ=255)
                                                                               aTGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCg                                                                        <  2:2215933/150‑1 (MQ=255)
                                                                               aTGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCg                                                                        >  2:543025/1‑150 (MQ=255)
                                                                                         ggTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTgcg                                                              >  2:6526634/1‑150 (MQ=255)
                                                                                          gTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTgcgc                                                             <  1:7310852/150‑1 (MQ=255)
                                                                                           tACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTgcgcg                                                            >  1:1354472/1‑150 (MQ=255)
                                                                                           tACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTgcgcg                                                            >  2:6706189/1‑150 (MQ=255)
                                                                                              aGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCg                                                         >  2:1396664/1‑150 (MQ=255)
                                                                                                    tCGTTCACGCTCAGGAACAGCGATTTCATAGGCATCACCACCGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGAGTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTAtttt                                                   <  2:7086081/150‑1 (MQ=255)
                                                                                                     cgtccccGCTCCGGAAACGCGTTTCCATAGGCACCACCAACGCGGCAGGGCGGCGGGCGCGGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTAttttt                                                  <  1:1995720/144‑1 (MQ=255)
                                                                                                     cGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCTGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTAttttt                                                  <  2:3758206/150‑1 (MQ=255)
                                                                                                      gTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTc                                                 >  2:2301855/1‑150 (MQ=255)
                                                                                                           cGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTa                                            <  2:6751634/150‑1 (MQ=255)
                                                                                                             cTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTAcc                                          <  2:7348543/150‑1 (MQ=255)
                                                                                                              tCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATa                                                                                       >  1:6837592/1‑104 (MQ=255)
                                                                                                              tCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCg                                         >  1:1552642/1‑150 (MQ=255)
                                                                                                                aGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGaa                                       >  1:1743565/1‑150 (MQ=255)
                                                                                                                     aaGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCTAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGt                                  <  2:2377752/150‑1 (MQ=255)
                                                                                                                     aaGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGt                                  <  2:1683543/150‑1 (MQ=255)
                                                                                                                     aaGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGAGTGGGCAAGGGCGGCTGTTCCTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGt                                  >  1:811603/1‑150 (MQ=255)
                                                                                                                        cGATTTCCTAGGCATCACCAACGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCCTCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGc                               <  2:3833570/150‑1 (MQ=255)
                                                                                                                             tCATAGGCATCTCCAACGCGGCAGGGCGGCGGGCGCCGCACTGCCGATGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGACGTATTTTTCGCGTACCTAAGCCGTGGcgtcgt                          <  2:2673078/150‑1 (MQ=255)
                                                                                                                              cATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTa                         >  2:5927537/1‑150 (MQ=255)
                                                                                                                                tAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTAGa                       >  1:4892426/1‑150 (MQ=255)
                                                                                                                                tAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTAGa                       >  2:4306693/1‑150 (MQ=255)
                                                                                                                                 aGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTAGaa                      <  2:5227857/150‑1 (MQ=255)
                                                                                                                                  ggCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTAGAAc                     <  2:4714301/150‑1 (MQ=255)
                                                                                                                                    cATCACCAACGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTAGAACgg                   <  1:6924101/150‑1 (MQ=255)
                                                                                                                                     aTCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTAGAACGGc                  >  1:2345046/1‑150 (MQ=255)
                                                                                                                                     aTCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTAGAACGGc                  >  1:4055929/1‑150 (MQ=255)
                                                                                                                                      tCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTAGAACGGct                 >  1:1026694/1‑150 (MQ=255)
                                                                                                                                      tCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTAGAACGGct                 <  2:1235157/150‑1 (MQ=255)
                                                                                                                                      tCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTAGAACGGct                 >  2:4210764/1‑150 (MQ=255)
                                                                                                                                        aCCAACGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTAGAACGGcttc               <  2:5843678/150‑1 (MQ=255)
                                                                                                                                         ccAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTAGAACGGcttct              >  1:6454899/1‑150 (MQ=255)
                                                                                                                                         ccAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTAGAACGGcttct              >  2:5751880/1‑150 (MQ=255)
                                                                                                                                           aaaGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTGTTCGCGTACCGAAGCCGTGGCGTCGTAGAACGGCTTCTTg            >  2:3088099/1‑150 (MQ=255)
                                                                                                                                            aaGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTAGAACGGCTTCTTGt           <  1:3653150/150‑1 (MQ=255)
                                                                                                                                            aaGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTAGAACGGCTTCTTGt           <  2:4497602/150‑1 (MQ=255)
                                                                                                                                            aaGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTAGAACGGCTTCTTGt           <  2:6171250/150‑1 (MQ=255)
                                                                                                                                            aaGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTAGAACGGCTTCTTGt           <  2:5078276/150‑1 (MQ=255)
                                                                                                                                             aGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTAGAACGGCTTCTTGTc          <  2:2917260/150‑1 (MQ=255)
                                                                                                                                              tcGGCAGTGCGGCGCCCGCCGCCCTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTAGAACGGCTTCTTGTCg         >  1:3495836/2‑150 (MQ=255)
                                                                                                                                              tcGGCAGTGCGGCGCCCGCCGCCCTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTAGAACGGCTTCTTGTCg         >  2:3970662/2‑150 (MQ=255)
                                                                                                                                              tcGGCAGTGCGGCGCCCGCCGCCCTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTAGAACGGCTTCTTGTCg         >  2:109457/2‑150 (MQ=255)
                                                                                                                                                ggCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTAGAACGGCTTCTTGTCGAc       >  2:4377753/1‑150 (MQ=255)
                                                                                                                                                 gCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTAGAACGGCTTCTTGTCGACg      <  1:3733658/150‑1 (MQ=255)
                                                                                                                                                 gCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTAGAACGGCTTCTTGTCGACg      >  1:6529088/1‑150 (MQ=255)
                                                                                                                                                    gggCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTAGAACGGCTTCTTGTCGACGGTa   <  2:5558276/150‑1 (MQ=255)
                                                                                                                                                     ggcggcGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTAGAACGGCTTCTTGTCGACGGTaa  <  1:5092277/150‑1 (MQ=255)
                                                                                                                                                     ggcggcGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGCGGCGTCGTAGAACGGCTTCTTGTCGACGGTaa  <  2:4594510/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGTGCCCAGCACGTAGCGGGCGAAGATGCCCCACGGGATGATCAGCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTAGAACGGCTTCTTGTCGACGGTAA  >  NC_002947/1348906‑1349204

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: