Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 311,734 G→C 11.2% T215S (ACC→AGC) ‡ PP_0256 ← putative oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947311,7340GC11.2% 579.8 / 51.5 188T215S (ACC→AGC) ‡PP_0256putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (74/93);  new base C (6/15);  total (80/108)
Fisher's exact test for biased strand distribution p-value = 2.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.37e-01

Only 100 of 189 total aligned reads displayed.

GACCACGCTCCTTGAGCGGGTTGATGCACACCACCTGGGCGCCACGCTTTACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTGACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTACGCCGCTTCGTTGCTGGCACGTCC  >  NC_002947/311588‑311882
                                                                                                                                                   |                                                                                                                                                     
gACCACGCTCCTTGAGCGGGTTGATGCACACCACCTGGGCGCCACGCTTTACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTGa                                                                                                                                                   >  2:2031952/1‑150 (MQ=255)
 aCCACGCTCCTTGAGCGGGTTGATGCACACCACCTGGGCGCCACGCTTTACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTGAc                                                                                                                                                  >  1:3603045/1‑150 (MQ=255)
   cACGCTCCTTGAGCGGGTTGATGCACACCACCTGGGCGCCACGCTTTACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTGACTg                                                                                                                                                >  1:1208242/1‑150 (MQ=255)
    aCGCTCCTTGAGCGGGTTGATGCACACCACCTGGGCGCCACGCTTTACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTGACTGt                                                                                                                                               <  1:651250/150‑1 (MQ=255)
     cGCTCCTTGAGCGGGTTGATGCACACCACCTGGGCGCCACGCTTTACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTGACTGTa                                                                                                                                              >  1:1338981/1‑150 (MQ=255)
            tGAGCGGGTTGATGCACACCACCTGGGCGCCACGCTTTACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTGACTGTACCTTTGc                                                                                                                                       <  1:2934812/150‑1 (MQ=255)
             gAGCGGGTTGATGCACACCACCTGGGCGCCACGCTTTACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTGACTGTACCTTTGcc                                                                                                                                      <  2:4894692/150‑1 (MQ=255)
                cGGGTTGATGCACACCACCTGGGCGCCACGCTTTACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTGACTGTACCTTTGCCGAc                                                                                                                                   <  1:5400436/150‑1 (MQ=255)
                  ggTTGATGCACACCACCTGGGCGCCACGCTTTACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTGACTGTACCTTTGCCGAccc                                                                                                                                 >  2:2399458/1‑150 (MQ=255)
                       aTGCACACCACCTGGGCGCCACGCTTTACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTGACTGTACCTTTGCCGACCCCGACg                                                                                                                            >  2:1142252/1‑150 (MQ=255)
                        tGCACACCACCTGGGCGCCACGCTTTACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTGACTGTACCTTTGCCGACCCCGACGc                                                                                                                           >  2:537889/1‑150 (MQ=255)
                          cacaCCACCTGGGCGCCACGCTTTACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTGACTGTACCTTTGCCGACCCCGACGCtt                                                                                                                         >  2:5285238/1‑150 (MQ=255)
                          cacaCCACCTGGGCGCCACGCTTTACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTGACTGTACCTTTGCCGACCCCGACGCtt                                                                                                                         <  1:2618999/150‑1 (MQ=255)
                           acacCACCTGGGCGCCACGCTTTACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTGACTGTACCTTTGCCGACCCCGACGCttt                                                                                                                        >  1:1042701/1‑150 (MQ=255)
                           acacCACCTGGGCGCCACGCTTTACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTGACTGTACCTTTGCCGACCCCGACGCttt                                                                                                                        >  1:1092420/1‑150 (MQ=255)
                            caccacCTGGGCGCCACGCTTTACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTGACTGTACCTTTGCCGACCCCGACGCTTTg                                                                                                                       >  1:3241060/1‑150 (MQ=255)
                             accaccTGGGCGCCACGCTTTACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTGACTGTACCTTTGCCGACCCCGACGCTTTGc                                                                                                                      <  2:2306839/150‑1 (MQ=255)
                             accaccTGGGCGCCACGCTTTACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTGACTGTACCTTTGCCGACCCCGACGCTTTGc                                                                                                                      <  1:5233092/150‑1 (MQ=255)
                              ccaccTGGGCGCCACGCTTTACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTGACTGTACCTTTGCCGACCCCGACGCTTTGcc                                                                                                                     <  1:591406/150‑1 (MQ=255)
                                 ccTGGGCGCCACGCTTTACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTGACTGTACCTTTGCCGACCCCGACGCTTTGCCCCa                                                                                                                  >  1:5357177/1‑150 (MQ=255)
                                  cTGGGCGCCACGCTTTACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTGACTGTACCTTTGCCGACCCCGACGCTTTGCCCCaa                                                                                                                 <  1:903978/150‑1 (MQ=255)
                                   tGGGCGCCACGCTTTACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTGACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAac                                                                                                                >  2:694469/1‑150 (MQ=255)
                                   tGGGCGCCACGCTTTACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTGACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAac                                                                                                                <  2:3829754/150‑1 (MQ=255)
                                       tgcCACGCTTTACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTGACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAacgcca                                                                                                            <  2:2201134/149‑1 (MQ=255)
                                        gcCACGCTTTACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACCAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTGACTGTACCTTTGCCGACCCCGACGCTTTGCCCCCacgccac                                                                                                           >  2:4342832/1‑150 (MQ=255)
                                        gcCACGCTTTACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTGACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAacgccac                                                                                                           >  2:1879979/1‑150 (MQ=255)
                                          cACGCTTTACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTGACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAacgccacgc                                                                                                         >  2:5347513/1‑150 (MQ=255)
                                           aCGCTTTACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTGACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAacgccacgcc                                                                                                        >  1:3475558/1‑150 (MQ=255)
                                              cTTTACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTGACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCt                                                                                                     <  2:3756026/150‑1 (MQ=255)
                                               tttACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTGACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTg                                                                                                    <  1:5111275/150‑1 (MQ=255)
                                               tttACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTGACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTg                                                                                                    <  2:1284917/150‑1 (MQ=255)
                                                 tACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTGACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGc                                                                                                  <  1:3753003/150‑1 (MQ=255)
                                                     gCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTGACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCa                                                                                              >  2:2858569/1‑150 (MQ=255)
                                                      cATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTGACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCAt                                                                                             >  1:284610/1‑150 (MQ=255)
                                                        tCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTGACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATgg                                                                                           >  1:4394451/1‑150 (MQ=255)
                                                               gTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTGACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGt                                                                                    >  1:412278/1‑150 (MQ=255)
                                                                 gggTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTGACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTc                                                                                  <  2:4862958/150‑1 (MQ=255)
                                                                 gggTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTGACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTc                                                                                  <  2:860692/150‑1 (MQ=255)
                                                                   gTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTGACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGa                                                                                <  1:1185687/150‑1 (MQ=255)
                                                                   gTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTGACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGa                                                                                <  1:1916050/150‑1 (MQ=255)
                                                                   gTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTGACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGa                                                                                <  2:2790975/150‑1 (MQ=255)
                                                                    cccaacATGCGGGGATGATTGGTGCCAGGATTTTGCCCCACGACGAAAATCGCGTCAGCCTGCTCGAAATCCTCAAGCAGAACAGTTCCCTTGCCAACCCCTACCGATTGTTTGAGCGCCACACCGCTGGCCTCATGGCACATGTTCGaa                                                                               <  2:2088191/145‑2 (MQ=255)
                                                                      gAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTGACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCa                                                                             >  1:4515812/1‑150 (MQ=255)
                                                                       aacATGCGGGGATGATTGGTGCCAGGATTTTGCCCCACGACGAAAATCGCGTCAGCCTGCTCGAAATCCTCAAGCAGAACAGTTCCCTTGCCAACCCCTACCGATTGTTTGAGCGCCACACCGCTGGCCTCATGGCACATGTTCGAACAg                                                                            >  2:3384367/3‑150 (MQ=255)
                                                                       aGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTGACTGTACCTTTGCCGACCCCGACGCTTTGcc                                                                                                                     <  2:3589351/109‑1 (MQ=255)
                                                                       aGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTGACTGTACCTTTGCCGACCCCGACGCTTTGcc                                                                                                                     >  1:3589351/1‑109 (MQ=255)
                                                                        gCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTGACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGt                                                                           <  2:3819626/150‑1 (MQ=255)
                                                                        gCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTGACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGt                                                                           <  1:5361846/150‑1 (MQ=255)
                                                                        gCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTGACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGt                                                                           >  1:361604/1‑150 (MQ=255)
                                                                        acATGCGGGGATGATTGGTGCCAGGATTTTGCCCCACGACGAAAATCGCGTCAGCCTGCTCGAAATCCTCAAGCAGAACAGGTCCCTTGCCAACCCCTACCGATTGTTTGAGCGCCACACCGCTGGCCTCATGGCACATGTTCGAACAGt                                                                           <  2:4781816/149‑1 (MQ=255)
                                                                         cATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTGACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTc                                                                          <  2:361604/150‑1 (MQ=255)
                                                                          aTGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTGACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCa                                                                         <  2:1680453/150‑1 (MQ=255)
                                                                           tGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTGACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAg                                                                        <  1:2803231/150‑1 (MQ=255)
                                                                                ggatGATTGGTGCCAGGATTTTGCCCCACGACGAAAATCGCGTCAGCCTGCTCGAAATCCTCAAGCAGAACAGTTCCCTTGCCAACCCCTACCGATTGTTTGAGCGCCACACCGCTGGCCTCATGGCACATGTTCGAACAGTCAGGGaaa                                                                   <  1:2939829/147‑2 (MQ=255)
                                                                                 ggTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTGACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATgg                                                                                           <  1:2296565/125‑1 (MQ=255)
                                                                                 ggTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTGACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATgg                                                                                           >  2:2296565/1‑125 (MQ=255)
                                                                                 ggTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTGACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAgt                                                                  >  2:4015479/1‑150 (MQ=255)
                                                                                  atGATTGGTGCCAGGATTTTGCCCCACGACGAAAATCGCGTCAGCCTGCTCGAAATCCTCAAGCAGAACAGTTCCCTTGCCAACCCCTACCGATTGTTTGAGCGCCACACCGCTGGCCTCATGGCACATGTTCGAACAGTCAGGGAAatt                                                                 <  2:1112865/149‑1 (MQ=255)
                                                                                    gATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTGACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAgttgt                                                               <  2:1913405/150‑1 (MQ=255)
                                                                                      ttGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTGACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGttgttg                                                             <  1:4602288/150‑1 (MQ=255)
                                                                                      ttGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTGACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGttgttg                                                             <  1:379824/150‑1 (MQ=255)
                                                                                      ttGGTGCCAGGATTTTGCCCCACGACGAAAATCGCGTCAGCCTGCTCGAAATCCTCAAGCAGAACAGTTCCCTTGCCAACCCCTACCGATTGTTTGAGCGCCACACCGCTGGCCTCATGGCACATGTTCGAACAGTCAGGGAAattgttc                                                             <  1:2149738/150‑2 (MQ=255)
                                                                                       tGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTGACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTc                                                                                  <  2:187281/128‑1 (MQ=255)
                                                                                       tGGTGCCAGGATTTTGCCCCACGACGAAAATCGCGTCAGCCTGCTCGAAATCCTCAAGCAGAACAGTTCCCTTGCCAACCCCTACCGATTGTTTGAGCGCCACACCGCTGGCCTCATGGCACATGTTCGAACAGTCAGGGAAattgttcg                                                            <  1:1847330/150‑3 (MQ=255)
                                                                                       tGGTGCCAGGATTTTGCCCCACGACGAAAATCGCGTCAGCCTGCTCGAAATCCTCAAGCAGAACAGTTCCCTTGCCAACCCCTACCGATTGTTTGAGCGCCACACCGCTGGCCTCATGGCACATGTTCGAACAGTCAGGGAAattgttcg                                                            >  2:1251424/1‑148 (MQ=255)
                                                                                         gTGCCAGGATTTTGCCCCACGACGAAAATCGCGTCAGCCTGCTCGAAATCCTCAAGCAGAACAGTTCCCTTGCCAACCCCTACCGATTGTTTGAGCGCCACACCGCTGGCCTCATGGCACATGTTCGAACAGTCAGGGAAATTGTTCGTg                                                          >  1:3339559/1‑150 (MQ=255)
                                                                                          tGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTGACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGc                                                         <  2:2686471/150‑1 (MQ=255)
                                                                                           gCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTGACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGcc                                                        <  2:5064580/150‑1 (MQ=255)
                                                                                           gCCGGGGTTCCGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTGACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGcc                                                        <  2:3147046/150‑1 (MQ=255)
                                                                                            ccGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTGACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCg                                                       >  1:3681481/1‑150 (MQ=255)
                                                                                               gggTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTGACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAg                                                    >  1:1333434/1‑150 (MQ=255)
                                                                                                ggTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTGACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAgg                                                   >  2:3872207/1‑150 (MQ=255)
                                                                                                    tttgCCCCACGACGAAAATCGCGTCAGCCTGCTCGAAATCCTCAAGCAGAACAGTTCCCTTGCCAACCCCTACCGATTGTTTGAGCGCCACACCGCTGGCCTCATGGCACATGTTCGAACAGTCAGGGAAATTGTTCGTGCCATAAgcgc                                                <  1:402085/147‑1 (MQ=255)
                                                                                                     tGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTGACTGTACCTTTGTCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGgcgcgc                                              >  2:3096652/1‑150 (MQ=255)
                                                                                                      ggCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTGACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCa                                             >  1:3503488/1‑150 (MQ=255)
                                                                                                       gCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTGACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCAc                                            <  1:694469/150‑1 (MQ=255)
                                                                                                          ccAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTGACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGaa                                         >  1:4222233/1‑150 (MQ=255)
                                                                                                          ccAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTGACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGaa                                         >  2:4420022/1‑150 (MQ=255)
                                                                                                           cAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTGACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAAc                                        <  2:1494475/150‑1 (MQ=255)
                                                                                                               aCAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTGACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGt                                     >  1:2511125/1‑150 (MQ=255)
                                                                                                                cAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTGACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGtt                                    <  1:3196537/150‑1 (MQ=255)
                                                                                                                cAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTGACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGtt                                    >  2:1315288/1‑150 (MQ=255)
                                                                                                                cAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTGACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGtt                                    <  2:2308561/150‑1 (MQ=255)
                                                                                                                  aaaaTGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTGACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTgg                                  <  2:1092420/150‑1 (MQ=255)
                                                                                                                   aaaTGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTGACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGt                                 <  2:523996/150‑1 (MQ=255)
                                                                                                                   aaaTGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTGACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGt                                 <  1:4886850/150‑1 (MQ=255)
                                                                                                                     aTGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTGACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTAc                               <  1:3163121/150‑1 (MQ=255)
                                                                                                                       ggCGTCGGCGTGCTCGAAGTCGTCGAAGGTGACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAg                             <  1:2721429/150‑1 (MQ=255)
                                                                                                                               cGTGCTCGAAGTCGTCGAAGGTGACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTAcgccg                     >  1:484852/1‑150 (MQ=255)
                                                                                                                                 tGCTCGAAGTCGTCGAAGGTGACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTACGCCGCt                   <  1:2033693/150‑1 (MQ=255)
                                                                                                                                   cTCGAAGTCGTCGAAGGTGACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTACGCCGCTTc                 <  2:3456951/150‑1 (MQ=255)
                                                                                                                                    tCGAAGTCGTCGAAGGTGACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTACGCCGCTTCg                <  2:3241060/150‑1 (MQ=255)
                                                                                                                                    tCGAAGTCGTCGAAGGTGACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTACGCCGCTTCg                >  2:1109370/1‑150 (MQ=255)
                                                                                                                                        aGTCGTCGTAGGTGTCTTTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCGGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTACGCCGCTTCGTTGc            <  1:1538650/150‑1 (MQ=255)
                                                                                                                                         gtcgtcGAAGGTGACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTACGCCGCTTCGTTGCt           <  1:4420022/150‑1 (MQ=255)
                                                                                                                                         gtcgtcGAAGGTGACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTACGCCGCTTCGTTGCt           >  2:2783412/1‑150 (MQ=255)
                                                                                                                                          tcgtcgAAGGTGACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTACGCCGCTTcacccagc          <  1:4265612/150‑8 (MQ=255)
                                                                                                                                            gtcgAAGGTGACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTACGCCGCTTCGTTGCTGGc        >  1:4744636/1‑150 (MQ=255)
                                                                                                                                                 aGGTGACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTACGCCGCTTCGTTGCTGGCACGTc   >  1:1823669/1‑150 (MQ=255)
                                                                                                                                                  ggTGACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTACGCCGCTTCGTTGCTGGCACGTcc  >  2:4347772/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GACCACGCTCCTTGAGCGGGTTGATGCACACCACCTGGGCGCCACGCTTTACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCAGACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTGACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTACGCCGCTTCGTTGCTGGCACGTCC  >  NC_002947/311588‑311882

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: