Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,833,371 C→G 5.5% P228P (CCG→CCC PP_4246 ← substrate‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,833,3710CG 5.5% 233.7 / 2.3 72P228P (CCG→CCCPP_4246substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (34/34);  new base G (2/2);  total (36/36)
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.58e-01

ACTGCAGGGAACGTTGCCGCAACACGGGAATTTAGTCTTGCGGTGCAGCCACTTTGATTCGGTACGGCTCGAAGATGCGCGCCAGTTCACCGCTGTCACGCAGGCCACGCAGCAGGGCAGCGAACGACTGCTCGCCGATAGGCGCCCCCGGGCGCAGCAAGGCGTAGTGGTGATAGACCTGGTCAATGCGTTGTGAAGGCACCAGCTGCGCCCTGCTCGCCGGGTTGCGTGCCAGGAAATCGCTGAGGTAGGAACGCGTGACCAGGGCAATATCGGCGCGGCCGGCCTGCACCATC  >  NC_002947/4833223‑4833518
                                                                                                                                                     |                                                                                                                                                    
acTGCAGGGAACGTTGCCGCAACACGGGAATTTAGTCTTGCGGTGCAGCCACTTTGATTCGGTACGGCTCGAAGATGCGCGCCAGTTCACCGCTGTCACGCAGGCCACGCAGCAGGGCAGCGAACGACTGCTCGCCGATAGGCGccccc                                                                                                                                                      >  8:273666/1‑149 (MQ=255)
acTGCAGGGAACGTTGCCGCAACACGGGAATTTAGTCTTGCGGTGCAGCCACTTTGATTCGGTACGGCTCGAAGATGCGCGCCAGTTCACCGCTGTCACGCAGGCCACGCAGCAGGGCAGCGAACGACTGCTCGCCGATAGGCGccccc                                                                                                                                                      >  7:196509/1‑149 (MQ=255)
      gggAACGTTGCCGCAACACGGGAATTTAGTCTTGCGGTGCAGCCACTTTGATTCGGTACGGCTCGAAGATGCGCGCCAGTTCACCGCTGTCACGCAGGCCACGCAGCAGGGCAGCGAACGACTGCTCGCCGATAGGCGCCCCCGGgcgc                                                                                                                                                <  1:70078/149‑1 (MQ=255)
        gAACGTTGCCGCAACACGGGAATTTAGTCTTGCGGTGCAGCCACTTTGATTCGGTACGGCTCGAAGATGCGCGCCAGTTCACCGCTGTCACGCAGGCCACGCAGCAGGGCAGCGAACGACTGCTCGCCGATAGGCGCCCCCGGGCgcag                                                                                                                                              <  5:65454/149‑1 (MQ=255)
           cGTTGCCGCAACACGGGAATTTAGTCTTGCGGTGCAGCCACTTTGATTCGGTACGGCTCGAAGATGCGCGCCAGTTCACCGCTGTCACGCAGGCCACGCAGCAGGGCAGCGAACGACTGCTCGCCGATAGGCGCCCCCGGGCGCAGCaa                                                                                                                                           <  6:4454/149‑1 (MQ=255)
                 cGCAACACGGGAATTTAGTCTTGCGGTGCAGCCACTTTGATTCGGTACGGCTCGAAGATGCGCGCCAGTTCACCGCTGTCACGCAGGCCACGCAGCAGGGCAGCGAACGACTGCTCGCCGATAGGCGCCCCCGGGCGCAGCAAGGCg                                                                                                                                      >  5:21112/1‑147 (MQ=255)
                    aaCACGGGAATTTAGTCTTGCGGTGCAGCCACTTTGATTCGGTACGGCTCGAAGATGCGCGCCAGTTCACCGCTGTCACGCAGGCCACGCAGCAGGGCAGCGAACGACTGCTCGCCGATAGGCGCCCCCGGGCGCAGCAAGGCGTAgtg                                                                                                                                 <  7:29739/149‑1 (MQ=255)
                       acGGGAATTTAGTCTTGCGGTGCAGCCACTTTGATTCGGTACGGCTCGAAGATGCGCGCCAGTTCACCGCTGTCACGCAGGCCACGCAGCAGGGCAGCGAACGACTGCTCGCCGATAGGCGCCCCCGGGCGCAGCAAGGCGTAgtggt                                                                                                                               >  8:16273/1‑148 (MQ=255)
                             aTTTAGTCTTGCGGTGCAGCCACTTTGATTCGGTACGGCTCGAAGATGCGCGCCAGTTCACCGCTGTCACGCAGGCCACGCAGCAGGGCAGCGAACGACTGCTCGCCGATAGGCGCCCCCGGGCGCAGCAAGGCGTAGTGGTGATAGAc                                                                                                                        <  4:74028/149‑1 (MQ=255)
                              tttAGTCTTGCGGTGCAGCCACTTTGATTCGGTACGGCTCGAAGATGCGCGCCAGTTCACCGCTGTCACGCAGGCCACGCAGCAGGGCAGCGAACGACTGCTCGCCGATAGGCGCCCCCGGGCGCAGCAAGGc                                                                                                                                        <  1:165678/133‑1 (MQ=255)
                              tttAGTCTTGCGGTGCAGCCACTTTGATTCGGTACGGCTCGAAGATGCGCGCCAGTTCACCGCTGTCACGCAGGCCACGCAGCAGGGCAGCGAACGACTGCTCGCCGATAGGCGCCCCCGGGCGCAGCAAGGc                                                                                                                                        >  2:165678/1‑133 (MQ=255)
                                  gTCTTGCGGTGCAGCCACTTTGATGCGGTACGGCTCGAAGATGCGCGCCAGTTCACCGCTGTCACGCAGGCCACGCAGCAGGGCAGCGAACGACTGCTCGCCGATAGGCGCCCCCGGGCGCAGCAAGGCGTAGTGGTGATAGACCTGGt                                                                                                                   <  5:272047/149‑1 (MQ=255)
                                     ttGCGGTGCAGCCACTTTGATTCGGTACGGCTCGAAGATGCGCGCCAGTTCACCGCTGTCACGCAGGCCACGCAGCAGGGCAGCGAACGACTGCTCGCCGATAGGCGCCCCCGGGCGCAGCAAGGCGTAGTGGTGATAGACCTGGTCaa                                                                                                                >  4:204955/1‑149 (MQ=255)
                                                ccACTTTGATTCGGTACGGCTCGAAGATGCGCGCCAGTTCACCGCTGTCACGCAGGCCACGCAGCAGGGCAGCGAACGACTGCTCGCCGATAGGCGCCCCCGGGCGCAGCAAGGCGTAGTGGTGATAGACCTGGTCAATGCGTTGTGaa                                                                                                     <  1:355004/149‑1 (MQ=255)
                                                  aCTTTGATTCGGTACGGCTCGAAGATGCGCGCCAGTTCACCGCTGTCACGCAGGCCACGCAGCAGGGCAGCGAACGACTGCTCGCCGATAGGCGCCCCCGGGCGCAGCAAGGCGTAGTGGTGATAGACCTGGTCAATGCGTTGTGAAgg                                                                                                   <  5:89315/149‑1 (MQ=255)
                                                    tttGATTCGGTACGGCTCGAAGATGCGCGCCAGTTCACCGCTGTCACGCAGGCCACGCAGCAGGGCAGCGAACGACTGCTCGCCGATAGGCGCCCCCGGGCGCAGCAAGGCGTAGTGGTGATAGACCTGGTCAATGCGTTGTGAAGGCa                                                                                                 >  5:203556/1‑149 (MQ=255)
                                                        aTTCGGTACGGCTCGAAGATGCGCGCCAGTTCACCGCTGTCACGCAGGCCACGCAGCAGGGCAGCGAACGACTCCTCGCCGATAGGCGCCCCCGGGCGCAGCAAGGCGTAGTGGTGATAGACCTGGTCAATGCGTTGTGAAGGCACCAg                                                                                             <  2:263025/149‑1 (MQ=255)
                                                         ttCGGTACGGCTCGAAGATGCGCGCCAGTTCACCGCTGTCACGCAGGCCACGCAGCAGGGCAGCGAACGACTGCTCGCCGATAGGCGCCCCCGGGCGCAGCAAGGCGTAGTGGTGATAGACCTGGTCAATGCGTTGTGAAGGCACCAGc                                                                                            >  3:311493/1‑149 (MQ=255)
                                                            ggTACGGCTCGAAGATGCGCGCCAGTTCACCGCTGTCACGCAGGCCACGCAGCAGGGCAGCGAACGACTGCTCGCCGATAGGCGCCCCCGGGCGCAGCAAGGCGTAGTGGTGATAGACCTGGTCAATGCGTTGTGAAGGCACCAGCTgc                                                                                         <  8:196509/149‑1 (MQ=255)
                                                             gTACGGCTCGAAGATGCGCGCCAGTTCACCGCTGTACACGCAGGCCACGCAGCAGGGCAGCGAACGACTGCTCGCCGATAGGCGCCCCCGGGCGCAGCAAGGCGTAGTGGTGATAGACCTGGTCAATGCGTTGTGAAGGCACCAGCTgc                                                                                         >  1:54718/1‑149 (MQ=255)
                                                                 ggCTCGAAGATGCGCGCCAGTTCACCGCTGTCACGCAGGCCACGCAGCAGGGCAGCGAACGACTGCTCGCCGATAGGCGCCCCCGGGCGCAGCAAGGCGTAGTGGTGATAGACCTGGTCAATGCGTTGTGAAGGCACCAGCTGCGCCCt                                                                                    >  7:328282/1‑149 (MQ=255)
                                                                  gCTCGAAGATGCGCGCCAGTTCACCGCTGTCACGCAGGCCACGCAGCAGGGCAGCGAACGACTGCTCGCCGATAGGCGCCCCCGGGCGCAGCAAGGCGTAGTGGTGATAGACCTGGTCAATGCGTTGTGAAGGCACCa                                                                                              <  1:358493/138‑1 (MQ=255)
                                                                  gCTCGAAGATGCGCGCCAGTTCACCGCTGTCACGCAGGCCACGCAGCAGGGCAGCGAACGACTGCTCGCCGATAGGCGCCCCCGGGCGCAGCAAGGCGTAGTGGTGATAGACCTGGTCAATGCGTTGTGAAGGCACCa                                                                                              >  2:358493/1‑138 (MQ=255)
                                                                   cTCGAAGATGCGCGCCAGTTCACCGCTGTCACGCAGGCCACGCAGCAGGGCAGCGAACGACTGCTCGCCGATAGGCGCCCCCGGGCGCAGCAAGGCGTAGTGGTGATAGACCTGGTCAATGCGTTGTGAAGGCACCAGCTGCGCCCTGc                                                                                  >  7:156665/1‑149 (MQ=255)
                                                                    tCGAAGATGCGCGCCAGTTCACCGCTGTCACGCAGGCCACGCAGCAGGGCAGCGAACGACTGCTCGCCGATAGGCGCCCCCGGgcg                                                                                                                                                 <  6:315155/86‑1 (MQ=255)
                                                                    tCGAAGATGCGCGCCAGTTCACCGCTGTCACGCAGGCCACGCAGCAGGGCAGCGAACGACTGCTCGCCGATAGGCGCCCCCGGGCGCAGCAAGGCGTAGTGGTGATAGACCTGGTCAATGCGTTGTGAAGGCACCAGCTGCGCCCTGCt                                                                                 >  5:130866/1‑149 (MQ=255)
                                                                      gAAGATGCGCGCCAGTTCACCGCTGTCACGCAGGCCACGCAGCAGGGCAGCGAACGACTGCTCGCCGATAGGCGCCCCCGGGCGCAGCAAGGCGTAGTGGTGATAGACCTGGTCAATGCGTTGTGAAGGCACCAGCTGCGCCCTGCTCg                                                                               >  1:352165/1‑149 (MQ=255)
                                                                        aGATGCGCGCCAGTTCACCGCTGTCACGCAGGCCACGCAGCAGGGCAGCGAACGACTGCTCGCCGATAGGCGCCCCCGGGCGCAGCAAGGCGTAGTGGTGATAGACCTGGTCAATGCGTTGTGAAGGCACCAGCTGCGCCCTGCTCGcc                                                                             >  6:43192/1‑149 (MQ=255)
                                                                          aTGCGCGCCAGTTCACCGCTGTCACGCAGGCCACGCAGCAGGGCAGCGAACGACTGCTCGCCGATAGGCGCCCCCGGGCGCAGCAAGGCGTAGTGGTGATAGACCTGGTCAATGCGTTGTGAAGGCACCAGCTGCGCCCTGCTCGCCgg                                                                           <  6:21112/149‑1 (MQ=255)
                                                                                     ttCACCGCTGTCACGCAGGCCACGCAGCAGGGCAGCGAACGACTGCTCGCCGATAGGCGCCCCCGGGCGCAGCAAGGCGTAGTGGTGATAGACCTGGTCAATGCGTTGTGAAGGCACCAGCTGCGCCCTGCTCGCCGGGTTGCGTGCCa                                                                <  2:177065/149‑1 (MQ=255)
                                                                                     ttCACCGCGGTCACGCAGGCCACGCAGCAGGGCAGCGATCGACTGCTCGCCGATAGGCGCCCCCGGGCGCAGCAAGGCGTAGTGGTGATAGACCTGGTCAATGCGTTCTGAAGGCACCAGCTGCGCCCTGCTCGCCGGGTTGCGTGCCa                                                                <  7:273666/149‑1 (MQ=255)
                                                                                        aCCGCTGTCACGCAGGCCACGCAGCAGGGCAGCGAACGACTGCTCGCCGATAGGCGCCCCCGGGCGCAGCAAGGCGTAGTGGTGATAGACCTGGTCAATGCGTTGTGAAGGCACCAGCTGCGCCCTGCTCGCCGGGTTGCGTGCCAgg                                                              >  7:50366/1‑148 (MQ=255)
                                                                                                  aCGCAGGCCACGCAGCAGGGCAGCGAACGACTGCTCGCCGATAGGCGCCCCCGGGCGCAGCAAGGCGTAGTGGTGATAGACCTGGTCAATGCGTTGTGAAGGCACCAGCTGCGCCCTGCTCGCCGGGTTGCGTGCCAGGAAATCGCt                                                      >  3:254949/1‑147 (MQ=255)
                                                                                                  aCGCAGGCCACGCAGCAGGGCAGCGAACGACTGCTCGCCGATAGGCGCCCCCGGGCGCAGCAAGGCGTAGTGGTGATAGACCTGGTCAATGCGTTGTGAAGGCACCAGCTGCGCCCTGCTCGCCGGGTTGCGTGCCAGGAAATCGCTGa                                                    <  7:136156/149‑1 (MQ=255)
                                                                                                  aCGCAGGCCACGCAGCAGGGCAGCGAACGACTGCTCGCCGATAGGCGCCCCCGGGCGCAGCAAGGCGTAGTGGTGATAGACCTGGTCAATGCGTTGTGAAGGCACCAGCTGCGCCCTGCTCGCCGGGTTGCGTGCCAGGAAATCGCTGa                                                    <  4:311493/149‑1 (MQ=255)
                                                                                                   cGCAGGCCACGCAGCAGGGCAGCGAACGACTGCTCGCCGATAGGCGCCCCCGGGCGCAGCAAGGCGTAGTGGTGATAGACCTGGTCAATGCGTTGTGAAGGCACCAGCTGCGCCCTGCTCGCCGGGTTGCGTGCCAGGAAATCGCTGAg                                                   <  2:54718/149‑1 (MQ=255)
                                                                                                    gCAGGCCACGCAGCAGGGCAGCGAACGACTGCTCGCCGATAGGCGCCCCCGGGCGCAGCAAGGCGTAGTGGTGATAGACCTGGTCAATGCGTTGTGAAGGCACCAGCTGCGCCCTGCTCGCCGGGTTGCGTGCCAGGAAATCGGTGAgg                                                  >  8:78658/1‑149 (MQ=255)
                                                                                                       ggCCACGCAGCAGGGCAGCGAACGACTGCTCGCCGATAGGCGCCCCCGGGCGCAGCAAGGCGTAGTGGTGATAGACCTGGTCAATGCGTTGTGAAGGCACCAGCTGCGCCCTGCTCGCCGGGTTGCGTGCCAGGAAATCGCTGAGGTAg                                               <  5:43192/149‑1 (MQ=255)
                                                                                                        gCCACGCAGCAGGGCAGCGAACGACTGCTCGCCGATAGGCGCCCCCGGGCGCAGCAAGGCGTAGTGGTGATAGACCTGGTCAATGCGTTGTGAAGGCACCAGCTGCGCCCTGCTCGCCGGGTTGCGTGCCAGGAAATCGCTGAGGTAgg                                              <  5:86813/149‑1 (MQ=255)
                                                                                                             gcagcaGGGCAGCGAACGACTGCTCGCCGATAGGCGCCCCCGGGCGCAGCAAGGCGTAGTGGTGATAGACCTGGTCAATGCGTTGTGAAGGCACCAGCTGCGCCCTGCTCGCCGGGTTGCGTGCCAGGAAATCGCTGAGGTAGGAAcgc                                         >  4:221173/1‑149 (MQ=255)
                                                                                                             gcagcaGGGCAGCGAACGACTGCTCGCCGATAGGCGCCCCCGGGCGCAGCAAGGCGTAGTGGTGATAGACCTGGTCAATGCGTTGTGAAGGCACCAGCTGCGCCCTGCTCGCCGGGTTGCGTGCCAGGAAATCGCTGAGGTAGGAAcgc                                         <  3:313505/149‑1 (MQ=255)
                                                                                                             gcagcaGGGCAGCGAACGACTGCTCGCCGATAGGCGCCCCCGGGCGCAGCAAGGCGTAGTGGTGATAGACCTGGTCAATGCGTTGTGAAGGCACCAGCTGCGCCCTGCTCGACGGGTTGCGTGCCAGGAAATCGCTGAGGTAGGAAcgc                                         >  5:343224/1‑149 (MQ=255)
                                                                                                               agcagGGCAGCGAACGACTGCTCGCCGATAGGCGCCCCCGGGCGCAGCAAGGCGTAGTGGTGATAGACCTGGTCAATGCGTTGTGAAGGCACCAGCTGCGCCCTGCTCGCCGGGTTGCGTGCCAGGAAATCGCTGAGGTAGGAACGCGt                                       <  3:204955/149‑1 (MQ=255)
                                                                                                               agcagGGCAGCGAACGACTGCTCGCCGATAGGCGCCCCCGGGCGCAGCAAGGCGTAGTGGTGATAGACCTGGTCAATGCGTTGTGAAGGCACCAGCTGCGCCCTGCTCGCCGGGTTGCGTGCCAGGAAATCGCTGAGGTAGGAACGCGt                                       <  8:156665/149‑1 (MQ=255)
                                                                                                               agcagGGCAGCGAACGACTGCTCGCCGATAGGCGCCCCCGGGCGCAGCAAGGCGTAGTGGTGATAGACCTGGTCAATGCGTTGTGAAGGCACCAGCTGCGCCCTGCTCGCCGGGTTGCGTGCCAGGAAATCGCTGAGGTAGGAACGCGt                                       <  2:352165/149‑1 (MQ=255)
                                                                                                                    ggCAGCGAACGACTGCTCGCCGATAGGCGCCCCCGGGCGCAGCAAGGCGTAGTGGTGATAGACCTGGTCAATGCGTTGTGAAGGCACCAGCTGCGCCCTGCTCGCCGGGTTGCGTGCCAGGAAATCGCTGAGGTGGGAACGCGTGAc                                    >  6:84904/1‑147 (MQ=255)
                                                                                                                      cAGCGAACGACTGCTCGCCGATAGGCGCCCCCGGGCGCAGCAAGGCGTAGTGGTGATAGACCTGGTCAATGCGTTGTGAAGGCACCAGCTGCGCCCTGCTCGCCGGGTTGCGTGCCAGGAAATCGCTGAGGTAGGAACGCGTGACCAgg                                <  2:208150/149‑1 (MQ=255)
                                                                                                                       aGCGAACGACTGCTCGCCGATAGGCGCCCCCGGGCGCAGCAAGGCGTAGTGGTGATAGACCTGGTCAATGCGTTGTGAAGGCACCAGCTGCGCCCTGCTCGCCGGGTTGCGTGCCAGGAAATCGCTGAGGTAGGAACGCGTGACCAggg                               <  3:266209/149‑1 (MQ=255)
                                                                                                                       aGCGAACGACTGCTCGCCGATAGGCGCCCCCGGGCGCAGCAAGGCGTAGTGGTGATAGACCTGGTCAATGCGTTGTGAAGGCACCAGCTGCGCCCTGCTCGCCGGGTTGCGTGCCAGGAAATCGCTGAGGTAGGAACGCGTGACCAggg                               <  8:102974/149‑1 (MQ=255)
                                                                                                                         cgaccgacTGCTCGCCGATAGGCGCCCCCGGGCGCAGCAAGGCGTAGTGGGGATAGACCTGGTCAATGCGTTGTGAAGGCACCAGCTGCGCCCTGCTCGCCGGGTTGCGTGCCAGGAAATCGCTGAGGTAGGAACGCGTGACCAGGGc                              >  4:108202/1‑148 (MQ=255)
                                                                                                                         cGAACGACTGCTCGCCGATAGGCGCCCCCGGGCGCAGCAAGGCGTAGTGGTGATAGACCTGGTCAATGCGTTGTGAAGGCACCAGCTGCGCCCTGCTCGCCGGGTTGCGTGCCAGGAAATCGCTGAGGTAGGAACGCGTGACCAGGGc                              >  1:215067/1‑148 (MQ=255)
                                                                                                                         cGAACGACTGCTCGCCGATAGGCGCCCCCGGGCGCAGCAAGGCGTAGTGGTGATAGACCTGGTCAATGCGTTGTGAAGGCACCAGCTGCGCCCTGCTCGCCGGGTTGCGTGCCAGGAAATCGCTGAGGTAGGAACGCGTGACCAGGGc                              >  1:220594/1‑148 (MQ=255)
                                                                                                                         cGAACGACTGCTCGCCGATAGGCGCCCCCGGGCGCAGCAAGGCGTAGTGGTGATAGACCTGGTCAATGCGTTGTGAAGGCACCAGCTGCGCCCTGCTCGCCGGGTTGCGTGCCAGGAAATCGCTGAGGTAGGAACGCGTGACCAGGGc                              >  4:346057/1‑148 (MQ=255)
                                                                                                                          gAACGACTGCTCGCCGATAGGCGCCCCCGGGCGCAGCAAGGCGTAGTGGTGATAGACCTGGTCAATGCGTTGTGAAGGCACCAGCTGCGCCCTGCTCGCCGGGTTGCGTGCCAGGAAATCGCTGAGGTAGGAACGCGTGACCAGGGCaa                            >  1:323686/1‑149 (MQ=255)
                                                                                                                            aCGACTGCTCGCCGATAGGCGCCCCCGGGCGCAGCAAGGCGTAGTGGTGATAGACCTGGTCAATGCGTTGTGAAGGCACCAGCTGCGCCCTGCTCGCCGGGTTGCGTGCCAGGAAATCGCTGAGGTAGGAACGCGTGACCAGGGCAata                          >  8:264339/1‑149 (MQ=255)
                                                                                                                               aCTGCTCGCCGATAGGCGCCCCCGGGCGCAGCAAGGCGTAGTGGTGATAGACCTGGTCAATGCGTTGTGAAGGCACCAGCTGCGCCCTGCTCGCCGGGTTGCGTGCCAGGAAATCGCTGAGGTAGGAACGCGTGACCAGGGCAATATCg                       <  2:323686/149‑1 (MQ=255)
                                                                                                                                   cTCGCCGATAGGCGCCCCCGGGCGCAGCAAGGCGTAGTGGTGATAGACCTGGTCAATGCGTTGTGAAGGCACCAGCTGCGCCCTGCTCGCCGGGTTGCGTGCCAGGAAATCGCTGAGGTAGGAACGCGTGACCAGGGCAATATCGgcgc                   >  3:8145/1‑149 (MQ=255)
                                                                                                                                   cTCGCCGATAGGCGCCCCCGGGCGCAGCAAGGCGTAGTGGTGATAGACCTGGTCAATGCGTTGTGAAGGCACCAGCTGCGCCCTGCTCGCCGGGTTGCGTGCCAGGAAATCGCTGAGGTAGGAACGCGTGACCAGGGCAATATCGgc                     <  8:50366/147‑1 (MQ=255)
                                                                                                                                    tCGCCGATAGGCGCCCCCGGGCGCAGCAAGGCGGTAGTGGTGATAGACCTGGTCAATGCGTTGTGAAGGCACCAGCTGCGCCCTGCTCGCCGGGTTGCGTGCCAGGAAATCGCTGAGGTAGGAACGCGTGACCAGGGCAATATCGgcgc                   >  1:44481/1‑149 (MQ=255)
                                                                                                                                    tCGCCGATAGGCGCCCCCGGGCGCAGCAAGGCGTAGTGGTGATAGACCTGGTCAATGCGTTGTGAAGGCACCAGCTGCGCCCTGCTCGCCGGGTTGCGTGCCAGGAAATCGCTGAGGTAGGAACGCGTGACCAGGGCAATATCGGcgcg                  >  3:210835/1‑149 (MQ=255)
                                                                                                                                    tCGCCGATAGGCGCCCCCGGGCGCAGCAAGGCGTAGTGGTGATAGACCTGGTCAATGCGTTGTGAAGGCACCAGCTGCGCCCTGCTCGCCGGGTTGCGTGCCAGGAAATCGCTGAGGTAGGAACGCGTGACCAGGGCAATATCGGcgcg                  >  7:65385/1‑149 (MQ=255)
                                                                                                                                    tCGCCGATAGGCGCCCCCGGGCGCAGCAAGGCGTAGTGGTGATAGACCTGGTCAATGCGTTGTGAAGGCACCAGCTGCGCCCTGCTCGCCGGGTTGCGTGCCAGGAAATCGCTGAGGTAGGAACGCGTGACCAGGGCAATATCGGcgcg                  >  5:276509/1‑149 (MQ=255)
                                                                                                                                        cGATAGGCGCCCCCGGGCGCAGCAAGGCGTAGTGGTGATAGACCTGGTCAATGCGTTGTGAAGGCACCAGCTGCGCCCTGCTCGCCGGGTTGCGTGCCAGGAAATCGCTGAGGTAGGAACGCGTGACCAGGGCAATAACGGCGcggccg              >  3:254019/1‑149 (MQ=255)
                                                                                                                                          tgctgcgcCCGGGGGCGCAGCAAGGCGTAGTGGTGATAGACCTGGTCAATGCGTTGTGAAGGCACCAGCTGCGCCCTGCTCGCCGGGTTGCGTGCCAGGAAATCGCTGAGGTAGGAACGCGTGACCAGGGCAATATCGGCGcggccggc            >  1:118913/5‑149 (MQ=255)
                                                                                                                                          tgctgcgcCCGGGGGCGCAGCAAGGCGTAGTGGTGATAGACCTGGTCAATGCGTTGTGAAGGCACCAGCTGCGCCCTGCTCGCCGGGTTGCGTGCCAGGAAATCGCTGAGGTAGGAACGCGTGACCAGGGCAATATCGGCGcggccggc            >  1:114329/5‑149 (MQ=255)
                                                                                                                                           tAGGCGCCCCCGGGCGCAGCAAGGCGTAGTGGTGATAGACCTGGTCAATGCGTTGTGAAGGCACCAGCTGCGCCCTGCTCGCCGGGTTGCGTGCCAGGAAATCGCTGAGGTAGGAACGCGTGACCa                                  <  1:65149/126‑1 (MQ=255)
                                                                                                                                           tAGGCGCCCCCGGGCGCAGCAAGGCGTAGTGGTGATAGACCTGGTCAATGCGTTGTGAAGGCACCAGCTGCGCCCTGCTCGCCGGGTTGCGTGCCAGGAAATCGCTGAGGTAGGAACGCGTGACCa                                  >  2:65149/1‑126 (MQ=255)
                                                                                                                                            aGGCGCCCCCGGGCGCAGCAAGGCGTAGTGGTGATAGACCTGGTCAATGCGTTGTGAAGGCACCAGCTGCGCCATGCTCGCCGGGTTGCGTGCCAGGAAATCGGGGAGGTAGGACCGCGTGACCAGGGCAATATCGGCGCGGCCGGCCt          >  3:147025/1‑149 (MQ=255)
                                                                                                                                              gcgcCCCCGGGCGCAGCAAGGCGTAGTGGTGATAGACCTGGTCAATGCGTTGTGAAGGCACCAGCTGCGCCCTGCTCGCCGGGTTGCGTGCCAGGAAATCGCTGAGGTAGGAACGCGTGACCAGGGCAATATCGGCGcggccg              <  3:182463/143‑1 (MQ=255)
                                                                                                                                              gcgcCCCCGGGCGCAGCAAGGCGTAGTGGTGATAGACCTGGTCAATGCGTTGTGAAGGCACCAGCTGCGCCCTGCTCGCCGGGTTGCGTGCCAGGAAATCGCTGAGGTAGGAACGCGTGACCAGGGCAATATCGGCGcggccg              >  4:182463/1‑143 (MQ=255)
                                                                                                                                               cgcCCGGGGGCGCAGCAAGGCGTAGTGGTGATAGACCTGGTCAATGCGTTGTGAAGGCACCAGCTGCGCCCTGCTCGCCGGGTTGCGTGCCAGGAAATCGCTGAGGTAGGAACGCGTGACCAGGGCAATATCGGCGCGGCCGGCCTg         <  2:114329/147‑1 (MQ=255)
                                                                                                                                               cgcCCGGGGGCGCAGCAAGGCGTAGTGGTGATAGACCTGGTCAATGCGTTGTGAAGGCACCAGCTGCGCCCTGCTCGCCGGGTTGCGTGCCAGGAAATCGCTGAGGTAGGAACGCGTGACCAGGGCAATATCGGCGCGGCCGGCCTg         <  2:118913/147‑1 (MQ=255)
                                                                                                                                               cgcCCCCGGGCGCAGCAAGGCGTAGTGGTGATAGACCTGGTCAATGCGTTGTGAAGGCACCAGCTGCGCCCTGCTCGCCGGGTTGCGTGCCAGGAAATCGCTGAGGTAGGAACGCGTGACCAGGGCAATATCGGCGCGGCCGGCCTGCa       <  7:16273/149‑1 (MQ=255)
                                                                                                                                                gcCCCCGGGCGCAGCAAGGCGTAGTGGTGATAGACCTGGTCAATGCGTTGTGAAGGCACCAGCTGCGCCCTGCTCGCCGGGTTGCGTGCCAGGAAATCGCTGAGGTAGGAACGCGTGACCAGGGCAATATCGGCGCGGCCGGCCTGCAc      <  4:210835/149‑1 (MQ=255)
                                                                                                                                                   cccGGGCGCAGCAAGGCGTAGTGGTGATAGACCTGGTCAATGCGTTGTGAAGGCACCAGCTGCGCCCTGCTCGCCGGGTTGCGTGCCAGGAAATCGCTGAGGTAGGACCGCGTGACCAGGGCAATATCGGCGCGGCCGGCCTGCACCAt   >  8:297040/1‑149 (MQ=255)
                                                                                                                                                   cccGGGCGCAGCAAGGCGTAGTGGTGATAGACCTGGTCAATGCGTTGTGAAGGCACCAGCTGCGCCCTGCTCGCCGGGTTGCGTGCCAGGAAATCGCTGAGGTAGGAACGCGTGACCAGGGCAATATCGGCGCGGCCGGCCTGCACCAt   >  6:88049/1‑149 (MQ=255)
                                                                                                                                                    ccGGGCGCAGCAAGGCGTAGTGGTGATAGACCTGGTCAATGCGTTGTGAAGGCACCAGCTGCGCCCTGCTCGCCGGGTTGCGTGCCAGGAAATCGCTGAGGTAGGAACGCGTGACCAGGGCAATATCGGCGCGGCCGGCCTGCACCAt   >  8:314347/1‑148 (MQ=255)
                                                                                                                                                    ccGGGCGCAGCAAGGCGTAGTGGTGATAGACCTGGTCAATGCGTTGTGAAGGCACCAGCTGCGCCCTGCTCGCCGGGTTGCGTGCCAGGAAATCGCTGAGGTAGGAACGCGTGACCAGGGCAATATCGGCGCGGCCGGCCTGCACCATc  <  6:130866/149‑1 (MQ=255)
                                                                                                                                                    ccGGGCGCAGCAAGGCGTAGTGGTGATAGACCTGGTCAATGCGTTGTGAAGGCACCAGCTGCGCCCTGCTCGCCGGGTTGCGTGCCAGGAAATCGCTGAGGTAGGAACGCGTGACCAGGGCAATATCGGCGCGGCCAGCCTGCACCATc  <  8:328282/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
ACTGCAGGGAACGTTGCCGCAACACGGGAATTTAGTCTTGCGGTGCAGCCACTTTGATTCGGTACGGCTCGAAGATGCGCGCCAGTTCACCGCTGTCACGCAGGCCACGCAGCAGGGCAGCGAACGACTGCTCGCCGATAGGCGCCCCCGGGCGCAGCAAGGCGTAGTGGTGATAGACCTGGTCAATGCGTTGTGAAGGCACCAGCTGCGCCCTGCTCGCCGGGTTGCGTGCCAGGAAATCGCTGAGGTAGGAACGCGTGACCAGGGCAATATCGGCGCGGCCGGCCTGCACCATC  >  NC_002947/4833223‑4833518

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: