Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,790,151 C→G 7.2% intergenic (+30/+38) ynfM → / ← PP_5073 efflux transporter/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,790,1510CG 7.2% 217.3 / 5.9 69intergenic (+30/+38)ynfM/PP_5073efflux transporter/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (27/37);  new base G (2/3);  total (29/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.14e-01

TGTTCTGGCACCAGTGGGGCTGGAACGGGATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCACCGGTCAGCTTCACCTTGGTCTTGTCGTTGAACGGTGTCGGCGGC  >  NC_002947/5790009‑5790297
                                                                                                                                               |                                                                                                                                                  
tgtTCTGGCACCAGTGGGGCTGGAACGGGATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTgg                                                                                                                                              <  4:296245/149‑1 (MQ=255)
  ttCTGGCACCAGTGGGGCTGGAACGGGATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTgg                                                                                                                                              >  1:324169/1‑147 (MQ=255)
     tGGCACCAGTGGGGCTGGAACGGGATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCg                                                                                                                                         <  7:119331/149‑1 (MQ=255)
       gCACCAGTGGGGCTGGAACGGGATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCgg                                                                                                                                        <  6:266499/148‑1 (MQ=255)
       gCACCAGTGGGGCTGGAACGGGATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCgg                                                                                                                                        <  2:292112/148‑1 (MQ=255)
       gCACCAGTGGGGCGGGAACGGGATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACACCCGCTCCCACAGGTTGTGGGTGCgg                                                                                                                                        <  8:299230/148‑1 (MQ=39)
        cACCAGTGGGGCTGGAACGGGATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGgt                                                                                                                                      <  3:2237/149‑1 (MQ=255)
         aCCAGTGGGGCTGGAACGGGATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGTCACACCCGCTCCCACAGGCTGTGGGTGCgggcg                                                                                                                                     <  8:290756/149‑3 (MQ=38)
          ccAGTGGGGCTGGAACGGGATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGgtgt                                                                                                                                    <  6:100688/149‑1 (MQ=255)
          ccAGTGGGGCTGGAACGGGATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGgtgt                                                                                                                                    <  4:234842/149‑1 (MQ=255)
              tGGGGCTGGAACGGGATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGGGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGAc                                                                                                                                 <  4:166202/148‑1 (MQ=255)
                    tGGAACGGGATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACg                                                                                                                             >  2:183945/1‑146 (MQ=255)
                      gAACGGGATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAagag                                                                                                                        >  6:43832/1‑149 (MQ=255)
                      gAACGGGATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGGTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAagag                                                                                                                        <  7:90621/149‑1 (MQ=255)
                       aaCGGGATCGGGCTGTTCATTGGCAGCTTGCTGGCTGTTGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCTAACCCACAACCTGTGGGAGCGGGTGTGACAACGAAGAgg                                                                                                                       <  4:162216/149‑1 (MQ=255)
                          gggATGGGGCTGTTCATTGGGAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCTCGAAGAGGCCa                                                                                                                    <  7:5494/149‑1 (MQ=255)
                           ggATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCat                                                                                                                   <  5:43832/149‑1 (MQ=255)
                                 ggCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCatcacat                                                                                                              >  4:604/1‑148 (MQ=255)
                                     gTTCATTGGCCGCTTGCTGCGGGTGCCCTTCTTGGTGGCGTGCAACCTCCGCATTTTGCAACCAAAAACGGGGTCATCCATTCGTGGCCAACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGt                                                                                                          <  1:285387/149‑1 (MQ=255)
                                        cATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAg                                                                                                                                           <  6:289536/112‑1 (MQ=255)
                                        cATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAg                                                                                                                                           >  5:289536/1‑112 (MQ=255)
                                         aTTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGtt                                                                                                         >  5:33487/1‑145 (MQ=255)
                                         aTTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGtt                                                                                                         <  6:33487/145‑1 (MQ=255)
                                         aTTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTtgatg                                                                                                     <  6:182146/149‑1 (MQ=255)
                                         aTTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTtgat                                                                                                      >  6:208442/1‑148 (MQ=255)
                                         aTTGGCAGCTTGCTGGCGGTGGCCTTGAAGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTtgat                                                                                                      >  1:203904/1‑148 (MQ=255)
                                            ggCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTgg                                                                                                                                              <  3:90178/105‑1 (MQ=255)
                                            ggCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTgg                                                                                                                                              >  4:90178/1‑105 (MQ=255)
                                             gCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGAtct                                                                                                 >  6:318175/1‑149 (MQ=255)
                                             ccAGCTTGCTGGCGGTGGCCTCGAGGGTGCCGTTGCACCGGAGCAATTTGCCACCCAAAAGGGTGGCATGCATTCGTGGGCACGCCGGCACCCACAACCTGGGGGTGGGGGGGTGACCACGAAGAGGCCATCAGATCAGTTGATGAtct                                                                                                 <  8:119023/148‑1 (MQ=255)
                                                 cTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGa                                                                                              >  8:314418/1‑148 (MQ=255)
                                                        gCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATatcc                                                                                      <  6:70668/149‑1 (MQ=255)
                                                        gCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATatc                                                                                       <  2:324169/148‑1 (MQ=255)
                                                              gCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGAtctc                                                                                                <  7:199956/133‑1 (MQ=255)
                                                              gCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGAtctc                                                                                                >  8:199956/1‑133 (MQ=255)
                                                                cTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATatccacatccac                                                                              >  3:136329/1‑149 (MQ=255)
                                                                    aTGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTcc                                                                          <  1:183945/149‑1 (MQ=255)
                                                                             ttGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGtt                                                                                                         <  8:251599/109‑1 (MQ=255)
                                                                             ttGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGtt                                                                                                         >  7:251599/1‑109 (MQ=255)
                                                                             ttGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGa                                                                                              >  7:47956/1‑120 (MQ=255)
                                                                             ttGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGa                                                                                              <  8:47956/120‑1 (MQ=255)
                                                                                 aCCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATatccacatcca                                                                               <  6:250024/131‑1 (MQ=255)
                                                                                 aCCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATatccacatcca                                                                               >  5:250024/1‑131 (MQ=255)
                                                                                        cAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTc                                                      <  5:208442/149‑1 (MQ=255)
                                                                                            ttGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGtt                                                                                                         <  7:261546/94‑1 (MQ=255)
                                                                                            ttGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGtt                                                                                                         >  8:261546/1‑94 (MQ=255)
                                                                                                  cccAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGa                                                                                              >  5:368651/1‑99 (MQ=255)
                                                                                                  cccAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGa                                                                                              <  6:368651/99‑1 (MQ=255)
                                                                                                   ccAAAACGGTGTGATGCATTCGTGGGCACGGCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCACCgg                                           >  1:28478/1‑149 (MQ=255)
                                                                                                   ccAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACCAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCACCgg                                           >  5:32649/1‑149 (MQ=255)
                                                                                                              tgATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCACCGGTCAGCTTCAc                                 <  2:28478/148‑1 (MQ=255)
                                                                                                                aTGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACAACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCACCGGTCAGCTTCACCtt                              >  2:208766/1‑149 (MQ=255)
                                                                                                                  gCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCACATCAGGTGCTTTTCAACCTCACCTGTCAGCTTCACCTTgt                            >  4:199071/1‑148 (MQ=255)
                                                                                                                   tcttCGTGGTCACACCCGCTCCCACAGGTTGTGGGTGCGGGCGTGCCCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCACCGGTCAGCTTCACCTTGGt                           >  5:118594/3‑149 (MQ=25)
                                                                                                                    aTTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCACCGGTCAGCTTCACCTTGGTc                          <  5:318175/149‑1 (MQ=255)
                                                                                                                     ttCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCACCGGTCAGCTTCACCTTGGTCt                         >  3:166269/1‑149 (MQ=255)
                                                                                                                        gTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCACCGGTCAGCTTCACCTTGGTCTTg                       <  3:604/148‑1 (MQ=255)
                                                                                                                          gggCACGCCCGCACCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCACCGGTCAGCTTCACCTAGGTCTTGTCGt                   >  1:291667/1‑149 (MQ=255)
                                                                                                                          gggCACGCCCGCACCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGATCATCAGGTGCTTTTCAACCTCACCGGTCAGCTTCACCTTGGTCTTGTCGt                   >  1:287348/1‑149 (MQ=255)
                                                                                                                            tcacaCCCGCTCCCACAGGTTGTGGGTGCGGGCGTGCCCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCACCGGTCAGCTTCACCTTGGTCTTGTCGtt                  <  6:118594/148‑1 (MQ=38)
                                                                                                                            gCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCACCGGTCAGCTTCACCTTGGTCTTGTCGtt                  <  2:203904/149‑1 (MQ=255)
                                                                                                                                    cGCTCCCACAGGTTGTGGGTGCGGGCGTGCCCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCACCGGTCAGCTTCACCTTGGTCTTGTCGTTGAACGgt           >  7:246502/1‑149 (MQ=38)
                                                                                                                                      cACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCACCGGTCAGCTTCACCTTGGTCTTGTCGTTGAACGgtgt         <  1:208766/149‑1 (MQ=255)
                                                                                                                                         ccACAACCTGTGGTAGCGGGTGTTACCACGAATAGGCCATCAGATCAGTTTATTATCTCGACAATATCCACATCCACTTCCCGTCTCATCAGGTGCTTTTCAACCTCACCGGt                                          <  1:27823/113‑1 (MQ=255)
                                                                                                                                         ccACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCACCTGTCAGCTTCACCTTGGTCTTGTCGTTGAACGGTGTcgg      >  1:248566/1‑149 (MQ=255)
                                                                                                                                           acaACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCACCGGTCAGCtt                                    >  7:5128/1‑117 (MQ=255)
                                                                                                                                           acaACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCACCGGTCAGCtt                                    <  8:5128/117‑1 (MQ=255)
                                                                                                                                             aaCCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCACCGGTCAGCTTCACCTTGGTCTTGTCGTTGAACGGTGTCggcggc  <  3:199071/149‑1 (MQ=255)
                                                                                                                                             aaCCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCACCGGTCAGCTTCACCTTGGTCTTGTCGTTGAACGGTGTCggcggc  <  7:314418/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                  
TGTTCTGGCACCAGTGGGGCTGGAACGGGATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCACCGGTCAGCTTCACCTTGGTCTTGTCGTTGAACGGTGTCGGCGGC  >  NC_002947/5790009‑5790297

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: