Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,010,934 A→T 8.2% R43R (CGT→CGA)  PP_0871 ← glycine betaine/carnitine/choline ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,010,9340AT 8.2% 158.9 / 5.2 61R43R (CGT→CGA) PP_0871glycine betaine/carnitine/choline ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base A (24/32);  new base T (2/3);  total (26/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.48e-01

AACTGCATGAAGCGCTCGGCTTTGTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAG‑CCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCG‑CGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAGGATGACGAGCACGGCGGCCCATCCCAACAGCCCCTTTCCGTAGCGTTTAGCCACAGGTGGACTCCTTGTGTAGTCGGCGAGCACACTTCTTGGTTGCCGGC  >  NC_002947/1010788‑1011079
                                                                                                                                                   |                                                                                                                                                  
aaCTGCATGAAGCGCTCGGCTTTGTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCAATCGAAG‑CCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAAcg                                                                                                                                                  >  8:28943/1‑149 (MQ=255)
  cTGCATGAAGCGCTCGGCTTTGTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAG‑CCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAAcgcg                                                                                                                                                <  6:34105/149‑1 (MQ=255)
      aTGAAGCGGTCGGATTTTTCTACCCGGTGCGTCCTACTTAACGAAATGCCGGCGGGTATGCCAACGGCAAGCGCCGCCGCCATCTAAGTCCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCG‑CGCGTTATcgcgacga                                                                                                                                            <  1:99103/149‑5 (MQ=255)
        gAAGCGCTCGGCTTTGTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAG‑CCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCG‑CGATCg                                                                                                                                         >  1:187906/1‑149 (MQ=255)
         aaGCGCTCGGCTTTGTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAG‑CCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCG‑CGCGTTATCGCGACGATctt                                                                                                                                        <  2:69247/149‑3 (MQ=255)
          aGCGCTCGGCTTTGTCGACCCGGCGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAG‑CCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCG‑CGATCgtg                                                                                                                                       >  5:228521/1‑149 (MQ=255)
           gcgcTCGGCTTTGTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAG‑CCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCG‑CGATCgtgt                                                                                                                                      >  5:336287/1‑149 (MQ=255)
           gcgcTCGGCTTTGTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAG‑CCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCG‑CGATCgtgt                                                                                                                                      <  6:63230/149‑1 (MQ=255)
           gcgcTCGGCTTTGTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAG‑CCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCG‑CGATCgtgt                                                                                                                                      <  7:30448/149‑1 (MQ=255)
                 ggCTTTGTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAG‑CCAGCACCAATACCAGGTGCGCTTGCAGGGAGAACCCAAGATCGTCGCGATAACGCG‑CGATCGTGTCGAt                                                                                                                                  <  2:39202/147‑1 (MQ=255)
                   cTTTGTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAG‑CCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCG‑CGATCGTGTCGATGCCg                                                                                                                              <  5:343944/149‑1 (MQ=255)
                   cTTTGTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAG‑CCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCG‑CGATCGTGTCGATGCCg                                                                                                                              <  7:28943/149‑1 (MQ=255)
                      tGTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAG‑CCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCG‑CGATCGTGTCGATGCCGATc                                                                                                                           >  3:281436/1‑149 (MQ=255)
                       gTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACTGCCAGCGCCGCCGCCATCGAAG‑CCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGTTCGTCGCGATAACGCG‑CGATCGTGTCGATGCCGATc                                                                                                                           <  1:203709/148‑1 (MQ=255)
                         cGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAG‑CCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCG‑CGCGTTATCGCGACGATCttg                                                                                                                                       >  2:75361/1‑134 (MQ=255)
                         cGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAG‑CCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCG‑CGCGTTATCGCGACGATCttg                                                                                                                                       <  1:75361/134‑1 (MQ=255)
                           aCCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAG‑CCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCG‑CGATCGTGTCGATGCCGATCCAGTg                                                                                                                      <  1:118393/149‑1 (MQ=255)
                                gTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAG‑CCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCG‑CGATCGTGTCGATGCCGATCCAGTGGAc                                                                                                                   >  1:178045/1‑147 (MQ=255)
                                    ggTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAG‑CCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCG‑CGATCGTGTCGATGCCGATCCAGTGGACcagcag                                                                                                             <  2:187906/149‑1 (MQ=255)
                                    ggTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAG‑CCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCG‑CGATCGTGTCGATGCCGATCCAGTGGACcagcag                                                                                                             <  6:336287/149‑1 (MQ=255)
                                        gACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAG‑CCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCG‑CGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGc                                                                                                          <  6:191276/148‑1 (MQ=255)
                                          cTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAG‑CCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCG‑CGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAg                                                                                                       <  4:281436/149‑1 (MQ=255)
                                             aaCGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAG‑CCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCG‑CGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAgg                                                                                                      <  1:300216/147‑1 (MQ=255)
                                              aCGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAG‑CCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCG‑CGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAGGAt                                                                                                    <  2:103406/148‑1 (MQ=255)
                                              aCGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAG‑CCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCG‑CGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAGGAt                                                                                                    <  1:21547/148‑1 (MQ=255)
                                                    tGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAG‑CCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCG‑CGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAGGATGACGAGc                                                                                             <  6:228521/149‑1 (MQ=255)
                                                                   ccACGGCCAGCGCCGCCGCCATCGAAG‑CCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCG‑CGATCgt                                                                                                                                        <  3:322621/91‑1 (MQ=255)
                                                                   ccACGGCCAGCGCCGCCGCCATCGAAG‑CCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCG‑CGATCgt                                                                                                                                        >  4:322621/1‑91 (MQ=255)
                                                                   ccACGGCCAGCGCCGCCGCCATCGAAG‑CCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCG‑CGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAGGATGACGAGCACGGCGGCCCATccc                                                                              >  6:184047/1‑149 (MQ=255)
                                                                   ccACGGCCAGCGCCGCCGCCATCGAAG‑CCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCG‑CGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAGGATGACGAGCACGGCGGCCCATcc                                                                               >  1:77276/1‑148 (MQ=255)
                                                                    cACGGCCAGCGCCGCCGCCATCGAAG‑CCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCG‑CGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAGGATGACGAGCACGGCGGCCCATccc                                                                              >  7:64887/1‑148 (MQ=255)
                                                                          cAGCGCCGCCGCCATCGAAG‑CCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCG‑CGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAGGATGACGAGCACGGCGGCCCATCCCAACAGcc                                                                       >  1:151757/1‑149 (MQ=255)
                                                                                  ccgccATCGAAG‑CCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCG‑CGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAGGATGACGAGCACGGCGGCCCATCCCAACAGCCCCTTTCCg                                                               <  5:215074/149‑1 (MQ=255)
                                                                                   cgccATCGAAG‑CCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCG‑CGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAGGATGACGAGCACGGCGGCGCTTCCGAACAGCCCCTTTCCg                                                               >  8:70971/1‑148 (MQ=255)
                                                                                           aaG‑CCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCG‑CGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAGGATGACGAGCACGGCGGCCCATCCCAACAGCCCCTTTCCGTAGCGTTTa                                                      <  2:155692/149‑1 (MQ=255)
                                                                                             g‑CCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCG‑CGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAGGATGACGAGCACGGCGGCCCATCCCAACAGCCCCTTTCCGTAGCGTTTAGc                                                    <  5:184047/149‑1 (MQ=255)
                                                                                                cAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCG‑CGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAGGATGACGAGCACGGCGGCCCATCCCAACAGCCCCTTTCCGTAGCGTTTAGCca                                                  <  8:64887/149‑1 (MQ=255)
                                                                                                      cAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCG‑CGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAGGATGACGAGCACGGCGGCCCATCCCAACAGCCCCTTTCCGTAGCGTTTAGCCACAGGTg                                            >  4:308015/1‑149 (MQ=255)
                                                                                                      cAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCG‑CGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAGGATGACGAGCACGGCGGCCCATCCCAACAGCCCCTTTCCGTAGCGTTTAGCCACAGGTg                                            <  2:151757/149‑1 (MQ=255)
                                                                                                      cAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCG‑CGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAGGATGACGAGCACGGCGGCCCATCCCAACAGCCCCTTTCCGTAGCGTTTAGCCACAGGTg                                            >  4:271262/1‑149 (MQ=255)
                                                                                                            cAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCG‑CGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAGGATGACGAGCACGGCGGCCCATCCCAACAGCCCCTTTCCGTAGCGTTTAGCCACAGGTGGACTcc                                      >  3:286201/1‑149 (MQ=255)
                                                                                                              ggTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCG‑CGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAGGATGACGAGCACGGCGGCCCATCCCAACAGCCCCTTTCCGTAGCGTTTAGCCACAGGTGGACTCCtt                                    <  2:178045/149‑1 (MQ=255)
                                                                                                                tccgcTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCG‑CGATCGGGTCGATGCCGATCCAGTGGACCAGCAGGGCCAGGATGACGAGCACGGCGGCCCATCCCAACAGCCCCTTTCCGTAGCGTTTAGCCACAGGTGGACTCCTtgt                                  <  7:70971/147‑1 (MQ=255)
                                                                                                                tGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCG‑CGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAGGATGACGAGCACGGCGGCCCATCCCAACAGCCCCTTTCCGTAGCGTTTAGCCACAGGTGGACTCCTtgt                                  <  2:77276/149‑1 (MQ=255)
                                                                                                                         aGGTAGAACCCAAGATCGTCGCGATAACGCG‑CGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAGGATGACGAGCACGGCGGCCCATCCCAACAGCCCCTTTCCGTAGCGTTTAGCCACAGGTGGACTCCTTGTGTAGTCGGc                         >  2:125512/1‑149 (MQ=255)
                                                                                                                             aGAACCCAAGATCGTCGCGATAACGCG‑CGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAGGATGACGAGCACGGCGGCCCATCCCAACAGCCCCTTTCCGTAGCGTTTAGCCACAGGTg                                            >  5:240600/1‑126 (MQ=255)
                                                                                                                             aGAACCCAAGATCGTCGCGATAACGCG‑CGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAGGATGACGAGCACGGCGGCCCATCCCAACAGCCCCTTTCCGTAGCGTTTAGCCACAGGTg                                            <  6:240600/126‑1 (MQ=255)
                                                                                                                             aGAACCCAAGATCGTCGCGATAACGCG‑CGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAGGATGACGAGCACGGCGGCCCATCCCAACAGCCCCTTTCCGTAGCGTTTAGCCACAGGTGGa                                          >  5:176009/1‑128 (MQ=255)
                                                                                                                             aGAACCCAAGATCGTCGCGATAACGCG‑CGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAGGATGACGAGCACGGCGGCCCATCCCAACAGCCCCTTTCCGTAGCGTTTAGCCACAGGTGGa                                          <  6:176009/128‑1 (MQ=255)
                                                                                                                                  ccAAGATCGTCGCGATAACGCG‑CGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAGGATGACGAGCACGGCGGCCCATCCCAACAGCCCCTTTCCGTAGCGTTTAGCCACAGGTGGACTCCTTGTGTAGTCGGCGAGCACAct                >  8:69675/1‑149 (MQ=255)
                                                                                                                                  ccAAGATCGTCGCGATAACGCG‑CGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAGGATGACGAGCACGGCGGCCCATCCCAACAGCCCCTTTCCGTAGCGTTTAGCCACAGGTGGACTCCTTGTGTAGTCGGCGAGCACAct                >  4:129130/1‑149 (MQ=255)
                                                                                                                                    aaGATCGTCGCGATAACGCG‑CGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAGGATGACGAGCACGGCGGCCCATCCCAACAGCCCCTTTCCGTAGCGTTTAGCCACAGGTGGACTCCTTGTGTAGTCGGCGAGCACActtc              >  5:312859/1‑149 (MQ=255)
                                                                                                                                      gATCGTCGCGATAACGCG‑CGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAGGATGACGAGCACGGCGGCCCATCCCAACAGCCCCTTTCCGTAGCGTTTAGCCACAGGTGGACTCCTTGTGTAGTCGGCGAGCACActtctt            <  4:286201/149‑1 (MQ=255)
                                                                                                                                       aTCGTCGCGATAACGCG‑CGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAGGATGACGAGCACGGCGGCCCATCCCAACAGCCCCTTTCCGTAGCGTTTAGCCACAGGTGGACTCCTTGTGTAGTCGGCGAGCACACTTCTTg           <  3:271262/149‑1 (MQ=255)
                                                                                                                                       aTCGTCGCGATAACGCG‑CGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAGGATGACGAGCACGGCGGCCCATCCCAACAGCCCCTTTCCGTAGCGTTTAGCCACAGGTGGACTCCTTGTGTAGTCGGCGAGCACACTTCTTg           >  6:157958/1‑149 (MQ=255)
                                                                                                                                       aTCGTCGCGATAACGCG‑CGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAGGATGACGAGCACGGCGGCCCATCCCAACAGCCCCTTTCCGTAGCGTTTAGCCACAGGTGGACTCCTTGTGTAGTCGGCGAGCACACTTCTTg           <  5:70405/149‑1 (MQ=255)
                                                                                                                                       aTCGTCACGATAACGCG‑CGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAGGATGACGAGCACGGCGGCCCATCCCAACAGCCCCTTTCCGTAGCGTTTAGCCACAGGTGGACTCCTTGTGTAGTCGGCGAGCACACTTCTTg           >  5:252869/1‑149 (MQ=255)
                                                                                                                                          cgcgcgTTATCGCGACGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAGGATGACGAGCACGGCGGCCCATCCCAACAGCCCCTTTCCGTAGCGTTTAGCCACAGGTGGACTCCTTGTGTAGTCGGCGAGCACACTTCTTGGt         >  3:210059/3‑149 (MQ=255)
                                                                                                                                            cgcgATAACGCG‑CGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAGGATGACGAGCACGGCGGCCCATCCCAACAGCCCCTTTCCGTAGCGTTTAGCCACAGGTGGACTCCTTGTGTAGTCGGCGAGCACACTTCTTGGTTGc      >  8:251149/1‑149 (MQ=255)
                                                                                                                                              cgATAACGCG‑CGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAGGATGACGAGCACGGCGGCCCATCCCAACAGCCCCTTTCCGTAGCGTTTAGCCACAGGTGGACTCCTTGTGTAGTCGGCGAGCACACTTCTTGGTTGCCg    <  1:125512/149‑1 (MQ=255)
                                                                                                                                                aTAACGCG‑CGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAGGATGACGAGCACGGCGGCCCATCCCAACAGCCCCTTTCCGTAGCGTTTAGCCACAGGTGGACTCCTTGTGTAGTCGGCGAGCACACTTCTTGGTTGCCGGc  <  5:157958/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
AACTGCATGAAGCGCTCGGCTTTGTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAG‑CCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCG‑CGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAGGATGACGAGCACGGCGGCCCATCCCAACAGCCCCTTTCCGTAGCGTTTAGCCACAGGTGGACTCCTTGTGTAGTCGGCGAGCACACTTCTTGGTTGCCGGC  >  NC_002947/1010788‑1011079

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