Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,716,648 G→C 6.0% V627L (GTG→CTG)  PP_4173 → two‑component system sensor histidine kinase/response regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,716,6480GC 6.0% 248.4 / 4.4 79V627L (GTG→CTG) PP_4173two‑component system sensor histidine kinase/response regulator
Reads supporting (aligned to +/- strand):  ref base G (35/39);  new base C (3/2);  total (38/41)
Fisher's exact test for biased strand distribution p-value = 6.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

GACGCATCCGCCAGTTGCCCAACGGCGCCGGGCTACCGATCATCGCCCTGACCGCCAACGCGTTGCAGGGCGACCGCGAACGTTGCATGGCCGCCGGTATGGACGATTACCTGAGCAAACCCTTGCGCCGCACCGAATTGCAGCAGGTGCTGCAGCGCTGGTTACCCGGTCAGACAACCGCGACTGGCGATAAATGCTAAAGTGCGGCAGTCTTAAACACTGGACAGGACCTTCATCGGACCTGAAAATAGACGTTTCAGTGCACACCTGTACTTCTTTCGCCAGGTGCGCT  >  NC_002947/4716502‑4716793
                                                                                                                                                  |                                                                                                                                                 
gACGCATCCGCCAGTTGCCCAACGGCGCCGGGCTACCGATCATCGCCCTGACCGCCAACGCGTTGCAGGGCGACCGCGAACGTTGCATGGCCGCCGGTATGGACGATTACCTGAGCAAACCCTTGCGCCGCACCGAATTGCAGCAGGtg                                                                                                                                                 <  3:259127/149‑1 (MQ=255)
   gCATCCGCCAGTTGCCCAACGGCGCCGGGCTACCGATCATCGCCCTGACCGCCAACGCGTTGCAGGGCGACCGCGAACGTTGCATGGCCGCCGGTATGGACGATTACCTGAGCAAACCCTTGCGCCGCACCGAATTGCAGCAGGtgctg                                                                                                                                              >  3:142382/1‑149 (MQ=255)
     aTCCGCCAGTTGCCCAACGGCGCCGGGCTACCGATCATCGCCCTGACCGCCAACGCGTTGCAGGGCGACCGCGAACGTTGCATGGCCGCCGGTATGGACGATTACCTGAGCAAACCCTTGCGCCGCCCCGAATTGCAGCAGGtgctg                                                                                                                                              >  4:108399/1‑147 (MQ=255)
     aTCCGCCAGTTGCCCAACGGCGCCGGGCTACCGATCATCGCCCTGACCGCCAACGCGTTGCAGGGCGACCGCGAACGTTGCATGGCCGCCGGTATGGACGATTACCTGAGCAAACCCTTGCGCCGCACCGAATTGCAGCAGGtgctg                                                                                                                                              >  2:209752/1‑147 (MQ=255)
     aTCCGCCAGTTGCCCAACGGCGCCGGGCTACCGATCATCGCCCTGACCGCCAACGCGTTGCAGGGCGACCGCGAACGTTGCATGGCCGCCGGTATGGACGATTACCTGAGCAAACCCTTGCGCCGCACCGAATTGCAGCAGGtgctg                                                                                                                                              >  4:25296/1‑147 (MQ=255)
     aTCCGCCAGTTGCCCAACGGCGCCGGGCTACCGATCATCGCCCTGACCGCCAACGCGTTGCAGGGCGACCGCGAACGTTGCATGGCCGCCGGTATGGACGATTACCTGAGCAAACCCTTGCGCCGCACCGAATTGCAGCAGGTGCTGCa                                                                                                                                            >  1:375406/1‑149 (MQ=255)
       ccGCCAGTTGCCCAACGGCGCCGGGCTACCGATCATCGCCCTGACCGCAAACGCGTTGCAGGGCGACCGCGAACGTTGCATGGCCGCCGGTATGGACGATTACCTGAGCAAACCATTGCGCCGAACCGAATTGCAGCAGCTGCTGCAgc                                                                                                                                          >  6:271221/1‑149 (MQ=255)
        cGCCAGTTGCCCAACGGCGCCGGGCTACCGATCATCGCCCTGACCGCCAACGCGTTGCAGGGCGACCGCGAACGTTGCATGGCCGCCGGTATGGACGATTACCTGAGCAAACCCTTGCGCCGCACCGAATTGCAGCAGGtgct                                                                                                                                               <  1:207672/143‑1 (MQ=255)
        cGCCAGTTGCCCAACGGCGCCGGGCTACCGATCATCGCCCTGACCGCCAACGCGTTGCAGGGCGACCGCGAACGTTGCATGGCCGCCGGTATGGACGATTACCTGAGCAAACCCTTGCGCCGCACCGAATTGCAGCAGGtgct                                                                                                                                               >  2:207672/1‑143 (MQ=255)
         gCCAGTTGCCCAACGGCGCCGGGCTACCGATCATCGCCCTGACCGCCAACGCGTTGCAGGGCGACCGCGAACGTTGCATGGCCGCCGGTATGGACGATTACCTGAGCAAACCCTTGCGCCGCACCGAATTGCAGCAGGTGCTGCAgcgc                                                                                                                                        <  4:3251/149‑1 (MQ=255)
           cAGTTGCCCAACGGCGCCGGGCTACCGATCATCGCCCTGACCGCCAACGCGTTGCAGGGCGACCGCGAACGTTGCATGGCCGCCGGTATGGACGATTACCTGAGCAAACCCTTGCGCCGCACCGAATTGCAGCAGGTGCTGCAGCGCTg                                                                                                                                      <  4:232456/149‑1 (MQ=255)
           cAGTTGCCCAACGGCGCCGGGCTACCGATCATCGCCCTGACCGCCAACGCGTTGCAGGGCGACCGCGAACGTTGCATGGCCGCCGGTATGGACGATTACCTGAGCAAACCCTTGCGCCGCACCGAATTGCAGCAGGTGCTGCAGCGCTg                                                                                                                                      <  2:187146/149‑1 (MQ=255)
           cAGTTGCCCAACGGCGCCGGGCTACCGATCATCGCCCTGACCGCCAACGCGTTGCAGGGCGACCGCGAACGTTGCATGGCCGCCGGTATGGACGATTACCTGAGCAAACCCTTGCGCCGCACCGAATTGCAGCAGGTGCTGCAGCGCTg                                                                                                                                      <  4:85333/149‑1 (MQ=255)
           cAGTTCCCCACCGGCGCCGGGCTACCTATCATCGCCCTGACCGCCAACGCGTTGCAGGGCGACCGCGAACGTTGCATGGCCGCCGGTATGGACGATTACCTGAGCAAACCCTTGCGCCGCACCGAATTGCAGCAGGTGCTGCAGCGCTg                                                                                                                                      <  7:341844/149‑1 (MQ=255)
            aGTTGCCCAACGGCGCCGGGCTACCGATCATCGCCCTGACCGCCAACGCGTTGCAGGGCGACCGCGAACGTTGCATGGCCGCCGGTATGGACGATTACCTGAGCAAACCCTTGCGCCGCACCGAATTGCAGCAGGTGCTGCAGCGCTgg                                                                                                                                     <  5:326045/149‑1 (MQ=255)
             gTTGCCCAACGGCGCCGGGCTACCGATCATCGCCCTGACCGCCAACGCGTTGCAGGGCGACCGCGAACGTTGCATGGCCGCCGGTATGGACGATTACCTGAGCAAACCCTTGCGCCGCACCGAATTGCAGCAGGTGCTGCAGCGCTGGt                                                                                                                                    <  1:8106/149‑1 (MQ=255)
                   cAACGGCGCCGGGCTACCGATCATCGCCCTGACCGCCAACGCGTTGCAGGGCGACCGCGAACGTTGCATGGCCGCCTGTTTGGACGATTACCTGAGCAAACCCTTGCGCCGCACCGAATTGCAGCAGGTGCTGCAGCGCTGGTTACCCg                                                                                                                              <  5:271221/149‑1 (MQ=255)
                   cAACGGCGCCGGGCTACCGATCATCGCCCTGACCGCCAACGCGTTGCAGGGCGACCGCGAACGTTGCATGGCCGCCGGTATGGACGATTACCTGAGCAAACCCTTGCGCCGCACCGAATTGCAGCAGGTGCTGCAGCGCTGGTTACCCg                                                                                                                              <  7:99355/149‑1 (MQ=255)
                        gcgcCGGGCTACCGATCATCGCCCTGACCGCCAACGCGTTGCAGGGCGACCGCGAACGTTGCATGGCCGCCGGTATGGACGATTACCTGAGCAAACCCTTGCGCCGCACCGAATTGCAGCAGGTGCTGCAGCGCTGGTTACCCGGTCAg                                                                                                                         <  3:40715/149‑1 (MQ=255)
                               gCTACCGATCATCGCCCTGACCGCCAACGCGTTGCAGGGCGACCGCGAACGTTGCATGGCCGCCGGTATGGACGATTACCTGAGCAAACCCTTGCGCCGCACCGAATTGCAGCAGGTGCTGCAGCGCTGGTTACCCGGTCAGACAACcg                                                                                                                  >  1:165145/1‑149 (MQ=255)
                                cTACCGATCATCGCCCTGACCGCCAACGCGTTGCAGGGCGACCGCGAACGTTGCATGGCCGCCGGTATGGACGATTACCTGAGCAAACCCTTGCGCCGCACCGAATTGCAGCACCTGCTGCaat                                                                                                                                          <  4:59088/124‑3 (MQ=255)
                                cTACCGATCATCGCCCTGACCGCCAACGCGTTGCAGGGCGACCGCGAACGTTGCATGGCCGCCGGTATGGACGATTACCTGAGCAAACCCTTGCGCCGCACCGAATTGCAGCACCTGCTGCaat                                                                                                                                          >  3:59088/1‑122 (MQ=255)
                                          tcGCCCTGACCGCCAACGCGTTGCAGGGCGACCGCGAACGTTGCATGGCCGCCGGTATGGACGATTACCTGAGCAAACCCTTGCGCCGCACCGAATTGCAGCAGGTGCTGCAGCGCTGGTTACCCGGTCAGACAACCGCGACTGGCGa                                                                                                        >  6:86344/1‑148 (MQ=255)
                                             cccTGACCGCCAACGCGTTGCAGGGCGACCGCGAACGTTGCATGGCCGCCGGTATGGACGATTACCTGAGCAAACCCTTGCGCCGCACCGAATTGCAGCAGGTGCTGCAGCGCTGGTTACCCGGTCAGACAACCGCGACTGGCGATaaa                                                                                                    >  2:209425/1‑149 (MQ=255)
                                              ccTGACCGCCAACGCGTTGCAGGGCGACCGCGAACGTTGCATGGCCGCCGGTATGGACGATTACCTGAGCAAACCCTTGCGCCGCACCGAATTGCAGCAGGTGCTGCAGCGCTGGTTACCCGGTCAAACAACCGCGACTGGCGATAAAt                                                                                                   <  4:104636/149‑1 (MQ=255)
                                                tGACCGCCAACGCGTTGCAGGGCGACCGCGAACGTTGCATGGCCGCCGGTATGGACGATTACCTGAGCAAACCCTTGCGCCGCACCGAATTGCAGCAGGTGCTGCAGCGCTGGTTACCCGGTCAGACAACCGCGACTGGCGATAAATGc                                                                                                 >  7:180890/1‑149 (MQ=255)
                                                  aCCGCCAACGCGTTGCAGGGCGACCGCGAACGTTGCATGGCCGCCGGTATGGACGATTACCTGAGCAAACCCTTGCGCCGCACCGAATTGCAGCAGGTGCTGCAGCGCTGGTTACCCGGTCAGACAACCGCGACTGGCGATAAATGCt                                                                                                <  8:180890/148‑1 (MQ=255)
                                                  aCCGCCAACGCGTTGCAGGGCGACCGCGAACGTTGCATGGCCGCCGGTATGGACGATTACCTGAGCAAACCCTTGCGCCGCACCGAATTGCAGCACCTGCTGCaattcg                                                                                                                                       <  8:250066/109‑6 (MQ=255)
                                                  aCCGCCAACGCGTTGCAGGGCGACCGCGAACGTTGCATGGCCGCCGGTATGGACGATTACCTGAGCAAACCCTTGCGCCGAACCGAATTGCAGCACCTGCTGCaattcg                                                                                                                                       >  7:250066/1‑104 (MQ=255)
                                                   ccGCCAACGCGTTGCAGGGCGACCGCGAACGTTGCATGGCCGCCGGTATGGACGATTACCTGAGCAAACCCTTGCGCCGCACCGAATTGCAGCAGGTGCTGCAGCGCTGGTTACCCGGTCAGACAACCGCGACTGGCGATAAATGCTa                                                                                               >  4:255884/1‑148 (MQ=255)
                                                    cGCCAACGCGTTGCAGGGCGACCGCGAACGTTGCATGGCCGCCGGTATGGACGATTACCTGAGCAAACCCTTGCGCCGCACCGAATTGCAGCAGGTGCTGCAGCGCTGGTTACCCGGTCAGACAACCGCGACTGGCGATAAATGCTaaa                                                                                             >  3:94889/1‑149 (MQ=255)
                                                    cGCCAACGCGTTGCAGGGCGACCGCGAACGTTGCATGGCCGCCGGTATGGACGATTACCTGAGCAAACCCTTGCGCCGCACCGAATTGCAGCAGGTGCTGCAGCGCTGGTTACCCGGTCAGACAACCGCGACTGGCGATAAATGCTa                                                                                               >  5:276300/1‑147 (MQ=255)
                                                      ccAACGCGTTGCAGGGCGACCGCGAACGTTGCATGGCCGCCGGTATGGACGATTACCTGAGCAAACCCTTGCGCCGCACCGAATTGCAGCAGGTGCTGCAGCGCTGGTTACCCGGTCAGACAACCGCGACTGGCGATAAATGCTAAAGt                                                                                           >  4:119338/1‑149 (MQ=255)
                                                       cAACGCGTTGCAGGGCGACCGCGAACGTTGCATGGCCGCCGGTATGGACGATTACCTGAGCAAACCCTTGCGCCGCACCGAATTGCAGCAGGTGCTGCAGCGCTGGTTACCCGGTCAGACAACCGCGACTGGCGATAAATGCTAAAGTg                                                                                          >  8:63216/1‑149 (MQ=255)
                                                       cAACGCGTTGCAGGGCGACCGCGAACGTTGCATGGCCGCCGGTATGGACGATTACCTGAGCAAACCCTTGCGCCGCACCGAATTGCAGCAGGTGCTGCAGCGCTGGTTACCCGGTCAGACAACCGCGACTGGCGATAAATGCTAAAGTg                                                                                          >  1:334706/1‑149 (MQ=255)
                                                         aCGCGTTGCAGGGCGACCGCGAACGTTGCATGGCCGCCGGTATGGACGATTACCTGAGCAAACCCTTGCGCCGCACCGAATTGCAGCAGGTGCTGCAGCGCTGGTTACCCGGTCAGACAACCGCGACTGGCGATAAATGCTAAAGTGCg                                                                                        <  1:209425/149‑1 (MQ=255)
                                                             gTTGCAGGGCGACCGCGAACGTTGCATGGCCGCCGGTATGGACGATTACCTGAGCAAACCCTTGCGCCGCACCGAATTGCAGCAGGTGCTGCAGCGCTGGTTACCCGGTCAGACAACCGCGACTGGCGATAAATGCTAAAGTGCGGCAg                                                                                    <  2:375406/149‑1 (MQ=255)
                                                                 cAGGGCGACCGCGAACGTTGCATGGCCGCCGGTATGGACGATTACCTGAGCAAACCCTTGCGCCGCACCGAATTGCAGCAGGTGCTGCAGCGCTGGTTACCCGGTCAGACAACCGCGACTGGCGATAAATGCTAAAGTGCGGCAGTCt                                                                                 <  2:334706/148‑1 (MQ=255)
                                                                    ggCGACCGCGAACGTTGCATGGCCGCCGGTATGGACGATTACCTGAGCAAACCCTTGCGCCGCACCGAATTGCAGCAGGTGCTGCAGCGCTGGTTACCCGGTCAGACAACCGCGACTGGCGATAAATGCTAAAGTGCGGCAGTCTTaaa                                                                             <  3:108399/149‑1 (MQ=255)
                                                                         ccGCGAACGTTGCATGGCCGCCGGTATGGACGATTACCTGAGCAAACCCTTGCGCCGCACCGAATTGCAGCAGGTGCTGCAGCGCTGGTTACCCGGTCAGACAACCGCGACTGGCGATAAATGCTAAAGTGCGGCAGTCTTAAACACTg                                                                        <  8:13091/149‑1 (MQ=255)
                                                                           gcgAACGTTGCATGGCCGCCGGTATGGACGATTACCTGAGCAAACCCTTGCGCCGCACCGAATTGCAGCAGGTGCTGCAGCGCTGGTTACCCGGTCAGACAACCGCGACTGGCGATAAATGCTAAAGTGCGGCAGTCTTAAACACTgg                                                                       >  2:325012/1‑148 (MQ=255)
                                                                             gAACGTTGCATGGCCGCCGGTATGGACGATTACCTGAGCAAACCCTTGCGCCGCACCGAATTGCAGCAGGTGCTGCAGCGCTGGTTACCCGGTCAGACAACCGCGACTGGCGATAAATGCTAAAGTGCGGCAGTCTTAAACACTGGa                                                                      >  2:364515/1‑147 (MQ=255)
                                                                             gAACGTTGCATGGCCGCCGGTATGGACGATTACCTGAGCAAACCCTTGCGCCGCACCGAATTGCAGCAGGTGCTGCAGCGCTGGTTACCCGGTCAGACAACCGCGACTGGCGATAAATGCTAAAGTGCGGCAGTCTTAAACACTGGACa                                                                    <  2:99798/149‑1 (MQ=255)
                                                                                 gTTGCATGGCCGCCGGTATGGAGGATTACCTGAGCAAACCCTTGCGCCGCACCGAATTGCAGCAGGTGCTGCAGCGCTGGTTACCCGGTCAGACAACCGCGACTGGCGATAAATGCTAAAGTGCGGCAGTCTTAAACACTGGACAGGAc                                                                >  7:59507/1‑149 (MQ=255)
                                                                                             ccgGTATGGACGATTACCTGAGCAAACCCTTGCGCCGCACCGAATTGCAGCAGGTGCTGCAGCGCTGGTTACCCGGTCAGACAACCGCGACTGGCGATAAATGCTAAAGTGCGGCAGTCTTAAACACTGGACAGGACCTTCATCGGAcc                                                    <  1:209752/149‑1 (MQ=255)
                                                                                                 tATGGACGATTACCTGAGCAAACCCTTGCGCCGCACCGAATTGCAGGAGGTGCTGCAGCGCTGGTTACCCGGTCAGCCAACCGCGACTGGCGATAAATGCTAAAGTGCGGCAGTCTTAAACACTGGACAGGACCTTCATCGGACCTGa                                                 >  1:281476/1‑148 (MQ=255)
                                                                                                     gACGATTACCTGAGCAAACCCTTGCGCCGCACCGAATTGCAGCAGGTGCTGCAGCGCTGGTTACCCGGTCAGACAACCGCGACTGGCGATAAATGCTAAAGTGCGGCAGTCTTAAACACTGGACAGGACCTTCATCGGACCTGaaaa                                              >  4:36861/1‑147 (MQ=255)
                                                                                                      aCGATTACCTGAGCAAACCCTTGCGCCGGACCGAATTGCAGCAGGTGCTGCAGCGCTGGTTACCCGGTCAGACAACCGCGACTGGCGATAAATGCTAAAGTGCGGCAGTCTTAAACACTGGACAGGCCCTTCATCGGACCTGaaaa                                              >  1:218836/1‑146 (MQ=255)
                                                                                                      aCGATTACCTGAGCAAACCCTTGCGCCGCACCGAATTGCAGCAGGTGCTGCAGCGCTGGTTACCCGGTCAGACAACCGCGACTGGCGATAAATGCTAAAGTGCGGCAGTCTTAAACACTGGACAGGACCTTCATCGGACCTGaaaa                                              >  6:297581/1‑146 (MQ=255)
                                                                                                       cGATTACCTGAGCAAACCCTTGCGCCGCACCGAATTGCAGCAGGTGCTGCAGCGCTGGTTACCCGGTCAGACAACCGCGACTGGCGATAAATGCTAAAGTGCGGCAGTCTTAAACACTGGACAGGACCTTCATCGGACCTGAAAATAGa                                          <  5:86344/149‑1 (MQ=255)
                                                                                                        gATTACCTGAGCAAACCCTTGCGCTGCACCGAATTGCAGCAGGTGCTGCAGCGCTGGTTACCCGGTCAGACAACCGCGACTGGCGATAAATGCTAAAGTGCGGCAGTCTTAAACACTGGACAGGACCTTCATCGGACCTGAAAATAGAc                                         >  6:109433/1‑149 (MQ=255)
                                                                                                              cTGAGCAAACCCTTGCGCCGCACCGAATTGCAGCAGGTGCTGCAGCGCTGGTTACCCGGTCAGACAACCGCGACTGGCGATAAATGCTAAAGTGCGGCAGTCTTAAACACTGGACAGGACCTTCATCGGACCTGAAAATAGACGTTTCa                                   <  4:311318/149‑1 (MQ=255)
                                                                                                              cTGAGCAAACCATTGCGCCGCACCGAATTGCAGCAGGTGCTGCAGCGCTGGTTACCCGGTCAGACAACCGCGACTGGCGATAAATGCTAAAGTGCGGCAGTCTTAAACACTGGACAGGACCTTCATCGGACCTGAAAATAGACGttt                                     >  1:167314/1‑147 (MQ=255)
                                                                                                                 aGCAAACCCTTGCGCCGCACCGAATTGCAGCAGGTGCTGCAGCGCTGGTTACCCGGTCAGACAACCGCGACTGGCGATAAATGCTAAAGTGCGGCAGTCTTAAACACTGGACAGGACCTTCATCGGACCTGAAAATAGACGTTTCAGTg                                <  1:364515/149‑1 (MQ=255)
                                                                                                                 aGCAAACCCTTGCGCCGCACCGAATTGCAGCAGGTGCTGCAGCGCTGGTTACCCGGTCAGACAACCGCGACTGGCGATAAATGCTAAAGTGCGGCAGTCTTAAACACTGGACAGGACCTTCATCGGACCTGAAAATAGACGTTTCAGTg                                <  6:357935/149‑1 (MQ=255)
                                                                                                                   cAAACCCTTGCGCCGCACCGAATTGCAGCAGGTGCTGCAGCGCTGGTTACCCGGTCAGACAACCGCGACTGGCGATAAATGCTAAAGTGCGGCAGTCTTAAACACTGGACAGGACCTTCATCGGACCTGAAAATAGACGTTTCAGTGca                              >  5:27947/1‑149 (MQ=255)
                                                                                                                    aaaCCCTTGCGCCGCACCGAATTGCAGCAGGTGCTGCAGCGCTGGTTACCCGGTCAGACAACCGCGACTGGCGATAAATGCTAAAGTGCGGCAGTCTTAAACACTGGACAGGACCTTCATCGGACCTGAAAATAGACGTTTCAGTGca                              >  6:376414/1‑148 (MQ=255)
                                                                                                                     aaCCCTTGCGCCGCACCGAATTGCAGCAGGTGCTGCAGCGCTGGTTACCCGGTCAGACAACCGCGACTGGCGATAAATGCTAAAGTGCGGCAGTCTTAAACACTGGACAGGACCTTCATCGGACCTGAAAATAGACGTTTCAGTGcaca                            >  6:85134/1‑149 (MQ=255)
                                                                                                                      aCCCTTGCGCCGCACCGAATTGCAGCAGGTGCTGCAGCGCTGGTTACCCGGTCAGACAACCGCGACTGGCGATAAATGCTAAAGTGCGGCAGTCTTAAACACTGGACAGGACCTTCATCGNACCTGAAAATAGACGTTTCAGTGCacac                           <  1:222563/149‑1 (MQ=255)
                                                                                                                      aCCCTTGCGCCGCACCGAATTGCAGCAGGTGCTGCAGCGCTGGTTACCCGGTCAGACAACCGCGACTGGCGATAAATGCTAAAGTGCGGCAGTCTTAAACACTGGACAGGACCTTCATCGGACCTGAAAATAGACGTTTCAGTGCacac                           <  6:27947/149‑1 (MQ=255)
                                                                                                                      aCCCTTGCGCCGCACCGAATTGCAGCAGGTGCTGCAGCGCTGGTTACCCGGTCAGACAACCGCGACTGGCGATAAATGCTAAAGTGCGGCAGTCTTAAACACTGGACAGGACCTTCATCGGACCTGAAAATAGACGTTTCAGTGCacac                           <  4:142382/149‑1 (MQ=255)
                                                                                                                        ccTTGCGCCGCACCGAATTGCAGCAGGTGCTGCAGCGCTGGTTACCCGGTCAGACAACCGCGACTGGCGATAAATGCTAAAGTGCGGCAGTCTTAAACACTGGACAGGACCTTCATCGGACCTGAAAATAGACGTTTCAGTGCACACCt                         >  4:30447/1‑149 (MQ=255)
                                                                                                                            gcgcCGCACCGAATTGCAGCAGGTGCTGCAGCGCTGGTTACCCGGTCAGACAACCGCGACTGGCGATAAATGCTAAAGTGCGGCAGTCTTAAACACTGGACAGGACCTTCATCGGACCTGAAAATAGACGTTTCAGTGCACACCTGt                       <  3:25296/147‑1 (MQ=255)
                                                                                                                               ccgcACCGAATTGCAGCAGGTGCTGCAGCGCTGGTTACCCGGTCAGACAACCGCGACTGGCGATAAATGCTAAAGTGCGGCAGTCTTAAACACTGGACAGGACCTTCATCGGAcc                                                    <  8:228695/115‑1 (MQ=255)
                                                                                                                                cgcACCGAATTGCAGCAGGTGCTGCAGCGCTGGTTACCCGGTCAGACAACCGCGACTGGCGATAAATGCTAAAGTGCGGCAGTCTTAAACACTGGACAGGACCTTCATCg                                                        <  7:141306/110‑1 (MQ=255)
                                                                                                                                cgcACCGAATTGCAGCAGGTGCTGCAGCGCTGGTTACCCGGTCAGACAACCGCGACTGGCGATAAATGCTAAAGTGCGGCAGTCTTAAACACTGGACAGGACCTTCATCg                                                        >  8:141306/1‑110 (MQ=255)
                                                                                                                                cgcACCGAATTGCAGCAGGTGCTGCAGCGCTGGTTACCCGGTCAGACAACCGCGACTGGCGATAAATGCTAAAGTGCGGCAGTCTTAAACACTGGACAGGACCTTCATCGGACCTGAAAATAGACGTTTCAGTGCACACCTGTActtct                 >  7:104349/1‑149 (MQ=255)
                                                                                                                                cgcACCGAATTGCAGCAGGTGCTGCAGCGCTGGTTACCCGGTCAGACAACCGCGACTGGCGATAAATGCTAAAGTGCGGCAGTCTTAAACACTGGACAGGACCTTCATCGGACCTGAAAATAGACGTTTCAGTGCACACCTGTActt                   <  7:212949/147‑1 (MQ=255)
                                                                                                                                cgcACCGAATTGCAGCAGGTGCTGCAGCGCTGGTTACCCGGTCAGACAACCGCGACTGGCGATAAATGCTAAAGTGCGGCAGTCTTAAACACTGGACAGGACCTTCATCGGACCTGAAAATAGACGTTTCAGTGCACACCTGTActt                   >  8:212949/1‑147 (MQ=255)
                                                                                                                                   aCCGAATTGCAGCAGGTGCTGCAGCGCTGGTTACCCGGTCAGACAACCGCGACTGGCGATAAATGCTAAAGTGCGGCAGTCTTAAACACTGGACAGGACCTTCATCGGACCTGAAAATAGACGTTTCAGTGCACACCTGTACTTCTTTc              <  2:165145/149‑1 (MQ=255)
                                                                                                                                     cGAATTGCAGCAGGTGCTGCAGCGCTGGTTACCCGGTCAGACAACCGCGACTGGCGATAAATGCTAAAGTGCGGCAGTCTTAAACACTGGACAGGACCTTCATCGGACCTGAAAATAGACGTTTCAGTGCAGACCTGTACTTCTTTCg             >  5:115718/1‑148 (MQ=255)
                                                                                                                                      gAATTGCAGCAGGTGCTGCAGCGCTGGTTACCCGGTCAGACAACCGCGACTGGCGATAAATGCTAAAGTGCGGCAGTCTTAAACACTGGACAGGACCTTCATCGGACCTGAAAATAGACGTTTCAGTGCACACCTGTACTTCTTTCGcc           <  3:255884/149‑1 (MQ=255)
                                                                                                                                       aaTTGCAGCAGGTGCTGCAGCGCTGGTTACCCGGTCAGACAACCGCGACTGGCGATAAATGCTAAAGTGCGGCAGTCTTAAACACTGGACAGGACCTTCATCgg                                                       <  2:205012/104‑1 (MQ=255)
                                                                                                                                       aaTTGCAGCAGGTGCTGCAGCGCTGGTTACCCGGTCAGACAACCGCGACTGGCGATAAATGCTAAAGTGCGGCAGTCTTAAACACTGGACAGGACCTTCATCgg                                                       >  1:205012/1‑104 (MQ=255)
                                                                                                                                        actGCAGCACGCGCCGCAGCGCTGGCTACCCCGCCACACAACCCCGACTCCCGATAAATGCTAAAGTGCGGCAGTCTTAAACACTGGACAGGACCTTCATCGGACCTGAAAATAGACGCTTCAGTGCACACCTGTACTTCTTTCGCCAg         <  8:59507/147‑1 (MQ=255)
                                                                                                                                        aTTGCAGCAGGTGCTGCAGCGCTGGTTACCCGGTCAGACAACCGCGACTGGCGATAAATGCTAAAGTGCGGCAGTCTTAAACACTGGACAGGACCTTCATCGGACCTGAAAATAGACGTTTCAGTGCACACCTGTACTTCTTTCGCCAg         <  7:63216/149‑1 (MQ=255)
                                                                                                                                           gcagcaGGTGCTGCAGCGCTGGTTACCCGGTCAGACAACCGCGACTGGCGATAAATGCTAAAGTGCGGCAGTCTTAAACACTGGACAGGACCTTCATCGGACCTGAAAATAGACGTTTCAGTGCACACCTGTACTTCTTTCGCCAGGTg      >  2:15914/1‑149 (MQ=255)
                                                                                                                                            cagcagGTGCTGCAGCGCTGGTTACCCGGTCAGACAACCGCGACTGGCGATAAATGCTAAAGTGCGGCAGTCTTAAACACTGGACAGGACCTTCATCGGACCTGAAAATAGACGTTTCAGTGCACACCTGTACTTCTTTCGCCAGGTg      <  3:119338/148‑1 (MQ=255)
                                                                                                                                               cagGTGCTGCAGCGCTGGTTACCCGGTCAGACAACCGCGACTGGCGATAAATGCTAAAGTGCGGCAGTCTTAAACACTGGACAGGACCTTCATCGGACCTGAAAATAGACGTTTCAGTGCACACCTGTACTTCTTTCGCCAGGTGCGCt  <  6:276300/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
GACGCATCCGCCAGTTGCCCAACGGCGCCGGGCTACCGATCATCGCCCTGACCGCCAACGCGTTGCAGGGCGACCGCGAACGTTGCATGGCCGCCGGTATGGACGATTACCTGAGCAAACCCTTGCGCCGCACCGAATTGCAGCAGGTGCTGCAGCGCTGGTTACCCGGTCAGACAACCGCGACTGGCGATAAATGCTAAAGTGCGGCAGTCTTAAACACTGGACAGGACCTTCATCGGACCTGAAAATAGACGTTTCAGTGCACACCTGTACTTCTTTCGCCAGGTGCGCT  >  NC_002947/4716502‑4716793

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