Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_002947 4,713,585 A→C T591P (ACC→CCC)  PP_4171 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,713,5850AC98.6% 745.1 / ‑0.5 216T591P (ACC→CCC) PP_4171hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (2/1);  new base C (94/117);  total (97/119)
Fisher's exact test for biased strand distribution p-value = 5.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 220 total aligned reads displayed.

CCTCGACCTGAGCCTGGATGTGAAGCTGTTCTACGACGAGATCACCGGCATGATCAGCGAGCCTCTGAGGAACAATCAGTACATCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCACCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAGCAACCCCGACGACCGTCGCCTCAGCGCCCGCAACAGCGGCTCGGCCGGGTGGATGCGCGAATGGGGTGCCGGCTGGTCGAGCGCGCTGTTCTAT  >  NC_002947/4713440‑4713731
                                                                                                                                                 |                                                                                                                                                  
cctcgacctGAGCCTGGATGTGAAGCTGTTCTACGACGAGATCACCGGCATGATCAGCGAGCCTCTGAGGAACAATCAGTACATCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCt                                                                                                                                                <  2:560924/150‑1 (MQ=255)
  tcgacctGAGCCTGGATGTGAAGCTGTTCTACGACGAGATCACCGGCATGATCAGCGAGCCTCTGAGGAACAATCAGTACATCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTgc                                                                                                                                              <  1:2345271/150‑1 (MQ=255)
  tcgacctGAGCCTGGATGTGAAGCTGTTCTACGACGAGATCACCGGCATGATCAGCGAGCCTCTGAGGAACAATCAGTACATCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTgc                                                                                                                                              >  2:3874632/1‑150 (MQ=255)
  tcgacctGAGCCTGGATGTGAAGCTGTTCTACGACGAGATCACCGGCATGATCAGCGAGCCTCTGAGGAACAATCAGTACATCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTgc                                                                                                                                              >  2:2131400/1‑150 (MQ=255)
  tcgacctGAGCCTGGATGTGAAGCTGTTCTACGACGAGATCACCGGCATGATCAGCGAGCCTCTGAGGAACAATCAGTACATCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTgc                                                                                                                                              >  2:2098507/1‑150 (MQ=255)
  tcgacctGAGCCTGGATGTGAAGCTGTTCTACGACGAGATCACCGGCATGATCAGCGAGCCTCTGAGGAACAATCAGTACATCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTgc                                                                                                                                              <  1:1291763/150‑1 (MQ=255)
       ctGAGCCTGGATGTGAAGCTGTTCTACGACGAGATCACCGGCATGATCAGCGAGCCTCTGAGGAACAATCAGTACATCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGAc                                                                                                                                         >  1:3855824/1‑150 (MQ=255)
       ctGAGCCTGGATGTGAAGCTGTTCTACGACGAGATCACCGGCATGATCAGCGAGCCTCTGAGGAACAATCAGTACATCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGAc                                                                                                                                         <  2:3268600/150‑1 (MQ=255)
            ccTGGATGTGAAGCTGTTCTACGACGAGATCACCGGCATGATCAGCGAGCCTCTGAGGAACAATCAGTACATCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACcgcct                                                                                                                                    <  2:930886/150‑1 (MQ=255)
              tGGATGTGAAGCTGTTCTACGACGAGATCACCGGCATGATCAGCGAGCCTCTGAGGAACAATCAGTACATCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTgc                                                                                                                                              >  2:1801852/1‑138 (MQ=255)
              tGGATGTGAAGCTGTTCTACGACGAGATCACCGGCATGATCAGCGAGCCTCTGAGGAACAATCAGTACATCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACcgcctgc                                                                                                                                  >  2:1981995/1‑150 (MQ=255)
              tGGATGTGAAGCTGTTCTACGACGAGATCACCGGCATGATCAGCGAGCCTCTGAGGAACAATCAGTACATCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACcgcctgc                                                                                                                                  >  1:115903/1‑150 (MQ=255)
               ggATGTGAAGCTGTTCTACGACGAGATCACCGGCATGATCAGCGAGCCTCTGAGGAACAATCAGTACATCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACcgcctgcg                                                                                                                                 >  2:2335669/1‑150 (MQ=255)
                   gtgAAGCTGTTCTACGACGAGATCACCGGCATGATCAGCGAGCCTCTGAGGAACAATCAGTACATCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACcgcctgcgcctg                                                                                                                             <  1:2878909/150‑1 (MQ=255)
                    tgAAGCTGTTCTACGACGAGATCACCGCCATGATCAGCGAGCCTCTGAGGAACAATCAGTACATCGCCAGCAACGCCAACACGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGa                                                                                                                            <  2:3431900/150‑1 (MQ=255)
                     gAAGCTGTTCTACGACGAGATCACCGGCATGATCAGCGAGCCTCTGAGGAACAATCAGTACATCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCGAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGAc                                                                                                                           >  1:1607396/1‑150 (MQ=255)
                     gAAGCTGTTCTACGACGAGATCACCGGCATGATCAGCGAGCCTCTGAGGAACAATCAGTACATCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGAc                                                                                                                           <  2:185009/150‑1 (MQ=255)
                         cTGTTCTACGACGAGATCACCGGCATGATCAGCGAGCCTCTGAGGAACAATCAGTACATCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTAt                                                                                                                       <  1:1240980/150‑1 (MQ=255)
                         cTGTTCTACGACGAGATCACCGGCATGATCAGCGAGCCTCTGAGGAACAATCAGTACATCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTAt                                                                                                                       <  1:2104221/150‑1 (MQ=255)
                          tGTTCTACGACGAGATCACCGGCATGATCAGCGAGCCTCTGAGGAACAATCAGTACATCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATg                                                                                                                      >  2:3761575/1‑150 (MQ=255)
                          tGTTCTACGACGAGATCACCGGCATGATCAGCGAGCCTCTGAGGAACAATCAGTACATCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATg                                                                                                                      >  1:3750243/1‑150 (MQ=255)
                           gTTCTACGACGAGATCACCGGCATGATCAGCGAGCCTCTGAGGAACAATCAGTACATCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGc                                                                                                                     <  1:1466060/150‑1 (MQ=255)
                              cTACGACGAGATCACCGGCATGATCAGCGAGCCTCTGAGGAACAATCAGTACATCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTa                                                                                                                  >  2:1044327/1‑150 (MQ=255)
                               tACGACGAGATCACCGGCATGATCAGCGAGCCTCTGAGGAACAATCAGTACATCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTAc                                                                                                                 <  1:1044327/150‑1 (MQ=255)
                                 cgacgaGATCACCGGCATGATCAGCGAGCCTCTGAGGAACAATCAGTACATCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGt                                                                                                               <  2:3386351/150‑1 (MQ=255)
                                  gacgaGATCACCGGCATGATCAGCGAGCCTCTGAGGAACAATCAGTACATCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTc                                                                                                              >  1:162327/1‑150 (MQ=255)
                                     gagaTCACCGGCATGATCAGCGAGCCTCTGAGGAACAATCAGTACATCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGAt                                                                                                           <  2:3008680/150‑1 (MQ=255)
                                       gaTCACCGGCATGATCAGCGAGCCTCTGAGGACCAATCAGTCCATCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGc                                                                                                         <  2:3560798/150‑1 (MQ=255)
                                       gaTCACCGGCATGATCAGCGAGCCTCTGAGGAACAATCAGTACATCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGc                                                                                                         >  1:3241321/1‑150 (MQ=255)
                                         tCACCGGCATGATCAGCGAGCCTCTGAGGAACAATCAGTACATCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGTGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCt                                                                                                       <  2:811443/150‑1 (MQ=255)
                                          cACCGGCATGATCAGCGAGCCTCTGAGGAACAATCAGTACATCGCCAGCAACGCCAACACGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTg                                                                                                      <  2:2845050/150‑1 (MQ=255)
                                           aCCGGCATGATCAGCGAGCCTCTGAGGAACAATCAGTACATCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTgg                                                                                                     >  2:4267017/1‑150 (MQ=255)
                                           aCCGGCATGATCAGCGAGCCTCTGAGGAACAATCAGTACATCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTgg                                                                                                     <  1:4339984/150‑1 (MQ=255)
                                           aCCGGCATGATCAGCGAGCCTCTGAGGAACAATCAGTACATCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTgg                                                                                                     <  2:4211722/150‑1 (MQ=255)
                                            ccGGCATGATCAGCGAGCCTCTGAGGAACAATCAGTACATCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTggg                                                                                                    <  1:2933202/150‑1 (MQ=255)
                                                   gATCAGCGAGCCTCTGAGGAACAATCAGTACATCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAGCaa                                                                                             >  1:152441/1‑150 (MQ=255)
                                                    aTCAGCGAGCCTCTGAGGAACAATCAGTACATCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAGCAAc                                                                                            <  1:4150956/150‑1 (MQ=255)
                                                     tCAGCGAGCCTCTGAGGAACAATCAGTACATCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAGCAAcc                                                                                           <  2:115903/150‑1 (MQ=255)
                                                     tCAGCGAGCCTCAGAGGAACAAACAGTACATCGCAAGCAACGCCAACAAGGCCCGTTTCAGCGGAAGCGAGGCTCAAAAGGACTCGCGCAACCCCCAGCGCGACCGCCTGCCCCAGACCTATGCCTACGTCGATGCCAGGGCCCGCAAcc                                                                                           >  2:4254654/1‑150 (MQ=255)
                                                      cAGCGAGCCTCTGAGGAACAATCAGTACATCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAGCAAccc                                                                                          <  1:3115590/150‑1 (MQ=255)
                                                        gCGAGCCTCTGAGGAACAATCAGTACATCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAGCAACCCcg                                                                                        >  1:4575497/1‑150 (MQ=255)
                                                         cGAGCCTCTGAGGAACAATCAGTACATCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAGCAACCCcga                                                                                       <  1:3268947/150‑1 (MQ=255)
                                                              ctctGAGGAACAATCAGTACATCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAGCAACCCCGACGAcc                                                                                  <  1:836753/150‑1 (MQ=255)
                                                                 tGAGGAACAATCAGTACATCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAGCAACCCCGACGACCGTc                                                                               <  1:172653/150‑1 (MQ=255)
                                                                 tGAGGAACAATCAGTACATCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAGCAACCCCGACGACCGTc                                                                               >  2:1333644/1‑150 (MQ=255)
                                                                 tGAGGAACAATCAGTACATCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAGCAACCCCGACGACCGTc                                                                               >  2:4015570/1‑150 (MQ=255)
                                                                 tGAGGAACAATCAGTACATCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAGCAACCCCGACGACCGTc                                                                               <  1:3441083/150‑1 (MQ=255)
                                                                 tGAGGAACAATCAGTACATCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAGCAACCCCGACGACCGTc                                                                               <  1:2521226/150‑1 (MQ=255)
                                                                   aGGAACAATCAGTACATCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAGCAACCCCGACGACCGTCGc                                                                             <  2:951914/150‑1 (MQ=255)
                                                                     gAACAATCAGTACATCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAGCAACCCCGACGACCATCGCCt                                                                           >  2:4104939/1‑150 (MQ=255)
                                                                       aCAATCAGTACATCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAGCAACCCCGACGACCGTCGCCTCa                                                                         <  1:2028425/150‑1 (MQ=255)
                                                                       aCAATCAGTACATCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAGCAACCCCGACGACCGTCGCCTCa                                                                         <  1:3529485/150‑1 (MQ=255)
                                                                       aCAATCAGTACATCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAGCAACCCCGACGACCGTCGCCTCa                                                                         <  1:3803026/150‑1 (MQ=255)
                                                                         aaTCAGTACATCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAGCAACCCCGACGACCGTCGCCTCAgc                                                                       >  1:2036010/1‑150 (MQ=255)
                                                                         aaTCAGTACATCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAGCAACCCCGACGACCGTCGCCTCAgc                                                                       <  1:4406566/150‑1 (MQ=255)
                                                                         aaTCAGTACATCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAGCAACCCCGACGACCGTCGCCTCAgc                                                                       >  1:4040673/1‑150 (MQ=255)
                                                                          aTCAGTACATCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAGCAACCCCGACGACCGTCGCCTCAgcg                                                                      >  2:1031067/1‑150 (MQ=255)
                                                                           tCAGTACATCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAGCAACCCCGACGACCGTCGCCTCAgcgc                                                                     >  1:2109903/1‑150 (MQ=255)
                                                                                aCATCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAGCAACCCCGACGACCGTCGCCTCAGCGCCCGCa                                                                >  1:3193985/1‑150 (MQ=255)
                                                                                   tCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAGCAACCCCGACGACCGTCGCCTCAGCGCCCGCAACa                                                             >  1:2639615/1‑150 (MQ=255)
                                                                                   tCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAGCAACCCCGACGACCGTCGCCTCAGCGCCCGCAACa                                                             >  1:143719/1‑150 (MQ=255)
                                                                                     gCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAGCAACCCCGACGACCGTCGCCTCAGCGCCCGCAACAGc                                                           >  2:1372898/1‑150 (MQ=255)
                                                                                     gCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAGCAACCCCGACGACCGTCGCCTCAGCGCCCGCAACAGc                                                           >  2:74282/1‑150 (MQ=255)
                                                                                      ccAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAGCAACCCCGACGACCGTCGCCTCAGCGCCCGCAACAGCg                                                          >  1:4635931/1‑150 (MQ=255)
                                                                                         gCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAGCAACCCCGACGACCGTCGCCTCAGCGCCCGCAACAGCGGCt                                                       >  1:1701699/1‑150 (MQ=255)
                                                                                           aaCGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAGCAACCCCGACGACCGTCGCCTCAGCGCCCGCAACAGCGGCTCg                                                     <  2:906944/150‑1 (MQ=255)
                                                                                              gCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAGCAACCCCGACGACCGTCGCCTCAGCGCCCGCAACAGCGGCTCGGcc                                                  >  2:637913/1‑150 (MQ=255)
                                                                                                cgcaaaGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAGCAACCCCGACGACCGTCGCCTCAGCGCCCGCAACAGCGGCTCGGCCgg                                                >  2:3665588/5‑150 (MQ=255)
                                                                                                 aacaaGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAGCAACCCCGACGACCGTCGCCTCAGCGCCCGCAACAGCGGCTCGGCCggg                                               <  2:1600943/150‑1 (MQ=255)
                                                                                                   caaGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAGCAACCCCGACGACCGTCGCCTCAGCGCCCGCAACAGCGGCTCGGCCGGGTg                                             >  1:3168166/1‑150 (MQ=255)
                                                                                                   caaGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAGCAACCCCGACGACCGTCGCCTCAGCGCCCGCAACAGCGGCTCGGCCGGGTg                                             >  1:2591471/1‑150 (MQ=255)
                                                                                                    aaGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAGCAACCCCGACGACCGTCGCCTCAGCGCCCGCAACAGCGGCTCGGCCGGGTgg                                            >  1:1707611/1‑150 (MQ=255)
                                                                                                        cccGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAGCAACCCCGACGACCGTCGCCTCAGCGCCCGCAACAGCGGCTCGGCCGGGTGGATgc                                        <  1:2861826/150‑1 (MQ=255)
                                                                                                          cGTTTCAGCGGCAGCGAGGCGCAACTGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAGCAACCCCGACGACCGTCGCCTCAGCGCCCGCAACAGCGGCTCGGCCGGGTGGATgcgc                                      <  1:1265331/150‑1 (MQ=255)
                                                                                                          cGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAGCAACCCCGACGACCGTCGCCTCAGCGCCCGCAACAGCGGCTCGGCCGGGTGGATgcgc                                      <  2:1481052/150‑1 (MQ=255)
                                                                                                          cGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCCGGGCCAGCAACCCCGACGACCGTCGCCTCAGCGCCCGCAACAGCGGCTCGGCCGGGTGGATgcgc                                      <  2:2484137/150‑1 (MQ=255)
                                                                                                             ttCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAGCAACCCCGACGACCGTCGCCTCAGCGCCCGCAACAGCGGCTCGGCCGGGTGGATGCGCGaa                                   <  1:3480182/150‑1 (MQ=255)
                                                                                                               cAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAGCAACCCCGACGACCGTCGCCTCAGCGCCCGCAACAGCGGCTCGGCCGGGTGGATGCGCGAATg                                 <  2:1831304/150‑1 (MQ=255)
                                                                                                                aGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAGCAACCCCGACGACCGTCGCCTCAGCGCCCGCAACAGCGGCTCGGCCTGGTGGATGCGCGAATgg                                <  2:152441/150‑1 (MQ=255)
                                                                                                                aGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAGCAACCCCGACGACCGTCGCCTCAGCGCCCGCAACAGCGGCTCGGCCGGGTGGATGCGCGAATgg                                <  2:161627/150‑1 (MQ=255)
                                                                                                                aGCGGCAGCGAGGCACAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAGCAACCCCGACGACCGTCGCCTCAGCGCCCGCAACAGCGGCTCGGCCGGGTGGATGCGCGAATgg                                <  2:2216522/150‑1 (MQ=255)
                                                                                                                 gCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAGCAACCCCGACGACCGTCGCCTCAGCGCCCGCAACAGCGGCTCGGCCGGGTGGATGCGCGAATggg                               >  1:2223212/1‑150 (MQ=255)
                                                                                                                 gCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAGCAACCCCGACGACCGTCGCCTCAGCGCCCGCAACAGCGGCTCGGCCGGGTGGATGCGCGAATggg                               <  2:1272785/150‑1 (MQ=255)
                                                                                                                   ggCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAGCAACCCCGACGACCGTCGCCTCAGCGCCCGCAACAGCGGCTCGGCCGGGTGGATGCGCGAATGGGGt                             <  2:2411951/150‑1 (MQ=255)
                                                                                                                   ggCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAGCAACCCCGACGACCGTCGCCTCAGCGCCCGCAACAGCGGCTCGGCCGGGTGGATGCGCGAATGGGGt                             <  2:66244/150‑1 (MQ=255)
                                                                                                                      aGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAGCAACCCCGACGACCGTCGCCTCAGCGCCCGCAACAGCGGCTCGGCCGGGTGGATGCGCGAATGGGGTGcc                          >  2:3079097/1‑150 (MQ=255)
                                                                                                                        cGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAGCAAccc                                                                                          <  2:4458955/84‑1 (MQ=255)
                                                                                                                           ggCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAGCAACCCCGACGACCGTCGCCTCAGCGCCCGCAACAGCGGCTCGGCCGGGTGGATGCGCGAATGGGGTGCCGGCTg                     <  2:4306680/150‑1 (MQ=255)
                                                                                                                              gcAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAGCAACCCCGACGACCGTCGCCTCAGCGCCCGCAACAGCGGCTCGGCCGGGTGGATGCGCGAATGGGGTGCCGGCTGGTc                  >  2:3858918/1‑150 (MQ=255)
                                                                                                                              gcAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAGCAACCCCGACGACCGTCGCCTCAGCGCCCGCAACAGCGGCTCGGCCGGGTGGATGCGCGAATGGGGTGCCGGCTGGTc                  <  1:1141267/150‑1 (MQ=255)
                                                                                                                              gcAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCCATGCCTGGGCCAGCAACCCCGACGACCGTCGCCTCAGCGCCCGCAACAGCGGCTCGGCCGGGTGGATGCGCGAATGGGGTGCCGGCTGGTc                  >  1:4521244/1‑150 (MQ=255)
                                                                                                                                aaaTGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGAt                                                                                                           <  2:3153788/59‑1 (MQ=255)
                                                                                                                                  aTGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAGCAACCCCGACGACCGTCGCCTCAGCGCCCGCAACAGCGGCTCGGCCGGGTGGATGCGCGAATGGGGTGCCGGCTGGTCGAgc              <  2:1701699/150‑1 (MQ=255)
                                                                                                                                          gcgcCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAGCAACCCCGACGACCGTCGCCTCAGCGCCCGCAACAGCGGCTCGGCCGGGTGGATGCGCGAATGGGGTGCCGGCTGGTCGAGCGCGCTGtt      >  1:2849368/1‑150 (MQ=255)
                                                                                                                                          gcgcCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAGCAACCCCGACGACCGTCGCCTCAGCGCCCGCAACAGCGGCTCGGCCGGGTGGATGCGCGAATGGGGTGCCGGCTGGTCGAGCGCGCTGtt      >  2:2036057/1‑150 (MQ=255)
                                                                                                                                            gcCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAGCAACCCCGACGACCGTCGCCTCAGCGCCCGCAACAGCGGCTCGGCCGGGTGGATGCGCGAATGGGGTGCCGGCTGGTCGAGCGCGCTGTTCt    >  2:372879/1‑150 (MQ=255)
                                                                                                                                            gcCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAGCAACCCCGACGACCGTCGCCTCAGCGCCCGCAACAGCGGCTCGGCCGGGTGGATGCGCGAATGGGGTGCCGGCTGGTCGAGCGCGCTGTTCt    >  1:1053374/1‑150 (MQ=255)
                                                                                                                                             ccccccccTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAGCAACCCCGACGACCGTCGCCTCAGCGCCCGCAACAGCGGCTCGGCCGGGTGGATGCGCGAATGGGGTGCCGGCTGGTCGAGCGCGCTGTTCTa   >  1:1899172/1‑150 (MQ=255)
                                                                                                                                              ccccccCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAGCAACCCCGACGACCGTCGCCTCAGCTCCCGCAACAGCGGCTCGTCCGGGTGGATGCGCGAATGGTGTGCCGGCTGGTCGAGCGCGCTGTTCTAt  <  1:4279410/150‑1 (MQ=255)
                                                                                                                                              ccccccCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGGCGATGCCTGGGCCAGCAACCCCGACGACCGTCGCCTCAGCGCCCGCAACAGCGGCTCGGCCGGGTGGATGCGCGAATGGGGTGCCGGCTGGTCGAGCGCGCGGTTCTAt  >  2:2054850/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
CCTCGACCTGAGCCTGGATGTGAAGCTGTTCTACGACGAGATCACCGGCATGATCAGCGAGCCTCTGAGGAACAATCAGTACATCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCACCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAGCAACCCCGACGACCGTCGCCTCAGCGCCCGCAACAGCGGCTCGGCCGGGTGGATGCGCGAATGGGGTGCCGGCTGGTCGAGCGCGCTGTTCTAT  >  NC_002947/4713440‑4713731

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