Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 461,037 T→G 7.6% P476P (CCA→CCC pqqF ← coenzyme PQQ synthesis protein F

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947461,0370TG 7.6% 273.8 / 11.0 91P476P (CCA→CCCpqqFcoenzyme PQQ synthesis protein F
Reads supporting (aligned to +/- strand):  ref base T (40/44);  new base G (3/4);  total (43/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.98e-01

CCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGGCGCCAGGGCCTGTTCCAGCACACGCTGTAACGGCTGGCGCATAGGTGAAGGTACATGC  >  NC_002947/460894‑461184
                                                                                                                                                 |                                                                                                                                                   
ccccGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCa                                                                                                                                                 <  6:232590/148‑1 (MQ=255)
ccccGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCAGCAGCAGTGCTGCGGGCa                                                                                                                                                 <  6:137725/148‑1 (MQ=255)
ccccGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCAGCAGCAGTGCTGCGGGCa                                                                                                                                                 <  1:917/148‑1 (MQ=255)
 cccGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAgg                                                                                                                                               <  2:190644/149‑1 (MQ=255)
  ccGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGc                                                                                                                                              <  5:359095/149‑1 (MQ=255)
  ccGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGc                                                                                                                                              >  6:359095/1‑149 (MQ=255)
  ccGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGc                                                                                                                                              >  1:217701/1‑149 (MQ=255)
      aCCGCGTCTGGCAGTTGCTTGAGCAGGGGCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCTTCAACGCCAAGGCTTGTTCCACTGCCAGCAGCAGTGCTGCGGGCAgtggaa                                                                                                                                           <  2:168440/148‑6 (MQ=255)
        cgcgTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGgcg                                                                                                                                         >  8:313022/1‑148 (MQ=255)
         gcgTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGgcg                                                                                                                                         >  7:95829/1‑147 (MQ=255)
           gTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGt                                                                                                                                     <  8:187715/149‑1 (MQ=255)
            tCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTg                                                                                                                                    >  3:264521/1‑149 (MQ=255)
                gCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAAc                                                                                                                                >  8:17713/1‑149 (MQ=255)
                gCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAAc                                                                                                                                >  8:350678/1‑149 (MQ=255)
                 cAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACt                                                                                                                               <  4:235220/149‑1 (MQ=255)
                  aGTTGCTTGAGCAGGGCCCGGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTg                                                                                                                              <  6:247838/149‑1 (MQ=255)
                   gTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGc                                                                                                                             <  4:185403/149‑1 (MQ=255)
                    ttGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGc                                                                                                                             >  3:208209/1‑148 (MQ=255)
                      gCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCa                                                                                                                           <  6:232654/148‑1 (MQ=255)
                        ttGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTa                                                                                                                        >  3:224435/1‑149 (MQ=255)
                         tGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTAt                                                                                                                       >  3:256003/1‑149 (MQ=255)
                               gggCGCTGATGGGGATGAGTGCCGGGGGGAGCGCCGTGCAGGGCAGCGAGCAACTGGCAGGTGGCTTGAGCATCTAGGCCTAGGCTTGTTCAACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCGCCg                                                                                                                 <  8:159340/149‑1 (MQ=255)
                                ggCCCTGATCCCCATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCgg                                                                                                                >  5:80197/1‑149 (MQ=255)
                                  cccTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCa                                                                                                              <  6:50335/149‑1 (MQ=255)
                                  cccTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCa                                                                                                              >  2:57638/1‑149 (MQ=255)
                                   ccTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCAc                                                                                                             >  5:279109/1‑149 (MQ=255)
                                     tGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCActg                                                                                                           >  8:220542/1‑149 (MQ=255)
                                      gATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAAc                                                                                                                                <  4:191431/127‑1 (MQ=255)
                                      gATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAAc                                                                                                                                >  3:191431/1‑127 (MQ=255)
                                       aTGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCActg                                                                                                           <  4:256003/147‑1 (MQ=255)
                                         ggggATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCTCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCActgct                                                                                                         <  1:108443/148‑1 (MQ=255)
                                             aTCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTg                                                                                                   >  2:235933/1‑149 (MQ=255)
                                               cAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTa                                                                                                                        <  3:213113/126‑1 (MQ=255)
                                               cAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTa                                                                                                                        >  4:213113/1‑126 (MQ=255)
                                                 gggCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGt                                                                                               >  3:370772/1‑149 (MQ=255)
                                                  ggCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCGGGCCGGCGCAGTGCAACt                                                                                                                               >  2:63356/1‑116 (MQ=255)
                                                  ggCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCGGGCCGGCGCAGTGCAACt                                                                                                                               <  1:63356/116‑1 (MQ=255)
                                                  ggCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGtt                                                                                              <  4:300998/149‑1 (MQ=255)
                                                    ccGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGc                                                                                            <  8:31106/149‑1 (MQ=255)
                                                       ggggCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACg                                                                                         <  2:179748/149‑1 (MQ=255)
                                                          gCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAg                                                                                      <  7:220542/149‑1 (MQ=255)
                                                                     cAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTcaca                                                                           <  5:360864/149‑1 (MQ=255)
                                                                      aGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACaa                                                                          <  8:95829/149‑1 (MQ=255)
                                                                       gggCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAAc                                                                         <  6:80197/149‑1 (MQ=255)
                                                                        ggCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACg                                                                        >  1:329634/1‑149 (MQ=255)
                                                                          cagccagcAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCActg                                                                                                           <  3:36391/112‑1 (MQ=255)
                                                                          cagccagcAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCActg                                                                                                           >  4:36391/1‑112 (MQ=255)
                                                                           agccagcAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTc                                                                     <  8:261446/149‑1 (MQ=255)
                                                                             ccagcAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCt                                                                   <  2:76953/149‑1 (MQ=255)
                                                                             ccagcAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCt                                                                   <  7:17713/149‑1 (MQ=255)
                                                                             ccagcAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCt                                                                   >  3:233420/1‑149 (MQ=255)
                                                                             ccagcAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCt                                                                   <  7:350678/149‑1 (MQ=255)
                                                                               agcAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGc                                                                 <  6:279109/149‑1 (MQ=255)
                                                                                 cAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCTGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAg                                                               <  4:208209/149‑1 (MQ=255)
                                                                                 cAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCa                                                                                                 <  8:57251/115‑1 (MQ=255)
                                                                                 cAACTGGCAGGTGGCTTGAGCAACAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCa                                                                                                 >  7:57251/1‑115 (MQ=255)
                                                                                     tGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGgc                                                           >  2:349534/1‑149 (MQ=255)
                                                                                       gCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGgc                                                           <  1:57638/147‑1 (MQ=255)
                                                                                       gCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGgc                                                           <  6:87786/147‑1 (MQ=255)
                                                                                           gTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCaa                                                                                                                                 >  2:22773/1‑73 (MQ=255)
                                                                                           gTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCaa                                                                                                                                 <  1:22773/73‑1 (MQ=255)
                                                                                           gTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGGCGCCAgg                                                     >  5:160919/1‑149 (MQ=255)
                                                                                           gTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGGCGCCAgg                                                     <  2:217701/149‑1 (MQ=255)
                                                                                               cTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGGCGCCAGGGCCt                                                 >  3:195325/1‑149 (MQ=255)
                                                                                               cTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGGCGCCAGGGCCt                                                 >  6:164982/1‑149 (MQ=255)
                                                                                                ttGAGCATCAAGGCCAAGGCTTGTTCCACGTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGGCGCCAGGGCCt                                                 >  1:28981/1‑149 (MQ=255)
                                                                                                ttGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGGCGCCAGGGCCTg                                                >  5:161420/1‑149 (MQ=255)
                                                                                                  gAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCAGTACAGCAGtgccg                                                                                           >  2:53347/1‑100 (MQ=255)
                                                                                                  gAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCAGTACAGCAGtgccg                                                                                           <  1:53347/104‑5 (MQ=255)
                                                                                                    gCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGGCGCCAGGGCCTGTTc                                             >  3:166639/1‑148 (MQ=255)
                                                                                                     cATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGGCGCCAGGGCCTGTTCCa                                           <  1:205263/149‑1 (MQ=255)
                                                                                                     cATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGGCGCCAGGGCCTGTTCCa                                           <  4:233420/149‑1 (MQ=255)
                                                                                                     cATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGGCGCCAGGACCTGTTCCa                                           <  2:329634/149‑1 (MQ=255)
                                                                                                      aTCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGGCGCCAGGGCCTGTTCCAg                                          >  1:89488/1‑149 (MQ=255)
                                                                                                        caaggccaaggcTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGGCGCCAGGGCCTGTTCCAGca                                        >  5:143557/1‑149 (MQ=255)
                                                                                                              caaggcTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGGCGCCAGGGCCTGTTCCAGCACACGCt                                  >  2:275622/1‑149 (MQ=255)
                                                                                                               aaggcTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGGCGCCAGGGCCTGTTCCAGCACACGCTg                                 >  6:326054/1‑149 (MQ=255)
                                                                                                                  gcTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGGCGCCAGGGCCTGTTCCAGCACACGCTGTaa                              >  3:308185/1‑149 (MQ=255)
                                                                                                                    ttGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGGCGCCAGGGCCTGTTCCAGCACACGCTGTAACg                            >  5:126947/1‑149 (MQ=255)
                                                                                                                      gTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGGCGCCAGGGCCTGTTCCAGCACACGCTGTAACgg                           <  2:28981/148‑1 (MQ=255)
                                                                                                                       ttCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGGCGCCAGGGCCTGTTCCAGCACGCGCTGTACCGGCt                         >  4:272732/1‑149 (MQ=255)
                                                                                                                       ttCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGGCGCCAGGGCCTGTTCCAGCACACGCTGTAACGGCt                         <  5:326054/149‑1 (MQ=255)
                                                                                                                        cggccTGCCAGCAGCAGTGCTGCGGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGGCGCCAGGGCCTGTTCCAGCACACGCTGTAACGGCTg                        >  4:85117/5‑149 (MQ=255)
                                                                                                                        cggccTGCCAGCAGCAGTGCTGCGGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGGCGCCAGGGCCTGTTCCAGCACACGCTGTAACGGCTg                        >  7:192044/5‑149 (MQ=255)
                                                                                                                           aCTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGGCGCCAGGGCCTGTTCCAGCACACGCTGTAACGGCTGgcg                     <  4:166639/149‑1 (MQ=255)
                                                                                                                              tGCCAGCAGCAGTGCTGCGGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGGCGCCAGGGCCTGTTCCAGCACACGCTGTAACgg                           >  4:73682/1‑141 (MQ=255)
                                                                                                                              tGCCAGCAGCAGTGCTGCGGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGGCGCCAGGGCCTGTTCCAGCACACGCTGTAACgg                           <  3:73682/141‑1 (MQ=255)
                                                                                                                                       gcaCTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGGCGCCAGGGCCTGTTCCAGCACACGCTGTAACGGCTGGCGCATAGGTGaa           <  4:308185/149‑1 (MQ=255)
                                                                                                                                        caCTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGGCGCCAGGGCCTGTTCCAGCACACGCTGTAACGGCTGGCGCATAGGTGAAg          <  7:13434/149‑1 (MQ=255)
                                                                                                                                        caCTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGGCGCCAGGGCCTGTTCCAGCACACGCTGTAACGGCTGGCGCATAGGTGAAg          <  4:195325/149‑1 (MQ=255)
                                                                                                                                         aCTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGGCGCCCGGGCCTGTTCCAGCACACGCTGTAACGGCTGGCGCATAGGTGAAgg         >  4:37845/1‑149 (MQ=255)
                                                                                                                                          ctgctgctgGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGGCGCCAGGGCCTGTTCCAGCACACGCTGTAACGGCTGGCGCATAGGTGAAGGt        <  1:275622/149‑1 (MQ=255)
                                                                                                                                             ctgctgGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGGCGCCAGGGCCTGTTCCAGCACACGCTGTAACGGCTGGCGCATAGGTGAAGGTAc      <  2:89488/148‑1 (MQ=255)
                                                                                                                                              tgctgGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGGCGCCAGGGCCTGTTCCAGCACACGCTGTAACGGCTGGCGCATAGGTGAAGGTACAt    <  7:313022/149‑1 (MQ=255)
                                                                                                                                                ctgGCAGGCCGGCGGAGTGCAACTGCCGGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGGCGCCAGGGCCTGTTCCAGCACACGCTGTAACGGCTGGCGCATAGGTGAAGGTACATGc  >  7:319216/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
CCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGGCGCCAGGGCCTGTTCCAGCACACGCTGTAACGGCTGGCGCATAGGTGAAGGTACATGC  >  NC_002947/460894‑461184

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: