Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,882,415 A→G 5.0% Q353R (CAG→CGG) ‡ PP_4289 → membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,882,4150AG 5.0% 516.0 / 9.5 158Q353R (CAG→CGG) ‡PP_4289membrane protein
Reads supporting (aligned to +/- strand):  ref base A (70/80);  new base G (6/2);  total (76/82)
Fisher's exact test for biased strand distribution p-value = 1.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.41e-01

Only 100 of 210 total aligned reads displayed.

TGAGCAAGCCGCAGCGAAGATCGAGTACCAGCCGCTGCCGGCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGCCCCTGCCGGCGTGATGTTCGACACTCCACAGCAGAT  >  NC_002947/4882266‑4882562
                                                                                                                                                     |                                                                                                                                                   
  aGCAAGCCGCAGCGAAGATCGAGTACCAGCCGCTGCCGGCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGgcccaggccgctgcccagg                                                                                                                                                   <  1:4516709/150‑1 (MQ=255)
   gcctcacGCAGCGAAGATCGAGTACCAGCCGCTGCCGGCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGgcccaggccgctgcccaggc                                                                                                                                                  >  2:835546/7‑150 (MQ=255)
    cAAGCCGCAGCGAAGATCGAGTACCAGCCGCTGCCGGCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGgcccaggccgctgcccaggcc                                                                                                                                                 >  2:1691312/1‑150 (MQ=255)
      aGCCGCAGCGAAGATCGAGTACCAGCCGCTGCCGGCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGgcccaggccgctgcccaggccgc                                                                                                                                               <  1:1333404/150‑1 (MQ=255)
         cGCAGCGAAGATCGAGTACCAGCCGCTGCCGGCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCctgggcagcggcctgggccggtgc                                                                                                                                            <  1:1588684/150‑1 (MQ=255)
          gCAGCGAAGATCGAGTACCAGCCGCTGCCGGCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCcaggccgctgcccaggccgctgcc                                                                                                                                           <  2:3628293/150‑1 (MQ=255)
          gCAGCGAAGATCGAGTACCAGCCGCTGCCGGCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCcaggccgctgcccaggccgctgcc                                                                                                                                           >  2:476701/1‑150 (MQ=255)
            aGCGAAGATCGAGTACCAGCCGCTGCCGGCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAggccgctgcccaggccgctgccca                                                                                                                                         <  2:1795160/150‑1 (MQ=255)
            aGCGAAGATCGAGTACCAGCCGCTGCCGGCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAggccgctgcccaggccgctgccca                                                                                                                                         <  2:70631/150‑1 (MQ=255)
             gCGAAGATCGAGTACCAGCCGCTGCCGGCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCAGCGGCctgggcagcggcctgggccggtgcttgc                                                                                                                                        <  1:1802508/150‑5 (MQ=38)
             gCGAAGATCGAGTACCAGCCGCTGCCGGCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCcaggccgcggcctgggcagcggcctgg                                                                                                                                        <  1:4289907/150‑4 (MQ=255)
               gAAGATCGAGTACCAGCCGCTGCCGGCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCcgctgcccaggccgctgcccaggc                                                                                                                                      >  1:3708160/1‑150 (MQ=255)
                aaGATCGAGTACCAGCCGCTGCCGGCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCgctgcccaggccgctgcccaggcc                                                                                                                                     <  1:5432370/150‑1 (MQ=255)
                   aTCGAGTACCAGCCGCTGCCGGCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGa                                                                                                                                  <  1:2519594/150‑1 (MQ=255)
                    tCGAGTACCAGCCGCTGCCGGCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAg                                                                                                                                 >  1:4681353/1‑150 (MQ=255)
                    tCGAGTACCAGCCGCTGCCGGCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAg                                                                                                                                 >  1:2805618/1‑150 (MQ=255)
                       aGTACCAGCCGCTGCCGGCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTg                                                                                                                              >  1:3218709/1‑150 (MQ=255)
                           ccAGCCGCTGCCGGCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGa                                                                                                                          >  1:1522143/1‑150 (MQ=255)
                           ccAGCCGCTGCCGGCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGa                                                                                                                          >  1:1310700/1‑150 (MQ=255)
                            cAGCCGCTGCCGGCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGgcccaggccgctgccca                                                                                                                                                     >  2:3045411/1‑122 (MQ=255)
                            cAGCCGCTGCCGGCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGgcccaggccgctgccca                                                                                                                                                     <  1:3045411/122‑1 (MQ=255)
                            cAGCCGCTGCCGGCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGaa                                                                                                                         >  1:1265133/1‑150 (MQ=255)
                            cAGCCGCTGCCGGCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGaa                                                                                                                         >  1:2900960/1‑150 (MQ=255)
                             aGCCGCTGCCGGCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAAt                                                                                                                        <  2:239408/150‑1 (MQ=255)
                              gCCGCTGCCGGCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAAtt                                                                                                                       <  2:1310472/150‑1 (MQ=255)
                                   tGCCGGCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACa                                                                                                                  <  1:232949/150‑1 (MQ=255)
                                   tGCCGGCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACa                                                                                                                  >  2:4831957/1‑150 (MQ=255)
                                   tGCCGGCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACa                                                                                                                  <  2:5262056/150‑1 (MQ=255)
                                    gCCGGCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACaa                                                                                                                 <  2:2318507/150‑1 (MQ=255)
                                    gCCGGCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACaa                                                                                                                 <  1:263027/150‑1 (MQ=255)
                                       ggCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGACTGCGACAAGag                                                                                                              >  2:4390410/1‑149 (MQ=255)
                                       ggCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGAt                                                                                                              >  1:3828215/1‑150 (MQ=255)
                                          aaCTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCca                                                                                                           <  2:1421109/150‑1 (MQ=255)
                                          aaCTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCca                                                                                                           <  1:1489327/150‑1 (MQ=255)
                                           aCTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCcac                                                                                                          <  1:585176/150‑1 (MQ=255)
                                                      cggcAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGTCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCa                                                                                               <  1:792420/150‑1 (MQ=255)
                                                       ggcAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAg                                                                                              >  2:2110102/1‑150 (MQ=255)
                                                              cGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCt                                                                                       <  1:965769/150‑1 (MQ=255)
                                                              cGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCt                                                                                       <  2:879663/150‑1 (MQ=255)
                                                                 cTGCCGAGCAGCGTGGCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCa                                                                                    >  1:4462612/1‑150 (MQ=255)
                                                                     cGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACtgt                                                                                <  1:2482455/150‑1 (MQ=255)
                                                                      gAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACtgtg                                                                               >  2:3312784/1‑150 (MQ=255)
                                                                       agcagcGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACtgtgt                                                                              <  2:1123522/150‑1 (MQ=255)
                                                                           gcGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCa                                                                          <  2:2759322/150‑1 (MQ=255)
                                                                              tGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTc                                                                       <  2:3708160/150‑1 (MQ=255)
                                                                                 cGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCTAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTc                                                                    <  1:835546/150‑1 (MQ=255)
                                                                                 cGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTc                                                                    <  2:3897231/150‑1 (MQ=255)
                                                                                 cGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGGGTGCCACTCGTc                                                                    >  2:367596/1‑150 (MQ=255)
                                                                                  gAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTcc                                                                   >  2:4366945/1‑150 (MQ=255)
                                                                                        gccgccAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCg                                                             >  2:5150091/1‑150 (MQ=255)
                                                                                        gccgccAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCg                                                             >  2:4242487/1‑150 (MQ=255)
                                                                                        gccgccAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCg                                                             >  2:5376148/1‑150 (MQ=255)
                                                                                        gccgccAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCg                                                             >  2:1117188/1‑150 (MQ=255)
                                                                                        gccgccAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGATGCCCAGGACGGAGGAGGCGAATTCGACAAGATCCACAAAGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCg                                                             >  1:2792360/1‑150 (MQ=255)
                                                                                         ccgccAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGa                                                            <  1:5296674/150‑1 (MQ=255)
                                                                                         ccgccAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGa                                                            <  2:2805618/150‑1 (MQ=255)
                                                                                           gccAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCATTGTGTGCCACTCGTCCAAGCCGAcc                                                          >  1:5058930/1‑150 (MQ=255)
                                                                                             cAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAg                                                        >  1:2381062/1‑150 (MQ=255)
                                                                                                 gCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCa                                                    <  1:4762396/150‑1 (MQ=255)
                                                                                                  cTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCAc                                                   <  1:5223136/150‑1 (MQ=255)
                                                                                                   tGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACt                                                  >  2:2953229/1‑150 (MQ=255)
                                                                                                   tGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACt                                                  <  2:3800155/150‑1 (MQ=255)
                                                                                                     ccGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCTGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGACGACCAGCCCACTGt                                                >  2:1182337/1‑150 (MQ=255)
                                                                                                      cGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACtgt                                                                                >  2:1537917/1‑117 (MQ=255)
                                                                                                      cGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGtt                                               <  1:4834882/150‑1 (MQ=255)
                                                                                                      cGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCAAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGtt                                               >  1:1663886/1‑150 (MQ=255)
                                                                                                       gAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTc                                              <  2:4378325/150‑1 (MQ=255)
                                                                                                       gAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTc                                              <  2:258439/150‑1 (MQ=255)
                                                                                                       gAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTc                                              >  1:674979/1‑150 (MQ=255)
                                                                                                        aaGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCa                                             <  1:2617199/150‑1 (MQ=255)
                                                                                                           cACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCa                                          >  1:3132807/1‑150 (MQ=255)
                                                                                                               cgccgcACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGc                                      <  2:5061078/150‑1 (MQ=255)
                                                                                                                 ccgcACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCa                                                                                               >  1:2601678/1‑91 (MQ=255)
                                                                                                                 ccgcACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCa                                                                                               <  2:2601678/91‑1 (MQ=255)
                                                                                                                  cgcACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGcccc                                   >  1:2334625/1‑150 (MQ=255)
                                                                                                                  cgcACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGcccc                                   <  2:1864530/150‑1 (MQ=255)
                                                                                                                   gcACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGCCCCt                                  >  2:1329979/1‑150 (MQ=255)
                                                                                                                        tGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGCCCCTGCCgg                             >  2:3910162/1‑150 (MQ=255)
                                                                                                                           gccacctccGGCCTGGGCAGCGGCCTGGGCAGCGGCCTGGGCCGGTGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGCCCCTGCCGGCGt                          >  2:2206550/8‑150 (MQ=25)
                                                                                                                           gcAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGCCCCTGCCGGCGt                          <  1:1079255/150‑1 (MQ=255)
                                                                                                                             cacctccGGCCTGGGCAGCGGCCTGGGCAGCGGCCTGGGCCGGTGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGCCCCTGCCGGCGTGa                        >  1:4391007/6‑150 (MQ=25)
                                                                                                                             cacctccGGCCTGGGCAGCGGCCTGGGCAGCGGCCTGGGCCGGTGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGCCCCTGCCGGCGTGa                        >  2:3458151/6‑150 (MQ=25)
                                                                                                                             cacctccGGCCTGGGCAGCGGCCTGGGCAGCGGCCTGGGCCGGTGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGGTCAGCAcc                                        <  2:2975912/129‑1 (MQ=21)
                                                                                                                             cacctccGGCCTGGGCAGCGGCCTGGGCAGCGGCCTGGGCCGGTGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGCTCAGCAcc                                        >  1:2975912/6‑134 (MQ=21)
                                                                                                                             aaGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGCCCCTGCCGGCGTGa                        <  2:4634319/150‑1 (MQ=255)
                                                                                                                                cACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGCCCCTGCCGGCGTGATGt                     <  1:5090722/150‑1 (MQ=255)
                                                                                                                                cACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGCCCCTGCCGGCGTGATGt                     >  1:5152836/1‑150 (MQ=255)
                                                                                                                                 aCCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGCCCCTGCCGGCGTGATGtt                    >  1:1673641/1‑150 (MQ=255)
                                                                                                                                  ccGGCCTGGGCAGCGGCCTGGGCAGCGGCCTGGGCCGGTGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCg                                                             >  2:2562005/1‑108 (MQ=21)
                                                                                                                                  ccGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGCCCCTGCCGGCGTGATGTTc                   <  1:1309089/150‑1 (MQ=255)
                                                                                                                                  ccGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGCCCCTGCCGGCGTGATGTTc                   <  1:4613591/150‑1 (MQ=255)
                                                                                                                                    ggCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGCCCCTGCCGGCGTGATGTTCGa                 <  2:3864165/150‑1 (MQ=255)
                                                                                                                                    ggCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGCCCCTGCCGGCGTGATGTTCGa                 <  2:1265133/150‑1 (MQ=255)
                                                                                                                                         aggccgctgcccaggccgctgcccAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGCCCCTGCCGGCGTGATGTTCGACACTc            >  1:4140410/1‑150 (MQ=255)
                                                                                                                                           gccgctgcccaggccgctgcccagGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGCCCCTGCCGGCGTGATGTTCGACACTCca          <  1:4911963/150‑1 (MQ=255)
                                                                                                                                             cgctgcccaggccgctgcccaggcCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGCCCCTGCCGGCGTGATGTTCGACACTCcaca        >  1:213197/1‑150 (MQ=255)
                                                                                                                                                tgcccaggccgctgcccaggccgGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGCCCCTGCCGGCGTGATGTTCGACACTCCAcagca     <  2:2624137/150‑1 (MQ=255)
                                                                                                                                                  cccaggccgctgcccaggccgGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGc                                      >  1:663732/1‑115 (MQ=255)
                                                                                                                                                   ccaggccgctgcccaggccgGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGCCCCTGCCGGCGTGATGTTCGACACTCCACAGCAGAt  >  1:3373920/1‑150 (MQ=255)
                                                                                                                                                    caggccgctgcccaggccgGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAg                                                        >  2:592391/1‑95 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
TGAGCAAGCCGCAGCGAAGATCGAGTACCAGCCGCTGCCGGCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGCCCCTGCCGGCGTGATGTTCGACACTCCACAGCAGAT  >  NC_002947/4882266‑4882562

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑