Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,860,738 A→T 5.1% I166N (ATT→AAT)  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,7380AT 5.1% 271.3 / 5.1 78I166N (ATT→AAT) ygbIDeoR family putative transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (38/36);  new base T (2/2);  total (40/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.24e-01

GCCACCTGACCATATTTCGTCGCATCACAGACCAAAACTCGCTGGCGGCTGGCACTGGCAATCGCCCGTTTCACCGTGACTTTATCTTCCGCTGGCGTAGAAATCCCCCGCACACTCCATGACGATGCAGAAATAAAAGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACGGTTTTCCCGACACACTGCACCGCCAGTGTGAATAATTGTGCAATTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTTGTGACCACAGT  >  NC_000913/2860591‑2860883
                                                                                                                                                   |                                                                                                                                                 
gCCACCTGACCATATTTCGTCGCATCACAGACCAAAACTCGCTGGCGGGTGGCAGTGGCAATCGCCCGTTTCACCGTGACTTTAGCTTCCGCTGGCGTAGAAATCCCCCGCACACTCCATGACGATGCAGAAATAAAAGCCTGatcaat                                                                                                                                                  <  2:8355/149‑1 (MQ=255)
gCCACCTGACCATATTTCGTCGCATCACAGACCAAAACTCGCTGGCGGCTGGCACTGGCAATCGCCCGTTTCACCGTGACTTTATCTTCCGCTGGCGTAGAAATCCCCCGCACACTCCATGACGATGCAGAAATAAAAGCCTGatcaat                                                                                                                                                  <  4:275539/149‑1 (MQ=255)
gCCACCTGACCATATTTCGTCGCATCACAGACCAAAACTCGCTGGCGGCTGGCACTGGCAATCGCCCGTTTCACCGTGACTTTATCTTCCGCTGGCGTAGAAATCCCCCGCACACTCCATGACGATGCAGAAATAAAAGCCTGatcaat                                                                                                                                                  <  3:549375/149‑1 (MQ=255)
   aCCTGACCATATTTCGTCGCATCACAGACCAAAACTCGCTGGCGGCTGGCACTGGCAATCGCCCGTTTCACCGTGACTTTATCTTCCGCTGGCGTAGAAATCCCCCGCACACTCCATGACGATGCAGAAATAAAAGCCTGATCAatcat                                                                                                                                               <  4:593503/149‑1 (MQ=255)
   aCCTGACCATATTTCGTCGCATCACAGACCAAAACTCGCTGGCGGCTGGCACTGGCAATCGCCCGTTTCACCGTGACTTTATCTTCCGCTGGCGTAGAAATCCCCCGCACACTCCATGACGATGCAGAAATAAAAGCCTGATCAatcat                                                                                                                                               <  1:289571/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                                                                                                                                        >  3:510217/1‑149 (MQ=255)
          cATATTTCGTCGCATCACAGACCAAAACTCGCTGGCGGCTGGCACTGGCAATCGCCCGTTTCACCGTGACTTTATCTTCCGCTGGCGTAGAAATCCCCCGCACACTCCATGACGATGCAGAAATAAAAGCCTGATCAATCATCAGGCTg                                                                                                                                        >  1:192782/1‑149 (MQ=255)
              tttCGTCGCATCACAGACCAAAACTCGCTGGCGGCTGGCACTGGCAATCGCCCGTTTCACCGTGACTTTATCTTCCGCTGGCGTAGAAATCCCCCGCACACTCCATGACGATGCAGAAATAAAAGCCTGATCAATCATCAGGCTGCgca                                                                                                                                    >  4:482248/1‑149 (MQ=255)
              tttCGTCGCATCACAGACCAAAACTCGCTGGCGGCTGGCACTGGCAATCGCCCGTTTCACCGTGACTTTATCTTCCGCTGGCGTAGAAATCCCCCGCACACTCCATGACGATGCAGAAATAAAAGCCTGATCAATCATCAGGCTGCgca                                                                                                                                    >  3:227244/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:12024/1‑149 (MQ=255)
                              aCCAAAACTCGCTGGCGGCTGGCACTGGCAATCGCCCGTTTCACCGTGACTTTATCTTCCGCTGGCGTAGAAATCCCCCGCACACTCCATGACGATGCAGAAATAAAAGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGc                                                                                                                    >  1:55592/1‑149 (MQ=255)
                                cAAAACTCGCTGGCGGCTGGCACTGGCAATCGCCCGTTTCACCGTGACTTTATCTTCCGCTGGCGTAGAAATCCCCCGCACACTCCATGACGATGCAGAAATAAAAGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCtt                                                                                                                  <  4:178400/149‑1 (MQ=255)
                                cAAAACTCGCTGGCGGCTGGCACTGGCAATCGCCCGTTTCACCGTGACTTTATCTTCCGCTGGCGTAGAAATCCCCCGCACACTCCATGACGATGCAGAAATAAAAGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCtt                                                                                                                  >  2:608520/1‑149 (MQ=255)
                                 aaaaCTCGCTGGCGGCTGGCACTGGCAATCGCCCGTTTCACCGTGACTTTATCTTCCGCTGGCGTAGAAATCCCCCGCACACTCCATGACGATGCAGAAATAAAAGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTc                                                                                                                 <  2:207095/149‑1 (MQ=255)
                                 aaaaCTCGCTGGCGGCTGGCACTGGCAATCGCCCGTTTCACCGTGACTTTATCTTCCGCTGGCGTAGAAATCCCCCGCACACTCCATGACGATGCAGAAATAAAAGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTc                                                                                                                 <  2:215590/149‑1 (MQ=255)
                                  aaaCTCGCTGGCGGCTGGCACTGGCAATCGCCCGTTTCACCGTGACTTTATCTTCCGCTGGCGTAGAAATCCCCCGCACACTCCATGACGATGCAGAAATAAAAGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTcc                                                                                                                >  4:626994/1‑149 (MQ=255)
                                      tCGCTGGCGGCTGGCACTGGCAATCGCCCGTTTCACCGTGACTTTATCTTCCGCTGGCGTAGAAATCCCCCGCACACTCCATGACGATGCAGAAATAAAAGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGa                                                                                                            <  1:75985/149‑1 (MQ=255)
                                      tCGCTGGCGGCTGGCACTGGCAATCGCCCGTTTCACCGTGACTTTATCTTCCGCTGGCGTAGAAATCCCCCGCACACTCCATGACGATGCAGAAATAAAAGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGa                                                                                                            <  2:572960/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                                                                                                <  4:1824/149‑1 (MQ=255)
                                                  ggcactggcaATCGCCCGTTTCACCGTGACTTTATCTTCCGCTGGCGTAGAAATCCCCCGCACACTCCATGACGATGCAGAAATAAAAGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACGGt                                                                                                >  2:767315/1‑149 (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:698782/149‑1 (MQ=255)
                                                               gCCCGTTTCACCGTGACTTTATCTTCCGCTGGCGTAGAAATCCCCCGCACACTCCATGACGATGCAGAAATAAAAGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACGGTTTTCCCGacacac                                                                                   >  2:416170/1‑149 (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                                      <  4:698420/149‑1 (MQ=255)
                                                                                                            cGCACACTCCATGACGATGCAGAAATAAAAGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACGGTTTTCCCGACACACTGCACCGCCAGTGTGAATAATTGTGCAATTACTGTTGTCGAGCaa                                      <  1:737790/149‑1 (MQ=255)
                                                                                                             gCACACTCCATGACGATGCAGAAATAAAAGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACGGTTTTCCCGACACACTGCACCGCCAGTGTGAATAATTGTGCAATTACTGTTGTCGAGCAag                                     <  2:55592/149‑1 (MQ=255)
                                                                                                                     cATGACGATGCAGAAATAAAAGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACGGTTTTCCCGACACACTGCACCGCCAGTGTGAATAATTGTGCAATTACTGTTGTCGAGCAAGTAGTCCGc                             >  4:521568/1‑149 (MQ=255)
                                                                                                                     cATGACGATGCAGAAATAAAAGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACGGTTTTCCCGACACACTGCACCGCCAGTGTGAATAATTGTGCAATTACTGTTGTCGAGCAAGTAGTCCGc                             >  1:702869/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)
                                                                                                                                                   |                                                                                                                                                 
GCCACCTGACCATATTTCGTCGCATCACAGACCAAAACTCGCTGGCGGCTGGCACTGGCAATCGCCCGTTTCACCGTGACTTTATCTTCCGCTGGCGTAGAAATCCCCCGCACACTCCATGACGATGCAGAAATAAAAGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACGGTTTTCCCGACACACTGCACCGCCAGTGTGAATAATTGTGCAATTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTTGTGACCACAGT  >  NC_000913/2860591‑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: ‑