Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,882,426 C→T 5.6% Q357* (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,4260CT 5.6% 507.9 / 11.4 160Q357* (CAG→TAG) PP_4289membrane protein
Reads supporting (aligned to +/- strand):  ref base C (71/80);  new base T (6/3);  total (77/83)
Fisher's exact test for biased strand distribution p-value = 3.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 202 total aligned reads displayed.

CAGCGAAGATCGAGTACCAGCCGCTGCCGGCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGCCCCTGCCGGCGTGATGTTCGACACTCCACAGCAGATCCAGGCCCAGGCC  >  NC_002947/4882277‑4882575
                                                                                                                                                     |                                                                                                                                                     
  gCGAAGATCGAGTACCAGCCGCTGCCGGCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGgccgctgcccaggccgctgcccag                                                                                                                                                     >  2:4613591/1‑150 (MQ=255)
    gAAGATCGAGTACCAGCCGCTGCCGGCACCTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGACGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCcgctgcccaggccgctgcccaggc                                                                                                                                                   >  2:5440537/1‑150 (MQ=255)
    gAAGATCGAGTACCAGCCGCTGCCGGCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCcgctgcccaggccgctgcccaggc                                                                                                                                                   <  1:570933/150‑1 (MQ=255)
    gAAGATCGAGTACCAGCCGCTGCCGGCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCcgctgcccaggccgctgcccaggc                                                                                                                                                   >  1:3708070/1‑150 (MQ=255)
        aTCGAGTACCAGCCGCTGCCGGCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGa                                                                                                                                               >  2:1820694/1‑150 (MQ=255)
            aGTACCAGCCGCTGCCGGCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTg                                                                                                                                           >  1:3218709/1‑150 (MQ=255)
                ccAGCCGCTGCCGGCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGa                                                                                                                                       >  1:1310472/1‑150 (MQ=255)
                ccAGCCGCTGCCGGCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGa                                                                                                                                       >  1:1522143/1‑150 (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)
                      gCTGCCGGCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGa                                                                                                                                 <  1:2967005/150‑1 (MQ=255)
                      gCTGCCGGCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGa                                                                                                                                 <  2:2567596/150‑1 (MQ=255)
                        tGCCGGCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACa                                                                                                                               <  2:5262056/150‑1 (MQ=255)
                        tGCCGGCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACa                                                                                                                               <  1:232953/150‑1 (MQ=255)
                         gCCGGCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACaa                                                                                                                              <  2:2318507/150‑1 (MQ=255)
                         gCCGGCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACaa                                                                                                                              <  1:263252/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                                                                                                                        <  1:1489327/150‑1 (MQ=255)
                               aaCTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCca                                                                                                                        <  2:1421109/150‑1 (MQ=255)
                                  tGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACaa                                                                                                                     >  1:79156/1‑150 (MQ=255)
                                  tGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACaa                                                                                                                     >  1:79246/1‑150 (MQ=255)
                                  tGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACaa                                                                                                                     >  1:17079/1‑150 (MQ=255)
                                         ggcggcAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATc                                                                                                              <  2:3104682/150‑1 (MQ=255)
                                          gcggcAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATcc                                                                                                             >  1:4634319/1‑150 (MQ=255)
                                            ggcAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAg                                                                                                           >  2:2110102/1‑150 (MQ=255)
                                                   cGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCt                                                                                                    <  1:965769/150‑1 (MQ=255)
                                                      cTGCCGAGCAGCGTGGCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCa                                                                                                 >  1:4462612/1‑150 (MQ=255)
                                                            agcagcGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACtgtgt                                                                                           <  2:1123522/150‑1 (MQ=255)
                                                              cagcGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGc                                                                                         <  2:2613186/150‑1 (MQ=255)
                                                              cagcGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGc                                                                                         <  2:4288137/150‑1 (MQ=255)
                                                                gcGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCa                                                                                       <  2:2759322/150‑1 (MQ=255)
                                                                  gTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACt                                                                                     <  2:821028/150‑1 (MQ=255)
                                                                   tGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTc                                                                                    <  2:3708160/150‑1 (MQ=255)
                                                                     ccGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGt                                                                                  <  1:3034914/150‑1 (MQ=255)
                                                                      cGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCTAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTc                                                                                 <  1:835546/150‑1 (MQ=255)
                                                                      cGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTc                                                                                 <  2:3897231/150‑1 (MQ=255)
                                                                      cGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTc                                                                                 <  2:4972793/150‑1 (MQ=255)
                                                                      cGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGGGTGCCACTCGTc                                                                                 >  2:367596/1‑150 (MQ=255)
                                                                       gAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTcc                                                                                >  2:3339756/1‑150 (MQ=255)
                                                                        aGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCa                                                                               >  1:5211106/1‑150 (MQ=255)
                                                                             gccgccAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCg                                                                          >  2:5150091/1‑150 (MQ=255)
                                                                             gccgccAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCg                                                                          >  2:1117188/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)
                                                                                      gCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCa                                                                 >  1:4558929/1‑150 (MQ=255)
                                                                                        tGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACt                                                               >  2:2953229/1‑150 (MQ=255)
                                                                                        tGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACt                                                               <  2:3800155/150‑1 (MQ=255)
                                                                                         gCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAAATCGACAAGATCCACAATGTCAACCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTg                                                              >  2:2588068/1‑150 (MQ=255)
                                                                                          ccGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCTGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGACGACCAGCCCACTGt                                                             >  2:1182337/1‑150 (MQ=255)
                                                                                          ccGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGACTTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGt                                                             >  2:5072244/1‑150 (MQ=255)
                                                                                           cGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACtgt                                                                                             >  2:1537917/1‑117 (MQ=255)
                                                                                           cGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGtt                                                            <  1:4866514/150‑1 (MQ=255)
                                                                                           cGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGtt                                                            <  1:2511482/150‑1 (MQ=255)
                                                                                           cGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCAAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGtt                                                            >  1:1663886/1‑150 (MQ=255)
                                                                                            gAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTc                                                           <  2:4378325/150‑1 (MQ=255)
                                                                                            gAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTc                                                           <  2:3625912/150‑1 (MQ=255)
                                                                                             aaGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCa                                                          <  2:1671629/150‑1 (MQ=255)
                                                                                             aaGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCa                                                          <  1:2617199/150‑1 (MQ=255)
                                                                                                cACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCa                                                       >  1:3132807/1‑150 (MQ=255)
                                                                                                   ccgccgCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCg                                                    >  1:4233710/1‑150 (MQ=255)
                                                                                                       cgcACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGcccc                                                >  1:2803021/1‑150 (MQ=255)
                                                                                                        gcACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGCCCCt                                               >  2:1329979/1‑150 (MQ=255)
                                                                                                             tGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGCCCCTGCCgg                                          >  2:3910162/1‑150 (MQ=255)
                                                                                                                gccacctccGGCCTGGGCAGCGGCCTGGGCAGCGGCCTGGGCCGGTGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGCCCCTGCCGGCGt                                       >  2:2206550/8‑150 (MQ=25)
                                                                                                                  cacctccGGCCTGGGCAGCGGCCTGGGCAGCGGCCTGGGCCGGTGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGCCCCTGCCGGCGTGa                                     >  1:4391007/6‑150 (MQ=25)
                                                                                                                  cacctccGGCCTGGGCAGCGGCCTGGGCAGCGGCCTGGGCCGGTGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGGTCAGCAcc                                                     <  2:2975912/129‑1 (MQ=21)
                                                                                                                  aaGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGCCCCTGCCGGCGTGa                                     >  2:4788489/1‑150 (MQ=255)
                                                                                                                     cACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGCCCCTGCCGGCGTGATGt                                  >  1:3722770/1‑150 (MQ=255)
                                                                                                                     cACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGCCCCTGCCGGCGTGATGt                                  <  1:5090722/150‑1 (MQ=255)
                                                                                                                      aCCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGCCCCTGCCGGCGTGATGtt                                 >  1:1673641/1‑150 (MQ=255)
                                                                                                                      aCCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGCCCCTGCCGGCGTGATGtt                                 >  1:2811484/1‑150 (MQ=255)
                                                                                                                       ccGGCCTGGGCAGCGGCCTGGGCAGCGGCCTGGGCCGGTGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCg                                                                          <  1:2562005/108‑1 (MQ=21)
                                                                                                                       ccGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGCCCCTGCCGGCGTGATGTTc                                <  1:2336661/150‑1 (MQ=255)
                                                                                                                         ggCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGCCCCTGCCGGCGTGATGTTCGa                              <  2:3864165/150‑1 (MQ=255)
                                                                                                                         ggCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGCCCCTGCCGGCGTGATGTTCGa                              <  2:4681353/150‑1 (MQ=255)
                                                                                                                         ggCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGCCCCTGCCGGCGTGATGTTCGa                              <  1:4831957/150‑1 (MQ=255)
                                                                                                                         ggCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGCCCCTGCCGGCGTGATGTTCGa                              <  2:1265133/150‑1 (MQ=255)
                                                                                                                             caggccgctgcccaggccgctgccCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGCCCCTGCCGGCGTGATGTTCGACACt                          <  1:1329979/150‑1 (MQ=255)
                                                                                                                                  cgctgcccaggccgctgcccaggcCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGCCCCTGCCGGCGTGATGTTCGACACTCcaca                     >  1:213197/1‑150 (MQ=255)
                                                                                                                                    ctgcccaggccgctgcccaggccgGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGCCCCTGCCGGCGTGATGTTCGACACTCCAcagc                   <  1:3948304/150‑1 (MQ=255)
                                                                                                                                     tgcccaggccgctgcccaggccgGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGCCCCTGCCGGCGTGATGTTCGACACTCCAcagca                  <  2:2624137/150‑1 (MQ=255)
                                                                                                                                       cccaggccgctgcccaggccgGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGc                                                   <  2:663732/115‑1 (MQ=255)
                                                                                                                                        ccaggccgctgcccaggccgGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGCCCCTGCCGGCGTGATGTTCGACACTCCACAGCAGAt               >  1:3373920/1‑150 (MQ=255)
                                                                                                                                        ccaggccgctgcccaggccgGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGCCCCTGCCGGCGTGATGTTCGACACTCCACAGCAGAt               >  1:665034/1‑150 (MQ=255)
                                                                                                                                         caggccgctgcccaggccgGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAg                                                                     <  1:592391/95‑1 (MQ=255)
                                                                                                                                          aggccgctgcccaggccgGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGCCCCTGCCGGCGTGATGTTCGa                              >  1:4748909/1‑133 (MQ=255)
                                                                                                                                               gctgcccaggccgGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGCCCCTGCCGGCGTGATGTTCGACACTCCACAGCAGATccaggcc        <  2:4139865/150‑1 (MQ=255)
                                                                                                                                                 tgcccaggccgGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGc                                                                    <  1:2435034/88‑1 (MQ=255)
                                                                                                                                                 tgcccaggccgGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGc                                                                    >  2:2435034/1‑88 (MQ=255)
                                                                                                                                                 tgcccaggccgGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGCCCCTGCCGGCGTGATGTTCGACACTCCACAGCAGATccaggccca      >  1:1036508/1‑150 (MQ=255)
                                                                                                                                                 tgcccaggccgGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGCCCCTGCCGGCGTGATGTTCGACACTCCACAGCAGATccaggccca      >  1:1595877/1‑150 (MQ=255)
                                                                                                                                                   cccaggccgGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGCCCCTGCCGGCGTGATGTTCGACACTCCACAGCAGATccaggcccagg    >  2:1789469/1‑150 (MQ=255)
                                                                                                                                                    ccaggccgGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGCCCCTGCCGGCGTGATGTTCGACACTCCACAGCAGATccaggcccaggc   >  2:2039533/1‑150 (MQ=255)
                                                                                                                                                    ccaggccgGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGCCCCTGCCGGCGTGATGTTCGACACTCCACAGCAGATccaggcccaggc   <  1:1253431/150‑1 (MQ=255)
                                                                                                                                                     caggccgGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGCCCCTGCCGGCGTGATGTTCGACACTCCACAGCAGATCcaggcccaggcc  >  1:2195371/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CAGCGAAGATCGAGTACCAGCCGCTGCCGGCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGCCCCTGCCGGCGTGATGTTCGACACTCCACAGCAGATCCAGGCCCAGGCC  >  NC_002947/4882277‑4882575

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: