Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 461,030 C→G 7.7% V479L (GTG→CTG)  pqqF ← coenzyme PQQ synthesis protein F

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947461,0300CG 7.7% 280.1 / 11.3 91V479L (GTG→CTG) pqqFcoenzyme PQQ synthesis protein F
Reads supporting (aligned to +/- strand):  ref base C (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 = 4.81e-01

GCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGGCGCCAGGGCCTGTTCCAGCACACGCTGTAACGGCTGGCGCATAGGTGAAGGT  >  NC_002947/460890‑461178
                                                                                                                                              |                                                                                                                                                    
gCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCACGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGATTGTTCCACTGACCGCAGCActgctgcag                                                                                                                                              >  6:365454/1‑147 (MQ=255)
    ccccGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCa                                                                                                                                           <  6:232590/148‑1 (MQ=255)
    ccccGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCAGCAGCAGTGCTGCGGGCa                                                                                                                                           <  1:917/148‑1 (MQ=255)
    ccccGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCAGCAGCAGTGCTGCGGGCa                                                                                                                                           <  6:137725/148‑1 (MQ=255)
     cccGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAgg                                                                                                                                         <  2:190644/149‑1 (MQ=255)
      ccGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGc                                                                                                                                        >  6:359095/1‑149 (MQ=255)
      ccGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGc                                                                                                                                        <  5:359095/149‑1 (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                                                                                                                          >  3:191431/1‑127 (MQ=255)
                                          gATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAAc                                                                                                                          <  4:191431/127‑1 (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                                                                                                     >  4:36391/1‑112 (MQ=255)
                                                                              cagccagcAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCActg                                                                                                     <  3:36391/112‑1 (MQ=255)
                                                                               agccagcAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTc                                                               <  8:261446/149‑1 (MQ=255)
                                                                                 ccagcAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCt                                                             <  2:76953/149‑1 (MQ=255)
                                                                                 ccagcAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCt                                                             >  3:233420/1‑149 (MQ=255)
                                                                                 ccagcAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCt                                                             <  7:17713/149‑1 (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                                                     <  6:87786/147‑1 (MQ=255)
                                                                                           gCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGgc                                                     <  1:57638/147‑1 (MQ=255)
                                                                                               gTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCaa                                                                                                                           >  2:22773/1‑73 (MQ=255)
                                                                                               gTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCaa                                                                                                                           <  1:22773/73‑1 (MQ=255)
                                                                                               gTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGGCGCCAgg                                               <  2:217701/149‑1 (MQ=255)
                                                                                               gTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGGCGCCAgg                                               >  5:160919/1‑149 (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                  >  7:192044/5‑149 (MQ=255)
                                                                                                                            cggccTGCCAGCAGCAGTGCTGCGGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGGCGCCAGGGCCTGTTCCAGCACACGCTGTAACGGCTg                  >  4:85117/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)
                                                                                                                                              |                                                                                                                                                    
GCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGGCGCCAGGGCCTGTTCCAGCACACGCTGTAACGGCTGGCGCATAGGTGAAGGT  >  NC_002947/460890‑461178

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: