Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,539,134 G→C 6.1% intergenic (‑96/+55) mrpAB ← / ← PP_2231 K(+)/H(+) antiporter subunit A/B/DMT superfamily permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,539,1340GC 6.1% 203.5 / 3.7 65intergenic (‑96/+55)mrpAB/PP_2231K(+)/H(+) antiporter subunit A/B/DMT superfamily permease
Reads supporting (aligned to +/- strand):  ref base G (35/26);  new base C (2/2);  total (37/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.97e-01

CCAGGCAACTGCCCACGAACGGCAGAAGCAATAGCACTATCAAGGACATAGCGGCCTATTCTGGAGAAGTTTGCAAAGAATCATACGTGCCGAGGTGTTTATCACCAACACGCAAGCTGTCGCAGGATCCTACAAACCGGCTGTGACAAGCGGTTTTTTTATAACAGTTTTTTACAAAGAATAGCCGCTGTTGACCGTTTTACAGCCCCGCCTGGCGTTCGCCCGGGTTCTGCACATCGCGTGTTTCGACAGTTTCGCGAGGGCGCCGCACCTTGAGTTCACTGAC  >  NC_002947/2538990‑2539275
                                                                                                                                                |                                                                                                                                              
ccAGGCAACTGCCCACGAACGGCAGAAGCAATAGCACTATCAAGGACATAGCGGCCTATTCTGGAGAAGTTTGCAAAGAATCATACGTGCCGAGGTGTTTATCACCAACACGCAAGCTGTCGCAGGATCCTACAAACCGGCTGTGAc                                                                                                                                              <  7:364728/147‑1 (MQ=255)
      aaCTGCCCACGAACGGCAGAAGCAATAGCACTATCAAGGACATAGCGGCCTATTCTGGAGAAGTTTGCAAAGAATCATACGTGCCGAGGTGTTTATCACCAACACGCAAGCTGTCGCAGGATCCTACAAACCGCTTGTCACAGCCGGt                                                                                                                                       <  6:152379/148‑1 (MQ=255)
      aaCTGCCCACGAACGGCAGAAGCAATAGCACTATCAAGGACATAGCGGCCTATTCTGGAGAAGTTTGCAAAGAATCATACGTGCCGAGGTGTTTATCACCAACACGCAAGCTGTCGCAGGATCCTACAAACCGCTTGTCACAGCCGGt                                                                                                                                       >  5:152379/1‑148 (MQ=255)
       aCTGCCCACGAACGGCAGAAGCAATAGCACTATCAAGGACATAGCGGCCTATTCTGGAGAAGTTTGCAAAGAATCATACGTGCCGAGGTGTTTATCACCAACACGCAAGCTGTCGCAGGATCCTACAAACCGGCTGTGACAAGCGGttt                                                                                                                                    <  7:354534/149‑1 (MQ=255)
       aCTGCCCACGAACGGCAGAAGCAATAGCACTATCAAGGACATAGCGGCCTATTCTGGAGAAGTTTGCAAAGAATCATACGTGCCGAGGTGTTTATCACCAACACGCAAGCTGTCGCAGGATCCTACAAACCGGCTGTGACAAGCGGttt                                                                                                                                    <  2:155331/149‑1 (MQ=255)
           cccACGAACGGCAGAAGCAATAGCACTATCAAGGACATAGCGGCCTATTCTGGAGAAGTTTGCAAAGAATCATACGTGCCGAGGTGTTTATCACCAACACGCAAGCTGTCGCAGGATCCTACAAACCGGCTGTGACAAGCGGtttttt                                                                                                                                 <  1:113639/148‑1 (MQ=255)
           cccACGAACGGCAGAAGCAATAGCACTATCAAGGACATAGCGGCCTATTCTGGAGAAGTTTGCAAAGAATCATACGTGCCGAGGTGTTTATCACAAACACGCAAGCTGTCGCAGGATCCTACAAAACGGCTGTGACAAGCGGttttttt                                                                                                                                >  4:33585/1‑149 (MQ=255)
              aCGAACGGCAGAAGCAATAGCACTATCAAGGACATAGCGGCCTATTCTGGAGAAGTTTGCAAAGAATCATACGTGCCGAGGTGTTTATCACCAACACGCAAGCTGTCGCAGGATCCTACAAACCGGCTGTGACAAGCGGTTTTTTtat                                                                                                                              >  4:101234/1‑148 (MQ=255)
              aCGAACGGCAGAAGCAATAGCACTATCAAGGACATAGCGGCCTATTCTGGAGAAGTTTGCAAAGAATCATACGTGCCGAGGTGTTTATCACCAACACGCAAGCTGTCGCAGGATCCTACAAACCGGCTGTGACAAGCGGTTTTTTtat                                                                                                                              >  6:411602/1‑148 (MQ=255)
                 aaCGGCAGAAGCAATAGCACTATCAAGGACATAGCGGCCTATTCTGGAGAAGTTTGCAAAGAATCATACGTGCCGAGGTGTTTATCACCAACACGCAAGCTGTCGCAGGATCCTACAAACCGGCTGTGACAAGCGGTTTTTTTATAACa                                                                                                                          <  3:101234/149‑1 (MQ=255)
                      cAGAAGCAATAGCACTATCAAGGACATAGCGGCCTATTCTGGAGAAGTTTGCAAAGAATCATACGTGCCGAGGTGTTTATCACCAACACGCAAGCTGTCGCAGGATCCTACAAACCGGCTGTGACAAGCGGTTTTTTTATAACAGtttt                                                                                                                     <  5:76625/149‑1 (MQ=255)
                      cAGAAGCAATAGCACTATCAAGGACATAGCGGCCTATTCTGGAGAAGTTTGCAAAGAATCATACGTGCCGAGGTGTTTATCACCAACACGCAAGCTGTCGCAGGATCCTACAAACCGGCTGTGACAAGCGGTTTTTTTATAACAGtttt                                                                                                                     >  1:341510/1‑149 (MQ=255)
                          aGCAATAGCACTATCAAGGAAATAGCGGCCTATTCTGGATAATTTTGCAAAGAATCATACGTGCCTAGTTGTTTATCACCAACACGCAAGCTGTCGCAGGATCCTACAAACCGGCTGTGACAAGCGGTTTTTTTATAACAGTTTTTTa                                                                                                                  <  1:218273/148‑1 (MQ=255)
                                        cAAGGACATAGCGGCCTATTCTGGAGAAGTTTGCAAAGAATCATACGTGCCGAGGTGTTTATCACCAACACGCAAGCTGTCGCAGGATCCTACAAACCGGCTGTGACAAGCGGTTTTTTTATAACAGTTTTTTACAAAGAATAGCCGCt                                                                                                   >  8:382978/1‑149 (MQ=255)
                                               aTAGCGGCCTATTCTGGAGAAGTTTGCAAAGAATCATACGTGCCGAGGTGTTTATCACCAACACGCAAGCTGTCGCAGGATCCTACAAACCGGCTGTGACAAGCGGTTTTTTTATAACAGTTTTTTACAAATAATAGCCGCTGTTGAcc                                                                                            <  2:341510/149‑1 (MQ=255)
                                                       cTATTCTGGAGAAGTTTGCAAAGAATCATACGTGCCGAGGTGTTTATCACCAACACGCAAGCTGTCGCAGGATCCTACAAACCGGCTGTGACAAGCGGTTTTTTTATAACAGTTTTTTACAAAGAATAGCCGCTGTTGACCGTTTTa                                                                                      >  8:117788/1‑147 (MQ=255)
                                                        tATTCTGGAGAAGTTTGCAAAGAATCATACGTGCCGAGGTGTTTATCACCAACACGCAAGCTGTCGCAGGATCCTACAAACCGGCTGTGACAAGCGGTTTTTTTATAACAGTTTTTTACAAAGAATAGCCGCTGTTGACCGTTTTACAg                                                                                   >  2:202655/1‑149 (MQ=255)
                                                        tATTCTGGAGAAGTTTGCAAAGAATCATACGTGCCGAGGTGTTTATCACCAACACGCAAGCTGTCGCAGGATCCTACAAACCGGCTGTGACAAGCGGTTTTTTTATAACAGTTTTTTACAAAGAATAGCCGCTGTTGACCGTTTTACAg                                                                                   >  6:298660/1‑149 (MQ=255)
                                                        tATTCTGGAGAAGTTTGCAAAGAATCATACGTGCCGAGGTGTTTATCACCAACACGCAAGCTGTCGCAGGATCCTACAAACCGGCTGTGACAAGCGGTTTTTTTATAACAGTTTTTTACAAAGAATAGCCGCTGTTGACCGTTTTACAg                                                                                   >  7:116621/1‑149 (MQ=255)
                                                          ttCTGGAGAAGTTTGCAAAGAATCATACGTGCCGAGGTGTTTATCACCAACACGCAAGCTGTCGCAGGATCCTACAAACCGGCTGTGACAAGCGGTTTTTTTATAACAGTTTTTTACAAAGAATAGCCGCTGTTGACCGTTTTACAGcc                                                                                 >  3:178817/1‑149 (MQ=255)
                                                             tGGAGAAGTTTGCAAAGAATCATACGTGCCGAGGTGTTTATCACCAACACGCAAGCTGTCGCAGGATCCTACAAACCGGCTGTGACAAGCGGTTTTTTTATAACAGTTTTTTACAAAGAATAGCCGCTGTTGACCGTTTTACAGCCCCg                                                                              >  5:347339/1‑149 (MQ=255)
                                                             tGGAGAAGTTTGCAAAGAATCATACGTGCCGAGGTGTTTATCACCAACACGCAAGCTGTCGCAGGATCCTACAAACCGGCTGTGACAAGCGGTTTTTTTATAACAGTTTTTTACAAAGAATAGCCGCTGTTGACCGTTTTACAGCCCCg                                                                              >  1:196763/1‑149 (MQ=255)
                                                             tGGAGAAGTTTCCAAAGAATCATACGTGCCGAGGTGTTTATCACCAACACGCAACCTGACGCAGGATCCTACAAACCCGCTGTGACAAGCGGATATATAATAACAGTTTTTTACAAAGAAAAGCCGCTGAAGACCCTTTTACAGcccag                                                                              >  5:367875/1‑147 (MQ=255)
                                                                          aaaGAATCATACGTGCCGAGGTGTTTATCACCAACACGCAAGCTGTCGCAGGATCCTACAAACCGGCTGTGACAAGCGGTTTTTTTATAACAGTTTTTTACAAAGAATAGCCGCTGTTGACCGTTTTACAGCCCCGCCTGGCGTTCGcc                                                                 >  3:93416/1‑149 (MQ=255)
                                                                          aaaGAATCATACGTGCCGAGGTGTTTATCACCAACACGCAAGCTGTCGCAGGATCCTACAAACCGGCTGTGACAAGCGGTTTTTTTATAACAGTTTTTTACAAAGAATAGCCGCTGTTGACCGTTTTACAGCCCCGCCTGGCGTTCGcc                                                                 >  1:269943/1‑149 (MQ=255)
                                                                          aaaGAATCATACGTGCCGAGGTGTTTATCACCAACACGCAAGCTGTCGCAGGATCCTACAAACCGGCTGTGACAAGCGGTTTTTTTATAACAGTTTTTTACAAAGAATAGCCGCTGTTGACCGTTTTACAGCCCCACCTGGCGtt                                                                     >  8:234739/1‑145 (MQ=255)
                                                                          aaaGAATCATACGTGCCGAGGTGTTTATCACCAACACGCAAGCTGTCGCAGGATCCTACAAACCGGCTGTGACAAGCGGTTTTTTTATAACAGTTTTTTACAAAGAATAGCCGCTGTTGACCGTTTTACAGCCCCACCTGGCGtt                                                                     <  7:234739/145‑1 (MQ=255)
                                                                           aaGAATCATACGTGCCGAGGTGTTTATCACCAACACGCAAGCTGTCGCAGGATCCTACAAACCGGCTGTGACAAGCGGTTTTTTTATAACAGTTTTTTACAAAGAATAGCCGCTGTTGACCGTTTTACAGCCCCGCCTGGCGTTCGccc                                                                >  2:255375/1‑149 (MQ=255)
                                                                           aaGAATCATACGTGCCGAGGTGTTTATCACCAACACGCAAGCTGTCGCAGGATCCTACAAACCGGCTGTGACAAGCGGTTTTTTTATAACAGTTTTTTACAAAGAATAGCCGCTGTTGACCGTTTTACAGCCCCGCCTGGCGTTCGccc                                                                <  5:298660/149‑1 (MQ=255)
                                                                             gAATCATACGTGCCGAGGTGTTTATCACCAACACGCAAGCTGTCGCAGGATCCTACAAACCGGCTGTGACAAGCGGTTTTTTTATAACAGTTTTTTACAAAGAATAGCCGCTGTTGACCGTTTTACAGCCCCGCCTGGCGTTCGCCCgg                                                              <  1:202655/149‑1 (MQ=255)
                                                                             gAATCATACGTGCCGAGGTGTTTATCACCAACACGCAAGCTGTCGCAGGATCCTACAAACCGGCTGTGACAAGCGGTTTTTTTATAACAGTTTTTTACAAAGAATAGCCGCTGTTGACCGTTTTACAGCCCCGCCTGGCGTTCGCCCgg                                                              <  2:196763/149‑1 (MQ=255)
                                                                                tCATACGTGCCGAGGTGTTTATCACCAACACGCAAGCTGTCGCAGGATCCTACAAACCGGCTGTGACAAGCGGTTTTTTTATAACAGTTTTTTACAAAGAATAGCCGCTGTTGACCGTTTTACAGCCCCGCCTGGCGTTCGCCCGGGtt                                                           >  1:291919/1‑149 (MQ=255)
                                                                                   tACGTGCCGAGGTGTTTATCACCAACACGCAAGCTGTCGCAGGATCCTACAAACCGGCTGTGACAAGCGGTTTTTTTATAACAGTTTTTTACAAAGAATAGCCGCTGTTGACCGTTTTACAGCCCCGCCTGGCGTTCGCCCGGGTTCTg                                                        <  3:33585/149‑1 (MQ=255)
                                                                                     cGTGCCGAGGTGTTTATCACCAACACGCAAGCTGTCGCAGGATCCTACAAACCGGCTGTGACAAGCGGTTTTTTTATAACAGTTTTTTACAAAGAATAGCCGCTGTTGACCGTTTTACAGCCCCGCCTGGCGTTCGCCCGGGTTCTGc                                                       <  7:382978/148‑1 (MQ=255)
                                                                                      gTGCCGAGGTGTTTATCACCAACACGCAAGCTGTCGCAGGATCCTACAAACCGGCTGTGACAAGCGGTTTTTTTATAACAGTTTTTTACAAAGAATAGCCGCTGTTGACCGTTTTACAGCCCCGCCTGGCGt                                                                      <  5:362074/132‑1 (MQ=255)
                                                                                      gTGCCGAGGTGTTTATCACCAACACGCAAGCTGTCGCAGGATCCTACAAACCGGCTGTGACAAGCGGTTTTTTTATAACAGTTTTTTACAAAGAATAGCCGCTGTTGACCGTTTTACAGCCCCGCCTGGCGt                                                                      >  6:362074/1‑132 (MQ=255)
                                                                                         ccGAGGTGTTTATCACCAACACGCAAGCTGTCGCAGGATCCTACAAACCGGCTGTGACAAGCGGTTTTTTTATAACAGTTTTTTACAAAGAATAGCCGCTGTTGACCGTTTTACAGCCCCGCCTGGCGTTCGCCCGGGTTCTGCACATc                                                  >  8:221310/1‑149 (MQ=255)
                                                                                             ggTGTTTATCACCAACACGCAAGCTGTCGCAGGATCCTACAAACCGGCTGTGACAAGCGGTTTTTTTATAACAGTTTTTTACAAAGAATAGCCGCTGTTGACCGTTTTACAGCCCCGCCTGGCGTTCGCCCGGGTTCTGCACATCGCgt                                              <  5:411602/149‑1 (MQ=255)
                                                                                                gtTTATCACCAACACGCAAGCTGTCGCAGGATCCTACAAACCGGCTGTGACAAGCGGTTTTTTTATAACAGt                                                                                                                        >  7:203470/1‑72 (MQ=255)
                                                                                                gtTTATCACCAACACGCAAGCTGTCGCAGGATCCTACAAACCGGCTGTGACAAGCGGTTTTTTTATAACAGt                                                                                                                        <  8:203470/72‑1 (MQ=255)
                                                                                                 tttATCACCAACACGCAAGCTGTCGCAGGATCCTACAAACCGGCTGTGACAAGCGGTTTTTTTATAACAGTTTTTTACAAAGAATAGCCGCTGTTGACCGTTTTACAGCCCCGCCTGGCGTTCGCCCGGGTTCTGCACATCGCGTGttt                                          >  4:304639/1‑149 (MQ=255)
                                                                                                 tttATCACCAACACGCAAGCTGTCGCAGGATCCTACAAACCGGCTGTGACAAGCGGTTTTTTTATAACAGTTTTTTACAAAGAATAGCCGCTGTTGACCGTTTTACAGCACCGCCTGGCGTTCGACCGGGTTATGCACATCTCGTGttt                                          >  7:179616/1‑149 (MQ=255)
                                                                                                    aTCACCAACACGCAAGCTGTCGCAGGATCCTACAAACCGGCTGTGACAAGCGGTTGTTTTTATAACAGTTTTTTACAAAGAATAGCCGCTGTTGACCGTTTTACAGCCCCGCCTGGCGTTCGCCCGGGTTCTGCACATCGCGTGTTTc                                         >  1:314940/1‑148 (MQ=255)
                                                                                                         cAACACGCAAGCTGTCGCAGGATCCTACAAACCGGCTGTGACAAGCGGTTTTTTTATAACAGTTTTTTACAAAGAATAGCCGCTGTTGACCGTTTTACAGCCCCGCCTGGCGTTCGCCCGGGTTCTg                                                        <  3:129344/127‑1 (MQ=255)
                                                                                                         cAACACGCAAGCTGTCGCAGGATCCTACAAACCGGCTGTGACAAGCGGTTTTTTTATAACAGTTTTTTACAAAGAATAGCCGCTGTTGACCGTTTTACAGCCCCGCCTGGCGTTCGCCCGGGTTCTg                                                        >  4:129344/1‑127 (MQ=255)
                                                                                                             acGCAAGCTGTCGCAGGATCCTACAAACCGGCTGTGACAAGCGGTTTTTTTATAACAGTTTTTTACAAAGAATAGCCGCTGTTGACCGTTTTACAGCCCCGCCTGGCGTTCGCCCGGGTTCTGCACATCGCGTGTTTCGACAGTTTcgc                              <  3:304639/149‑1 (MQ=255)
                                                                                                             acGCAAGCTGTCGCAGGATCCTACAAACCGGCTGTGACAAGCGGTTTTTTTATAACAGTTTTTTACAAAGAATAGCCGCTGTTGACCGTTTTACAGCCCCGCCTGGCGTTCGCCCGGGTTCTGCACATCGCGTGTTTCGACAGTTTcgc                              >  1:108557/1‑149 (MQ=255)
                                                                                                               gCAAGCTGTCGCAGGATCCTACAAACCGGCTGTGACAAGCGGTTTTTTTATAACAGTTTTTTACAAAGAATAGCCGCTGTTGACCGTTTTACAGCCCCGCCTGGCGTTCGCCCGGGTTCTGCACATCGCGTGTTTCGACAGTTTCGCGa                            <  7:117788/149‑1 (MQ=255)
                                                                                                                   gCTGTCGCAGGATCCTACAAACCGGCTGTGACAAGCGGTTTTTTTATAACAGTTTTTTACAAAGAATAGCCGCTGTTGACCGTTTTACAGCCCCGCCTGGCGTTCGCCCGGGTTCTGCACATCGCGTGTTTCGACAGTTTCGCGAGGgc                        >  8:263062/1‑149 (MQ=255)
                                                                                                                       tCGCAGGATCCTACAAACCGGCTGTGACAAGCGGTTTTTTTATAACAGTTTTTTACAAAGAATAGCCGCTGTTGACCGTTTTACAGCCCCGCCTGGCGTTCGCCCGGGTTCTGCACATCGCGTGTTTCGACAGTTTCGCGAGGGcgccg                    >  7:369050/1‑149 (MQ=255)
                                                                                                                           aGGATCCTACAAACCGGCTGTGACAAGCGGTTTTTTTATAACAGTTTTTTACAAAGAATAGCCGCTGTTGACCGTTTTACAGCCCCGCCTGGCGTTCGCCCGGGTTCTGCACATCGCGTGTTTCGACAGTTTCGCGAGGGCGCCGCAcc                >  7:369984/1‑149 (MQ=255)
                                                                                                                              aTCCTACAAACCGGCTGTGACAAGCGGTTTTTTTATAACAGTTTTTTACAAAGAATAGCCGCTGTTGACCGTTTTACAGCCCCGCCTGGCGTTCGCCCGGGTTCTGCACATCGCGTGTTTCGACAGTTTCGCGAGGGCGCCGCACCTTg             <  8:369984/149‑1 (MQ=255)
                                                                                                                               ataaaaaaaCCGCTTGTCACAGCCGGTTTTTTTATAACAGTTTTTTACAAAGAATAGCCGCTGTTGACCGTTTTACAGCCCCGCCTGGCGTTCGCCCGGGTTCTGCACATCGCGTGTTTCGACAGTTTCGCGa                            >  1:32498/7‑133 (MQ=255)
                                                                                                                               ataaaaaaaCCGCTTGTCACAGCCGGTTTTTTTATAACAGTTTTTTACAAAGAATAGCCGCTGTTGACCGTTTTACAGCCCCGCCTGGCGTTCGCCCGGGTTCTGCACATCGCGTGTTTCGACAGTTTCGCGa                            <  2:32498/127‑1 (MQ=255)
                                                                                                                                  tACAAACCGGCTGTGACAAGCGGTTTTTTTATAACAGTTTTTTACAAAGAATAGCCGCTGTTGACCGTTTTACAGCCCCGCCTGGCGTTCGCCCGGGTTCTGCACATCGCGTGTTTCGACAGTTTCGCGAGGGCGCCGCACCTTGAGtt         >  1:87819/1‑149 (MQ=255)
                                                                                                                                  tACAAACCGGCTGTGACAAGCGGTTTTTTTATAACAGTTTTTTACAAAGAATAGCCGCTGTTGACCGTTTTACAGCCCCGCCTGGCGTTCGCCCGGGTTCTGCACATCGCGTGTTTCGACAGTTTCGCGAGGGCGCCGCACCTTGAGt          >  7:260306/1‑148 (MQ=255)
                                                                                                                                   aCAAACCGGCTGTGACAAGCGGTTTTTTTATAACAGTTTTTTACAAAGAATAGCCGCTGTTGACCGTTTTACAGCCCCGCCTGGCGTTCGCCCGGGTTCTGCACATCGCGTGTTTCGACAGTTTCGCGAGGGCGCCGCACCTTGAGtt         >  3:223979/1‑148 (MQ=255)
                                                                                                                                     aaaCCGGCTGTGACAAGCGGTTTTTTTATAACAGTTTTTTACAAAGAATAGCCGCTGTTGACCGTTTTACAGCCCCGCCTGGCGTTCGCCCGGGTTCTGCACATCGCGTGTTTCGACAGTTTCGCGAGGGCGCCGCAcc                <  8:312531/139‑1 (MQ=255)
                                                                                                                                     aaaCCGGCTGTGACAAGCGGTTTTTTTATAACAGTTTTTTACAAAGAATAGCCGCTGTTGACCGTTTTACAGCCCCGCCTGGCGTTCGCCCGGGTTCTGCACATCGCGTGTTTCGACAGTTTCGCGAGGGCGCCGCAcc                >  7:312531/1‑139 (MQ=255)
                                                                                                                                      aaCCGGCTGTGACAAGCGGTTTTTTTATAACAGTTTTTTACAAAGAATAGCCGCTGTTGACCGTTTTACAGCCCCGCCTGGCGTTCGCCCGGGTTCTGCACATCGCGTGTTTCGACAGTTTCGCGAGGGCGCCGCACCTTGAGTTCACt     <  4:235468/149‑1 (MQ=255)
                                                                                                                                      aaCCGGCTGTGACAAGCGGTTTTTTTATAACAGTTTTTTACAAAGAATAGCCGCTGTTGACCGTTTTACAGCCCCGCCTGGCGTTCGCCCGGGTTCTGCACATCGCGTGTTTCGACAGTTTCGCGAGGGCGCCGCACCTTGAGTTCACt     <  2:269943/149‑1 (MQ=255)
                                                                                                                                      aaCCGGCTGTGACAAGCGGTTTTTTTATAACAGTTTTTTACAAAGAATAGCCGCTGTTGACCGTTTTACAGCCCCGCCTGGCGTTCGCCCGGGTTCTGCACATCGCGTGTTACGACAGTTTCGCGAGGGAGCCGCACCTTGAGTTCACt     >  7:254434/1‑149 (MQ=255)
                                                                                                                                        ccGGCTGTGACAAGCGGTTTTTTTATAACAGTTTTTTACAAAGAATAGCCGCTGTTGACCGTTTTACAGCCCCGCCTGGCGTTCGCCCGGGTTCTGCACATCGCGTGTTTCGACAGCTTCGCGAGGGCGCCGCACCTTGAGTTCACTGa   <  1:255375/149‑1 (MQ=255)
                                                                                                                                        ccGGCTGTGACAAGCGGTTTTTTTATAACAGTTTTTTACAAAGAATAGCCGCTGTTGACCCTTTTACAGCCCCGCCTGTCGTTCACCCGGGTTCTGCACATCGCTTGTTTCGACAGTTACGCGAGGGAGCCGCACCTTGAGTTCACTGa   >  7:281051/1‑149 (MQ=255)
                                                                                                                                         cGGCTGTGACAAGCGGTTTTTTTATAACAGTTTTTTACAAAGAATAGCCGCTGTTGACCGTTTTACAGCCCCGCCTGGCGTTCGCCCGGGTTCTGCACATCGCGTGTTTCGACAGTTTCGCGAGGGCGCCGCACCTTGAGTTCACTGac  <  4:178817/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                              
CCAGGCAACTGCCCACGAACGGCAGAAGCAATAGCACTATCAAGGACATAGCGGCCTATTCTGGAGAAGTTTGCAAAGAATCATACGTGCCGAGGTGTTTATCACCAACACGCAAGCTGTCGCAGGATCCTACAAACCGGCTGTGACAAGCGGTTTTTTTATAACAGTTTTTTACAAAGAATAGCCGCTGTTGACCGTTTTACAGCCCCGCCTGGCGTTCGCCCGGGTTCTGCACATCGCGTGTTTCGACAGTTTCGCGAGGGCGCCGCACCTTGAGTTCACTGAC  >  NC_002947/2538990‑2539275

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