Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,882,419 C→A 5.1% A354A (GCC→GCA)  PP_4289 → membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,882,4190CA 5.1% 531.7 / 12.7 158A354A (GCC→GCA) PP_4289membrane protein
Reads supporting (aligned to +/- strand):  ref base C (70/80);  new base A (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 = 9.67e-01

Only 100 of 201 total aligned reads displayed.

CAAGCCGCAGCGAAGATCGAGTACCAGCCGCTGCCGGCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGCCCCTGCCGGCGTGATGTTCGACACTCCACAGCAGATCCAG  >  NC_002947/4882270‑4882566
                                                                                                                                                     |                                                                                                                                                   
  aGCCGCAGCGAAGATCGAGTACCAGCCGCTGCCGGCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGgcccaggccgctgcccaggccgc                                                                                                                                                   <  1:2603602/150‑1 (MQ=255)
         gCGAAGATCGAGTACCAGCCGCTGCCGGCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCAGCGGCctgggcagcggcctgggccggtgcttgc                                                                                                                                            <  2:3702180/150‑5 (MQ=38)
         gCGAAGATCGAGTACCAGCCGCTGCCGGCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCAGCGGCctgggcagcggcctgggccggtgcttgc                                                                                                                                            <  1:1802508/150‑5 (MQ=38)
         gCGAAGATCGAGTACCAGCCGCTGCCGGCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCcaggccgcggcctgggcagcggcctgg                                                                                                                                            <  1:4289907/150‑4 (MQ=255)
          cGAAGATCGAGTACCAGCCGCTGCCGGCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGccgctgcccaggccgctgcccagg                                                                                                                                           >  1:1658473/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)
               aTCGAGTACCAGCCGCTGCCGGCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCAAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGa                                                                                                                                      <  1:3908085/150‑1 (MQ=255)
                tCGAGTACCAGCCGCTGCCGGCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAg                                                                                                                                     >  1:2805618/1‑150 (MQ=255)
                   aGTACCAGCCGCTGCCGGCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTg                                                                                                                                  <  1:1517711/150‑1 (MQ=255)
                       ccAGCCGCTGCCGGCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGa                                                                                                                              >  1:1310472/1‑150 (MQ=255)
                        cAGCCGCTGCCGGCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGac                                                                                                                             >  2:242952/1‑149 (MQ=255)
                        cAGCCGCTGCCGGCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGaa                                                                                                                             >  1:1265133/1‑150 (MQ=255)
                        cAGCCGCTGCCGGCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGaa                                                                                                                             >  1:2900960/1‑150 (MQ=255)
                          gCCGCTGCCGGCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGACtt                                                                                                                           >  1:4969079/1‑150 (MQ=255)
                          gCCGCTGCCGGCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAAtt                                                                                                                           >  1:1671629/1‑150 (MQ=255)
                          gCCGCTGCCGGCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAAtt                                                                                                                           <  2:1310472/150‑1 (MQ=255)
                            cGCTGCCGGCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCg                                                                                                                         >  1:2632279/1‑150 (MQ=255)
                             gCTGCCGGCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGa                                                                                                                        <  2:2233519/150‑1 (MQ=255)
                             gCTGCCGGCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGa                                                                                                                        <  2:2567596/150‑1 (MQ=255)
                               tGCCGGCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACa                                                                                                                      <  1:232949/150‑1 (MQ=255)
                               tGCCGGCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACa                                                                                                                      <  2:5262056/150‑1 (MQ=255)
                                      aaCTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCca                                                                                                               <  2:1421109/150‑1 (MQ=255)
                                         tGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACaa                                                                                                            >  1:79246/1‑150 (MQ=255)
                                         tGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACaa                                                                                                            >  1:79156/1‑150 (MQ=255)
                                          gCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAAt                                                                                                           <  2:1699330/150‑1 (MQ=255)
                                                ggcggcAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATc                                                                                                     <  2:3104682/150‑1 (MQ=255)
                                                ggcggcAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATc                                                                                                     <  1:3092580/150‑1 (MQ=255)
                                                 gcggcAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATcc                                                                                                    >  1:4634319/1‑150 (MQ=255)
                                                   ggcAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAg                                                                                                  >  2:2110044/1‑150 (MQ=255)
                                                        gACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACg                                                                                             >  1:1480570/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)
                                                                   agcagcGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACtgtgt                                                                                  <  2:1123522/150‑1 (MQ=255)
                                                                     cagcGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGc                                                                                <  2:2613186/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)
                                                                             cGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCTAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTc                                                                        <  1:835546/150‑1 (MQ=255)
                                                                             cGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTc                                                                        <  2:3897231/150‑1 (MQ=255)
                                                                             cGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTc                                                                        <  2:4972793/150‑1 (MQ=255)
                                                                              gAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTcc                                                                       >  2:3339756/1‑150 (MQ=255)
                                                                              gAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTcc                                                                       >  2:4366945/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)
                                                                                    gccgccAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCg                                                                 >  2:4242487/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)
                                                                                       gccAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGAcc                                                              >  1:3325586/1‑150 (MQ=255)
                                                                                         cAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAg                                                            >  1:2381062/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)
                                                                                                  cGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACtgt                                                                                    <  1:1537917/117‑1 (MQ=255)
                                                                                                  cGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGtt                                                   <  1:4866514/150‑1 (MQ=255)
                                                                                                  cGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGtt                                                   <  1:2511482/150‑1 (MQ=255)
                                                                                                   gAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTc                                                  <  2:4378418/150‑1 (MQ=255)
                                                                                                   gAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTc                                                  <  2:3625912/150‑1 (MQ=255)
                                                                                                    aaGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCa                                                 <  2:1671629/150‑1 (MQ=255)
                                                                                                      gCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGc                                               <  1:3191571/150‑1 (MQ=255)
                                                                                                              cgcACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGcccc                                       >  1:2334625/1‑150 (MQ=255)
                                                                                                               gcACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGCCCCt                                      >  2:1329979/1‑150 (MQ=255)
                                                                                                                    tGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGCCCCTGCCgg                                 >  2:3910162/1‑150 (MQ=255)
                                                                                                                       gccacctccGGCCTGGGCAGCGGCCTGGGCAGCGGCCTGGGCCGGTGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGCCCCTGCCGGCGt                              >  2:2206393/8‑150 (MQ=25)
                                                                                                                       gcAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGCCCCTGCCGGCGt                              <  1:1079255/150‑1 (MQ=255)
                                                                                                                         cacctccGGCCTGGGCAGCGGCCTGGGCAGCGGCCTGGGCCGGTGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGCCCCTGCCGGCGTGa                            >  1:4391007/6‑150 (MQ=25)
                                                                                                                         cacctccGGCCTGGGCAGCGGCCTGGGCAGCGGCCTGGGCCGGTGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGGTCAGCAcc                                            <  2:2975912/129‑1 (MQ=21)
                                                                                                                         cacctccGGCCTGGGCAGCGGCCTGGGCAGCGGCCTGGGCCGGTGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGCTCAGCAcc                                            >  1:2975912/6‑134 (MQ=21)
                                                                                                                         aaGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGCCCCTGCCGGCGTGa                            >  2:4788489/1‑150 (MQ=255)
                                                                                                                         aaGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGCCCCTGCCGGCGTGa                            <  2:4634319/150‑1 (MQ=255)
                                                                                                                            cACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGCCCCTGCCGGCGTGATGt                         <  1:5090722/150‑1 (MQ=255)
                                                                                                                            cACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGCCCCTGCCGGCGTGATGt                         >  1:3722770/1‑150 (MQ=255)
                                                                                                                              ccGGCCTGGGCAGCGGCCTGGGCAGCGGCCTGGGCCGGTGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCg                                                                 <  1:2562005/108‑1 (MQ=21)
                                                                                                                                ggCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGCCCCTGCCGGCGTGATGTTCGa                     <  1:4831957/150‑1 (MQ=255)
                                                                                                                                ggCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGCCCCTGCCGGCGTGATGTTCGa                     <  2:3864165/150‑1 (MQ=255)
                                                                                                                                ggCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGCCCCTGCCGGCGTGATGTTCGa                     <  2:1265133/150‑1 (MQ=255)
                                                                                                                                    caggccgctgcccaggccgctgccCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGCCCCTGCCGGCGTGATGTTCGACACt                 <  1:1329979/150‑1 (MQ=255)
                                                                                                                                    caggccgctgcccaggccgctgccCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGCCCCTGCCGGCGTGATGTTCGACACt                 >  1:2624137/1‑150 (MQ=255)
                                                                                                                                     aggccgctgcccaggccgctgcccAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGCCCCTGCCGGCGTGATGTTCGACACTc                >  1:4140410/1‑150 (MQ=255)
                                                                                                                                     aggccgctgcccaggccgctgcccAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGCCCCTGCCGGCGTGATGTTCGACACTc                >  2:2342147/1‑150 (MQ=255)
                                                                                                                                     aggccgctgcccaggccgctgcccAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGCCCCTGCCGGCGTGATGTTCGACACTc                <  1:1691312/150‑1 (MQ=255)
                                                                                                                                         cgctgcccaggccgctgcccaggcCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGCCCCTGCCGGCGTGATGTTCGACACTCcaca            >  1:213197/1‑150 (MQ=255)
                                                                                                                                         cgctgcccaggccgctgcccaggcCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGCCCCTGCCGGCGTGATGTTCGACACTCcaca            >  2:3671714/1‑150 (MQ=255)
                                                                                                                                            tgcccaggccgctgcccaggccgGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGCCCCTGCCGGCGTGATGTTCGACACTCCAcagca         <  2:2624137/150‑1 (MQ=255)
                                                                                                                                              cccaggccgctgcccaggccgGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGc                                          <  2:663732/115‑1 (MQ=255)
                                                                                                                                               ccaggccgctgcccaggccgGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGCCCCTGCCGGCGTGATGTTCGACACTCCACAGCAGAt      >  1:665034/1‑150 (MQ=255)
                                                                                                                                                caggccgctgcccaggccgGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAg                                                            <  1:592391/95‑1 (MQ=255)
                                                                                                                                                 aggccgctgcccaggccgGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGCCCCTGCCGGCGTGATGTTCGa                     <  2:4748909/133‑1 (MQ=255)
                                                                                                                                                 aggccgctgcccaggccgGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGCCCCTGCCGGCGTGATGTTCGa                     >  1:4748909/1‑133 (MQ=255)
                                                                                                                                                  ggccgctgcccaggccgGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGCCCCTGCCGGCGTGATGTTCGACACTCCACAGCAGATcca   <  2:4558929/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CAAGCCGCAGCGAAGATCGAGTACCAGCCGCTGCCGGCAACTGCGGTAGGCGGCAAGACGGCTGCCGAGCAGCGTGCCGAGGACGCCGCCAAGGCTGCCGAAGCACCCGCCGCACCTGCGCAAGCACCGGCCCAGGCCGCTGCCCAGGCCGCTGCCCAGGCCGGAGGTGGCGAATTCGACAAGATCCACAATGTCATCCAGGAACGCTGCACTGTGTGCCACTCGTCCAAGCCGACCAGCCCACTGTTCAGCACCGCCCCTGCCGGCGTGATGTTCGACACTCCACAGCAGATCCAG  >  NC_002947/4882270‑4882566

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: ‑