Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,080,449 A→G 6.0% I396T (ATC→ACC)  tktA ← transketolase 1, thiamine triphosphate‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,080,4490AG 6.0% 182.5 / 2.5 67I396T (ATC→ACC) tktAtransketolase 1, thiamine triphosphate‑binding
Reads supporting (aligned to +/- strand):  ref base A (42/21);  new base G (2/2);  total (44/23)
Fisher's exact test for biased strand distribution p-value = 6.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.58e-01

GGCGTTACGTGCGTATTCCACGAACATCAGGAAGGTGGAGGTGTACGGCAGGAAGCCACCGTGCAGGGAGATACCGTTAGCAATCGCGGTCATACCGAACTCGCGAACACCGTAGTGGATGTAGTTACCCGCAGCATCTTCGTTGATTGCTTTAGAACCAGACCACAGGGTCAGGTTAGACGGCGCCAGGTCAGCAGAACCGCCGAGGAATTCCGGCAACAGCGGACCGAACGCTTCGATAGCATTCTGAGACGCTTTACGGCTGGCGATTTTCGCCGGATTAGCCTGCA  >  NC_000913/3080304‑3080593
                                                                                                                                                 |                                                                                                                                                
ggCGTTACGTGCGTATTCCACGAACATCAGGAAGGTGGAGGTGTACGGCAGGAAGCCACCGTGCAGGGAGATACCGTTAGCAATCGCGGTCATACCGAACTCGCGAACACCGTAGTGGATGTAGTTACCCGCAGCATCTTCGTTGATTg                                                                                                                                               >  2:316724/1‑149 (MQ=255)
     tACGTGCGTATTCCACGAACATCAGGAAGGTGGAGGTGTACGGCAGGAAGCCACCGTGCAGGGAGATACCGTTAGCAATCGCGGTCATACCGAACTCGCGAACACCGTAGTGGATGTAGTTACCCGCAGCATCTTCGTTGATTGCTTTa                                                                                                                                          <  2:1053492/149‑1 (MQ=255)
     tACGTGCGTATTCCACGAACATCAGGAAGGTGGAGGTGTACGGCAGGAAGCCACCGTGCAGGGAGATACCGTTAGCAATCGCGGTCATACCGAACTCGCGAACACCGTAGTGGATGTAGTTACCCGCAGCATCTTCGTTGATTGCTTTa                                                                                                                                          <  1:32603/149‑1 (MQ=255)
        gTGCGTATTCCACGAACATCAGGAAGGTGGAGGTGTACGGCAGGAAGCCACCGTGCAGGGAGAAACCGATAGCAATCGCGGTCATACCTAACTCGCGAACACCGTAGTGGATGTAGTTACCCGCAGCATATTCGTTGATTGCTTTTGaa                                                                                                                                       >  2:808312/1‑149 (MQ=255)
          gCGTATTCCACGAACATCAGGAAGGTGGAGGTGTACGGCAGGAAGCCACCGTGCAGGGAGATACCGTTAGCAATCGCGGTCATACCGAACTCGCGAACACCGTAGTGGATGTAGTTACCCGCAGCATCTTCGTTGATTGCTTTAGAAcc                                                                                                                                     <  1:272863/149‑1 (MQ=255)
           agTATTCCACGAACATAAGGAAGGTGGAGGTGTACGGCAGGAAGCCCCCGTGCAGGGAGATACCGTTAGCAATCGCGGTCATACCGAACTCGCGAACCCCGTAGTGGATGTATTTACCAGCAGCATCTTCGAGGATTGCTTTAGACCCa                                                                                                                                    >  1:921172/2‑149 (MQ=255)
              aTTCCACGAACATCAGGAAGGTGGAGGTGTACGGCAGGAAGCCACCGTGCAGGGAGATACCGTTAGCAATCGCGGTCATACCGAACTCGCGAACACCGTAGTGGATGTAGTTACCCGCAGCATCTTCGTTGATTGCTTTAGAACCAGAc                                                                                                                                 >  1:692211/1‑149 (MQ=255)
                   aCGAACATCAGGAAGGTGGAGGTGTACGGCAGGAAGCCACCGTGCAGGGAGATACCGTTAGCAATCGCGGTCATACCGAACTCGCGAACCCCGTAGTGGATGTAGTTACCCGCAGCATCTTCGTTGATTGCTTTAGAACCAGACCACAg                                                                                                                            >  1:39796/1‑149 (MQ=255)
                   aCGAACATCAGGAAGGTGGAGGTGTACGGCAGGAAGCCACCGTGCAGGGAGATACCGTTAGCAATCGCGGTCATACCGAACTCGCGAACACCGTAGTGGATGTAGTTACCCGCAGCATCTTCGTTGATTGCTTTAGAACCAGACCACAg                                                                                                                            <  2:649159/149‑1 (MQ=255)
                    cGAACATCAGGAAGGTGGAGGTGTACGGCAGGAAGCCACCGTGCAGGGAGATACCGTTAGCAATCGCGGTCATACCGAACTCGCGAACACCGTAGTGGATGTAGTTACCCGCAGCATCTTCGTTGATTGCTTTAGAACCAGACCACAgg                                                                                                                           <  2:214464/149‑1 (MQ=255)
                    cGAACATCAGGAAGGTGGAGGTGTACGGCAGGAAGCCACCGTGCAGGGAGATACCGTTAGCAATCGCGGTCATACCGAACTCGCGAACACCGTAGTGGATGTAGTTACCCGCAGCATCTTCGTTGATTGCTTTAGAACCAGACCACAgg                                                                                                                           <  1:709938/149‑1 (MQ=255)
                       aCATGAGGAAGGTGGAGGTGTACGGCAGGAAGCCACCGTGCAGGGAGATACCGTTTGCAATCGCGGTCATACCGAACTCGCGAACACCGTAGTGGATGTAGTTACCCGCAGCATCTTCGTTGAATGCTTTAGAACCAGACCACAGGGTc                                                                                                                        >  2:9957/1‑149 (MQ=255)
                         aTCAGGAAGGTGGAGGTGTACGGCAGGAAGCCACCGTGCAGGGAGATACCGTTAGCAATCGCGGTCATACCGAACTCGCGAACACCGTAGTGGATGTAGTTACCCGCAGCATCTTCGTTGATTGCTTTAGAACCAGACCACAGGGTCAg                                                                                                                      >  2:508140/1‑149 (MQ=255)
                            aGGAAGGTGGAGGTGTACGGCAGGAAGCCACCGTGCAGGGAGATACCGTTAGCAATCGCGGTCATACCGAACTCGCGAACACCGTAGTGGATGTAGTTACCCGCAGCATCTTCGTTGATTGCTTTAGAACCAGACCACAGGGTCAGGtt                                                                                                                   >  1:554322/1‑149 (MQ=255)
                               aaGTTGGAGGTGTACGGGAGGAAGCCACCGTGCAGGGTGATACCGTTAGCAATCGCGGTCATACGGTACTCGCGAACACCGTAGTGGATGTTTTTACCCGCAGCTTCTTCGTTGATTGCTTTAGAACCAGTCCACAGGGTCAGGTTAGa                                                                                                                <  2:921172/149‑1 (MQ=255)
                               aaGGTTGAGGTGGACGGCAGGAAGCCACCGGGCAGGGAGATACCGTTAGCAATCGCGGGCATACCGAACTCGCGAACACCGTAGTGGATGTAGTTACCCGCAGCATCTTCGTTGATTGCTTTAGAACCAGACCACAGGGTCAGGTTAGa                                                                                                                <  1:1036196/149‑1 (MQ=255)
                                   tGGAGGTGTACGGCAGGAAGCCACCGTGCAGGGAGATACCGTTAGCAATCGCGGTCATACCGAACTCGCGAACACCGTAGTGGATGTAGTTACCCGCAGCATCTTCGTTGATTGCTTTAGAACCAGCCCACAGGGTCAGGTTAGCCGgc                                                                                                            >  1:504041/1‑149 (MQ=255)
                                   tGGAGGTGTACGGCAGGAAGCCACCGTGCAGGGAGATACCGTTAGCAATCGCGGTCATACCGAACTCGCGAACACCGTAGTGGATGTAGTTACCCGCAGCATCTTCGTTGATTGCTTTAGAACCAGACCACAGGGTCAGGTTAGACGgc                                                                                                            >  2:245417/1‑149 (MQ=255)
                                        gtgtACGGCAGGAAGCCACCGTGCAGGGAGATACCGTTAGCAATCGCGGTCATACCGAACTCGCGAACACCGTAGTGGATGTAGTTACCCGCAGCATCTTCGTTGATTGCTTTAGAACCAGACCACAGGGTCAGGTTAGACGGCGCCAg                                                                                                       >  2:330928/1‑149 (MQ=255)
                                        gggtACGGCAGGAAGCCACCGTGCAGGGAGATACCGTTAGCAATCGCGGTCATACCGAACTCGCGAACACCGTAGTGGATGTAGTTACCCGCAGCATCTTCGTTGATTGCGTTAGAACCAGACCACAGGGTCAGGTTAGACGGCGCCAg                                                                                                       <  2:790958/147‑1 (MQ=255)
                                           tACGGCAGGAAGCCACCGTGCAGGGAGATACCGTTAGCAATCGCGGTCATACCGAACTCGCGAACACCGTAGTGGATGTAGTTACCCGCAGCATCTTCGTTGGTTGCTTTAGAACCAGACCACAGGGTCAGGTTAGACGGCGCCAGGTc                                                                                                    >  1:879434/1‑149 (MQ=255)
                                           tACGGCAGGAAGCCACCGTGCAGGGAGATACCGTTAGCAATCGCGGTCATACCGAACTCGCGAACACCGGAGTGGATGTAGTTACCCGCCGCATCTTCGTTGGTTGCTTTAGAACCAGACCACAGGGTCAGGTTAGACGGCGCCAGGTc                                                                                                    >  1:879438/1‑149 (MQ=255)
                                              ggCAGGAAGCCACCGTGCAGGGAGATACCGTTAGCAATCGCGGTCATACCGAACTCGCGAACACCGTAGTGGATGTAGTTACCCGCAGCATCTTCGTTGATTGCTTTAGAACCAGACCACAGGGTCAGGTTAGACGGCGCCAGGTcagc                                                                                                 >  1:477415/1‑149 (MQ=255)
                                              ggCAGGAAGCCACCGTGCAGGGAGATACCGTTAGCAATCGCGGGCCTACCGAACTCGCGAACACCGTAGTGGATGTAGTTACCCGCAGCATCTTCGTTGATTGCTTTAGAACCAGACCACAGGGGCAGGTTAGACGGCGCCAGGTcagc                                                                                                 <  2:691958/149‑1 (MQ=255)
                                                 aGGAAGCCACCGTGCAGGGAGATACCGTTAGCAATCGCGGTCATACCGAACTAGCGAACACCGTAGTGGAGGTCGTTACCCGCAGCATCTTCGTTGATTGCTTTAGAACCCGACCACAGGGTCAGGTTAGACGGCGCCCGGTcagcagc                                                                                              >  1:1111332/1‑148 (MQ=255)
                                                   gAAGCCACCGTGCAGGGAGATACCGTTAGCAATCGCGGTCATACCGAACTCGCGAACACCGTAGTGGATGTAGTTACCCGCAGCATCTTCGTTGATTGCTTTAGAACCAGACCACAGGGTCAGGTTAGACGGCGCCAGGTCAGCAGAAc                                                                                            >  1:274562/1‑149 (MQ=255)
                                                   gAAGCCACCGTGCAGGGAGAAACCGTTAGCAATCGCGGTCATACGAAAAGCGCAAACACGGAAGTGTATGAAGTAACCCGAAGAATGTACTTTGATGGCTTAAGAACAACACCCAAGGGGCCGGTAACAAGGCGCCAGTTCAGCAGAAc                                                                                            >  2:898928/1‑149 (MQ=255)
                                                        cACCGTGCAGGGAGATACCGTTAGCAATCGCGGTCATACCGAACTCGCGAACACCGTAGTGGATGTAGTTACCCGCAGCATCTTCGTTGATTGCTTTAGAACCAGACCACAGGGTCAGGTTAGACGGCGCCAGGTCAGCAGAAccgccg                                                                                       >  2:639161/1‑149 (MQ=255)
                                                                aGGGAGATACCGTTAGCAATCGCGGTCATACCGAAATCGCGCACACCGTAGTGGATGTAGATAACCGCAGCATCTTCTTTCAGAGCTTTAGAACCACACCACAGCGTCAGTTTTGACGGCGGCAGGTCATCAGAACCGCCGAGGAATTc                                                                               >  2:245826/1‑149 (MQ=255)
                                                                    agaTACCGTTATCAATCGCGGTCATACCGAACTCGCGAACACCGTAGTGGATGTAGTTACCCGCAGCATCTTCGTTGATTGCATTAGAACCAGACCACAGGGTCAGGTTAGACGGCGCCAGGTCAGCAGAACCGCCGAGGAATTCCGGc                                                                           <  2:95263/149‑1 (MQ=255)
                                                                       tACCGTTAGCAATCGCGGTCATACCGAACTCGCGAACACCGTAGTGGATGTAGTTACCCGCAGCATCTTCGTTGATTGCTTTAGAACCAGACCACAGGGTCAGGTTAGACGGCGCCAGGTCAGCAGAACCGCCGAGGAATTCCGGCAAc                                                                        <  1:245417/149‑1 (MQ=255)
                                                                        aCCGTTAGCAATCGCGGTCATACCGAACTCGCGAACACCGTAGTGGATGTAGTTACCCGCAGCATCTTCGTTGATTGCTTTAGAACCAGACCACAGGGTCAGGTTAGACGGCGCCAGGTCAGCAGAACCGCCGa                                                                                      <  2:290579/134‑1 (MQ=255)
                                                                        aCCGTTAGCAATCGCGGTCATACCGAACTCGCGAACACCGTAGTGGATGTAGTTACCCGCAGCATCTTCGTTGATTGCTTTAGAACCAGACCACAGGGTCAGGTTAGACGGCGCCAGGTCAGCAGAACCGCCGa                                                                                      >  1:290579/1‑134 (MQ=255)
                                                                         ccGTTAGCAATCGCGGTCATACCGAACTCGCGAACACCGTAGTGGATGTAGTTACCCGCAGCATCTTCGTTGATTGCTTTAGAACCAGACCACAGGGTCAGGTTAGACGGCGCCAGGTCAGCAGAACCGCCGAGGAATTCCGGCAACAg                                                                      <  1:508140/149‑1 (MQ=255)
                                                                           gTTAGCAATCGCGGTCATACCGAACTCGCGAACACCGTAGTGGATGTAGTTACCCGCAGCATCTTCGTTGATTGCTTTAGAACCAGACCACAGGGTCAGGTTTGACGGCGCCAGGTCAGCAGAACCGCCGAGGAATTCCGGCAACAGCg                                                                    >  2:75399/1‑149 (MQ=255)
                                                                           gTTAGCAATCGCGGTCATACCGAACTCGCGAACACCGTAGTGGATGTAGTTACCCGCAGCATCTTCGTTGATTGCTTTAGAACCAGACCACAGGGTCAGGTTAGACGGCGCCAGGTCAGCAGAACCGCCGAGGAATTCCGGCAACAGCg                                                                    >  1:25256/1‑149 (MQ=255)
                                                                            ttAGCAATCGCGGTCATACCGAACTCGCGAACACCGTAGTGGATGTAGTTACCCGCAGCATCTTCGTTGATTGCTTTAGAACCAGACCACAGGGTCAGGTTAGACGGCGCCAGGTCAGCAGAACCGCCGAGGAATTCCGGCAACAGCgg                                                                   >  1:618860/1‑149 (MQ=255)
                                                                             tAGCAATCGCGGTCATACCGACCTCGCGAACACCGTAGTGGATGTAGTTCCCCGCAGCATCTTCGGTGAATGCTTTAGAACCATAACACAGGGTCAGGTTAGACGGCGCCAGGTCAGCAGCCCCGCCGAGGAATTCCTGCAAAAGCggc                                                                  >  2:994471/1‑148 (MQ=255)
                                                                             aaGCAATCGCGGTCATACCGAACTCGCGAACACCGTAGTGGATGTAGTTACCCGCAGCATCTTCGGTGATTGCTTTAGACCCAGACCACAGGGTCAGGTTAGACGGCGCCAGGTCAGCAGAACCGCCGAGGAATTCCGGCAACAGCGGa                                                                  <  1:153400/148‑1 (MQ=255)
                                                                              aGCAATCGCGGTCATACCGAACTCGCGAACACCGTAGTGGATGTAGTTACCCGCAGCATCTTCGTTGATTGCTTTAGAACCAGACCACAGGGTCAGGTTAGACGGCGCCAGGTCAGCAGAACCGCCGAGGAATTCCGGCAACAGCGGAc                                                                 <  1:787501/149‑1 (MQ=255)
                                                                               gCAATCGCGGTCATACCGAACTCGCGAACACCGTAGTGGATGTAGTTACCCGCAGCATCTTCGTTGATTGCTTTAGAACCAGACCACAGGGTCAGGTTAGACGGCGCCAGGTCAGCAGAACCGCCGAGGAATTCCGGCAACAGCGGAcc                                                                <  2:692211/149‑1 (MQ=255)
                                                                                  atagcgGTCATACCGCACTCGCGAACACCGTAGTGGATGTAGATACCCGCAGCATCTTCCTTGATTGCTTTAGAACCAGACCACAGGGTCAGGTTAGACGGCGCCAGGTCAGCAGAACCGCCGAGGAATTCCGGCAACAGCGGACCGaa                                                             <  2:1175959/146‑1 (MQ=255)
                                                                                   tCGCGGTCATACCGAACTCGCGAACACCGTAGTGGATGTAGTTACCCGCAGCATCTTCGTTGATTGCTTTAGAACCAGACCACAGGGTCAGGTTAGACGGCGCCAGGTCAGCAGAAccg                                                                                          >  1:440320/1‑119 (MQ=255)
                                                                                   tCGCGGTCATACCGAACTCGCGAACACCGTAGTGGATGTAGTTACCCGCAGCATCTTCGTTGATTGCTTTAGAACCAGACCACAGGGTCAGGTTAGACGGCGCCAGGTCAGCAGAAccg                                                                                          <  2:440320/119‑1 (MQ=255)
                                                                                             aCCGAACTCGCGAACACCGTAGTGGATGTAGTTACCCGCAGCATCTTCGTTGATTGCTTTAGAACCAGACCACAGGGTCAGGTTAGACGGCGCCAGGTCAGCAGAACCGCCGAGGAATTCCGGCAACAGCGGACCGTACGCTTCGATAg                                                  >  1:1007915/1‑149 (MQ=255)
                                                                                              ccGAACTCGCGAACACCGTAGTGGATGTAGTTACCCGCAGCATCTTCGTTGATTGCTTTAGAACCAGACCACAGGGTCAGGTTAGACGGCGCCAGGTCAGCAGAACCGCCGAGGAATTCCGGCAACAGCGGACCGAACGCTTCGATAGc                                                 >  1:986781/1‑149 (MQ=255)
                                                                                                    tCGCGAACACCGTAGTGGATGTAGTTACCCGCAGCATCTTCGTTGATTGCTTTAGAACCAGACCACAGGGTCAGGTTAGACGGCGCCAGGTCAGCAGAACCGCCGAGGAATTCCGGCAACAGCGGTCCGCACGCTTCGATAGCATTCTg                                           >  2:391467/1‑149 (MQ=255)
                                                                                                    tCGCGAACACCGTAGTGGATGTAGTTACCCGCAGCATCTTCGTTGATTGCTTTAGAACCAGACCACAGGGTCAGGTTAGACGGCGCCAGGTCAGCAGAACCGCCGAGGAATTCCGGCAACAGCGGACCGAACGCTTCGATAGCATTCTg                                           >  2:1084584/1‑149 (MQ=255)
                                                                                                    tCGCGAACACCGTAGTGGATGTAGTTACCCGCAGCATCTTCGTTGATTGCTTTAGAACCAGACCACAGGGTCAGGTTAGACGGCGCCAGGTCAGCAGAACCGCCGAGGAATTCCGGCAACAGCGGACCGAACGCTTCGATAGCATTCTg                                           >  1:13406/1‑149 (MQ=255)
                                                                                                      gcgAACACCGTAGTGGATGTAGTTACCCGCAGCATCTTCGTTGATTGCTTTAGAACCAGACCACAGGGTCAGGTTAGACGGCGCCAGGTCAGCAGAACCGCCGAGGAATTCCGGCAACAGCGGACCGAACGCTTCGATAGCATTCTgag                                         >  1:890571/1‑149 (MQ=255)
                                                                                                      gcgAACACCGTAGTGGATGTAGTTACCCGCAGCATCTTCGTTGATTGCTTTAGAACCAGACCACAGGGTCAGGTTAGACGGCGCCAGGTCAGCAGAACCGCCGAGGAATTCCGGCAACAGCGGACCGAACGCTTCGATAGCATTCTgag                                         >  1:285315/1‑149 (MQ=255)
                                                                                                              cGTAGTGGATGTAGTTACCCGCAGCATCTTCGTTGATTGCTTTAGAACCAGACCACAGGGTCAGGTTAGACGGCGCCAGGTCAGCAGAACCGCCGAGGAATTCCGGCAACAGCGGACCGAACGCTTCGATAGCATTCTGAGACGCTTTa                                 >  2:1166130/1‑149 (MQ=255)
                                                                                                               gTAGTGGATGGAGTTACCCGCAGCATGTTCGTGGGTTGCTTTAGCACCAGACCACCGGTTCAGGTTAGACGGCGCCAGGTCAGCAGAACCGCCGAGGAATTCCGGCAACAGCGGACCGAACGCTTCGATAGCATTCTGAGACGCTTTAc                                <  2:879438/149‑1 (MQ=255)
                                                                                                               gTAGTGGAGTTAGTTACCCGCAGCAGCTTCGTGTGGTGCTTTAGAACCAGACCACAGGGTCAGGTTAGACGGCGCCAGGTCAGCAGAACCGCCGAGGAATTCCGGCAACAGCGGACCGAACGCTTCGATAGCATTCTGAGACGCTTTAc                                <  2:879434/149‑1 (MQ=255)
                                                                                                                   tGGATGTAGTTACCCGCAGCATCTTCGTTGATTGCTTTAGAACCAGACCACAGGGTCAGGTTAGACGGCGCCAGGTCAGCAGAACCGCCGAGGAATTCCGGCAACAGCGGACCGAACGCTTCGATAGCATTCTGAGACGCTTTACGGCt                            >  1:454470/1‑149 (MQ=255)
                                                                                                                        gTAGTTACCCGCAGCATCATCGTTGATTGCTTTAGAACCAGACCACAGGGTCAGGTTAGACGGCGCCAGGTCAGCAGAACCGCCGAGGAATTCCGGCAACAGCGGACCGAACGCTTCGAAAGCATTCTGAGACGCTTTACGGCTGGCGa                       >  1:573026/1‑149 (MQ=255)
                                                                                                                           gTTACCCGCAGCATCTTCGTTGATCGCTTTAGAACCAGACCACAGGGTCAGGTTAGACGGCGCCAGGTCAGCAGAACCGCCGAGGAATTCCGGCAACAGCGGACCGAACGCTGCGATAGCATTCTGAGACGCTTTACGGCTGGCGAttt                    >  1:747964/1‑149 (MQ=255)
                                                                                                                            ttACCCGCAGCATCTTCGTTGATTGCTTTAGAACCAGACCACAGGGTCAGGTTAGACGGCGCCAGGTCAGCAGAACCGCCGAGGAATTCCGGCAACAGCGGACCGCACGCTTCGATAGCAGTCTGAGACGCTTTACGGCTGGCGAtttt                   >  2:663261/1‑149 (MQ=255)
                                                                                                                            ttACCCGCAGCATCTTCGTTGATTGCTTTAGAACCAGACCACAGGGTCAGGTTAGACGGCGCCAGGTCAGCAGAACCGCCGAGGAATTCCGGCAACAGCGGACCGAACGCTTCGATAGCATTCTGAGACGCTTTACGGCTGGCGAtttt                   >  1:1135520/1‑149 (MQ=255)
                                                                                                                             tACCCGCAGCATCTTCGTTGATTGCTTTAGAACCAGACCACCGGGGCAGGTTAGACGGCGCCAGGTCAGCAGAACCGCCGAGGAATTCCGGCAACAGCGGACCGAACGCTTCGATAGCATTCTGAGACGCTTTACGGCTGGCGATTTTc                  <  2:554322/149‑1 (MQ=255)
                                                                                                                             tACCCGCAGCATCTTCGTTGATTGCTTTAGAACCAGACCACAGGGTCAGGTTAGACGGCGCCAGGTCAGCAGAACCGCCGAGGAATTCCGGCAACAGCGGACCGAACGCTTCGATAGCATTCTGAGACGCTTTACGGCTGGCGATTTTc                  >  2:430754/1‑149 (MQ=255)
                                                                                                                                ccGCAGCATCTTCGTTGATTGCTTTAGAACCAGACCACAGGGTCAGGTTAGACGGCGCCAGGTCAGCAGAACCGCCGAGGAATTCCGGCAACAGCGGACCGAACGCTTCGATAGCATTCTGAGACGCTTTACGGCTGGCGCTTTTCGcc               >  1:518665/1‑149 (MQ=255)
                                                                                                                                   cagcaTCTTCGTTGATTGCTTTAGAACCAGACCACAGGGTCAGGTTAGACGGCGCCAGGTCAGCAGAACCGCCGAGGAATTCCGGCTACAGCGGACCGAACGCTTCGATAGCATTCTGAGACGCTTTACGGCTGGCGATTTTCGCCGGa            >  2:367198/1‑149 (MQ=255)
                                                                                                                                    agcaTCTTCGTTGATTGCTTTAGAACCAGACCACAGGGTCAGGTTAGACGGCGCCAGGTCAGCAGAACCGCCGAGGAATTCCGGCAACAGCGGACCGAACGCTTCGATAGCATTCTGAGACGCTTTACGGCTGGCGATTTTCGCCGGAt           >  2:345895/1‑149 (MQ=255)
                                                                                                                                           tCGTTGATTGCTTTAGAACCAGACCACAGGGTCAGGTTAGACGGCGCCAGGTCAGCAGAACCGCCGAGGAATTCCGGCAACAGCGGACCGAACGCTTCGATAGCATTCTGAGACGCTTTACGGCTGGCGATTTTCGCCGGATTAGCCTg    >  1:130215/1‑149 (MQ=255)
                                                                                                                                           tCGTTGATTGCTTTAGAACCAGACCACAGGGTCAGGTTAGACGGCGCCAGGTCAGCAGAACCGCCGAGGAATTCCGGCAACAGCGGACCGAACGCTTCGATAGCATTCTGAGACGCTTTACGGCTGGCGATTTTCGCCGGATTAGCCTg    >  1:576453/1‑149 (MQ=255)
                                                                                                                                             gTTGAGTGCTTTAGAACCAGTCCCGAGGGGCAGCGAAGACGGCGCCCGGGCAGCAGAACCTCCGAGGAAGTCCGGCAACAGCTGACCGACCGCTTCGATAGCATTCTGCGAAGCTTTACGGCTGGCGTTTTTCGCCGGATTAGCCTGCa  <  1:994471/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
GGCGTTACGTGCGTATTCCACGAACATCAGGAAGGTGGAGGTGTACGGCAGGAAGCCACCGTGCAGGGAGATACCGTTAGCAATCGCGGTCATACCGAACTCGCGAACACCGTAGTGGATGTAGTTACCCGCAGCATCTTCGTTGATTGCTTTAGAACCAGACCACAGGGTCAGGTTAGACGGCGCCAGGTCAGCAGAACCGCCGAGGAATTCCGGCAACAGCGGACCGAACGCTTCGATAGCATTCTGAGACGCTTTACGGCTGGCGATTTTCGCCGGATTAGCCTGCA  >  NC_000913/3080304‑3080593

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑