Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,649,301 T→C 5.5% H631R (CAT→CGT)  PP_4959 ← diguanylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,649,3010TC 5.5% 214.5 / 2.2 72H631R (CAT→CGT) PP_4959diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base T (31/37);  new base C (2/2);  total (33/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.60e-01

GCAACCATTGGCCTTTAGCCAGGCGGCCTGCGCCGGGGTCTCCACACCCTCGGCGATCACCGTGATACGGTAGTCGGCGCACAGGCCGATCACGCTGCGAACCAACGCAGCATCTACCGTCGAATCGGGCAGACGTGCCACCAGATGACGGTCAAGCTTGAGGGTGTCGATAGGCAGGTCGC‑‑GCAGCATGCGCAACGAACAGTCACCCGCACCGAAGTCATCCAGGGCCACCCTCACGCCCAGCTCACGCAAGCGGTGTATCTGTTTGACGGCCGCGTCGATGTTATA  >  NC_002947/5649156‑5649443
                                                                                                                                                 |                                                                                                                                                
gCAACCATTGGCCTTTAGCCAGGCGGCCTGCGCCGGGGTCTCCACACCCTCGGCGATCACCGTGATACGGTAGTCGGCGCACAGGCCGATCACGCTGCGAACCAACGCAGCATCTACCGTCGAATCGGGCAGACGTGCCACCAGATGAc                                                                                                                                               >  5:163744/1‑149 (MQ=255)
gCAACCATTGGCCTTTAGCCAGGCGGCCTGCGCCGGGGTCTCCACACCCTCGGCGATCACCGTGATACGGTAGTCGGCGCACAGGCCGATCACGCTGCGAACCAACGCAGCATCTACCGTCGAATCGGGCAGACGTGCCACCAGATGAc                                                                                                                                               >  6:27861/1‑149 (MQ=255)
  aaCCATTGGCCTTTAGCCAGGCGGCCTGCGCCGGGGTCTCCACACCCTCGGCGATCACCGTGATACGGTAGTCGGCGCACAGGCCGATCACGCTGCGAACCAACGCAGCATCTACCGTCGAATCGGGCAGACGTGCCACCAGATGACgg                                                                                                                                             >  7:246936/1‑149 (MQ=255)
      aTTGGCCTTTAGCCAGGCGGCCTGCGCCGGGGTCTCCACACCCTCGGCGATCACCGTGATACGGTAGTCGGCGCACAGGCCGATCACGCTGCGAACCAACGCAGCATCTACCGTCGAATCGGGCAGACGTGCCACCAGATGACGGTCaa                                                                                                                                         <  1:80203/149‑1 (MQ=255)
          gCCTTTAGCCAGGCGGCCTGCGCCGGGGTCTCCACACCCTCGGCGATCACCGTGATACGGTAGTCGGCGCACAGGCCGATCACGCTGCGAACCAACGCAGCATCTACCGTCGAATCGGGCAGACGTGCCACCAGATGACGGTCAAGCtt                                                                                                                                     >  2:68077/1‑149 (MQ=255)
            cTTTAGCCAGGCGGCCTGCGCCGGGGTCTCCACACCCTCGGCGATCACCGTGATACGGTAGTCGGCGCACAGGCCGATCACGCTGCGAACCAACGCAGCATCTACCGTCGAATCGGGCAGACGTGCCACCAGATGACGGTCAAGCTTGa                                                                                                                                   <  2:79598/149‑1 (MQ=255)
             tttAGCCAGGCGGCCTGCGCCGGGGTCTCCACACCCTCGGCGATCACCGTGATACGGTAGTCGGCGCACAGGCCGATCACGCTGCGAACCAACGCAGCATCTACCGTCGAATCGGGCAGACGTGCCACCAGATGACGGTCAAGCTTGAg                                                                                                                                  >  5:98594/1‑149 (MQ=255)
               tAGCCAGGCGGCCTGCGCCGGGGTCTCCACACCCTCGGCGATCACCGTGATACGGTAGTCGGCGCACAGGCCGATCACGCTGCGAACCAACGCAGCATCTACCGTCGAATCGGGCAGACGTGCCACCAGATGACGGTCAAGCTTGAg                                                                                                                                  <  7:57360/147‑1 (MQ=255)
                  ccAGGCGGCCTGCGCCGGGGTCTCCACACCCTCGGCGATCACCGTGATACGGTAGTCGGCGCACAGGCCGATCACGCTGCGAACCAACGCAGCATCTACCGTCGAATCGGGCAGACGTGCCACCAGATGACGGTCAAGCTTGAGGgtgt                                                                                                                             <  7:101713/149‑1 (MQ=255)
                    aGGCGGCCTGCGCCGGGGTCTCCACACCCTCGGCGATCACCGTGATACGGTAGTCGGCGCACAGGCCGATCACGCTGCGAACCAACGCAGCATCTACCGTCGAATCGGGCAGACGTGCCACCAGATGACGGTCAAGCTTGAGGGTGTc                                                                                                                            >  4:242320/1‑148 (MQ=255)
                    aGGCGGCCTGCGCCGGGGTCTCCACACCCTCGGCGATCACCGTGATACGGTAGTCGGCGCACAGGCCGATCACGCTGCGAACCAACGCAGCATCTACCGTCGAATCGGGCAGACGTGCCACCAGATGACGGTCAAGCTTGAGGGTGTCg                                                                                                                           >  3:220276/1‑149 (MQ=255)
                      gcggcCTGCGCCGGGGTCTCCACACCCTCGGCGAGCACCGTGATACGGTAGTCGGCGCACAGGCCGATCACGCTGCGAACCAACGCAGCATCTACCGTCGAATCGGGCAGACGTGCCACCAGATGACGGTCAAGCTTGAGGGTGTCGAt                                                                                                                         <  8:131351/149‑1 (MQ=255)
                          ccTGCGCCGGGGTCTCCACACCCTCGGCGATCACCGTGATACGGTAGTCGGCGCACAGGCCGATCACGCTGCGAACCAACGCAGCATCTACCGTCGAATCGGGCAGACGTGCCACCAGATGACGGTCAAGCTTGAGGGTGTCGATAGGc                                                                                                                     <  5:242820/149‑1 (MQ=255)
                          ccTGCGCCGGGGGCTCCACACCCTCGGCGATCACCGTGATACGGTAGTCGGCGCACAGGCCGATCACGCTGCGAACCAACGCAGCATCTACCGTCGAAGCGGGCAGACGTGCCACGAGATGACGGTCAAGCTTGAGGGTGTCGATAGGc                                                                                                                     <  2:247298/149‑1 (MQ=255)
                              cgcCGGGGTCTCCACACCCTCGGCGATCACCGTGATACGGTAGTCGGCGCACAGGCCGATCACGCTGCGAACCAACGCAGCATCTACCGTCGAATCGGGCAGACGTGCCACCAGATGACGGTCAAGCTTGAGGGTGTCGATAGGCAg                                                                                                                   <  4:94969/147‑1 (MQ=255)
                              cgcCGGGGTCTCCACACCCTCGGCGATCACCGTGATACGGTAGTCGGCGCACAGGCCGATCACGCTGCGAACCAACGCAGCATCTACCGTCGAATCGGGCAGACGTGCCACCAGATGACGGTCAAGCTTGAGGGTGTCGATAGGCAg                                                                                                                   <  4:309577/147‑1 (MQ=255)
                              cgcCGGGGTCTCCACACCCTCGGCGATCACCGTGATACGGTAGTCGGCGCACAGGCCGATCACGCTGCGAACCAACGCAGCATCTACCGTCGAATCGGGCAGACGTGCCACCAGATGACGGTCAAGCTTGAGGGTGTCGATAGGCAGGt                                                                                                                 <  1:177425/149‑1 (MQ=255)
                               gcCGGGGTCTCCACACCCTCGGCGATCACCGTGATACGGTAGTCGGCGCACAGGCCGATCACGCTGCGAACCAACGCAGCATCTACCGTCGAATCGGGCAGACGTGCCACCAGATGACGGTCAAGCTTGAGGGTGTCGATAGGCAGGTc                                                                                                                >  2:116396/1‑149 (MQ=255)
                               gcCGGGGTCTCCACACCCTCGGCGATCACCGTGATACGGTAGTCGGCGCACAGGCCGATCACGCTGCGAACCAACGCAGCATCTACCGTCGAATCGGGCAGACGTGCCACCAGATGACGGTCAAGCTTGAGGGTGTCGATAGGCAGGTc                                                                                                                <  8:246936/149‑1 (MQ=255)
                                 cGGGGTCTCCACACCCTCGGCGATCACCGTGATACGGTAGTCGGCGCACAGGCCGATCACGCTGCGAACCAACGCAGCATCTACCGTCGAATCGGGCAGACGTGCCACCAGATGACGGTCAAGCTTGAGGgtgt                                                                                                                             >  6:140429/1‑134 (MQ=255)
                                 cGGGGTCTCCACACCCTCGGCGATCACCGTGATACGGTAGTCGGCGCACAGGCCGATCACGCTGCGAACCAACGCAGCATCTACCGTCGAATCGGGCAGACGTGCCACCAGATGACGGTCAAGCTTGAGGgtgt                                                                                                                             <  5:140429/134‑1 (MQ=255)
                                   gggTCTCCACACCCTCGGCGATCACCGTGATACGGTAGTCGGCGCACAGGCCGATCACGCTGCGAACCAACGCAGCATCTACCGTCGAATCGGGCAGACGTGCCACCAGATGACGGTCAAGCTTGAGGGTGTCGATAGGCAGGTcgc‑‑g                                                                                                           >  6:257625/1‑148 (MQ=255)
                                   gggTCTCCACACCCTCGGCGATCACCGTGATACGGTAGTCGGCGCACAGGCCGATCACGCTGCGAACCAACGCAGCATCTACCGTCGAATCGGGCAGACGTGCCACCAGATGACGGTCAAGCTTGAGGGTGTCGATAGGCAGGTCgc‑‑gc                                                                                                          <  6:332441/149‑1 (MQ=255)
                                          cacaCCCTCGGCGATCACCGTGATACGGTAGTCGGCGCACAGGCCGATCACGCTGCGAACCAACGCAGCATCTACCGTCGAATCGGGCAGACGTGCCACCAGATGACGGTCAAGCTTGAGGGTGTCGATAGGCAGGTCGC‑‑GCAGCAt                                                                                                     <  5:150553/147‑1 (MQ=255)
                                          cacaCCCTCGGCGATCACCGTGATACGGTAGTCGGCGCACAGGCCGATCACGCTGCGAACCAACGCAGCATCTACCGTCGAATCGGGCAGACGTGCCACCAGATGACGGTCAAGCTTGAGGGTGTCGATAGGCAGGTCGC‑‑GCAGCATgc                                                                                                   >  1:139303/1‑149 (MQ=255)
                                              cccTCGGCGATCACCGTGATACGGTAGTCGGCGCACAGGCCGATCACGCTGCGAACCAACGCAGCATCTACCGTCGAATCGGGCAGACGTGCCACCAGATGACGGTCAAGCTTGAGGGTGTCGATAGGCAGGTCGC‑‑GCAGCATGCGCaa                                                                                               <  5:27861/149‑1 (MQ=255)
                                               ccTCGGCGATCACCGTGATACGGTAGTCGGCGCACAGGCCGATCACGCTGCGAACCAACGCAGCATCTACCGTCGAATCGGGCAGACGTGCCACCAGATGACGGTCAAGCTTGAGGGTGTCGATAGGCAGGTCGC‑‑GCAGCATgc                                                                                                   >  4:114110/1‑144 (MQ=255)
                                               ccTCGGCGATCACCGTGATACGGTAGTCGGCGCACAGGCCGATCACGCTGCGAACCAACGCAGCATCTACCGTCGAATCGGGCAGACGTGCCACCAGATGACGGTCAAGCTTGAGGGTGTCGATAGGCAGGTCGC‑‑GCAGCATgc                                                                                                   <  3:114110/144‑1 (MQ=255)
                                                cTCGGCGATCACCGTGATACGGTAGTCGGCGCACAGGCCGATCACGCTGCGAACCAACGCAGCATCTACCGTCGAATCGGGCAGACGTGCCACCAGATGACGGTCAAGCTTGAGGGTGTCGATAGGCAGGTCGC‑‑GCAGCATGCGCAAc                                                                                              >  6:4704/1‑148 (MQ=255)
                                                 tCGGCGATCACCGTGATACGGTAGTCGGCGCACAGGCCGATCACGCTGCGAACCAACGCAGCATCTACCGTCGAATCGGGCAGACGTGCCACCAGATGACGGTCAAGCTTGAGGGTGTCGATAGGCAGGTCGC‑‑GCAGCATGCGCAACg                                                                                             >  1:79710/1‑148 (MQ=255)
                                                 tCGGCGATCACCGTGATACGGTAGTCGGCGCACAGGCCGATCACGCTGCGAACCAACGCAGCATCTACCGTCGAATCGGGCAGACGTGCCACCAGATGACGGTCAAGCTTGAGGGTGTCGATAGGCAGGTCGC‑‑GCAGCATGCGCAACg                                                                                             >  5:267216/1‑148 (MQ=255)
                                                      gATCACCGTGATACGGTAGTCGGCGCACAGGCCGATCACGCTGCGAACCAACGCAGCATCTACCGTCGAATCGGGCAGACGTGCCACCAGATGACGGTCAAGCTTGAGGGTGTCGATAGGCAGGTCGC‑‑gcagca                                                                                                      >  1:48850/1‑134 (MQ=255)
                                                      gATCACCGTGATACGGTAGTCGGCGCACAGGCCGATCACGCTGCGAACCAACGCAGCATCTACCGTCGAATCGGGCAGACGTGCCACCAGATGACGGTCAAGCTTGAGGGTGTCGATAGGCAGGTCGC‑‑GCAGCATGCGCAACGAACAGt                                                                                       >  7:237495/1‑149 (MQ=255)
                                                      gATCACCGTGATACGGTAGTCGGCGCACAGGCCGATCACGCTGCGAACCAACGCAGCATCTACCGTCGAATAGGGCAGACGTGCCACCAGATGACTGTCAAGCTTGAGGGTGTCGATAGGCAGGTCGC‑‑gcagca                                                                                                      <  2:48850/134‑1 (MQ=255)
                                                           ccGTGATACGGTAGTCGGCGCACAGGCCGATCACGCTGCGAACCAACGCAGCATCTACCGTCGAATCGGGCAGACGTGCCACCAGATGACGGTCAAGCTTGAGGGTGTCGATAGGCAGGTCGC‑‑GCAGCATGCGCAACGAACAGTCAccc                                                                                  >  1:241621/1‑149 (MQ=255)
                                                            cGTGATACGGTAGTCGGCGCACAGGCCGATCACGCTGCGAACCAACGCAGCATCTACCGTCGAATCGGGCAGACGTGCCACCAGATGACGGTCAAGCTTGAGGGTGTCGATAGGCAGGTCGC‑‑GCAGCATGCGCAACGAACAGTCACCCg                                                                                 >  4:44871/1‑149 (MQ=255)
                                                                  aCGGTAGTCGGCGCACAGGCCGATCACGCTGCGAACCAACGCAGCATCTACCGTCGAATCGGGCAGACGTGCCACCAGATGACGGTCAAGCTTGAGGGTGTCGATAGGCAGGTCGC‑‑GCAGCATGCGCAACGAACAGTCACCCGCACCg                                                                            >  2:99278/1‑148 (MQ=255)
                                                                     gTAGTCGGCGCACAGGCCGATCACGCTGCGAACCAACGCAGCATCTACCGTCGAATCGGGCAGACGTGCCACCAGATGACGGTCAAGCTTGAGGGTGTCGATAGGCAGGTCGC‑‑GCAGCATGCGCAACGAACAGTCACCCGCACCGAAGt                                                                        >  7:19272/1‑149 (MQ=255)
                                                                          cGGCGCACAGGCCGATCACGCTGCGAACCAACGCAGCATCTACCGTCGAATCGGGCAGACGTGCCACCAGATGACGGTCAAGCTTGAGGGTGTCGATAGGCAGGTCGC‑‑GCAGCATGCGCAACGAACAGTCACCCGCACCGaa                                                                          <  5:307734/142‑1 (MQ=255)
                                                                          cGGCGCACAGGCCGATCACGCTGCGAACCAACGCAGCATCTACCGTCGAATCGGGCAGACGTGCCACCAGATGACGGTCAAGCTTGAGGGTGTCGATAGGCAGGTCGC‑‑GCAGCATGCGCAACGAACAGTCACCCGCACCGaa                                                                          >  6:307734/1‑142 (MQ=255)
                                                                          cGGCGCACAGGCCGATCACGCTGCGAACCAACGCAGCATCTACCGTCGAATCGGGCAGACGTGCCACCAGATGACGGTCAAGCTTGAGGGTGTCGATAGGCAGGTCGC‑‑GCAGCATGCGCAACGAACAGTCACCCGCACCGAAGTCATc                                                                    <  5:4704/148‑1 (MQ=255)
                                                                             cgcACAGGCCGATCACGCTGCGAACCAACGCAGCATCTACCGTCGAATCGGGCAGACGTGCCACCAGATGACGGTCAAGCTTGAGGGTGTCGATAGGCAGGTCGC‑‑GCAGCATGCGCAACGAACAGTCACCCGCACCGAAGTCATCCAgg                                                                <  3:277655/149‑1 (MQ=255)
                                                                             cgcACAGGCCGATCACGCTGCGAACCAACGCAGCATCTACCGTCGAATCGGGCAGACGTGCCACCAGATGACGGTCAAGCTTGAGGGTGTCGATAGGCAGGTCGC‑‑GCAGCATGCGCAACGAACAGTCACCCGCACCGAAGTCATCCAgg                                                                <  5:107164/149‑1 (MQ=255)
                                                                              gcACAGGCCGATCACGCTGCGAACCAACGCAGCATCTACCGTCGAATCGGGCAGACGTGCCACCAGATGACGGTCAAGCTTGAGGGTGTCGATAGGCAGGTCGC‑‑GCAGCATGCGCAACGAACAGTCACCCGCACCGAAGTCATCCAggg                                                               <  8:19272/149‑1 (MQ=255)
                                                                              gcACAGGCCGATCACGCTGCGAACCAACGCAGCATCTACCGTCGAATCGGGCAGACGTGCCACCAGATGACGGTCAAGCTTGAGGGTGTCGATAGGCAGGTCGC‑‑GCAGCATGCGCAACGAACAGTCACCCGCACCGAAGTCATCCAggg                                                               <  4:220276/149‑1 (MQ=255)
                                                                                      cGATCACGCTGCGAACCAACGCAGCATCTACCGTCGAATCGGGCAGACGTGCCACCAGATGACGGTCAAGCTTGAGGGTGTCGATAGGCAGGTCGC‑‑GCAGCATGCGCAACGAACAGTCACCCGCACCGaa                                                                          <  8:45054/130‑1 (MQ=255)
                                                                                      cGATCACGCTGCGAACCAACGCAGCATCTACCGTCGAATCGGGCAGACGTGCCACCAGATGACGGTCAAGCTTGAGGGTGTCGATAGGCAGGTCGC‑‑GCAGCATGCGCAACGAACAGTCACCCGCACCGaa                                                                          >  7:45054/1‑130 (MQ=255)
                                                                                      cGATCACGCTGCGAACCAACGCAGCATCTACCGTCGAATCGGGCAGACGTGCCACCAGACGACGGTCAAGCTTGAGGGTGTCGATAGGCAGGTCGC‑‑GCAGCATGCGCAACGAACAGTCACCCGCACCGAAGTCATCCAGGGCCACCCt                                                        >  3:107201/1‑148 (MQ=255)
                                                                                      cGATCACGCTGCGAACCAACGCAGCATCTACCGTCGAATCGGGCAGACGTGCCACCAGACGACGGTCAAGCTTGAGGGTGTCGATAGGCAGGTCGC‑‑GCAGCATGCGCAACGAACAGTCACCCGCACCGAAGTCATCCAGGGCCACCCt                                                        >  3:261822/1‑148 (MQ=255)
                                                                                       gATCACGCTGCGAACCAACGCAGCATCTACCGTCGAATCGGGCAGACGTGCCACCAGATGACGGTCAAGCTTGAGGGTGTCGATAGGCAGGTCGC‑‑GCAGCATGCGCAACGAACAGTCCCCCGCACCGAAGTCATCCAGGGCCAACCTCa                                                      >  5:173292/1‑149 (MQ=255)
                                                                                           aCGCTGCGAACCAACGCAGCATCTACCGTCGAATCGGACAGACGTGCCACCAGATGACGGTCAAGCTTGAGGGTGTCGATAGGCAGGTCGC‑‑GCAGCATGCGCAACGAACAGTCACCCGCACCGAAGTCATCCAGGGCCACCCTCACGc                                                   <  6:163744/148‑1 (MQ=255)
                                                                                                   aaCCAACGCAGCATCTACCGTCGAATCGGGCAGACGTGCCACCAGATGACGGTCAAGCTTGAGGGTGTCGATAGGCAGGTCGC‑‑GCAGCATGCGCaa                                                                                               >  3:189072/1‑96 (MQ=255)
                                                                                                   aaCCAACGCAGCATCTACCGTCGAATCGGGCAGACGTGCCACCAGATGACGGTCAAGCTTGAGGGTGTCGATAGGCAGGTCGC‑‑GCAGCATGCGCaa                                                                                               <  4:189072/96‑1 (MQ=255)
                                                                                                    aCCAACGCAGCATCTACCGTCGAATCGGGCAGACGTGCCACCAGATGACGGTCAAGCTTGAGGGTGTCGATAGGCAGGTCGCCTGCAGCATGCGCAACGAACAGTCACCCGCACCGAAGTCATCCAGGGCCACCCTCACGCCCAGCTCa                                           <  1:99278/149‑1 (MQ=255)
                                                                                                         cGCAGCATCTACCGTCGAATCGGGCAGACGTGCCACCAGATGACGGTCAAGCTTGAGGGTGTCGATAGGCAGGTCGC‑‑GCAGCATGCGCAACGAACAGTCACCCGCACCGAAGTCATCCAGGGCCACCCTCACGCCCAGCTCACGCAAGc                                    <  3:44871/149‑1 (MQ=255)
                                                                                                          gcagcaTCTACCGTCGAATCGGGCAGACGTGCCACCAGATGACGGTCAAGCTTGAGGGTGTCGATAGGCAGGTCGC‑‑GCAGCATGCGCAACGAACAGTCACCCGCACCGAAGTCATCCAGGGCCACCCTCACGCCCAGCTCACGCAAGCg                                   <  2:79710/149‑1 (MQ=255)
                                                                                                          gcagcaTCTACCGTCGAATCGGGCAGACGTGCCACCAGATGACGGTCAAGCTTGAGGGTGTCGATAGGCAGGTCGC‑‑GCAGCATGCGCAACGAACAGTCACCCGCACCGAAGTCATCCAGGGCCACCCTCACGCCCAGCTCACGCAAGCg                                   <  5:277505/149‑1 (MQ=255)
                                                                                                          gcagcaTATACCGTCGAATCGGGCAGACGTGCCACCAGATGACGGTCAAGCTTGAGGGTGTCGATAGGCAGGTCGC‑‑GCAGCATGCGCAACGAACTGTCACCCGCACCGAAGTCATCCAGGGCCACCCTCACGCCCAGCTCACGCAAGCg                                   <  7:115025/149‑1 (MQ=255)
                                                                                                              caTCTACCGTCGAATCGGGCAGACGTGCCACCAGATGACGGTCAAGCTTGAGGGTGTCGATAGGCAGGTCGC‑‑GCAGCATGCGCAACGAACAGTCACCCGCACCGAAGTCATCCAGGGCCACCCTCACGCCCAGCTCACGCAAGCGgtgt                               <  1:26296/149‑1 (MQ=255)
                                                                                                              caTCTACCGTCGAATCGGGCAGACGTGCCACCAGATGACGGTCAAGCTTGAGGGTGTCGATAGGCAGGTCGC‑‑GCAGCATGCGCAACGAACAGTCACCCGCACCGAAGTCATCCAGGGCCACCCTCACGCCCAGCTCACGCAAGCGgtg                                >  2:26296/1‑148 (MQ=255)
                                                                                                                  tACCGTCGAATCGGGCAGACGTGCCACCAGATGACGGTCAAGCTTGAGGGTGTCGATAGGCAGGTCGC‑‑GCAGCATGCACAACGAACAGTCACCCGCACCGAAGTCATCCAGGGCCACCCTCACGCCCAGCTCACGCAAGCGGTGTATCt                           >  1:151466/1‑149 (MQ=255)
                                                                                                                      gTCGAATCGGGCAGACGTGCCACCAGATGACGGTCAAGCTTGAGGGTGTCGATAGGCAGGTCGC‑‑GCAGCATGCGCAACGAACAGTCACCCGCACCGAAGTCATCCAGGGCCACCCTCACGCCCAGCTCACGCAAGCGGTGTATCTGttt                       <  6:267216/149‑1 (MQ=255)
                                                                                                                       tCGAATCGGGCAGACGTGCCACCAGACGACGGTCAAGCTTGAGGGTGTCGATAGGCAGGTCGC‑‑GCAGCATGCGCAACGAACAGTCACCCGCACCGAAGTCATCCAGGGCCACCCTCACGCCCAGCTCACGCAAGCGGTGTATCTGTTTg                      <  4:107201/149‑1 (MQ=255)
                                                                                                                       tCGAATCGGGCAGACGTGCCACCAGACGACGGTCAAGCTTGAGGGTGTCGATAGGCAGGTCGC‑‑GCAGCATGCGCAACGAACAGTCACCCGCACCGAAGTCATCCAGGGCCACCCTCACGCCCAGCTCACGCAAGCGGTGTATCTGTTTg                      <  4:261822/149‑1 (MQ=255)
                                                                                                                          aaTCGGGCAGACGTGCCACCAGATGACGGTCAAGCTTGAGGGTGTCGATAGGCAGGTCGC‑‑GCAGCATGCGCAACGAACAGTCACCCGCACCGAAGTCATCCAGGGCCACCCTCACGCCCAGCTCACGCCAGCGGTGTATCTGTTTGACg                   >  4:73808/1‑149 (MQ=255)
                                                                                                                          aaTCGGGCAGACGTGCCACCAGATGACGGTCAAGCTTGAGGGTGTCGATAGGCAGGTCGC‑‑GCAGCATGCGCAACGAACAGTCACCCGCACCGAAGTCATCCAGGGCCACCCTCACGCCCAGCTCACGCAAGCGGTGTATCTGTTTGACg                   >  4:69466/1‑149 (MQ=255)
                                                                                                                           aTCGGGCAGACGTGCCACCAGATGACGGTCAAGCTTGAGGGTGTCGATAGGCAGGTCGC‑‑GCAGCATGCGCAACGAACAGTCACCCGCACCGAAGTCATCCAGGG‑CACCCTCACGCCCAGCTCACGCAAGCGGTGTATCTGTTTGACGGc                 <  6:98594/149‑1 (MQ=255)
                                                                                                                                     cGTGCCACCAGATGACGGTCAAGCTTGAGGGTGTCGATAGGCAGGTCGC‑‑GCAGCATGCGCAACGAACAGTCACCCGCACCGAAGTCATCCAGGGCCACCCTCACGCCCAGCTCACGCAAGCGGTGTATCTGTTTGACGGCCGCGTCGAt        <  2:241621/149‑1 (MQ=255)
                                                                                                                                         ccaccaGATGACGGTCAAGCTTGAGGGTGTCGATAGGCAGGTCGC‑‑GCAGCATGCGCAACGAAc                                                                                          >  7:302033/1‑63 (MQ=255)
                                                                                                                                         ccaccaGATGACGGTCAAGCTTGAGGGTGTCGATAGGCAGGTCGC‑‑GCAGCATGCGCAACGAAc                                                                                          <  8:302033/63‑1 (MQ=255)
                                                                                                                                         ccaccaGATGACGGTCAAGCTTGAGGGTGTCGATAGGCAGGTCGC‑‑GCAGCATGCGCAACGAACAGTCACCCGCACCGAAGTCATCCAGGGCCACCCTCACGCCCAGCTCACGCAAGCGGTGTATCTGTTTGACGGCCGCGTCGATGTta    >  7:61638/1‑149 (MQ=255)
                                                                                                                                           accaGATGACGTTCAATCTTGAGGGTGTCGATAGGCAGGTCTC‑‑TCAGCATGCGCAACGAACAGTCACCCGCACCGAAGTCATCCAGGGCCACCCTCACGCCCAGCTCACGCAAGCGGTGTATCTGTTTGACGGCCGCGTCTATGTtata  <  2:189664/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
GCAACCATTGGCCTTTAGCCAGGCGGCCTGCGCCGGGGTCTCCACACCCTCGGCGATCACCGTGATACGGTAGTCGGCGCACAGGCCGATCACGCTGCGAACCAACGCAGCATCTACCGTCGAATCGGGCAGACGTGCCACCAGATGACGGTCAAGCTTGAGGGTGTCGATAGGCAGGTCGC‑‑GCAGCATGCGCAACGAACAGTCACCCGCACCGAAGTCATCCAGGGCCACCCTCACGCCCAGCTCACGCAAGCGGTGTATCTGTTTGACGGCCGCGTCGATGTTATA  >  NC_002947/5649156‑5649443

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