Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,467,909 A→G 5.4% S267G (AGC→GGC)  queF → NADPH‑dependent 7‑cyano‑7‑deazaguanine reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,467,9090AG 5.4% 277.0 / 4.1 93S267G (AGC→GGC) queFNADPH‑dependent 7‑cyano‑7‑deazaguanine reductase
Reads supporting (aligned to +/- strand):  ref base A (43/45);  new base G (2/3);  total (45/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 = 8.53e-01

GCATGCCGACTTCCATGAGCAGTGTGTGGAGCGGATTTACCTGGACTTGAAGAATCTGCTGCAGCCGGAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCAGGGTGGCAGTGAGGTCGGGGTGGTCGGCCTCGAAGCGCTCGATT  >  NC_002947/2467764‑2468057
                                                                                                                                                 |                                                                                                                                                    
gCATCCCGACTTCCATGAGCAGGGCGTGGAGCGCATTTACCTGGACTTGAACAATCTGCTGCACCCGGAGCACCTGACGCTCTATCCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTAGCGCAGCACCGGGCCGATCAGcc                                                                                                                                                   <  6:95191/149‑1 (MQ=255)
  aTGCCGACTTCCATGAGCAGTGTGTGGAGCGGATTTACCTGGACTTGAAGAATCTGCTGCAGCCGGAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCg                                                                                                                                                 <  3:49100/149‑1 (MQ=255)
  aTGCCGACTTCCATGAGCAGTGTGTGGAGCGGATTTACCTGGACTTGAAGAATCTGCTGCAGCCGGAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCg                                                                                                                                                 <  8:129130/149‑1 (MQ=255)
   tGCCGACTTCCATGAGCAGTGGGTGGAGCGGATTTACCTGGACTTGAAGAATCTGCTGCAGCCGGAGCACCTGACTGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCgg                                                                                                                                                <  3:177536/149‑1 (MQ=255)
     ccGACTTCCATGAGCAGTGTGTGGAGCGGATTTACCTGGACTTGAAGAATCTGCTGCAGCCGGAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGa                                                                                                                                               <  1:247992/148‑1 (MQ=255)
       gACTTCCATGAGCAGTGTGTGGAGCGGATTTACCTGGACTTGAAGAATCTGCTGCAGCCGGAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCTGATCGGCCCggtgct                                                                                                                                            <  7:334285/149‑5 (MQ=255)
        aCTTCCATGAGCAGTGTGTGGAGCGGATTTACCTGGACTTGAAGAATCTGCTGCAGCCGGAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAAc                                                                                                                                           <  4:262909/149‑1 (MQ=255)
        aCTTCCATGAGCAGTGTGTGGAGCGGATTTACCTGGACTTGAAGAATCCGCTGCAGCCGGAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAAc                                                                                                                                           <  1:291287/149‑1 (MQ=255)
             cATGAGCAGTGTGTGGAGCGGATTTACCTGGACTTGAAGAATCTGCTGCAGCCGGAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAAc                                                                                                                                       <  6:7151/148‑1 (MQ=255)
               tGAGCAGTGTGTGGAGCGGATTTACCTGGACTTGAAGAATCTGCTGCAGCCGGAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGc                                                                                                                                    <  7:222520/149‑1 (MQ=255)
                        tgtgGAGCGGATGTACCTGGACTTGAAGAATCTGCTGCAGCCGGAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCgg                                                                                                                            <  8:24447/148‑1 (MQ=255)
                           ggtgCGTATTTACCTGGACTTGATGAATCTGCTGCAGCCGGAGCACCTGACGGTCTATGCGCGCTATGTGAGGCGCGTTGGGCTGGACATCAACCCGTATCGCTGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAg                                                                                                                         <  7:248286/145‑1 (MQ=255)
                           ggAGCGGATTTACCTGGACTTGAAGAATCTGCTGCAGCCGGAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGAt                                                                                                                                              >  8:41970/1‑127 (MQ=255)
                           ggAGCGGATTTACCTGGACTTGAAGAATCTGCTGCAGCCGGAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGAt                                                                                                                                              <  7:41970/127‑1 (MQ=255)
                           ggAGCGGATTTACCTGGACTTGAAGAATCTGCTGCAGCCGGAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAg                                                                                                                         >  1:273493/1‑148 (MQ=255)
                            gAGCGGATTTACCTGGACTTGAAGAATCTGCTGCAGCCGGAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGt                                                                                                                        >  5:306804/1‑148 (MQ=255)
                              gCGGATTTACCTGGACTTGAAGAATCTGCTGCAGCCGGAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAAt                                                                                                                     <  7:258123/149‑1 (MQ=255)
                                 gATTTACCTGGACTTGAAGAATCTGCTGCAGCCGGAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATa                                                                                                                    >  7:76429/1‑147 (MQ=255)
                                      aCCTGGACTTGAAGAATCTGCTGCAGCCGGAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATAc                                                                                                                   >  1:3256/1‑143 (MQ=255)
                                      aCCTGGACTTGAAGAATCTGCTGCAGCCGGAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATAc                                                                                                                   <  2:3256/143‑1 (MQ=255)
                                       ccTGGACTTGAAGAATCTGCTGCAGCCGGAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGAt                                                                                                                                              >  2:307885/1‑115 (MQ=255)
                                        cTGGACTTGAAGAATCTGCTGCAGCCGGAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAAc                                                                                                                                           <  7:179661/117‑1 (MQ=255)
                                        cTGGACTTGAAGAATCTGCTGCAGCCGGAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAAc                                                                                                                                           >  8:179661/1‑117 (MQ=255)
                                          ggACTTGAAGAATCTGCTGCAGCCGGAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGc                                                                                                                           >  1:248025/1‑131 (MQ=255)
                                          ggACTTGAAGAATCTGCTGCAGCCGGAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGc                                                                                                                           <  2:248025/131‑1 (MQ=255)
                                          ggACTTGAAGAATCTGCTGCAGCCGGAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGc                                                                                                         >  3:37732/1‑149 (MQ=255)
                                           gACTTGAAGAATCTGCTGCAGCCGGAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCACCCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCa                                                                                                                          >  5:157231/1‑131 (MQ=255)
                                           gACTTGAAGAATCTGCTGCAGCCGGAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCACCCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCa                                                                                                                          <  6:157231/131‑1 (MQ=255)
                                             cTTGAAGAATCTGCTGCAGCCGGAGCACCTGCCGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCgg                                                                                                                            <  5:369294/127‑1 (MQ=255)
                                             cTTGAAGAATCTGCTGCAGCCGGAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCgg                                                                                                                            >  6:369294/1‑127 (MQ=255)
                                              ttGAAGAATCTGCTGCAGCCGGAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAACCGGCTGGTCCGGCAGTAATACCGGGGGCCGCTTTg                                                                                                     >  5:188747/1‑149 (MQ=255)
                                              ttGAAGAATCTGCTGCAGCCGGAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTg                                                                                                     >  6:323213/1‑149 (MQ=255)
                                                gaagaaTCTGCTGCAGCCGGAGCACCTGACGGTCTATGTGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCg                                                                                                   >  3:339735/1‑149 (MQ=255)
                                                gaagaaTCTGCTGCAGCCGGAGCACCTGACGGTCTATGCGCGCTATGTGCACCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCg                                                                                                   >  1:89830/1‑149 (MQ=255)
                                                 aagaaTCTGCTGCAGCCGGAGCACCTGACGGTCTATGTGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCg                                                                                                   <  4:339735/148‑1 (MQ=255)
                                                   gaaTCTGCTGCAGCCGGAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGcc                                                                                                >  1:113916/1‑149 (MQ=255)
                                                      tCTGCTGCAGCCGGAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGccccc                                                                                             >  7:115000/1‑149 (MQ=255)
                                                          ctgcAGCCGTCGCCCCTGCCGGTCAATGCGCGCTATGTGCGCGGCGGTGGGCTGGACATCAACCCGTATCGCACCACCGGGCCGATCAGCCCGGATAAAAAACGGCTGGTCCGGCATTAATACAGGGGCCCGCTTTGCGGCCCCCGCTa                                                                                         <  6:202711/149‑1 (MQ=255)
                                                               gCCGGAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCa                                                                                    <  2:113916/149‑1 (MQ=255)
                                                                ccGGAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAg                                                                                   <  5:323213/149‑1 (MQ=255)
                                                                ccGGAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAg                                                                                   <  7:75887/149‑1 (MQ=255)
                                                                   gAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATc                                                                                >  3:45104/1‑149 (MQ=255)
                                                                       aCCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCa                                                                            <  4:45104/149‑1 (MQ=255)
                                                                        ccTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCAt                                                                           >  6:6768/1‑149 (MQ=255)
                                                                           gACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCt                                                                        >  1:10891/1‑149 (MQ=255)
                                                                           gACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCt                                                                        >  8:263028/1‑149 (MQ=255)
                                                                            aCGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTa                                                                                                                       <  2:211435/101‑1 (MQ=255)
                                                                            aCGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTa                                                                                                                       >  1:211435/1‑101 (MQ=255)
                                                                            aCGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCAt                                                                           <  5:135839/145‑1 (MQ=255)
                                                                            aCGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCAt                                                                           >  6:135839/1‑145 (MQ=255)
                                                                               gTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCGCCGCTATTTCAGATCCCCATGCtgtg                                                                     >  3:212106/1‑148 (MQ=255)
                                                                               gTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCtgtg                                                                     >  2:17469/1‑148 (MQ=255)
                                                                                tCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCt                                                                        >  1:134075/1‑144 (MQ=255)
                                                                                tCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCt                                                                        <  2:134075/144‑1 (MQ=255)
                                                                                tCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCa                                                                   >  2:244870/1‑149 (MQ=255)
                                                                                     gcgcgcTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGa                                                               <  8:76429/148‑1 (MQ=255)
                                                                                     gcgcgcTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAg                                                              >  8:164803/1‑149 (MQ=255)
                                                                                       gcgcTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGtt                                                            >  8:316748/1‑149 (MQ=255)
                                                                                         gcttttggcggcgCGGTGGGCTGGTAATGAAACCGTATTGGATGACCGGTGCGTAGAGGCCGGTTTAAAAAAGGCTTGTCCGGGATTAATACCCGGGGCCGGTTTGGGGGCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTgg                                                          <  4:43286/142‑1 (MQ=255)
                                                                                         gcTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTgg                                                          <  4:37732/149‑1 (MQ=255)
                                                                                             tgtgCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCa                                                        <  5:6768/147‑1 (MQ=255)
                                                                                                 cgccgcGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAAtgctg                                                  <  7:263028/149‑1 (MQ=255)
                                                                                                 cgccgcGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAAtgct                                                   <  6:306804/148‑1 (MQ=255)
                                                                                                  gccgcGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAAtgctgc                                                 >  1:38332/1‑149 (MQ=255)
                                                                                                      cgGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGc                                                         >  6:63678/1‑137 (MQ=255)
                                                                                                      cgGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGc                                                         <  5:63678/137‑1 (MQ=255)
                                                                                                       ggTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCAgg                                            <  1:17469/149‑1 (MQ=255)
                                                                                                          gggCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCAGGGTg                                         <  8:115000/149‑1 (MQ=255)
                                                                                                               ggACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAACCGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCAGGGTGGCAGt                                    <  6:160592/149‑1 (MQ=255)
                                                                                                                gACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCAGGGTGGCAGTg                                   <  4:329171/149‑1 (MQ=255)
                                                                                                                 aCATCACCCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCAGGGTGGCAGt                                    >  2:223588/1‑147 (MQ=255)
                                                                                                                     cAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCAGGGTGGCAGTGAGGTc                              >  4:128937/1‑149 (MQ=255)
                                                                                                                       aCCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCAGGGTGGCAGTGAGGTCgg                            <  3:128937/149‑1 (MQ=255)
                                                                                                                       aCCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCAGGGTGGCAGTGAGGTCgg                            >  6:304872/1‑149 (MQ=255)
                                                                                                                        cccGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCAGGGTGGCAGTGAGGTCggg                           <  7:316748/149‑1 (MQ=255)
                                                                                                                         ccGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAg                                                                  <  2:210356/109‑1 (MQ=255)
                                                                                                                         ccGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCAGGGTGGCAGTGAGGTCgggg                          <  7:83263/149‑1 (MQ=255)
                                                                                                                         ccGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCAGGGTGGCAGTGAGGTCgggg                          >  8:83263/1‑149 (MQ=255)
                                                                                                                         ccGTATCGCAGCACCGGGCCGATCAGCCCCGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAg                                                                  >  1:210356/1‑109 (MQ=255)
                                                                                                                          cGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCAGGGTGGCAGTGAGGTCGGggt                         >  7:334617/1‑149 (MQ=255)
                                                                                                                              tCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCAGGGTGGCAGTGAGGTCGGGGTGGTc                     >  8:120002/1‑149 (MQ=255)
                                                                                                                                gcagcaCCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCAGGGTGGCAGTGAGGTCGGGGTGGTCgg                   >  4:108518/1‑149 (MQ=255)
                                                                                                                                   gcaCCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCAGGGTGGCAGTGAGGTCGGGGTGGTCGGCCt                <  6:322548/149‑1 (MQ=255)
                                                                                                                                       cGGGCTGATCGGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCAGGGTGGCAGTGAGGTCGGGGTGGTCGGCCTCg              <  6:53642/147‑1 (MQ=255)
                                                                                                                                       cGGGCTGATCGGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCAGGGTGGCAGTGAGGTCGGGGTGGTCGGCCTCg              <  6:48276/147‑1 (MQ=255)
                                                                                                                                       cGGGCTGATCGGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCAGGGTGGCAGTGAGGTCGGGGTGGTCGGCCTCg              >  5:53642/1‑147 (MQ=255)
                                                                                                                                       cGGGCTGATCGGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCAGGGTGGCAGTGAGGTCGGGGTGGTCGGCCTCg              >  5:48276/1‑147 (MQ=255)
                                                                                                                                         ggCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGCGGGCCGCTTTGCGGCCCCCGCTATTTCAGAACCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCAGGGTGGCAGTGAGGTCGGGGTGGTCGGCCTCGAAgc          >  4:352192/1‑149 (MQ=255)
                                                                                                                                          gCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCAGGGTGGCAGTGAGGTCGGGGTGGTCGGCCTCGAAgcg         >  2:85695/1‑149 (MQ=255)
                                                                                                                                            cGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCAGGGTGGCAGTGAGGTCGGGGTGGTc                     <  6:318368/135‑1 (MQ=255)
                                                                                                                                            cGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCAGGGTGGCAGTGAGGTCGGGGTGGTc                     >  5:318368/1‑135 (MQ=255)
                                                                                                                                            cGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCAGGGTGGCAGTGAGGTCGGGGTGGTCGGCCTCGAAGCGCt       <  7:164803/149‑1 (MQ=255)
                                                                                                                                               tCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCAGGGTGGCAGTGAGGTCGGGGTGGTCGGCCTCGAAGCGCTCGa    >  4:293811/1‑149 (MQ=255)
                                                                                                                                                 aGCCCGGATAACAACCGGCTGGTCCGGCAGTAATACCGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCAGGGTGGCAGTGAGGTCGGGGTGGGCGGCCTCGAAGCGCTCGAtt  <  6:188747/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
GCATGCCGACTTCCATGAGCAGTGTGTGGAGCGGATTTACCTGGACTTGAAGAATCTGCTGCAGCCGGAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCAGGGTGGCAGTGAGGTCGGGGTGGTCGGCCTCGAAGCGCTCGATT  >  NC_002947/2467764‑2468057

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