Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,541,841 A→G 5.1% Q186Q (CAA→CAG) ‡ nhaB → Na(+)/H(+) antiporter NhaB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,541,8410AG 5.1% 235.4 / 2.5 78Q186Q (CAA→CAG) ‡nhaBNa(+)/H(+) antiporter NhaB
Reads supporting (aligned to +/- strand):  ref base A (36/38);  new base G (2/2);  total (38/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.01e-01

TCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATGGCGCTGTGGGTACCGCCCTGGGCGGCGTTTGCACGTTGGTGGGCGAGCCGCAGAACCTGCTGATCGGCCACGAGATGG  >  NC_002947/4541698‑4541986
                                                                                                                                                |                                                                                                                                                 
tCGGCGTTTCTATCCGCGTTTCTCGATGCCCTAACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGc                                                                                                                                              >  1:232703/1‑149 (MQ=255)
 cGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACTTGCTGGTCGCt                                                                                                                                             >  2:353390/1‑149 (MQ=255)
   gCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCa                                                                                                                                           >  4:165/1‑149 (MQ=255)
    cGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACTTGCTGGTCGCTctc                                                                                                                                          <  1:353390/149‑3 (MQ=255)
       ttCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTg                                                                                                                                       >  8:305552/1‑149 (MQ=255)
         cTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGTGAGCGACTTGCTGGTCGCTCTccagt                                                                                                                                      <  3:181641/148‑4 (MQ=255)
          ttgtggcgTTTCTGGATGCCCTGACCGTCACGGCGGTGATCTTCAGGGCTGGGGTCGGCTTATACGGGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCAc                                                                                                                                    <  2:178908/144‑1 (MQ=255)
           aTCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCaa                                                                                                                                                   <  5:150163/133‑1 (MQ=255)
           aTCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCaa                                                                                                                                                   >  6:150163/1‑133 (MQ=255)
            tCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCg                                                                                                                                  >  7:349093/1‑149 (MQ=255)
            tCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCg                                                                                                                                  >  2:64599/1‑149 (MQ=255)
              cgcgTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCggg                                                                                                                                <  1:297563/149‑1 (MQ=255)
              cgcgTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCgg                                                                                                                                 >  2:297563/1‑148 (MQ=255)
               gcgTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCAc                                                                                                                                    >  1:178908/1‑144 (MQ=255)
               gcgTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCggg                                                                                                                                >  6:29446/1‑148 (MQ=255)
                cgTTTCTCGATGCGCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACGACCGCGTCGCGTCCGGGGCGAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGaa                                                                                                                              <  2:342806/149‑1 (MQ=255)
                cgTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGaa                                                                                                                              <  1:113193/149‑1 (MQ=255)
                  tttCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGa                                                                                                                            >  2:301366/1‑149 (MQ=255)
                       cGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCt                                                                                                                                        <  4:69192/132‑1 (MQ=255)
                       cGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGAGACCAGCAAGTCGCTCACCt                                                                                                                                        >  3:69192/1‑132 (MQ=255)
                        gATGCCCTGACCGTCACGGCGGTGATCATCTGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGa                                                                                                                      >  1:34104/1‑149 (MQ=255)
                        gATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGaa                                                                                                                              <  8:134491/141‑1 (MQ=255)
                        gATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGaa                                                                                                                              >  7:134491/1‑141 (MQ=255)
                        gATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTc                                                                                                                        <  4:350317/147‑1 (MQ=255)
                        gATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGa                                                                                                                      <  8:355029/149‑1 (MQ=255)
                             ccTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGt                                                                                                                 >  7:334425/1‑149 (MQ=255)
                                  ccGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCcgc                                                                                                            <  2:34104/149‑1 (MQ=255)
                                  ccGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCcgc                                                                                                            <  7:305552/149‑1 (MQ=255)
                                    gTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGGGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCgcgc                                                                                                          <  1:66708/149‑1 (MQ=255)
                                         ggcggTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTcc                                                                                                     >  8:220997/1‑149 (MQ=255)
                                                   atcaGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCAGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAACTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTCCGTAGCct                                                                                           >  4:18032/1‑149 (MQ=255)
                                                           tGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCTACAAGTTCCGCGCCTTCATGCGTAGCCTGCTGATg                                                                                    >  6:205740/1‑148 (MQ=255)
                                                           tGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATg                                                                                    >  6:302224/1‑148 (MQ=255)
                                                             cGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGc                                                                                   <  2:82417/147‑1 (MQ=255)
                                                              ggTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATg                                                                                <  1:159129/149‑1 (MQ=255)
                                                               gTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATgg                                                                               <  8:334425/149‑1 (MQ=255)
                                                                 cGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATGgcg                                                                             >  5:218956/1‑149 (MQ=255)
                                                                    cTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATGGCGCtg                                                                          >  2:204512/1‑149 (MQ=255)
                                                                    cTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATGGCGCtg                                                                          >  2:305881/1‑149 (MQ=255)
                                                                     ttCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGTCGACCAGCAAGTCGCTCACCTGCCCCGGGAAGCCCTAGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATGgcgcgg                                                                          >  1:296275/1‑147 (MQ=255)
                                                                       cTACGCGGTGTACCACCGGGTCGCGTCGGGGGCCAATCGGCGCTAGGAAAGCGCACTGTATAGGGACCAGCAAGTCGCTTACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATGGCGCtgtg                                                                        <  1:301366/148‑1 (MQ=255)
                                                                        tACCCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATgg                                                                               <  8:39514/140‑1 (MQ=255)
                                                                        gaCGCGGTGGTACACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCctgct                                                                                        >  1:15133/2‑131 (MQ=255)
                                                                        gaCGCGGTGGTACACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCctgct                                                                                        <  2:15133/130‑1 (MQ=255)
                                                                        gaCGCGGTGGTACACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATGGCGCTGTggg                                                                      >  5:371266/2‑149 (MQ=255)
                                                                               tgtACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATGGCGCTGTGGGTACCGcc                                                               <  6:218956/149‑1 (MQ=255)
                                                                                  accaccGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATGGCGCTGTGGGTACCGCCCTg                                                            >  7:293045/1‑149 (MQ=255)
                                                                                   ccaccGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATGGCGCTGTGGGTACCGCCCTgg                                                           <  8:334834/149‑1 (MQ=255)
                                                                                   ccaccGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATGGCGCTGTGGGTACCGCCCTgg                                                           <  5:29446/149‑1 (MQ=255)
                                                                                     accGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATGGCGCTGTGGGTACCGCCCTGggc                                                         <  1:204512/149‑1 (MQ=255)
                                                                                     accGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACATCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATGGCGCTGTGGGTACCGCCCTGggc                                                         <  5:205740/149‑1 (MQ=255)
                                                                                         cgtcgcgtcCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATGGCGCTGTGGGTACCGCCCTGGgcggcg                                                     <  1:64599/149‑1 (MQ=255)
                                                                                           tcgcgtcCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCctg                                                                                          >  4:9111/1‑110 (MQ=255)
                                                                                           tcgcgtcCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCctg                                                                                          <  3:9111/110‑1 (MQ=255)
                                                                                             gcgtcCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATGGCGCTGTGGGTACCGCCCTGGGCGGCGttt                                                  >  8:126802/1‑148 (MQ=255)
                                                                                             gcgtcCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATGGCGCTGTGGGTACCGCCCTGGGCGGCGttt                                                  <  7:126802/148‑1 (MQ=255)
                                                                                             gcggcCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATGGCGCTGTGGGTACCGCCCTGGGCGGCGttt                                                  <  2:232703/148‑1 (MQ=255)
                                                                                              cgtcCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATGGCGCTGTGGGTACCGCCCTGGGCGGCGTTTg                                                 >  2:168836/1‑148 (MQ=255)
                                                                                                tcCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATGGCGCTGTGGGTACCGCCCTGGGCGGCGTTTGCAc                                              <  8:293045/149‑1 (MQ=255)
                                                                                                    gggCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATGGCGCTGTGGGTACCGCCCTGGGCGGCGTTTGCACGTtg                                          <  1:168836/149‑1 (MQ=255)
                                                                                                    gggCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATGGCGCTGTGGGTACCGCCCTGGGCGGCGTTTGCACGTtg                                          <  3:18032/149‑1 (MQ=255)
                                                                                                        cAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATGGCGCTGTGGGTACCGCCCTGGGCGGCGTTTGCACGTtggtgg                                      >  7:60322/1‑149 (MQ=255)
                                                                                                        cAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATGGCGCTGTGGGTACCGCCCTGGGCGGCGTTTGCACGTtggtgg                                      <  8:60322/149‑1 (MQ=255)
                                                                                                        cAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATGGCGCTGTGGGTACCGCCCTGGGCGGCGTTTGCACGTtggtgg                                      <  6:371266/149‑1 (MQ=255)
                                                                                                             cgcgcgAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATGGCGCTGTGGGTACCGCCCTGGGCGGCGTTTGCACGTTGGTGGGCGa                                  <  7:220997/148‑1 (MQ=255)
                                                                                                                         gcgcACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATGGCGCTGTGGGTACCGCCCTGGGCGGCGTTTGCACGTTGGTGGGCGAGCCGCAGAACctg                     <  1:305881/149‑1 (MQ=255)
                                                                                                                         gcgcACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATGGCGCTGTGGGTACCGCCCTGGGCGGCGTTTGCACGTTGGTGGGCGAGCCGCAGAACctg                     >  2:64258/1‑149 (MQ=255)
                                                                                                                          cgcACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATGGCGCTGTGGGTACCGCCCTGGGCGGCGTTTGCACGTTGGTGGGCGAGCCGCAGAACctgc                    <  2:296275/149‑1 (MQ=255)
                                                                                                                          cgcACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATGGCGCTGTGGGTACCGCCCTGGGCGGCGTTTGCACGTTGGTGGGCGAGCCGCAGAACctgc                    <  5:302224/149‑1 (MQ=255)
                                                                                                                            cACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCAGAAGCATGGCGCTGTCGGTACCGCCCTGGGCGGCGTTTGCACGTTGGTGGGCGAGCCGCAGACCctgct                   >  6:199641/1‑148 (MQ=255)
                                                                                                                                 gagaGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATGGCGCTGTGGGTACCGCCCTGGGCGGCGTTTGCACGTTGGTGGGCGAGCCGCAGAACCTGCTGATCgg             >  3:282817/1‑149 (MQ=255)
                                                                                                                                   gagCGACTTGCTGGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATGGCGCTGTGGGTACCGCCCTGGGCGGCGTTTGCACGTTGGTGGGCGAGCCGCAGAACCTGCTGATCGGc            >  2:196476/1‑148 (MQ=255)
                                                                                                                                   gagCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATGGCGCTGTGGGTACCGCCCTGGGCGGCGTTTGCACGTTGGTGGGCGAGCCGCAGAACCTGCTGATCGGcc           <  3:165/149‑1 (MQ=255)
                                                                                                                                    agCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATGGCGCTGTGGGTACCGCCCTGGGCGGCGttt                                                  <  8:13822/109‑1 (MQ=255)
                                                                                                                                    agCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATGGCGCTGTGGGTACCGCCCTGGGCGGCGttt                                                  >  7:13822/1‑109 (MQ=255)
                                                                                                                                        gACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATGGCGCTGTGGGTACCGCCCTGGGCGGCGTTTGCACGTTGGTGGGCGAGCCGCAGAACCTGCTGATCGGCCAc         >  5:210139/1‑147 (MQ=255)
                                                                                                                                          ccAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATGGCGCTGTGGGTACCGCCCTGGGCGGCGTTTGCACGTTGGTGGGCGAGCCGCAGAACCTGCTGATCGGCCAc         >  2:197652/1‑145 (MQ=255)
                                                                                                                                           cAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATGGCGCTGTGGGTACCGCCCTGggcggc                                                      >  1:45877/1‑99 (MQ=255)
                                                                                                                                           cAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATGGCGCTGTGGGTACCGCCCTGggcggc                                                      <  2:45877/99‑1 (MQ=255)
                                                                                                                                             gCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATGGCGCTGTGGGTACCGCCCTGGGCGGCGTTTGCACGTTGGTGGGCGAGCCGCAGAACCTGCTGATCGGCCACGAGATgg  >  6:13078/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
TCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATGGCGCTGTGGGTACCGCCCTGGGCGGCGTTTGCACGTTGGTGGGCGAGCCGCAGAACCTGCTGATCGGCCACGAGATGG  >  NC_002947/4541698‑4541986

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