Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,210,088 C→A 10.8% G213V (GGT→GTT)  PP_2814 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,210,0880CA10.8% 188.1 / 10.8 64G213V (GGT→GTT) PP_2814hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (25/32);  new base A (3/4);  total (28/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.29e-01

GGATCTTCTGTACCCGCGGAGAAAGCGTGCGCACCAGCTCGAGCGGCACGATCAACGGGTAGGCTTCTTCAGGCGTGATTATCGAAGCCTGCATCGGGATGGTCAGCCCGCGACGGTAATAGTCGGTGGTCACGTTGGCCCCGCCACCGTGAACGACCTTGCCGCCGAACAGCACGGCATCCCCGGCTTTGAGCAGGGCAGGAATCGTCATCTCTGGGGTGCCGACATCGTCGAAGTCATCCCAGTCCTGGCTCCCGGGGATCAAGCGCGTGGCGCCGTTTTCCTCGGT  >  NC_002947/3209942‑3210230
                                                                                                                                                  |                                                                                                                                              
ggATCTTCTGTACCCGCGGAGAAAGCGTGCGCACCAGCTCGAGCGGCACGATCAACGGGTAGGCTTCTTCAGGCGTGATTATCGAAGCCTGCATCGGGATCGTCAGCCCGCGACGGTAATAGTCGGTGGTCACGTTGGCCCCGCCACCg                                                                                                                                              <  8:224377/149‑1 (MQ=255)
     ttctGTACCCGCGGAGAAAGCGTGCGCACCAGCTCGAGCGGCACGATCAACGGGTAGGCTTCTTCAGGCGTGATTATCGAAGCCTGCATCGGGATGGTCAGCCCGCGACGGTAATAGTCGGTGGTCACGTTGGCCCCGCCACCGTGAAc                                                                                                                                         >  3:233671/1‑149 (MQ=255)
       ctGTACCCGCGGAGAAAGCGTGCGCACCAGCTCGAGCGGCACGATCAACGGGTAGGCTTCTTCAGGCGTGATTATCGAAGCCTGCATCGGGATGGTCAGCCCGCGACGGTAATAGTCGGTGGTCACGGTGGCGGGGCCAACGTGa                                                                                                                                           <  7:5285/145‑1 (MQ=255)
        tGTACCCGCGGAGAAAGCGTGCGCACCAGCTCGAGCGGCACGATCAACGGGTAGGCTTCTTCAGGCGTGATTATCGAAGCCTGCATCGGGATGGTCAGCCCGCGACGGTAATAGTCGGTGGTCACGTTGGCCCCGCCACCGTGAACGAc                                                                                                                                      >  6:326297/1‑149 (MQ=255)
          tACCCGCGGAGAAAGCGTGCGCACCAGCTCGAGCGGCACGATCAACGGGTAGGCTTCTTCAGGCGTGATTATCGAAGCCTGCATCGGGATGGTCAGCCCGCGACGGTAATAGTCGGTGGTCACGTTGGCCCCGCCACCGTGAACGAc                                                                                                                                      <  3:189261/147‑1 (MQ=255)
          tACCCGCGGAGAAAGCGTGCGCACCAGCTCGAGCGGCACGATCAACGGGTAGGCTTCTTCAGGCGTGATTATCGAAGCCTGCATCGGGATGGTCAGCCCGCGACGGTAATAGTCGGTGGTCACGGTGGCGGGGCCAACGTGaccacc                                                                                                                                      <  7:28674/147‑6 (MQ=255)
            cccGCGGAGAAAGCGTGCGCACCAGCTCGAGCGGCACGATCAACGGGTAGGCTTCTTCAGGCGTGATTATCGAAGCCTGCATCGGGATGGTCAGCCCGCGACGGTAATAGTCGGTGGTCACGTTGGCCCCGCCACCGTGAACGACCTTg                                                                                                                                  >  3:39701/1‑149 (MQ=255)
             ccGCGGAGAAAGCGTGCGCACCAGCTCGAGCGGCACGATCAACGGGTAGGCTTCTTCAGGCGTGATTATCGAAGCCTGCATCGGGATGGTCAGCCCGCGACGGTAATAGTCGGTGGTCACGTTGGCCCCGCCACCGTGAACGACCTTgc                                                                                                                                 >  3:284181/1‑149 (MQ=255)
                     aaaGCGTGCGCACCAGCTCGAGCGGCACGATCAACGGGTAGGCTTCTTCAGGCGTGATTATCGAAGCCTGCATCGGGATGGTCAGCCCGCGACGGTAATAGTCGGTGGTCACGTTGGCCCCGCCACCGTGAACGACCTTGCCGCCGAAc                                                                                                                         <  5:170821/149‑1 (MQ=255)
                      aaGCGTGCGCACCAGCTCGAGCGGCACGATCAACGGGTAGGCTTCTTCAGGCGTGATTATCGAAGCCTGCATCGGGATGGTCAGCCCGCGACGGTAATAGTCGGTGGTCACGTTGGCCCCGCCACCGTGAACGACCTTGCCGCCGAACa                                                                                                                        <  2:241631/149‑1 (MQ=255)
                           tGCGCACCAGCTCGAGCGGCACGATCAACGGGTAGGCTTCTTCAGGCGTGATTATCGAAGCCTGCATCGGGATGGTCAGCCCGCGACGGTAATAGTCGGTGGTCACGTTGGCCCCGCCACCGTGAACGACCTTGCCGCCGAACAGCACg                                                                                                                   <  4:39701/149‑1 (MQ=255)
                              gcACCAGCTCGAGCGGCACGATCAACGGGTAGGCTTCTTCAGGCGTGATTATCGAAGCCTGCATCGGGATGGTCAGCCCGCGACGGTAATAGTCGGTGGTCACGTTGGCCCCGCCACCGTGAACGACCTTGCCGCCGAACAGCACGGCa                                                                                                                >  6:118823/1‑149 (MQ=255)
                                        gAGCGGCACGATCAACGGGTAGGCTTCTTCAGGCGTGATTATCGAAGCCTGCATCGGGATGGTCAGCCCGCGACGGTAATAGTCGGTGGTCACGTTGGCCCCGCCACCGTGAACGACCTTgccgc                                                                                                                              >  5:87320/1‑125 (MQ=255)
                                        gAGCGGCACGATCAACGGGTAGGCTTCTTCAGGCGTGATTATCGAAGCCTGCATCGGGATGGTCAGCCCGCGACGGTAATAGTCGGTGGTCACGTTGGCCCCGCCACCGTGAACGACCTTgccgc                                                                                                                              <  6:87320/125‑1 (MQ=255)
                                        gAGCGGCACGATCAACGGGTAGGCTTCTACAGGCGTGATTATCGAAGCCTGCATCGGGATGGTCAGCCCGCGACGGTAATAGTCGGTGGTCACGTTGGCCACCCCACCCTGAACGACCTTGCCGCCGAACAGCACGTCATACCCGGCtt                                                                                                      >  8:249411/1‑149 (MQ=255)
                                           cGGCACGATCAACGGGTAGGCTTCTTCAGGCGTGATTATCGAAGCCTGCATCGGGATGGTCAGCCCGCGACGGTAATAGTCGGTGGTCACGTTGGCCCCGCCACCGTGAACGACCTTGCCGCCGAACAGCACGGCATCCCCGGCTTTGa                                                                                                   <  3:338292/149‑1 (MQ=255)
                                           cGGCACGATCAACGGGTAGGCTTCTTCAGGCGTGATTATCGAAGCCTGCATCGGGATGGTCAGCCCGCGACGGTAATAGTCGGTGGTCACGTTGGCCCCGCCACCGTGAACGACCTTGCCGCCGAACAGCACGGCATCCCCGGCTTTGa                                                                                                   <  5:326297/149‑1 (MQ=255)
                                                 gATCAACGGGTAGGCTTCTTCAGGCGTGATTATCGAAGCCTGCATCGGGATGGTCAGCCCGCGACGGTAATAGTCGGTGGTCACGTTGGCCCCGCCACCGTGAACGACCTTGCCGCCGAACAGCACGGCATCCCCGGCTTTGAgcaggg                                                                                             <  4:154489/149‑1 (MQ=255)
                                                   tCAACGGGTAGGCTTCTTCAGGCGTGATTATCGAAGCCTGCATCGGGATGGTCAGCCCGCGACGGTAATAGTCGGTGGTCACGTTGGCCCCGCCACCGTGAACGACCTTGCCGCCGAACAGAACGGCATCCCCGGCTTTGAgcagggca                                                                                           <  8:200063/149‑1 (MQ=255)
                                                     aaCGGGTAGGCTTCTTCAGGCGTGATTATCGAAGCCTGCATCGGGATGGTCAGCCCGCGACGGTAATAGTCGGTGGTCACGGTGGCGGGGCCAACGTGa                                                                                                                                           <  6:262845/99‑1 (MQ=255)
                                                     aaCGGGTAGGCTTCTTCAGGCGTGATTATCGAAGCCTGCATCGGGATGGTCAGCCCGCGACGGTAATAGTCGGTGGTCACGGTGGCGGGGCCAACGTGa                                                                                                                                           >  5:262845/1‑99 (MQ=255)
                                                          gTAGGCTTCTTCAGGCGTGATTATCGAAGCCTGCATCGGGATGGTCAGCCCGCGACGGTAATAGTCGGTGGTCACGTTGGCCCCGCCACCGTGAACGACCTTGCCGCCGAACAGCACGGCATCCCCGGCTTTGAGCAGGGCAGGAATCg                                                                                    <  4:233671/149‑1 (MQ=255)
                                                                ttcttcAGGCGTGATTATCGAAGCCTGCATCGGGATGGTCAGCCCGCGACGGTAATAGTCGGTGGTCACGTTGGCCCCGCCACCGTGAACGACCTTGCCGCCGAACAGCACGGCATCCCCGGCTTTGAGCAGGGCAGGAATCGTCAtc                                                                               >  3:325136/1‑148 (MQ=255)
                                                                ttcttcAGGCGTGATTATCGAAGCCTGCATCGGGATGGTCAGCCCGCGACGGTAATAGTCGGTGGTCACGTTGGCCCCGCCACCGTGAACGACCTTGCCGCCGAACAGCACGGCATCCCCGGCTTTGAGCAGGGCAGGAATCGTCAtc                                                                               >  1:279985/1‑148 (MQ=255)
                                                                ttcttcAGGCGTGATTATCGAAGCCTGCATCGGGATGGTCAGCCCGCGACGGTAATAGTCGGTGGTCACGTTGGCCCCGCCACCGTGAACGACCTTGCCGCCGAACAGCACGGCATCCCCGGCTTTGAGCAGGGCAGGAAACGTCAtct                                                                              >  2:28350/1‑149 (MQ=255)
                                                                   ttcAGGCGTGATTATCGAAGCCTGCATCGGGATGGTCAGCCCGCGACGGTAATAGTCGGTGGTCACGTTGGCCCCGCCACCGTGAACGACCTTGCCGCCGAACAGCACGGc                                                                                                                 <  8:39935/111‑1 (MQ=255)
                                                                   ttcAGGCGTGATTATCGAAGCCTGCATCGGGATGGTCAGCCCGCGACGGTAATAGTCGGTGGTCACGTTGGCCCCGCCACCGTGAACGACCTTGCCGCCGAACAGCACGGCATCCCCGGCTTTGAGCAGGGCAGGAATCGTCATCTCTg                                                                           >  4:269523/1‑149 (MQ=255)
                                                                             aTTATCGAAGCCTGCATCGGGATGGTCAGCCCGCGACGGTAATAGTCGGTGGTCACGTTGGCCCCGCCACCGTGAACGACCTTGCCGCCGAACAGCACGGCATCCCCGGCTTTGAGCAGGGCAGGAATCGTCATCTCTGGGGTGCCGa                                                                  <  3:198983/148‑1 (MQ=255)
                                                                              ttATCGAAGCCTGCATCGGGATGGTCAGCCCGCGACGGTAATAGTCGGTGGTCACGTTGGCCCCGCCACCGTGAACGACCTTGCCGCCGAACAGCACgg                                                                                                                  <  8:330188/99‑1 (MQ=255)
                                                                              ttATCGAAGCCTGCATCGGGATGGTCAGCCCGCGACGGTAATAGTCGGTGGTCACGTTGGCCCCGCCACAGTGAACGACCTTGCCGCCGAACAGCACgg                                                                                                                  >  7:330188/1‑99 (MQ=255)
                                                                                   gAAGCCTGCATCGGGATGGTCAGCCCGCGACGGTAATAGTCGGTGGTCACGTTGGCCCCGCCACCGTGAACGACCTTGCCGCCGAACAGCACGGCATCCCCGGCTTTGAGCAGGGCAGGAATCGTCATCTCTGGGGTGCCGACAtcgtc                                                           >  3:51534/1‑149 (MQ=255)
                                                                                   gAAGCCTGCATCGGGATGGTCAGCCCGCGACGGTAATAGTCGGTGGTCACGTTGGCCCCGCCACCGTGAACGACCTTGCCGCCGAACAGCACGGCATCCCCGGCTTTGAGCAGGGCAGGAATCGTCATCTCTGGGGTGCCGACAtcgt                                                            <  4:284181/148‑1 (MQ=255)
                                                                                      gCCTGCATCGGGATGGTCAGCCCGCGACGGTAATAGTCGGTGGTCACGTTGGCCCCGCCACCGTGAACGACCTTGCCGCCGAACAGCACGGCATCCCCGGCTTTGAGCAGGGCAGGAATCGTCATCTCTGGGGTGCCGACATCGTCGaa                                                        <  2:148002/149‑1 (MQ=255)
                                                                                      gCCTGCATCGGGATGGTCAGCCCGCGACGGTAATAGTCGGTGGTCACGTTGGCCCCGCCACCGTGAACGACCTTGCCGCCGAACAGCACGGCATCCCCGGCTTTGAGCAGGGCAGGAATCGTCATCTCTGGGGTGCCGACATCGTCGaa                                                        <  5:118823/149‑1 (MQ=255)
                                                                                      gCCTGCATCGGGATGGTCAGCCCGCGACGGTAATAGTCGGTGGGCACGTTGGCCCCGCCACCGTGAACGACCTTGCCGCCGAACAGCACGGCATCCCCGGCTTTGAGCAGGGCAGGAAACGTCATCTCTGGGGTGCCGACATCGTCGaa                                                        <  1:28350/149‑1 (MQ=255)
                                                                                        cTGCATCGGGATGGTCAGCCCGCGACGGTAATAGTCGGTGGTCACGTTGGCCCCGCCACCGTGAACGACCTTGCCGCCGAACAGCACGGCATCCCCGGCTTTGAGCAGGGCAGGAATCGTCATCTCTGGGGTGCCGACATCGTCGaa                                                        <  2:258142/147‑1 (MQ=255)
                                                                                        cTGCATCGGGATGGTCAGCCCGCGACGGTAATAGTCGGTGGTCACGTTGGCCCCGCCACCGTGAACGACCTTGCCGCCGAACAGCACGGCATCCCCGGCTTTGAGCAGGGCAGGAATCGTCATCTCTGGGGTGCCGACATCGTCGaa                                                        >  1:258142/1‑147 (MQ=255)
                                                                                        cTGCATCGGGATGGTCAGCCCGCGACGGTAATAGTCGGTGGTCACGTTGGCCCCGCCACCGTGAACGACCTTGCCGCCGAACAGCACGGCATCCCCGGCTTTGAGCAGGGCAGGAATCGTCATCTCTGGGGTGCCGACATCGTCGa                                                         <  6:175034/146‑1 (MQ=255)
                                                                                        cTGCATCGGGATGGTCAGCCCGCGACGGTAATAGTCGGTGGTCACGTTGGCCCCGCCACCGTGAACGACCTTGCCGCCGAACAGCACGGCATCCCCGGCTTTGAGCAGGGCAGGAATCGTCATCTCTGGGGTGCCGACATCGTCGa                                                         >  5:175034/1‑146 (MQ=255)
                                                                                             tCGGGATGGTCAGCCCGCGACGGTAATAGTCGGTGGTCACGTTGGCCCCGCCACCGTGAACGACCTTGCCGCCGAACAGCACGGCATCCCCGGCTTTGAgcag                                                                                               >  4:294685/1‑103 (MQ=255)
                                                                                             tCGGGATGGTCAGCCCGCGACGGTAATAGTCGGTGGTCACGTTGGCCCCGCCACCGTGAACGACCTTGCCGCCGAACAGCACGGCATCCCCGGCTTTGAgcag                                                                                               <  3:294685/103‑1 (MQ=255)
                                                                                                ggATGGTCAGCCCGCGACGGTAATAGTCGGTGGTCACGTTGGCCCCGCCACCGTGAACGACCTTGCCGCCGAACAGCACGGCATCCCCGGCTTTGAGCAGGGCAGGAATCGTCATCTCTGGGGTGCCGACATCGTCGAAGTCATCCCAg                                              >  5:64714/1‑149 (MQ=255)
                                                                                                       cAGCCCGCGACGGTAATAGTCGGTGGTCACGTTGGCCCCGCCACCGTGAACGACCTTGCCGCCGAACAGCACGGCATCCCCGGCTTTGAGCAGGGCAGGAATCGTCATCTCTGGGGTGCCGACATCGTCGAAg                                                       <  5:32221/133‑1 (MQ=255)
                                                                                                       cAGCCCGCGACGGTAATAGTCGGTGGTCACGTTGGCCCCGCCACCGTGAACGACCTTGCCGCCGAACAGCACGGCATCCCCGGCTTTGAGCAGGGCAGGAATCGTCATCTCTGGGGTGCCGACATCGTCGAAg                                                       >  6:32221/1‑133 (MQ=255)
                                                                                                           ccGCGACGGTAATAGTCGGTGGTCACGTTGGCCCCGCCACCGTGAACGACCTTGCCGCCGAACAGCACGGCATCCCCGGCTTTGAGCAGGGCAGGAATCGTCATCTCTGGGGTGCCGACATCGTCGAAGTCATCCCAGTCCTGGCTccc                                   <  2:279985/149‑1 (MQ=255)
                                                                                                            cgcgACGGTAATAGTCGGTGGTCACGTTGGCCCCGCCACCGTGAACGACCTTGCCGCCGAACAGCACGGCATCCCCGGCTTTGAGCAGGGCAGGAATCGTCATCTCTGGGGTGCCGACATCGTCGAAGTCATCCCAGTCCTGGCTCCCg                                  >  6:307723/1‑149 (MQ=255)
                                                                                                             gcgACGGTAATAGTCGGTGGTCACGTTGGCCCCGCCACCGTGAACGACCTTGCCGCCGAACAGCACGGCATCCCCGGCTTTGAGCAGGGCAGGAATCGTCATCTCTGGGGTGCCGACATCGTCGAAGTCATCCCAGTCCTGGCTCCCgg                                 <  4:51534/149‑1 (MQ=255)
                                                                                                                aCGGTAATAGTCGGTGGTCACGTTGGCCCCGCCACCGTGAACGACCTTGCCGCCGAACAGCACGGCATCCCCGGCTTTGAGCAGGGCAGGAATCGTCATCTCTGGGGTGCCGACATCGTCGAAGTCATCCCAGTCCTGGCTCCCGGGGa                              <  3:269523/149‑1 (MQ=255)
                                                                                                                     aaTAGTCGGTGGTCACGTTGGCCCCGCCACCGTGAACGACCTTGCCGCCGAACAGCACGGCATCCCCGGCTTTGAGCAGGGCAGGAATc                                                                                     >  6:98870/1‑89 (MQ=255)
                                                                                                                     aaTAGTCGGTGGTCACGTTGGCCCCGCCACCGTGAACGACCTTGCCGCCGAACAGCACGGCATCCCCGGCTTTGAGCAGGGCAGGAATCGTCATCTCTGGGGTGCCGACATCGTCGAAGTCATCCCAGTCCTGGCTCCCGGGGATCAAg                         >  3:171525/1‑149 (MQ=255)
                                                                                                                      aTAGTCGGTGGTCACGTTGGCCCCGCCACCGTGAACGACCTTGCCGCCGATCAGCACGGCATCCCCGGCTTTGAGCAGGGCAGGAATCGTCATCTCTGGGGTGCCGACATCGTCGAAGGCATCCCAGTCCTGGCTCCCGGGGATCAAg                         <  7:249411/148‑1 (MQ=255)
                                                                                                                               gttCACGGTGGCGGGGCCAACGTGACCGACCTTGCCGCCGAACAGCACGGCATCCCCGGCTTTGAGCAGGGCAGGAATCGt                                                                                   <  2:250658/79‑1 (MQ=25)
                                                                                                                               gttCACGGTGGCGGGGCCAACGTGACCGACCTTGCCGCCGAACAGCACGGCATCCCCGGCTTTGAGCAGGGCAGGAATCGt                                                                                   >  1:250658/3‑81 (MQ=25)
                                                                                                                               gttCACGGTGGCGGGGCCAACGTGACCGACCTTGCCGCCGAACAGCACGGCATCCCCGGCTTTGAGCAGGGCAGGAATCGTCATCTCTGGGGTGCCGACATCGTCGAAGTCATCCCAGTCCTGGCTCCCGGGGATCAAGCGCGTGgcgc               >  2:120173/3‑149 (MQ=255)
                                                                                                                                 tCACGTTGGCCCCGCCACCGTGAACGACCTTGCCGCCGAACAGCACGGCATCCCCGGCTTTGAGCAGGGCAGGAATCGTCATCTCTGGGGTGCCGACATCGTCGAAGTCATCCCAGTCCTGGCTCCCGGGGATCAAGCGCGTGGCGCCg             >  1:18894/1‑149 (MQ=255)
                                                                                                                                     gTTGGCCCCGCCACCGTGAACGACCTTGCCGCCGAACAGCACGGCATCCCCGGCTTTGAGCAGGGCAGGAATCGTCATCTCTGGGGTGCCGACATCGTCGAAGTCATCCCAGTCCTGGCTCCCGGGGATCAAGCGCGTGGCGCCGtttt         >  6:181879/1‑149 (MQ=255)
                                                                                                                                         gCCCCGCCACCGTGAACGACCTTGCCGCCGAACAGCACGGCATCCCCGGCTTTGAGCAGGGCAGGAATCGTCATCTCTGGGGTGCCGACATCGTCGAAGTCATCCCAGTCCTGGCTCCCGGGGATCAAGCGCGTGGCGcc              >  4:236/1‑140 (MQ=255)
                                                                                                                                         gCCCCGCCACCGTGAACGACCTTGCCGCCGAACAGCACGGCATCCCCGGCTTTGAGCAGGGCAGGAATCGTCATCTCTGGGGTGCCGACATCGTCGAAGTCATCCCAGTCCTGGCTCCCGGGGATCAAGCGCGTGGCGcc              <  3:236/140‑1 (MQ=255)
                                                                                                                                         gCCCCGCCACCGTGAACGACCTTGCCGCCGAACAGCACGGCATCCCCGGCTTTGAGCAGGGCAGGAATCGTCATCTCTGGGGTGCCGACATCGTCGAAGTCATCCCAGTCCTGGCTCCCGGGGATCAAGCGCGTGGCGCCGTTTTCCTc     >  7:70639/1‑149 (MQ=255)
                                                                                                                                           cccGCCACCGTGAACGACCTTGCCGCCGGACAGCACGGCATCCCCGGCTTTGAGCAGGGCAGGAATCGTCATCTCTGGGGTGCCGACATCGTCGAAGTCATCCCAGTCCTGGCTCCCGGGGATCAAGCGCGTGGCGCCGTTTTCCTCgg   <  3:148477/149‑1 (MQ=255)
                                                                                                                                           cccGCCACCGTGAACGACCTTGCCGCCGAACAGCACGGCATCCCCGGCTTTGAGCAGGGCAGGAATCGTCATCTCTGGGGTGCCGACATCGTCGAAGTCATCCCAGTCCTGGCTCCCGGGGATCAAGCGCGTGGCGCCGTTTTCCTCgg   <  2:18894/149‑1 (MQ=255)
                                                                                                                                            ccGCCACCGTGAACGACCTTGCCGCCGAACAGCACGGCATCCCCGGCTTTGAGCAGGGCAGGAATCGTCATCTCTGGGGTGCCGACATCGTCGAAGTCATCCCAGTCCTGGCTCCCGGGGATCAAGCGCGTGGCGCCGTTTTCCTCGGt  <  4:171525/149‑1 (MQ=255)
                                                                                                                                            ccGCCACCGTGAACGACCTTGCCGCCGAACAGCACGGCATCCCCGGCTTTGAGCAGGGCAGGAATCGTCATCTCTGGGGTGCCGACATCGTCGAAGTCATCCCAGTCCTGGCTCCCGGGGATCAAGCGCGTGGCGCCGTTTTCCTCGGt  <  5:181879/149‑1 (MQ=255)
                                                                                                                                            ccGCCACCGTGAACGACCTTGCCGCCGAACAGCACGGCATCCCCGGCTTTGAGCAGGGCAGGAATCGTCATCTCTGGGGTGCCGACATCGTCGAAGTCATCCCAGTCCTGGCTCCCGGGGATCAAGCGCGTGGCGCCGTTTTCCTCGGt  <  4:325136/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                              
GGATCTTCTGTACCCGCGGAGAAAGCGTGCGCACCAGCTCGAGCGGCACGATCAACGGGTAGGCTTCTTCAGGCGTGATTATCGAAGCCTGCATCGGGATGGTCAGCCCGCGACGGTAATAGTCGGTGGTCACGTTGGCCCCGCCACCGTGAACGACCTTGCCGCCGAACAGCACGGCATCCCCGGCTTTGAGCAGGGCAGGAATCGTCATCTCTGGGGTGCCGACATCGTCGAAGTCATCCCAGTCCTGGCTCCCGGGGATCAAGCGCGTGGCGCCGTTTTCCTCGGT  >  NC_002947/3209942‑3210230

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: