Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,261,067 A→G 5.9% D178D (GAT→GAC) ‡ PP_4637 ← 5‑methyltetrahydropteroyltriglutamate‑ homocysteine S‑methyltransferase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,261,0670AG 5.9% 247.0 / 4.1 84D178D (GAT→GAC) ‡PP_46375‑methyltetrahydropteroyltriglutamate‑ homocysteine S‑methyltransferase family protein
Reads supporting (aligned to +/- strand):  ref base A (40/39);  new base G (2/3);  total (42/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

GTGCAGCAGCCCCAGGGTCTTGCCCAAATCACAAGGCGAGGCGTTGCGACCGTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTCAACGGTGCGCGTTGAAGGATGTTGTAGACGCGTTCATAGGCGTTTTTCCAGTCCTGGGGCAGCGCCTGGGCGAGGA  >  NC_002947/5260920‑5261215
                                                                                                                                                    |                                                                                                                                                     
gTGCAGCAGCCCCAGGGTCTTGCCCAAATCACAAGGCGAGGCGTTGCGACCGTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCAt                                                                                                                                                      >  7:6653/1‑149 (MQ=255)
 tGCAGCAGCCCCAGGGTCTTGCCCAAATCACAAGGCGAGGCGTTGTGACCGTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATc                                                                                                                                                     <  3:200344/149‑1 (MQ=255)
    cgaagcCCCAGGGTCTTGCCCAAATCACAAGGCGAGGCGTTGCGACCGTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAAc                                                                                                                                                  <  4:102607/146‑1 (MQ=255)
    agcagcCCCAGGGTCTTGCCCAAATCACAAGGCGAGGCGTTGCGACCGTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAAc                                                                                                                                                  >  3:312195/1‑149 (MQ=255)
     gcagcCCCAGGGTCTTGCCCAAATCACAAGGCGAGGCGTTGCGACCGTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCGCCAAcc                                                                                                                                                 >  2:433324/1‑149 (MQ=255)
       agcCCCAGGGTCTTGCCCAAATCACAAGGCGAGGCGTTGCGACCGTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGATCCAGATCGACGTGCAGGCCATCAACCgg                                                                                                                                               >  4:320770/1‑149 (MQ=255)
        gcCCCAGGGTCTTGCCCAAATCACAAGGCGAGGCGTTGCGACCGTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGc                                                                                                                                              >  5:428171/1‑149 (MQ=255)
            cAGGGTCTTGCCCAAATCACAAGGCGAGGCGTTGCGACCGTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCGCCAACCTGCCGGt                                                                                                                                          >  5:333085/1‑149 (MQ=255)
               ggTCTTGCCCAAATCACAAGGCGAGGCGTTGCGACCGTCGACCAGGCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTggc                                                                                                                                      <  4:417188/149‑1 (MQ=255)
                 tCTTGCCCAAATCACAAGGCGAGGCGTTGCGACCGTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTggcg                                                                                                                                     <  2:371785/149‑1 (MQ=255)
                   ttGCCCAAATCACAAGGCGAGGCGTTGCGACCGTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCGCCAACCTGCCGGTTgatggc                                                                                                                                   <  6:333085/149‑6 (MQ=255)
                     gCCCAAATCACAAGGCGAGGCGTTGCGACCGTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCGCCAACCTGCCGGTTgatggc                                                                                                                                   <  1:433324/147‑6 (MQ=255)
                     gCCCAAATCACAAGGCGAGGCGTTGCGACCGTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTggcgg                                                                                                                                    <  4:312195/146‑1 (MQ=255)
                      cccAAATCACCAGGCGAGGCGTTGCGACCGTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCGCCAACCTGCCGGTTGATGGcct                                                                                                                                 <  8:83056/148‑2 (MQ=255)
                      cccAAATCACAAGGCGAGGCGTTGCGACCGTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCa                                                                                                                                 <  1:415744/148‑1 (MQ=255)
                       ccAAATCACAAGGCGAGGCGTTGCGACCGTCGACCAGTCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAgg                                                                                                                               >  2:313632/1‑149 (MQ=255)
                       ccAAATCACAAGGCGAGGCGTTGCGACCGTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAgg                                                                                                                               >  7:120314/1‑149 (MQ=255)
                       ccAAATCACAAGGCGAGGCGTTGCGACCGTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAgg                                                                                                                               >  2:292422/1‑149 (MQ=255)
                              acaAGGCGAGGCGTTGCGACCGTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCAGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAgg                                                                                                                         >  1:347621/1‑149 (MQ=255)
                              acaAGGCGAGGCGTTGCGACCGTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGt                                                                                                                        >  2:355239/1‑149 (MQ=255)
                              acaAGGCGAGGCGTTGCGACCGTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGt                                                                                                                        >  3:25366/1‑149 (MQ=255)
                                aaGGCGAGGCGTTGCGACCGTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTg                                                                                                                      <  2:176579/149‑1 (MQ=255)
                                 aggcgaggcgTTGCGACCGTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGc                                                                                                                    <  7:19370/149‑1 (MQ=255)
                                 aggcgaggcgTTGCGACCGTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGc                                                                                                                    <  8:39404/149‑1 (MQ=255)
                                      aggcgTTGCGACCGTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTcc                                                                                                               <  2:376772/149‑1 (MQ=255)
                                      aggcgTTGCGACCGTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTc                                                                                                                >  5:381793/1‑148 (MQ=255)
                                       ggcgTTGCGACCGTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCCATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCa                                                                                                              <  1:355239/149‑1 (MQ=255)
                                        gcgTTGCGACCGTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAg                                                                                                             <  8:19894/149‑1 (MQ=255)
                                           ttGCGACCGTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGc                                                                                                           >  1:185674/1‑148 (MQ=255)
                                               gACCGTCGACCAGTCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCAcc                                                                                                       >  1:430561/1‑148 (MQ=255)
                                                  cGTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGtt                                                                                                                       <  4:98398/130‑1 (MQ=255)
                                                  cGTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGtt                                                                                                                       >  3:98398/1‑130 (MQ=255)
                                                      gACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTACCAGGCCACCGAAGTACg                                                                                               >  4:100905/1‑149 (MQ=255)
                                                       aCCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGt                                                                                              >  2:101840/1‑149 (MQ=255)
                                                        ccAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTg                                                                                             >  8:93812/1‑149 (MQ=255)
                                                         cAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTgg                                                                                            >  6:216416/1‑149 (MQ=255)
                                                            gCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGa                                                                                         <  6:381793/149‑1 (MQ=255)
                                                             cccAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTgg                                                                                            >  6:246934/1‑145 (MQ=255)
                                                             cccAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTgg                                                                                            <  5:246934/145‑1 (MQ=255)
                                                              ccAGCGACAGTACCTTGTAGGCCGGCAGGCGATTGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGAt                                                                                        >  1:415906/1‑148 (MQ=255)
                                                              ccAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGAt                                                                                        >  3:424885/1‑148 (MQ=255)
                                                               cAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCCGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCa                                                                                      <  8:105409/149‑1 (MQ=255)
                                                                 gCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGc                                                                                    >  6:122302/1‑149 (MQ=255)
                                                                 gCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGc                                                                                    <  6:105546/149‑1 (MQ=255)
                                                                  cGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGct                                                                                   >  1:104713/1‑149 (MQ=255)
                                                                      aGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGGCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGcttctt                                                                               <  1:313632/149‑1 (MQ=255)
                                                                      aGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGcttctt                                                                               <  8:5334/149‑1 (MQ=255)
                                                                      aGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGcttctt                                                                               <  2:347621/149‑1 (MQ=255)
                                                                      aGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGcttct                                                                                >  2:240514/1‑148 (MQ=255)
                                                                            ttGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCATGAACCAGATCGACGTGCAGGCCATCAACCGGCAGCTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACCTGGCGATCAGCTTCTTCAACgg                                                                         >  8:333352/1‑149 (MQ=255)
                                                                             tGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTCAACGGt                                                                        >  3:76450/1‑149 (MQ=255)
                                                                             tGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTCAACGGt                                                                        >  3:424846/1‑149 (MQ=255)
                                                                              gTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTCAACGGt                                                                        <  3:258787/148‑1 (MQ=255)
                                                                              gTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTCAACGGt                                                                        <  4:76450/148‑1 (MQ=255)
                                                                                aGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGAGCAGCTTCTTCAAGGGTgcg                                                                     >  3:319005/1‑149 (MQ=255)
                                                                                 ggCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTCAACGGTgcg                                                                     <  1:101840/149‑1 (MQ=255)
                                                                                    cGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTCAACGGTGCGCGtt                                                                 >  7:256201/1‑149 (MQ=255)
                                                                                     ggCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTCAACGGTGCGCGTTg                                                                <  1:292422/149‑1 (MQ=255)
                                                                                        aGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTCAACGGTGCGCGTTGAAg                                                             <  7:93812/149‑1 (MQ=255)
                                                                                           cGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTCAACGGTGCGCGTTGAAGGAt                                                          <  5:247284/149‑1 (MQ=255)
                                                                                           cGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTCAACGGTGCGCGTTGAAGGAt                                                          <  2:49487/149‑1 (MQ=255)
                                                                                            gATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAAGCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTCAACGGTGCGCGTTGAAGGAtg                                                         >  1:290204/1‑149 (MQ=255)
                                                                                            gATAGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTCAACGGTGCGCGTTGAAGGAtg                                                         >  5:210364/1‑149 (MQ=255)
                                                                                                  agaATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTCAACGGTGCGCGTTGAAGGATGTTGTa                                                    <  5:122302/149‑1 (MQ=255)
                                                                                                        gTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGGGATCAGCTTCTTCAACGGTGCGCGTTGAAGGATGTTGTAGAcgc                                               <  4:95262/148‑1 (MQ=255)
                                                                                                                cTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTCAACGGTGCGCGTTGAAGGATGTTGTAGACGCGTTCATAgg                                      >  1:255803/1‑149 (MQ=255)
                                                                                                                   cTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGTGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTCAACGGTGCGCGTTGAAGGATGTTGTAGAGGCGTTCATAGGCGt                                   <  2:90185/149‑1 (MQ=255)
                                                                                                                   cTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTCAACGGTGCGCGTTGAAGGATGTTGTAGACGCGTTCATAGGc                                     <  2:281801/147‑1 (MQ=255)
                                                                                                                   cTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTCAACGGTGCGCGTTGAAGGATGTTGTAGACGCGTTCATAGGc                                     <  2:430561/147‑1 (MQ=255)
                                                                                                                   cTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTCAACGGTGCGCGTTGAAGGATGTTGTAGACGCGTTCATAGGCGt                                   >  4:148073/1‑149 (MQ=255)
                                                                                                                      gggCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTCAACGGTGCGCGTTGAAGGATGTTGTAGACGCGTTCATAGGCGtttt                                >  6:347402/1‑149 (MQ=255)
                                                                                                                        gcgcCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTCAACGGTGCGCGTTGAAGGATGTTGTAGACGCGTTCATAGGCGTTTTTc                              <  2:347999/149‑1 (MQ=255)
                                                                                                                          gcatGATCCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCAGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTCAACGGTGCGCGTTGAAGGATGTAGTAGACGCGTACATAGGCGTTTTTCCa                            >  6:152922/4‑149 (MQ=255)
                                                                                                                            cTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTCAACGGTGCGCGTTGAAGGATGTTGTAGACGCGTTCATAGGCGTTTTTCCAGt                          >  1:199408/1‑149 (MQ=255)
                                                                                                                            cTGAACCAGATCGACGTGCAGGCCATCAAACGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTCAACGGTGCGCGTTGAAGGATGTTGTAGACGCATTCATAGGCGTTTTTCCAGt                          >  1:32232/1‑149 (MQ=255)
                                                                                                                               aaCCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTCAACGGTGCGCGTTGAAGGATGTTGTAGACGCGTTCATAGGCGTTTTTCCAGTcc                        <  4:424885/148‑1 (MQ=255)
                                                                                                                                 ccAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTCAACGGTGCGCGTTGAAGGATGTTGTAGACGCGTTCATAGGCGTTTTTCCAGTCCTgg                     >  7:146843/1‑149 (MQ=255)
                                                                                                                                 ccAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTCAACGGTGCGCGTTGAAGGATGTTGTAGACGCGTTCATAGGCGTTTTTCCAGTCCTgg                     <  8:120314/149‑1 (MQ=255)
                                                                                                                                 ccAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTCAACGGTGCGCGTTGAAGGATGTTGTAGACGCATTCATAGGCGTTTTTCCAGTCCTgg                     <  2:32232/149‑1 (MQ=255)
                                                                                                                                       cGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTCAACGGTgcgc                                                                    <  6:228453/96‑1 (MQ=255)
                                                                                                                                       cGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTCAACGGTgcgc                                                                    >  5:228453/1‑96 (MQ=255)
                                                                                                                                              cAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTCAACGGTGCGCGTTGAAGGATGTTGTAGACGCGTTCATAGGCGTTTTTCCAGTCCTGGGGCAGCGCCTggg        >  5:442138/1‑149 (MQ=255)
                                                                                                                                                ggCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTCAACGGTGCGCGTTGAAGGATGTTGTAGACGCGTTCATAGGCGTTTTTCCAGTCCTGGGGCAGCGCCTGGGCg      <  2:270353/149‑1 (MQ=255)
                                                                                                                                                  ccATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTCAACGGTGCGCGTTGAAGGATGTTGTAGACGCGTTCATAGGCGTTTTTCCAGTCCTGGGGCAGCGCCTGGGCGAg    <  2:415906/149‑1 (MQ=255)
                                                                                                                                                    aTCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTCCTTCAACGGTGCGCGTTGAAGGATGTTGTAGACGCGTTCATAGGCGTTTTTCCAGTCCTGGGGCAGCGCCTGGGCGAGGa  >  6:96669/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GTGCAGCAGCCCCAGGGTCTTGCCCAAATCACAAGGCGAGGCGTTGCGACCGTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTCAACGGTGCGCGTTGAAGGATGTTGTAGACGCGTTCATAGGCGTTTTTCCAGTCCTGGGGCAGCGCCTGGGCGAGGA  >  NC_002947/5260920‑5261215

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: