Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,906,124 Δ1 bp 5.1% coding (437/567 nt) mntP → putative Mn(2+) efflux pump, mntR‑regulated

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,906,1240T. 5.1% 264.5 / 7.3 78coding (437/567 nt)mntPputative Mn(2+) efflux pump, mntR‑regulated
Reads supporting (aligned to +/- strand):  ref base T (43/31);  new base . (2/2);  total (45/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.21e-01

AAGAGCCGCGCCGTCGACACGGTTTCTGGCTACTGGTAACCACCGCGATTGCCACCAGCCTGGATGCCATGGCTGTGGGTGTTGGTCTTGCTTTCCTGCAGGTCAACATTATCGCGACCGCATTGGCCATTGGTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGG  >  NC_000913/1905977‑1906272
                                                                                                                                                   |                                                                                                                                                    
gagagCCGCGCCGTCGACACGGTTTCTGGCTACTGGTAACCACCGCGATTGCCACCAGCCTGGATGCCATGGCTGTGGGTGTTGGTCTTGCTTTCCTGCAGGTCAACATTATCGCGACCGCATTGGCCATTGGTTGTGCAACCTTGAtt                                                                                                                                                     <  4:394765/148‑1 (MQ=255)
   agCCGCGCCGTCGACACGGTTTCTGGCTACTGGTAACCACCGCGATTGCCACCAGCCTGGATGCCATGGCTGTGGGTGTTGGTCTTGCTTTCCTGCAGGTCAACATTATCGCGACCGCATTGGCCATTGGTTGTGCAACCTTGATTATg                                                                                                                                                  <  1:419792/149‑1 (MQ=255)
       gcgcCGTCGACACGGTTTCTGGCTACTGGTAACCACCGCGATTGCCACCAGCCTGGATGCCATGGCTGTGGGTGTTGGTCTTGCTTTCCTGCAGGTCAACATTATCGCGACCGCATTGGCCATTGGTTGTGCAACCTTGATTATGTCaa                                                                                                                                              >  1:596805/1‑149 (MQ=255)
            gTCGACACGGTTTCTGGCTACTGGTAACCACCGCGATTGCCACCAGCCTGGATGCCATGGCTGTGGGTGTTGGTCTTGCTTTCCTGCAGGTCAACATTATCGCGACCGCATTGGCCATTGGTTGTGCAACCTTGATTATGTCAACATTa                                                                                                                                         >  1:126630/1‑149 (MQ=255)
            gTCGACACGGTTTCTGGCTACTGGTAACCACCGCGATTGCCACCAGCCTGGATGCCATGGCTGTGGGTGTTGGTCTTGCTTTCCTGCAGGTCAACATTATCGCGACCGCATTGGCCATTGGTTGTGCAACCTTGATTATGTCAACATTa                                                                                                                                         <  2:419674/149‑1 (MQ=255)
             tCGACACGGTTTCTGGCTACTGGTAACCACCGCGATTGCCACCAGCCTGGATGCCATGGCTGTGGGTGTTGGTCTTGCTTTCCTGCAGGTCAACATTATCGCGACCGCATTGGCCATTGGTTGTGCAACCTTGATTATGTCAACATTAg                                                                                                                                        <  3:565034/149‑1 (MQ=255)
              cGACACGGTTTCTGGCTACTGGTAACCACCGCGATTGCCACCAGCCTGGATGCCATGGCTGTGGGTGTTGGTCTTGCTTTCCTGCAGGTCAACATTATCGCGACCGCATTGGCCATTGGTTGTGCAACCTTGATTATGTCAACATTAgg                                                                                                                                       <  3:37341/149‑1 (MQ=255)
                   cGGTTTCTGGCTACTGGTAACCACCGCGATTGCCACCAGCCTGGATGCCATGGCTGTGGGTGTTGGTCTTGCTTTCCTGCAGGTCAACATTATCGCGACCGCATTGGCCATTGGTTGTGCAACCTTGATTATGTCAACATTAGGgatga                                                                                                                                  >  1:539088/1‑149 (MQ=255)
                    ggTTTCTGGCTACTGGTAACCACCGCGATTGCCACCAGCCTGGATGCCATGGCTGTGGGTGTTGGTCTTGCTTTCCTGCAGGTCAACATTATCGCGACCGCATTGGCCATTGGTTGTGCAACCTTGATTATGTCAACATTAGGgatgat                                                                                                                                 >  2:382264/1‑149 (MQ=255)
                    ggTTTCTGGCTACTGGTAACCACCGCGATTGCCACCAGCCTGGATGCCATGGCTGTGGGTGTTGGTCTTGCTTTCCTGCAGGTCAACATTATCGCGACCGCATTGGCCATTGGTTGTGCAACCTTGATATGTCAACATTAGGGatgatg                                                                                                                                >  1:556876/1‑149 (MQ=255)
                    ggTTTCTGGCTACTGGTAACCACCGCGATTGCCACCAGCCTGGATGCCATGGCTGTGGGTGTTGGTCTTGCTTTCCTGCAGGTCAACATTATCGCGACCGCATTGGCCATTGGTTGTGCAACCTTGATATGTCAACATTAGGGatgatg                                                                                                                                >  1:485294/1‑149 (MQ=255)
                     gTTTCTGGCTACTGGTAACCACCGCGATTGCCACCAGCCTGGATGCCATGGCTGTGGGTGTTGGTCTTGCTTTCCTGCAGGTCAACATTATCGCGACCGCATTGGCCATTGGTTGTGCAACCTTGATTATGTCAACATTAGGGatgatg                                                                                                                                >  2:350697/1‑149 (MQ=255)
                             cTACTGGTAACCACCGCGATTGCCACCAGCCTGGATGCCATGGCTGTGGGTGTTGGTCTTGCTTTCCTGCAGGTCAACATTATCGCGACCGCATTGGCCATTGGTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCg                                                                                                                        >  4:293973/1‑149 (MQ=255)
                              tACTGGTAACCACCGCGATTGCCACCAGCCTGGATGCCATGGCTGTGGGTGTTGGTCTTGCTTTCCTGCAGGTCAACATTATCGCGACCGCATTGGCCATTGGTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGc                                                                                                                       >  4:504874/1‑149 (MQ=255)
                              tACTGGTAACCACCGCGATTGCCACCAGCCTGGATGCCATGGCTGTGGGTGTTGGTCTTGCTTTCCTGCAGGTCAACATTATCGCGACCGCATTGGCCATTGGTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGc                                                                                                                       >  3:174498/1‑149 (MQ=255)
                               aCTGGTAACCACCGCGATTGCCACCAGCCTGGATGCCATGGCTGTGGGTGTTGGTCTTGCTTTCCTGCAGGTCAACATTATCGCGACCGCATTGGCCATTGGTTGTGCAACCTTGATATGTCAACATTAGGGATGATGGTTGGTCGCtt                                                                                                                     <  2:485294/149‑1 (MQ=255)
                               aCTGGTAACCACCGCGATTGCCACCAGCCTGGATGCCATGGCTGTGGGTGTTGGTCTTGCTTTCCTGCAGGTCAACATTATCGCGACCGCATTGGCCATTGGTTGTGCAACCTTGATATGTCAACATTAGGGATGATGGTTGGTCGCtt                                                                                                                     <  2:556876/149‑1 (MQ=255)
                                 tGGTAACCACCGCGATTGCCACCAGCCTGGATGCCATGGCTGTGGGTGTTGGTCTTGCTTTCCTGCAGGTCAACATTATCGCGACCGCATTGGCCATTGGTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCttt                                                                                                                    >  3:202374/1‑149 (MQ=255)
                                   gTAACCACCGCGATTGCCACCAGCCTGGATGCCATGGCTGTGGGTGTTGGTCTTGCTTTCCTGCAGGTCAACATTATCGCGACCGCATTGGCCATTGGTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTAt                                                                                                                  >  3:446931/1‑149 (MQ=255)
                                     aaCCACCGCGATTGCCACCAGCCTGGATGCCATGGCTGTGGGTGTTGGTCTTGCTTTCCTGCAGGTCAACATTATCGCGACCGCATTGGCCATTGGTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCg                                                                                                                >  3:475046/1‑149 (MQ=255)
                                       ccaccGCGATTGCCACCAGCCTGGATGCCATGGCTGTGGGTGTTGGTCTTGCTTTCCTGCAGGTCAACATTATCGCGACCGCATTGGCCATTGGTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGc                                                                                                              >  1:118546/1‑149 (MQ=255)
                                       ccaccGCGATTGCCACCAGCCTGGATGCCATGGCTGTGGGTGTTGGTCTTGCTTTCCTGCAGGTCAACATTATCGCGACCGCATTGGCCATTGGTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGc                                                                                                              >  3:213406/1‑149 (MQ=255)
                                         accGCGATTGCCACCAGCCTGGATGCCATGGCTGTGGGTGTTGGTCTTGCTTTCCTGCAGGTCAACATTATCGCGACCGCATTGGCCATTGGTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATc                                                                                                                 >  1:233502/1‑144 (MQ=255)
                                         accGCGATTGCCACCAGCCTGGATGCCATGGCTGTGGGTGTTGGTCTTGCTTTCCTGCAGGTCAACATTATCGCGACCGCATTGGCCATTGGTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATc                                                                                                                 <  2:233502/144‑1 (MQ=255)
                                         accGCGATTGCCACCAGCCTGGATGCCATGGCTGTGGGTGTTGGTCTTGCTTTCCTGCAGGTCAACATTATCGCGACCGCATTGGCCATTGGTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTc                                                                                                            >  1:188393/1‑149 (MQ=255)
                                           cgcgATTGCCACCAGCCTGGATGCCATGGCTGTGGGTGTTGGTCTTGCTTTCCTGCAGGTCAACATTATCGCGACCGCATTGGCCATTGGTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCaa                                                                                                          >  3:248473/1‑149 (MQ=255)
                                                ttGCCACCAGCCTGGATGCCATGGCTGTGGGTGTTGGTCTTGCTTTCCTGCAGGTCAACATTATCGCGACCGCATTGGCCATTGGTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAattatt                                                                                                     <  1:382264/149‑1 (MQ=255)
                                                    caccaGCCTGGATGCCATGGCTGTGGGTGTTGGTCTTGCTTTCCTGCAGGTCAACATTATCGCGACCGCATTGGCCATTGGTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGa                                                                                                 >  2:163670/1‑149 (MQ=255)
                                                      ccaGCCTGGATGCCATGGCTGTGGGTGTTGGTCTTGCTTTCCTGCAGGTCAACATTATCGCGACCGCATTGGCCATTGGTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAatt                                                                                                        >  4:560637/1‑140 (MQ=255)
                                                      ccaGCCTGGATGCCATGGCTGTGGGTGTTGGTCTTGCTTTCCTGCAGGTCAACATTATCGCGACCGCATTGGCCATTGGTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAatt                                                                                                        <  3:560637/140‑1 (MQ=255)
                                                       caGCCTGGATGCCATGGCTGTGGGTGTTGGTCTTGCTTTCCTGCAGGTCAACATTATCGCGACCGCATTGGCCATTGGTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGaaaa                                                                                              >  1:583033/1‑149 (MQ=255)
                                                       caGCCTGGATGCCATGGCTGTGGGTGTTGGTCTTGCTTTCCTGCAGGTCAACATTATCGCGACCGCATTGGCCATTGGTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGaaaa                                                                                              >  1:417056/1‑149 (MQ=255)
                                                         gCCTGGATGCCATGGCTGTGGGTGTTGGTCTTGCTTTCCTGCAGGTCAACATTATCGCGACCGCATTGGCCATTGGTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGaaaaaa                                                                                            <  2:126630/149‑1 (MQ=255)
                                                                    aTGGCTGTGGGTGTTGGTCTTGCTTTCCTGCAGGTCAACATTATCGCGACCGCATTGGCCATTGGTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATtct                                                                                 >  4:182627/1‑149 (MQ=255)
                                                                    aTGGCTGTGGGTGTTGGTCTTGCTTTCCTGCAGGTCAACATTATCGCGACCGCATTGGCCATTGGTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATtct                                                                                 >  4:242322/1‑149 (MQ=255)
                                                                        cTGTGGGTGTTGGTCTTGCTTTCCTGCAGGTCAACATTATCGCGACCGCATTGGCCATTGGTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTcggc                                                                             <  4:213406/149‑1 (MQ=255)
                                                                        cTGTGGGTGTTGGTCTTGCTTTCCTGCAGGTCAACATTATCGCGACCGCATTGGCCATTGGTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTcggc                                                                             <  4:319365/149‑1 (MQ=255)
                                                                              gtgtTGGTCTTGCTTTCCTGCAGGTCAACATTATCGCGACCGCATTGGCCATTGGTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTg                                                                       >  1:537834/1‑149 (MQ=255)
                                                                                 ttGGTCTTGCTTTCCTGCAGGTCAACATTATCGCGACCGCATTGGCCATTGGTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTg                                                                    <  1:137521/149‑1 (MQ=255)
                                                                                     tCTTGCTTTCCTGCAGGTCAACATTATCGCGACCGCATTGGCCATTGGTTGTGCAACCTTGATTATGTCAACAt                                                                                                                                           >  1:91467/1‑74 (MQ=255)
                                                                                     tCTTGCTTTCCTGCAGGTCAACATTATCGCGACCGCATTGGCCATTGGTTGTGCAACCTTGATTATGTCAACAt                                                                                                                                           <  2:91467/74‑1 (MQ=255)
                                                                                     tCTTGCTTTCCTGCAGGTCAACATTATCGCGACCGCATTGGCCATTGGTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGa                                                                >  2:116226/1‑149 (MQ=255)
                                                                                       ttGCTTTCCTGCAGGTCAACATTATCGCGACCGCATTGGCCATTGGTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGatc                                                              >  3:101067/1‑149 (MQ=255)
                                                                                             tCCTGCAGGTCAACATTATCGCGACCGCATTGGCCATTGGTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGatcggcatc                                                        <  4:42905/149‑1 (MQ=255)
                                                                                             tCCTGCAGGTCAACATTATCGCGACCGCATTGGCCATTGGTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGatcggcatc                                                        <  4:291590/149‑1 (MQ=255)
                                                                                             tCCTGCAGGTCAACATTATCGCGACCGCATTGGCCATTGGTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGatcggcatc                                                        <  4:475046/149‑1 (MQ=255)
                                                                                                tGCAGGTCAACATTATCGCGACCGCATTGGCCATTGGTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGatcggcatcggc                                                     <  4:367255/149‑1 (MQ=255)
                                                                                                tGCAGGTCAACATTATCGCGACCGCATTGGCCATTGGTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGatcggcatcggc                                                     >  2:247868/1‑149 (MQ=255)
                                                                                                   aGGTCAACATTATCGCGACCGCATTGGCCATTGGTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTc                                                  >  2:321821/1‑149 (MQ=255)
                                                                                                     gTCAACATTATCGCGACCGCATTGGCCATTGGTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCa                                                >  1:526109/1‑149 (MQ=255)
                                                                                                      tCAACATTATCGCGACCGCATTGGCCATTGGTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAg                                               >  1:161605/1‑149 (MQ=255)
                                                                                                        aaCATTATCGCGACCGCATTGGCCATTGGTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGAt                                             >  2:23498/1‑149 (MQ=255)
                                                                                                        aaCATTATCGCGACCGCATTGGCCATTGGTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTAAATTATTGGGAAAAAAGCGGAAATTCTCGGCTGGCTGGTGCTGATCGGCATCGGCGTCCAGAt                                             <  1:350697/149‑1 (MQ=255)
                                                                                                         aCATTATCGCGACCGCATTGGCCATTGGTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATc                                            <  2:77250/149‑1 (MQ=255)
                                                                                                           aTTATCGCGACCGCATTGGCCATTGGTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCct                                          <  1:5019/149‑1 (MQ=255)
                                                                                                            ttATCGCGACCGCATTGGCCATTGGTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCctc                                         <  2:161605/149‑1 (MQ=255)
                                                                                                            ttATCGCGACCGCATTGGCCATTGGTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCctc                                         >  1:396316/1‑149 (MQ=255)
                                                                                                            ttATCGCGACCGCATTGGCCATTGGTTGTGCAACCTTGATTATGACAACATTAGGGATGATGGTTGGTAGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCctc                                         <  3:501003/149‑1 (MQ=255)
                                                                                                                cgcgACCGCATTGGCCATTGGTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGa                                     <  4:446931/149‑1 (MQ=255)
                                                                                                                   gACCGCATTGGCCATTGGTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGAGGATCTGGACGc                                  <  4:101067/149‑1 (MQ=255)
                                                                                                                          ttGGCCATTTGTTGTGCAACATTGCTTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGACTCAATTATTTGGAAAAAAGCGGAAATTATAGGCGGGCTTGTGCTGATCGGCATCGGCGTCCAGATCATCTGGACGAACTTCCa                           >  4:301783/1‑149 (MQ=255)
                                                                                                                          ttGGCCATTGGTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCa                           <  1:116226/149‑1 (MQ=255)
                                                                                                                             gCCATTGGTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACgg                        >  1:542383/1‑149 (MQ=255)
                                                                                                                              ccATTGGTTGTGCAACCTTGATTATGTCAACATTAGGGatgatg                                                                                                                                >  4:570226/1‑44 (MQ=255)
                                                                                                                              ccATTGGTTGTGCAACCTTGATTATGTCAACATTAGGGatgatg                                                                                                                                <  3:570226/44‑1 (MQ=255)
                                                                                                                              ccATTGGTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGt                       <  4:174498/149‑1 (MQ=255)
                                                                                                                              ccATTGGTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGt                       <  2:396316/149‑1 (MQ=255)
                                                                                                                                 ttGGTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTtaa                    >  3:375647/1‑149 (MQ=255)
                                                                                                                                   ggTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTtaata                  >  1:225770/1‑149 (MQ=255)
                                                                                                                                    gTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTtaataa                 <  4:202374/149‑1 (MQ=255)
                                                                                                                                    gTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTtaataa                 >  2:103060/1‑149 (MQ=255)
                                                                                                                                    gTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCCCAATTATTGGGAAAAAAGCGGAAATtctc                                                                                <  4:29070/86‑1 (MQ=255)
                                                                                                                                    gTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCCCAATTATTGGGAAAAAAGCGGAAATtctc                                                                                >  3:29070/1‑86 (MQ=255)
                                                                                                                                          cAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGc           >  2:298683/1‑149 (MQ=255)
                                                                                                                                           aaCCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCt          >  2:43644/1‑149 (MQ=255)
                                                                                                                                           aaCCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCt          >  1:280706/1‑149 (MQ=255)
                                                                                                                                            aCCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTg         >  1:589121/1‑149 (MQ=255)
                                                                                                                                                 gATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCaga    <  3:293973/149‑1 (MQ=255)
                                                                                                                                                  aTTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCagag   <  1:43644/149‑1 (MQ=255)
                                                                                                                                                   ttATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAgg  >  2:518573/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
AAGAGCCGCGCCGTCGACACGGTTTCTGGCTACTGGTAACCACCGCGATTGCCACCAGCCTGGATGCCATGGCTGTGGGTGTTGGTCTTGCTTTCCTGCAGGTCAACATTATCGCGACCGCATTGGCCATTGGTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGG  >  NC_000913/1905977‑1906272

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