Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,882,414 C→T 5.1% Q353* (CAG→TAG) ‡ PP_4289 → membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,882,4140CT 5.1% 505.6 / 9.6 158Q353* (CAG→TAG) ‡PP_4289membrane protein
Reads supporting (aligned to +/- strand):  ref base C (69/80);  new base T (6/2);  total (76/82)
Fisher's exact test for biased strand distribution p-value = 1.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 208 total aligned reads displayed.

TGAGCAAGCCGCAGCGAAGATCGAGTACCAGCCGCTGCCGGCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGCCCCTGCCGGCGTGATGTTCGACACTCCACAGCAGAT  >  NC_002947/4882266‑4882562
                                                                                                                                                    |                                                                                                                                                    
 gAGCAAGCCGCAGCGAAGATCGAGTACCAGCCGCTGCCGGCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGgcccaggccgctgcccag                                                                                                                                                    >  2:43776/1‑150 (MQ=255)
   gcctcacGCAGCGAAGATCGAGTACCAGCCGCTGCCGGCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGgcccaggccgctgcccaggc                                                                                                                                                  >  2:835546/7‑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)
           cAGCGAAGATCGAGTACCAGCCGCTGCCGGCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCaggccgctgcccaggccgctgccc                                                                                                                                          >  2:2779641/1‑150 (MQ=255)
            aGCGAAGATCGAGTACCAGCCGCTGCCGGCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAggccgctgcccaggccgctgccca                                                                                                                                         <  2:70631/150‑1 (MQ=255)
               gAAGATCGAGTACCAGCCGCTGCCGGCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCcgctgcccaggccgctgcccaggc                                                                                                                                      >  1:3708160/1‑150 (MQ=255)
                aaGATCGAGTACCAGCCGCTGCCGGCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCgctgcccaggccgctgcccaggcc                                                                                                                                     <  1:5432370/150‑1 (MQ=255)
                   aTCGAGTACCAGCCGCTGCCGGCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGa                                                                                                                                  >  2:1820694/1‑150 (MQ=255)
                   aTCGAGTACCAGCCGCTGCCGGCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGa                                                                                                                                  >  1:4496359/1‑150 (MQ=255)
                   aTCGAGTACCAGCCGCTGCCGGCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGa                                                                                                                                  <  1:2519594/150‑1 (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)
                            cAGCCGCTGCCGGCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGaa                                                                                                                         >  1:1265133/1‑150 (MQ=255)
                             aGCCGCTGCCGGCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAAt                                                                                                                        <  2:239408/150‑1 (MQ=255)
                              gCCGCTGCCGGCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGACtt                                                                                                                       >  1:4969079/1‑150 (MQ=255)
                              gCCGCTGCCGGCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAAtt                                                                                                                       <  2:1310472/150‑1 (MQ=255)
                               ccGCTGCCGGCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTc                                                                                                                      >  2:2967005/1‑150 (MQ=255)
                               ccGCTGCCGGCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTc                                                                                                                      <  2:3737378/150‑1 (MQ=255)
                                 gCTGCCGGCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGa                                                                                                                    <  1:2967005/150‑1 (MQ=255)
                                   tGCCGGCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACa                                                                                                                  <  1:2157997/150‑1 (MQ=255)
                                   tGCCGGCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACa                                                                                                                  >  2:4831957/1‑150 (MQ=255)
                                   tGCCGGCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACa                                                                                                                  <  2:5262056/150‑1 (MQ=255)
                                   tGCCGGCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGACGAATTCGACa                                                                                                                  <  2:2647509/150‑1 (MQ=255)
                                    gCCGGCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACaa                                                                                                                 <  1:263252/150‑1 (MQ=255)
                                    gCCGGCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACaa                                                                                                                 <  1:263027/150‑1 (MQ=255)
                                    gCCGGCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACaa                                                                                                                 <  2:2318507/150‑1 (MQ=255)
                                       ggCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGACTGCGACAAGag                                                                                                              >  2:4390410/1‑149 (MQ=255)
                                       ggCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGAt                                                                                                              >  1:3828215/1‑150 (MQ=255)
                                          aaCTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCca                                                                                                           <  1:1489327/150‑1 (MQ=255)
                                          aaCTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCca                                                                                                           <  2:4110233/150‑1 (MQ=255)
                                          aaCTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCca                                                                                                           <  2:1421109/150‑1 (MQ=255)
                                           aCTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCcac                                                                                                          <  1:585176/150‑1 (MQ=255)
                                             tGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACaa                                                                                                        >  1:17079/1‑150 (MQ=255)
                                              gCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAAt                                                                                                       <  2:1699330/150‑1 (MQ=255)
                                                  tAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCa                                                                                                   >  1:189377/1‑150 (MQ=255)
                                                     gcggcAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATcc                                                                                                >  1:4634319/1‑150 (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)
                                                                     cGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACtgt                                                                                <  1:2482455/150‑1 (MQ=255)
                                                                      gAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACtgtg                                                                               >  2:3312784/1‑150 (MQ=255)
                                                                       agcagcGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACtgtgt                                                                              <  2:1123522/150‑1 (MQ=255)
                                                                        gcagcGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACtgtgtg                                                                             <  1:4316512/150‑1 (MQ=255)
                                                                         cagcGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGc                                                                            <  2:4288137/150‑1 (MQ=255)
                                                                              tGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTc                                                                       <  2:3708070/150‑1 (MQ=255)
                                                                                ccGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGt                                                                     <  1:3034914/150‑1 (MQ=255)
                                                                                 cGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCTAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTc                                                                    <  1:835546/150‑1 (MQ=255)
                                                                                 cGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTc                                                                    <  2:3897017/150‑1 (MQ=255)
                                                                                 cGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGGGTGCCACTCGTc                                                                    >  2:367596/1‑150 (MQ=255)
                                                                                        gccgccAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCg                                                             >  2:1117188/1‑150 (MQ=255)
                                                                                        gccgccAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCg                                                             >  2:5150091/1‑150 (MQ=255)
                                                                                        gccgccAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCg                                                             >  2:5376148/1‑150 (MQ=255)
                                                                                        gccgccAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCg                                                             >  2:2617199/1‑150 (MQ=255)
                                                                                         ccgccAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGa                                                            <  1:5296674/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)
                                                                                                 gCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCa                                                    >  1:4558929/1‑150 (MQ=255)
                                                                                                  cTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCAc                                                   <  1:5223136/150‑1 (MQ=255)
                                                                                                     ccGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCTGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGACGACCAGCCCACTGt                                                >  2:1182337/1‑150 (MQ=255)
                                                                                                      cGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACtgt                                                                                >  2:1537917/1‑117 (MQ=255)
                                                                                                      cGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGtt                                               <  1:4866514/150‑1 (MQ=255)
                                                                                                      cGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCAAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGtt                                               >  1:1663886/1‑150 (MQ=255)
                                                                                                       gAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTc                                              >  1:674979/1‑150 (MQ=255)
                                                                                                       gAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTc                                              <  2:258439/150‑1 (MQ=255)
                                                                                                        aaGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCa                                             <  1:2617199/150‑1 (MQ=255)
                                                                                                        aaGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCa                                             <  2:1671629/150‑1 (MQ=255)
                                                                                                           cACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCa                                          >  1:3132807/1‑150 (MQ=255)
                                                                                                              ccgccgCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCg                                       >  1:4233710/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                                   >  1:2803021/1‑150 (MQ=255)
                                                                                                                   gcACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGCCCCt                                  >  2:1329979/1‑150 (MQ=255)
                                                                                                                        tGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGCCCCTGCCgg                             >  2:3909985/1‑150 (MQ=255)
                                                                                                                           gccacctccGGCCTGGGCAGCGGCCTGGGCAGCGGCCTGGGCCGGTGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGCCCCTGCCGGCGt                          >  2:2206550/8‑150 (MQ=25)
                                                                                                                           gcAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGCCCCTGCCGGCGt                          <  1:1079255/150‑1 (MQ=255)
                                                                                                                           gcAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGCCCCTGCCGGCGt                          <  1:4396325/150‑1 (MQ=255)
                                                                                                                             cacctccGGCCTGGGCAGCGGCCTGGGCAGCGGCCTGGGCCGGTGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGCCCCTGCCGGCGTGa                        >  2:3458151/6‑150 (MQ=25)
                                                                                                                             aaGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGCCCCTGCCGGCGTGa                        <  2:4634319/150‑1 (MQ=255)
                                                                                                                                cACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGCCCCTGCCGGCGTGATGt                     <  1:5090722/150‑1 (MQ=255)
                                                                                                                                cACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGCCCCTGCCGGCGTGATGt                     >  1:5152836/1‑150 (MQ=255)
                                                                                                                                 aCCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGCCCCTGCCGGCGTGATGtt                    >  1:2811484/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)
                                                                                                                                    ggCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGCCCCTGCCGGCGTGATGTTCGa                 <  2:3809560/150‑1 (MQ=255)
                                                                                                                                    ggCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGCCCCTGCCGGCGTGATGTTCGa                 <  2:1265133/150‑1 (MQ=255)
                                                                                                                                               ctgcccaggccgctgcccaggccgGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGCCCCTGCCGGCGTGATGTTCGACACTCCAcagc      <  1:3948304/150‑1 (MQ=255)
                                                                                                                                                tgcccaggccgctgcccaggccgGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGCCCCTGCCGGCGTGATGTTCGACACTCCAcagca     <  2:2811484/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: ‑