Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,541,840 A→T 6.6% Q186L (CAA→CTA) ‡ nhaB → Na(+)/H(+) antiporter NhaB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,541,8400AT 6.6% 180.3 / 2.7 61Q186L (CAA→CTA) ‡nhaBNa(+)/H(+) antiporter NhaB
Reads supporting (aligned to +/- strand):  ref base A (28/29);  new base T (2/2);  total (30/31)
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.37e-01

TCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATGGCGCTGTGGGTACCGCCCTGGGCGGCGTTTGCACGTTGGTGGGCGAGCCGCAGAACCTGCTGATCGGCCACGAGATGGG  >  NC_002947/4541698‑4541987
                                                                                                                                              |                                                                                                                                                   
tCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGc                                                                                                                                               >  6:189748/1‑149 (MQ=255)
   gCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCa                                                                                                                                            >  1:156839/1‑149 (MQ=255)
            tCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCg                                                                                                                                   <  4:225110/149‑1 (MQ=255)
            tCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCg                                                                                                                                   >  1:63795/1‑149 (MQ=255)
              cgcgTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCATAGCTGCGGTCGGCTTCTACGCGGTGTACCAGCTCGTCGCTTCAGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCgg                                                                                                                                  <  6:292526/148‑1 (MQ=255)
              cgcgTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGTGAGCGACTTGCTGGTCGCTCTccagtgc                                                                                                                                     <  8:98751/145‑6 (MQ=255)
              cgcgTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGTGAGCGACTTGCTGGTCGCTCTccagtgc                                                                                                                                     >  7:98751/1‑140 (MQ=255)
               gcgTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACTTGCTGGTCGCTctc                                                                                                                                           <  1:149115/138‑3 (MQ=255)
               gcgTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACTTGCTGGTCGCTctc                                                                                                                                           >  2:149115/1‑136 (MQ=255)
                cgTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGaa                                                                                                                               <  2:156839/149‑1 (MQ=255)
                cgTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGaa                                                                                                                               <  3:91740/149‑1 (MQ=255)
                  tttCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAg                                                                                                                              >  5:325744/1‑148 (MQ=255)
                  tttCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAg                                                                                                                              <  6:325744/148‑1 (MQ=255)
                     ctcGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCt                                                                                                                          <  6:290108/149‑1 (MQ=255)
                     ctcGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCt                                                                                                                          >  4:231472/1‑149 (MQ=255)
                        gATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCg                                                                                                                        <  1:235068/148‑1 (MQ=255)
                                     tCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTc                                                                                                                                                 <  3:39970/110‑1 (MQ=255)
                                     tCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTc                                                                                                                                                 >  4:39970/1‑110 (MQ=255)
                                     tCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGcc                                                                                                          >  3:198279/1‑149 (MQ=255)
                                       aCGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCtt                                                                                                        >  6:291646/1‑149 (MQ=255)
                                             gTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCg                                                                                                  >  5:146923/1‑149 (MQ=255)
                                             gTGAGCCTCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGCGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCg                                                                                                  <  8:155035/149‑1 (MQ=255)
                                              tGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGt                                                                                                 <  6:146923/149‑1 (MQ=255)
                                                atcatcAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTa                                                                                                <  5:291646/148‑1 (MQ=255)
                                                   atcaGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCAACGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCct                                                                                            >  8:85625/1‑149 (MQ=255)
                                                    tcaGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCctg                                                                                           >  1:85942/1‑149 (MQ=255)
                                                      aGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCctgc                                                                                          <  6:320084/148‑1 (MQ=255)
                                                       gcgcTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCctgctg                                                                                        <  4:291796/149‑1 (MQ=255)
                                                                tCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATGgc                                                                               <  2:261254/149‑1 (MQ=255)
                                                                     ttCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATGGCGCtgt                                                                          >  5:270462/1‑149 (MQ=255)
                                                                     ttCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATGGCGCtgt                                                                          >  2:50556/1‑149 (MQ=255)
                                                                       cgaCGCGGTGGTACACCGCGTAGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATGGCGCTGTgg                                                                        >  4:148501/3‑149 (MQ=255)
                                                                       cTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATGGCGCtgt                                                                          <  3:66556/147‑1 (MQ=255)
                                                                       cTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATGGCGCtgt                                                                          <  5:49291/147‑1 (MQ=255)
                                                                        tACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCCTGGCGCTGTggg                                                                       >  8:38500/1‑149 (MQ=255)
                                                                          cgcgGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATGGCGCTGTgg                                                                        <  4:198279/146‑1 (MQ=255)
                                                                          cgcgGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATGGCGCTGTGGGTa                                                                     >  3:51785/1‑149 (MQ=255)
                                                                                   ccaccGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTTGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGAGCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATGGCGCTGTGGGTACCGCCCTgg                                                            >  1:40074/1‑149 (MQ=255)
                                                                                   ccaccGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATGGCGCTGTGGGTACCGCCCTgg                                                            <  7:5428/149‑1 (MQ=255)
                                                                                    caccGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATGGCGCTGTGGGTACCGCCCTggg                                                           >  1:56798/1‑149 (MQ=255)
                                                                                     accGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATGGCGCTGTGGGTACCGCCCTGggc                                                          >  2:93287/1‑149 (MQ=255)
                                                                                      ccGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTCGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCAGCTCCTTCCTGCGTAGCCTGCTGATGCATGGCGCTGTGGGTACCGCCCTGggcg                                                         >  3:266810/1‑149 (MQ=255)
                                                                                         cgtcgcgtcCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATGGCGCTGTGGGTACCGCCCTGGgcggcg                                                      <  1:93287/149‑1 (MQ=255)
                                                                                          gtcgcgtcCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATGGCGCTGTGGGTACCGCCCTGGGCGGCGt                                                     <  2:63795/149‑1 (MQ=255)
                                                                                             gcgtcCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATGGCGCTGTGGGTACCGCCCTGGGCGGCGttt                                                   <  2:56798/148‑1 (MQ=255)
                                                                                                      gCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGTAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATGGCGCTGTGGGTACCGCCCTGGGCGGCGTTTGCACGTtggt                                         <  1:205299/149‑1 (MQ=255)
                                                                                                        cAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATGGCGCTGTGGGTACCGCCCTGGGCGGCGTTTGCACGTtggtgg                                       <  5:37834/149‑1 (MQ=255)
                                                                                                              gcgcgAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAACACCTCGACCAGTTCCGCGCCTTACTCAGTAGCCTGCTGATGCATGGAGCTGTGGGTACCGCCCTGCGCGGCGTTTGCACGTTGGTCGGCGAGc                                 >  2:321795/1‑149 (MQ=255)
                                                                                                                      aCAGCGCACTGGAGAGCGACCAGCAAGTCTCTCACCTGCACGGGGAAGACCTCGACCAGTTCCGCGCCTTGGTGCGTAGCCTGGTGATTCATGGCGCTGTGGGTACCGCCCTGGGCGGGGTTTGCACGTTGGTGGTCGAGCCGCAGAAc                         <  2:40074/149‑1 (MQ=255)
                                                                                                                             aCTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTa                                                                                                >  8:77285/1‑71 (MQ=255)
                                                                                                                             aCTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTa                                                                                                <  7:77285/71‑1 (MQ=255)
                                                                                                                               tGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGTCCTCGAGCAGTTCCGCGCCTTGCTGGGTAGCGTGCTGATGGATGGCGCTGTGGGTACGGCCCTGGGCGGCGTTTGCACGTTGGTGGGCGAGCCGCAGAACCTGCTGATc                <  1:321795/149‑1 (MQ=255)
                                                                                                                                 gagaGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATGGCGCTGTGGGTACCGCCCTGGGCGGCGTTTGCACGTTGGTGGGCGAGCCGCAGAACctgc                     >  4:14813/1‑142 (MQ=255)
                                                                                                                                 gagaGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATGGCGCTGTGGGTACCGCCCTGGGCGGCGTTTGCACGTTGGTGGGCGAGCCGCAGAACctgc                     <  3:14813/142‑1 (MQ=255)
                                                                                                                                 gagaGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATGGCGCTGTGGGTACCGCCCTGGGCGGCGTTTGCACGTTGGTGGGCGAGCCGCAGAACCTGCTGATCgg              >  6:51964/1‑149 (MQ=255)
                                                                                                                                   gagCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATGGCGCTGTGGGTACCGCCCTGGGCGGCGTTTGCACGTTGGTGGGCGAGCCGCAGAACCTGCTGATCGGc             >  2:47923/1‑148 (MQ=255)
                                                                                                                                        aCCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATGGCGCTGTGGGTACCGCCCTGGGCGGCGTTTGCACGTTGGTGGGCGAGCCGCAGAACCTGCTGATCGGCCa           <  6:46858/145‑1 (MQ=255)
                                                                                                                                        aCCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATGGCGCTGTGGGTACCGCCCTGGGCGGCGTTTGCACGTTGGTGGGCGAGCCGCAGAACCTGCTGATCGGCCa           >  5:46858/1‑145 (MQ=255)
                                                                                                                                        aCCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATGGCGCTGTGGGTACCGCCCTGGGCGGCGTTTGCACGTTGGTGGGCGAGCCGCAGAACCTGCTGATCGGCCAc          >  7:74627/1‑146 (MQ=255)
                                                                                                                                            gCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATGGCGCTGTGGGTACCGCCCTGGGCGGCGTTTGCACGTTGGTGGGCGAGCCGCAGAACCTGCTGATCGGCCACGAGATgg   >  6:132984/1‑149 (MQ=255)
                                                                                                                                             cAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATGGCGCTGTGGGTACCGCCCTGGGCGGCGTTTGCACGTTGGTGGGCGAGCCGCAGAACCTGCTGATCGGCCACGAGATggg  <  6:270462/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                   
TCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATGGCGCTGTGGGTACCGCCCTGGGCGGCGTTTGCACGTTGGTGGGCGAGCCGCAGAACCTGCTGATCGGCCACGAGATGGG  >  NC_002947/4541698‑4541987

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: ‑