Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,860,740 C→G 5.3% M165I (ATG→ATC)  ygbI ← DeoR family putative transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,860,7400CG 5.3% 270.6 / 5.7 76M165I (ATG→ATC) ygbIDeoR family putative transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (37/35);  new base G (2/2);  total (39/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.23e-01

ACCTGACCATATTTCGTCGCATCACAGACCAAAACTCGCTGGCGGCTGGCACTGGCAATCGCCCGTTTCACCGTGACTTTATCTTCCGCTGGCGTAGAAATCCCCCGCACACTCCATGACGATGCAGAAATAAAAGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACGGTTTTCCCGACACACTGCACCGCCAGTGTGAATAATTGTGCAATTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTTGTGACCACAGT  >  NC_000913/2860594‑2860883
                                                                                                                                                  |                                                                                                                                               
aCCTGACCATATTTCGTCGCATCACAGACCAAAACTCGCTGGCGGCTGGCACTGGCAATCGCCCGTTTCACCGTGACTTTATCTTCCGCTGGCGTAGAAATCCCCCGCACACTCCATGACGATGCAGAAATAAAAGCCTGATCAatcat                                                                                                                                               <  1:289571/149‑1 (MQ=255)
aCCTGACCATATTTCGTCGCATCACAGACCAAAACTCGCTGGCGGCTGGCACTGGCAATCGCCCGTTTCACCGTGACTTTATCTTCCGCTGGCGTAGAAATCCCCCGCACACTCCATGACGATGCAGAAATAAAAGCCTGATCAatcat                                                                                                                                               <  4:593503/149‑1 (MQ=255)
  cTGACCATATTTCGTCGCATCACAGACCAAAACTCGCTGGCGGCTGGCACTGGCAATCGCCCGTTTCACCGTGACTTTATCTTCCGCTGGCGTAGAAATCCCCCGCACACTCCATGACGATGCAGAAATAAAAGCCTGATCAAtcatca                                                                                                                                             >  1:536429/1‑149 (MQ=255)
  cTGACCATATTTCGTCGCATCACAGACCAAAACTCGCTGGCGGCTGGCACTGGCAATCGCCCGTTTCACCGTGACTTTATCTTCCGCTGGCGTAGAAATCCCCCGCACACTCCATGACGATGCAGAAATAAAAGCCTGATCAAtcatca                                                                                                                                             <  1:247366/149‑1 (MQ=255)
   tGACCATATTTCGTCGCATCACAGACCAAAACTCGCTGGCGGCTGGCACTGGCAATCGCCCGTTTCACCGTGACTTTATCTTCCGCTGGCGTAGAAATCCCCCGCACACTCCATGACGATGCAGAAATAAAAGCCTGATCAATCATCAg                                                                                                                                            <  2:70126/149‑1 (MQ=255)
    gACCATATTTCGTCGCATCACAGACCAAAACTCGCTGGCGGCTGGCACTGGCAATCGCCCGTTTCACCGTGACTTTATCTTCCGCTGGCGTAGAAATCCCCCGCACACTCCATGACGATGCAGAAATAAAAGCCTGATCAATCATCAgg                                                                                                                                           <  4:160672/149‑1 (MQ=255)
      ccATATTTCGTCGCATCACAGACCAAAACTCGCTGGCGGCTGGCACTGGCAATCGCCCGTTTCACCGTGACTTTATCTTCCGCTGGCGTAGAAATCCCCCGCACACTCCATGACGATGCAGAAATAAAAGCCTGATGATTGATCAGGCt                                                                                                                                         <  1:318820/149‑1 (MQ=255)
       cATATTTCGTCGCATCACAGACCAAAACTCGCTGGCGGCTGGCACTGGCAATCGCCCGTTTCACCGTGACTTTATCTTCCGCTGGCGTAGAAATCCCCCGCACACTCCATGACGATGCAGAAATAAAAGCCTGATGATTGATCAGGCtt                                                                                                                                        <  3:712378/149‑2 (MQ=255)
       cATATTTCGTCGCATCACAGACCAAAACTCGCTGGCGGCTGGCACTGGCAATCGCCCGTTTCACCGTGACTTTATCTTCCGCTGGCGTAGAAATCCCCCGCACACTCCATGACGATGCAGAAATAAAAGCCTGATCAATCATCAGGCTg                                                                                                                                        >  1:192782/1‑149 (MQ=255)
       cATATTTCGTCGCATCACAGACCAAAACTCGCTGGCGGCTGGCACTGGCAATCGCCCGTTTCACCGTGACTTTATCTTCCGCTGGCGTAGAAATCCCCCGCACACTCCATGACGATGCAGAAATAAAAGCCTGATCAATCATCAGGCTg                                                                                                                                        >  3:510217/1‑149 (MQ=255)
           tttCGTCGCATCACAGACCAAAACTCGCTGGCGGCTGGCACTGGCAATCGCCCGTTTCACCGTGACTTTATCTTCCGCTGGCGTAGAAATCCCCCGCACACTCCATGACGATGCAGAAATAAAAGCCTGATCAATCATCAGGCTGCgca                                                                                                                                    >  3:227244/1‑149 (MQ=255)
           tttCGTCGCATCACAGACCAAAACTCGCTGGCGGCTGGCACTGGCAATCGCCCGTTTCACCGTGACTTTATCTTCCGCTGGCGTAGAAATCCCCCGCACACTCCATGACGATGCAGAAATAAAAGCCTGATCAATCATCAGGCTGCgca                                                                                                                                    >  4:482248/1‑149 (MQ=255)
               gtcgCATCACAGACCAAAACTCGCTGGCGGCTGGCACTGGCAATCGCCCGTTTCACCGTGACTTTATCTTCCGCTGGCGTAGAAATCCCCCGCACACTCCATGACGATGCAGAAATAAAAGCCTGATCAATCATCAGGCTGCGCAGCAt                                                                                                                                >  3:64310/1‑149 (MQ=255)
                  gCATCACAGACCAAAACTCGCTGGCGGCTGGCACTGGCAATCGCCCGTTTCACCGTGACTTTATCTTCCGCTGGCGTAGAAATCCCCCGCACACTCCATGACGATGCAGAAATAAAAGCCTGATCAATCATCAGGCTGCGCAGCATGGt                                                                                                                             <  3:482248/149‑1 (MQ=255)
                   cATCACAGACCAAAACTCGCTGGCGGCTGGCACTGGCAATCGCCCGTTTCACCGTGACTTTATCTTCCGCTGGCGTAGAAATCCCCCGCACACTCCATGACGATGCAGAAATAAAAGCCTGATCAATCATCAGGCTGCGCAGCATGGTc                                                                                                                            <  4:57771/149‑1 (MQ=255)
                       acaGACCAAAACTCGCTGGCGGCTGGCACTGGCAATCGCCCGTTTCACCGTGACTTTATCTTCCGCTGGCGTAGAAATCCCCCGCACACTCCATGACGATGCAGAAATAAAAGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAg                                                                                                                        >  3:280714/1‑149 (MQ=255)
                         aGACCAAAACTCGCTGGCGGCTGGCACTGGCAATCGCCCGTTTCACCGTGACTTTATCTTCCGCTGGCGTAGAAATCCCCCGCACACTCCATGACGATGCAGAAATAAAAGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCg                                                                                                                      >  3:501380/1‑149 (MQ=255)
                           aCCAAAACTCGCTGGCGGCTGGCACTGGCAATCGCCCGTTTCACCGTGACTTTATCTTCCGCTGGCGTAGAAATCCCCCGCACACTCCATGACGATGCAGAAATAAAAGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGc                                                                                                                    >  1:55592/1‑149 (MQ=255)
                           aCCAAAACTCGCTGGCGGCTGGCACTGGCAATCGCCCGTTTCACCGTGACTTTATCTTCCGCTGGCGTAGAAATCCCCCGCACACTCCATGACGATGCAGAAATAAAAGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGc                                                                                                                    >  1:12024/1‑149 (MQ=255)
                             cAAAACTCGCTGGCGGCTGGCACTGGCAATCGCCCGTTTCACCGTGACTTTATCTTCCGCTGGCGTAGAAATCCCCCGCACACTCCATGACGATGCAGAAATAAAAGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCtt                                                                                                                  >  2:608520/1‑149 (MQ=255)
                             cAAAACTCGCTGGCGGCTGGCACTGGCAATCGCCCGTTTCACCGTGACTTTATCTTCCGCTGGCGTAGAAATCCCCCGCACACTCCATGACGATGCAGAAATAAAAGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCtt                                                                                                                  <  4:178400/149‑1 (MQ=255)
                              aaaaCTCGCTGGCGGCTGGCACTGGCAATCGCCCGTTTCACCGTGACTTTATCTTCCGCTGGCGTAGAAATCCCCCGCACACTCCATGACGATGCAGAAATAAAAGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTc                                                                                                                 <  2:215590/149‑1 (MQ=255)
                              aaaaCTCGCTGGCGGCTGGCACTGGCAATCGCCCGTTTCACCGTGACTTTATCTTCCGCTGGCGTAGAAATCCCCCGCACACTCCATGACGATGCAGAAATAAAAGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTc                                                                                                                 <  2:207095/149‑1 (MQ=255)
                               aaaCTCGCTGGCGGCTGGCACTGGCAATCGCCCGTTTCACCGTGACTTTATCTTCCGCTGGCGTAGAAATCCCCCGCACACTCCATGACGATGCAGAAATAAAAGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTcc                                                                                                                >  4:626994/1‑149 (MQ=255)
                                   tCGCTGGCGGCTGGCACTGGCAATCGCCCGTTTCACCGTGACTTTATCTTCCGCTGGCGTAGAAATCCCCCGCACACTCCATGACGATGCAGAAATAAAAGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGa                                                                                                            <  2:572960/149‑1 (MQ=255)
                                   tCGCTGGCGGCTGGCACTGGCAATCGCCCGTTTCACCGTGACTTTATCTTCCGCTGGCGTAGAAATCCCCCGCACACTCCATGACGATGCAGAAATAAAAGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGa                                                                                                            <  1:75985/149‑1 (MQ=255)
                                    cGCTGGCGGCTGGCACTGGCAATCGCCCGTTTCACCGTGACTTTATCTTCCGCTGGCGTAGAAATCCCCCGCACACTCCATGACGATGCAGAAATAAAAGCCTGATCAATCATCAGGCTGCGTAGCATGGTCGCAGCGGCTTCCCCGac                                                                                                           >  3:115455/1‑149 (MQ=255)
                                       tGGCGGCTGGCACTGGCAATCGCCCGTTTCACCGTGACTTTATCTTCCGCTGGCGTAGAAATCCCCCGCACACTCCATGACGATGCAGAAATAAAAGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGAcaca                                                                                                        >  1:517273/1‑149 (MQ=255)
                                           ggcTGGCACTGGCAATCGCCCGTTTCACCGTGACTTTATCTTCCGCTGGCGTAGAAATCCCCCGCACACTCCATGACGATGCAGAAATAAAAGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGaa                                                                                                    <  2:328968/149‑1 (MQ=255)
                                             ctggcactggcaATCGCCCGTTTCACCGTGACTTTATCTTCCGCTGGCGTAGAAATCCCCCGCACACTCCATGACGATGCAGAAATAAAAGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACg                                                                                                  <  1:427933/149‑1 (MQ=255)
                                              tggcactggcaATCGCCCGTTTCACCGTGACTTTATCTTCCGCTGGCGTAGAAATCCCCCGCACACTCCATGACGATGCAGAAATAAGAGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACgg                                                                                                 >  3:691944/1‑149 (MQ=255)
                                               ggcactggcaATCGCCCGTTTCACCGTGACTTTATCTTCCGCTGGCGTAGAAATCCCCCGCACACTCCATGACGATGCAGAAATAAAAGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACGGt                                                                                                >  2:767315/1‑149 (MQ=255)
                                               ggcactggcaATCGCCCGTTTCACCGTGACTTTATCTTCCGCTGGCGTAGAAATCCCCCGCACACTCCATGACGATGCAGAAATAAAAGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACGGt                                                                                                <  4:1824/149‑1 (MQ=255)
                                                  actggcaATCGCCCGTTTCACCGTGACTTTATCTTCCGCTGGCGTAGAAATCCCCCGCACACTCCATGACGATGCAGAAATAAAAGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACGGtttt                                                                                             >  1:394502/1‑149 (MQ=255)
                                                    tggcaATCGCCCGTTTCACCGTGACTTTATCTTCCGCTGGCGTAGAAATCCCCCGCACACTCCATGACGATGCAGAAATAAAAGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACGGTTTTcc                                                                                           <  4:357321/149‑1 (MQ=255)
                                                     ggcaATCGCCCGTTTCACCGTGACTTTATCTTCCGCTGGCGTAGAAATCCCCCGCACACTCCATGACGATGCAGAAATAAAAGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACGGTTTTccc                                                                                          <  4:510217/149‑1 (MQ=255)
                                                        aaTCGCCCGTTTCACCGTGACTTTATCTTCCGCTGGCGTAGAAATCCCCCGCACACTCCATGACGATGCAGAAATAAAAGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACGGTTTTCCCGac                                                                                       >  3:6689/1‑149 (MQ=255)
                                                            gCCCGTTTCACCGTGACTTTATCTTCCGCTGGCGTAGAAATCCCCCGCACACTCCATGACGATGCAGAAATAAAAGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACGGTTTTCCCGacacac                                                                                   >  2:416170/1‑149 (MQ=255)
                                                            gCCCGTTTCACCGTGACTTTATCTTCCGCTGGCGTAGAAATCCCCCGCACACTCCATGACGATGCAGAAATAAAAGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACGGTTTTCCCGacacac                                                                                   <  2:698782/149‑1 (MQ=255)
                                                             cccGTTTCACCGTGACTTTATCTTCCGCTGGCGTAGAAATCCCCCGCACACTCCATGACGATGCAGAAATAAAAGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACGGTTTTCCCGACACACt                                                                                  >  3:538714/1‑149 (MQ=255)
                                                              ccGTTTCACCGTGACTTTATCTTCCGCTGGCGTAGAAATCCCCCGCACACTCCATGACGATGCAGAAATAAAAGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACGGtttt                                                                                             >  1:304975/1‑137 (MQ=255)
                                                              ccGTTTCACCGTGACTTTATCTTCCGCTGGCGTAGAAATCCCCCGCACACTCCATGACGATGCAGAAATAAAAGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACGGtttt                                                                                             <  2:304975/137‑1 (MQ=255)
                                                              ccGTTTCACCGTGACTTTATCTTCCGCTGGCGTAGAAATCCCCCGCACACTCCATGACGATGCAGAAATAAAAGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACGGTTTTCCCGACACACTg                                                                                 >  4:505793/1‑149 (MQ=255)
                                                                gTTTCACCGTGACTTTATCTTCCGCTGGCGTAGAAATCCCCCGCACACTCCATGACGATGCAGAAATAAAAGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACGGTTTTCCCGACACACTGCa                                                                               <  3:44448/149‑1 (MQ=255)
                                                                               tATCTTCCGCTGGCGTAGAAATCCCCCGCACACTCCATGACGATGCAGAAATAAAAGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACGGTTTTCCCGACACACTGCACCGCCAGTGTGaata                                                                >  3:421342/1‑149 (MQ=255)
                                                                                      cGCTGGCGTAGAAATCCCCCGCACACTCCATGACGATGCAGAAATAAAAGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACGGTTTTCCCGACACACTGCACCGCCAGTGTGAATAATTGTGc                                                         <  1:33462/149‑1 (MQ=255)
                                                                                         tGGCGTAGAAATCCCCCGCACACTCCATGACGATGCAGAAATAAAAGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACGGTTTTCCCGACACACTGCACCGCCAGTGTGAATAATTGTGCAAt                                                      >  2:491287/1‑149 (MQ=255)
                                                                                             gTAGAAATCCCCCGCACACTCCATGACGATGCAGAAATAAAAGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACGGTTTTCCCGACACACTGCACCGCCAGTGTGAATAATTGTGCAATTACt                                                  <  2:71635/149‑1 (MQ=255)
                                                                                             gTAGAAATCCCCCGCACACTCCATGACGATGCAGAAATAAAAGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACGGTTTTCCCGACACACTGCACCGCCAGTGTGAATAATTGTGCAATTACt                                                  >  1:34176/1‑149 (MQ=255)
                                                                                              tAGAAATCCCCCGCACACTCCATGACGATGCAGAAATAAAAGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACGGTTTTCCCGACACACTGCACCGCCAGTGTGAATAATTGTGCAATTACtg                                                 <  1:327637/149‑1 (MQ=255)
                                                                                                  aaTCCCCCGCACACTCCATGACGATGCAGAAATAAAAGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACGGTTTTCCCGACACACTGCACCGCCAGTGTGAATAATTGTGCAATTACtgttgt                                             <  3:613690/149‑1 (MQ=255)
                                                                                                   aTCCCCCGCACACTCCATGACGATGCAGAAATAAAAGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACGGTTTTCCCGACACACTGCACCGCCAGTGTGAATAATTGTGCAATTACTGTTGTc                                            <  3:395022/149‑1 (MQ=255)
                                                                                                    tCCCCCGCACACTCCATGACGATGCAGAAATAAAAGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACGGTTTTCCCGACACACTGCACCGCCAGTGTGAATAATTGTGCAATTACTGTTGTCg                                           >  1:22881/1‑149 (MQ=255)
                                                                                                     cccccGCACACTCCATGACGATGCAGAAATAAAAGCCTGATCAATCATCAGGCTGCGTAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACGGTTTTCCCGACACACTGCACCGCCAGTGTGAATAATTGTGCAATTACTGTTGTCGa                                          <  4:115455/149‑1 (MQ=255)
                                                                                                     cccccGCACACTCCATGACGATGCAGAAATAAAAGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACGGTTTTCCCGACACACTGCACCGCCAGTGTGAATAATTGTGCAATTACTGTTGTCGa                                          <  1:702523/149‑1 (MQ=255)
                                                                                                        ccGCACACTCCATGACGATGCAGAAATAAAAGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACGGTTTTCCCGACACACTGCACCGCCAGTGTGAATAATTGTGCAATTAc                                                   >  3:157218/1‑137 (MQ=255)
                                                                                                        ccGCACACTCCATGACGATGCAGAAATAAAAGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACGGTTTTCCCGACACACTGCACCGCCAGTGTGAATAATTGTGCAATTAc                                                   <  4:157218/137‑1 (MQ=255)
                                                                                                        ccGCACACTCCATGACGATGCAGAAATAAAAGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACGGTTTTCCCGACACACTGCACCGCCAGTGTGAATAATTGTGCAATTACTGTTGTCGAGCa                                       >  2:516377/1‑149 (MQ=255)
                                                                                                         cGCACACTCCATGACGATGCAGAAATAAAAGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACGGTTTTCCCGACACACTGCACCGCCAGTGTGAATAATTGTGCAATTACTGTTGTCGAGCaa                                      <  1:737790/149‑1 (MQ=255)
                                                                                                         cGCACACTCCATGACGATGCAGAAATAAAAGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACGGTTTTCCCGACACACTGCACCGCCAGTGTGAATAATTGTGCAATTACTGTTGTCGAGCaa                                      <  4:698420/149‑1 (MQ=255)
                                                                                                          gCACACTCCATGACGATGCAGAAATAAAAGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACGGTTTTCCCGACACACTGCACCGCCAGTGTGAATAATTGTGCAATTACTGTTGTCGAGCAag                                     <  2:55592/149‑1 (MQ=255)
                                                                                                                  cATGACGATGCAGAAATAAAAGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACGGTTTTCCCGACACACTGCACCGCCAGTGTGAATAATTGTGCAATTACTGTTGTCGAGCAAGTAGTCCGc                             >  1:702869/1‑149 (MQ=255)
                                                                                                                  cATGACGATGCAGAAATAAAAGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACGGTTTTCCCGACACACTGCACCGCCAGTGTGAATAATTGTGCAATTACTGTTGTCGAGCAAGTAGTCCGc                             >  4:521568/1‑149 (MQ=255)
                                                                                                                      aCGATGCAGAAATAAAAGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACGGTTTTCCCGACACACTGCACCGCCAGTGTGAATAATTGTGCAATTACTGTTGTCGAGCAAGTAGTCCGCAATa                         >  2:545163/1‑149 (MQ=255)
                                                                                                                       cGATGCAGAAATAAAAGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACGGTTTTCCCGACACACTGCACCGCCAGTGTGAATAATTGTGCAATTACTGTTGTCGAGCAAGTAGTCCGCAATaa                        <  1:746078/149‑1 (MQ=255)
                                                                                                                        gATGCAGAAATAAAAGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACGGTTTTCCCGACACACTGCACCGCCAGTGTGAATAATTGTGCAATTACTGTTGTCGAGCAAGTAGTCCGCAATAAc                       >  1:159462/1‑149 (MQ=255)
                                                                                                                        gATGCAGAAATAAAAGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACGGTTTTCCCGACACACTGCACCGCCAGTGTGAATAATTGTGCAATTACTGTTGTCGAGCAAGTAGTCCGCAATAAc                       >  4:108385/1‑149 (MQ=255)
                                                                                                                             aGAAATAAAAGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACGGTTTTCCCGACACACTGCACCGCCAGTGTGAATAATTGTGCAATTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAAt                  >  1:212528/1‑149 (MQ=255)
                                                                                                                               aaaTAAAAGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACGGTTTTCCCGACACACTGCACCGCCAGTGTGAATAATTGTGCAATTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCg                >  2:1376/1‑149 (MQ=255)
                                                                                                                                 aTAAAAGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACGGTTTTCCCGACACACTGCACCGCCAGTGTGAATAATTGTGCAATTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGtt              <  4:6689/149‑1 (MQ=255)
                                                                                                                                  tAAAAGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACGGTTTTCCCGACACACTGCACCGCCAGTGTGAATAATTGTGCAATTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGttt             <  4:409388/149‑1 (MQ=255)
                                                                                                                                  tAAAAGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACGGTTTTCCCGACACACTGCACCGCCAGTGTGAATAATTGTGCAATTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGttt             <  4:93021/149‑1 (MQ=255)
                                                                                                                                  tAAAAGCCTGATCAACCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACGGTTTTCCCGACACACTGCACCGCCAGTGTGAATAATTGTGCAATTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGttt             <  3:523318/149‑1 (MQ=255)
                                                                                                                                    aaaGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACGGTTTTCCCGACACACTGCACCGCCAGTGTGAATAATTGTGCAATTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTtgt           <  2:113709/149‑1 (MQ=255)
                                                                                                                                    aaaGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACGGTTTTCCCGACACACTGCACCGCCAGTGTGAATAATTGTGCAATTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTtgt           <  1:545163/149‑1 (MQ=255)
                                                                                                                                    aaaGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACGGTTTTCCCGACACACTGCACCGCCAGTGTGAATAATTGTGCAATTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTtgt           >  4:754353/1‑149 (MQ=255)
                                                                                                                                     caGCCTGATGATTGATCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACGGTTTTCCCGACACACTGCACCGCCAGTGTGAATAATTGTGCAATTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTtgtg          >  1:373231/2‑149 (MQ=255)
                                                                                                                                       gCCTGATGATTGATCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACGGTTTTCCCGACACACTGCACCGCCAGTGTGAATAATTGTGCAATTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTTGTGAc        >  3:437750/1‑149 (MQ=255)
                                                                                                                                          tGATCAATCATCAGGCTGCGCAGAATGGTCGCAGCGGCTTCCCCGACACAGGAACGGTTTTCCCGACACACTGCACCGCCAGTGTGAATAATTGTGCAATTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTTGTGACcac     >  4:48097/1‑149 (MQ=255)
                                                                                                                                           gATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACGGTTTTCCCGACACACTGCACCGCCAGTGTGAATAATTGTGCAATTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTTGTGACcaca    >  1:204429/1‑149 (MQ=255)
                                                                                                                                             tcaatcaTCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACGGTTTTCCCGACACACTGCACCGCCAGTGTGAATAATTGTGCAATTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTTGTGACCAcagt  <  4:697444/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
ACCTGACCATATTTCGTCGCATCACAGACCAAAACTCGCTGGCGGCTGGCACTGGCAATCGCCCGTTTCACCGTGACTTTATCTTCCGCTGGCGTAGAAATCCCCCGCACACTCCATGACGATGCAGAAATAAAAGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACGGTTTTCCCGACACACTGCACCGCCAGTGTGAATAATTGTGCAATTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTTGTGACCACAGT  >  NC_000913/2860594‑2860883

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑