Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,294,288 T→G 14.5% L536R (CTG→CGG)  PP_2893 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,294,2880TG14.5% 450.9 / 16.0 187L536R (CTG→CGG) PP_2893hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (50/101);  new base G (30/2);  total (84/103)
Fisher's exact test for biased strand distribution p-value = 9.02e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 198 total aligned reads displayed.

GCTCAAGGCCGAGACCCAGGGGGCAGGGCAAAACTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAA‑GGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCCAGCGCTTGGCCCAGGCCGAAAACCGCCTGCAGCAGTTGGACCTGGATGCCCGCCTGGACA  >  NC_002947/3294139‑3294437
                                                                                                                                                     |                                                                                                                                                      
gCTCAAGGCCGAGACCCAGGGGGCAGGGCAAAACTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCt                                                                                                                                                        <  1:2937348/150‑1 (MQ=255)
 cTCAAGGCCGAGACCCAGGGGGCAGGGCAAAACTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCACGGCCAGCtg                                                                                                                                                       >  1:1285658/1‑150 (MQ=255)
  tCAAGGCCGAGACCCAGGGGGCAGGGCAAAACTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGGCTCCGGCAGCtgt                                                                                                                                                      >  1:16690/1‑150 (MQ=255)
  tCAAGGCCGAGACCCAGGGGGCAGGGCAAAACTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCtgt                                                                                                                                                      >  2:866081/1‑150 (MQ=255)
  tCAAGGCCGAGACCCAGGGGGCAGGGCAAAACTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCtgt                                                                                                                                                      >  1:1740228/1‑150 (MQ=255)
  tCAAGGCCGAGACCCAGGGGGCAGGGCAAAACTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCGCGGCCAGCtgt                                                                                                                                                      >  1:2208327/1‑150 (MQ=255)
   cAAGGCCGAGACCCAGGGGGCAGGGCAAAACTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCtgtg                                                                                                                                                     <  1:1072915/150‑1 (MQ=255)
    aaGGCCGAGACCCAGGGGGCAGGGCAAAACTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTgg                                                                                                                                                    <  1:2618927/150‑1 (MQ=255)
     aGGCCGAGACCCAGGGGGCAGGGCAAAACTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGGCGCGGCCCGCGGTggg                                                                                                                                                   >  1:1810123/1‑149 (MQ=255)
     aGGCCGAGACCCAGGGGGCAGGGCAAAACTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCGGTGGc                                                                                                                                                   >  1:4688073/1‑150 (MQ=255)
        ccGAGACCCAGGGGGCAGGGCAAAACTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGa                                                                                                                                                <  2:2918403/150‑1 (MQ=255)
         cGAGACCCAGGGGGCAGGGCAAAACTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCGCGGCCCGCAGTGGCCgcg                                                                                                                                               >  1:3679342/1‑148 (MQ=255)
          gagaCCCAGGGGGCAGGGCAAAACTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAgg                                                                                                                                              >  2:2985810/1‑150 (MQ=255)
          gagaCCCAGGGGGCAGGGCAAAACTGGACCCTGGGCACACGGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCACAGCTGGCCGAgg                                                                                                                                              <  1:2512137/150‑1 (MQ=255)
           agaCCCAGGGGGCAGGGCAAAACTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCACAGCTGGCCGAGGc                                                                                                                                             <  2:5022695/150‑1 (MQ=255)
               ccAGGGGGCAGGGCAAAACTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACa                                                                                                                                         >  2:3403521/1‑150 (MQ=255)
               ccAGGGGGCAGGGCAAAACTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCAGTGGCCGAGGCTACa                                                                                                                                         >  1:5042897/1‑150 (MQ=255)
                cAGGGGGCAGGGCAAAACTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCACGGCCCGCGGTGGCCGAGGCTACaa                                                                                                                                        >  1:168451/1‑150 (MQ=255)
                 aGGGGGCAGGGCAAAACTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAA‑g                                                                                                                                      >  2:3507950/1‑150 (MQ=255)
                 aGGGGGCAGGGCAAAACTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAA‑g                                                                                                                                      >  2:5407862/1‑150 (MQ=255)
                      gcagggcaAAACTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAA‑GGCCgg                                                                                                                                 <  2:4100938/150‑1 (MQ=255)
                      gcagggcaAAACTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCGGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAA‑GGCCgg                                                                                                                                 <  2:1643021/150‑1 (MQ=255)
                        agggcaAAACTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAA‑GGCCGGGt                                                                                                                               <  1:184438/150‑1 (MQ=255)
                           gcaAAACTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACCA‑GGCCGGGTCaa                                                                                                                            <  1:1276318/150‑1 (MQ=255)
                            caAAACTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAA‑GGCCGGGTCAAg                                                                                                                           <  2:1740228/150‑1 (MQ=255)
                             aaaaCTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAA‑GGCCGGGTCAAgg                                                                                                                          <  1:1828051/150‑1 (MQ=255)
                             aaaaCTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAA‑GGCCGGGTCAAgg                                                                                                                          <  2:431729/150‑1 (MQ=255)
                               aaCTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAA‑GGCCGGGTCAAGGGc                                                                                                                        <  2:4144042/150‑1 (MQ=255)
                                    gACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCCGCGGTGGCCGAGGCGACAA‑GGCCGGGTCAAGGGCCCGcg                                                                                                                   >  1:3689815/1‑149 (MQ=255)
                                    gACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAA‑GGCCGGGTCAAGGGCCGGCt                                                                                                                   <  2:4659452/150‑1 (MQ=255)
                                     aCCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAA‑GGCCGGGTCAAGGGCCGGCTg                                                                                                                  <  2:3846522/150‑1 (MQ=255)
                                        cTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAA‑GGCCGGGTCAAGGGCCGGCTGGAt                                                                                                               >  2:5025792/1‑150 (MQ=255)
                                         tGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCACCGCCAGCTGTGGCCCAGGCTAAAA‑AGCCCGGTCAAAGGCCCGCGGGAtt                                                                                                              >  1:4791669/1‑149 (MQ=255)
                                         tGGGCACACGGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAA‑GGCCGGGTCAAGGGCCGGCttgtag                                                                                                              <  1:5000904/150‑6 (MQ=255)
                                             cacaCTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAA‑GGCCGGGTCAAGGGCCGGCTGGATGTCGc                                                                                                          <  2:4759285/150‑1 (MQ=255)
                                             cacaCTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAA‑GGCCGGGTCAAGGGCCGGCTGGATGTCGc                                                                                                          <  2:5303501/150‑1 (MQ=255)
                                              acacTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAA‑GGCCGGGTCAAGGGCCGGCTGGATGTCGcc                                                                                                         <  1:680541/150‑1 (MQ=255)
                                                       aTCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAA‑GGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTg                                                                                                <  1:59192/150‑1 (MQ=255)
                                                       aTCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCGCGGCCAGCTGTGGCCGAGGCTACAA‑GGCCGGGTCAAGGGCCGGCTGGTTGTCGCCGGCCCCATg                                                                                                >  1:3042107/1‑150 (MQ=255)
                                                         ccAGCTGGGTGACCACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTAGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAC‑GGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCa                                                                                              <  2:4570304/150‑1 (MQ=255)
                                                         ccAGCTGGGTGACAATCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAA‑GGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCa                                                                                              <  2:959613/150‑1 (MQ=255)
                                                             cTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAA‑GGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGcc                                                                                          <  2:1540927/150‑1 (MQ=255)
                                                               gggTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAA‑GGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGcccc                                                                                        <  2:4403590/150‑1 (MQ=255)
                                                                ggTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTACCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAA‑GGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCg                                                                                       >  1:3698776/1‑150 (MQ=255)
                                                                      aaCCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTTACTTGCCGCGCCGCGGCCCGCGGGGGCCGCGGCTACAA‑GGCCGGGGCAAGGGCCGGCGGGATGTCGCCGGCCCCTTGCCGGCCCCCCAAGGt                                                                                 >  1:3518811/1‑150 (MQ=255)
                                                                      aaCCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAA‑GGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGt                                                                                 <  1:833685/150‑1 (MQ=255)
                                                                       aCCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAA‑GGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTa                                                                                <  2:2799301/150‑1 (MQ=255)
                                                                       aCCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCGCCGCCCGCGGTGGCCGAGGCGACAA‑GGCCCGGGCAAGGGCCGGCGGGATGTCGCCGGCCCCCTGCAGGCCCCCCAAGGTa                                                                                >  1:3147326/1‑150 (MQ=255)
                                                                        ccGCATCAACGGCAGAGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTAGGCCACCAGTGGCCGAGGCAACAA‑TGCCGCGTCAAGGGGCGGCTGGATGTCGCCGCCACCATTCAGGCCCCGCCAGGAAc                                                                               >  2:4264028/1‑150 (MQ=255)
                                                                         cGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAA‑GGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACg                                                                              >  1:1716196/1‑150 (MQ=255)
                                                                              cAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAA‑GGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCATGCAGGCCCCGCCAGGTACGCTCAc                                                                         >  1:4913740/1‑150 (MQ=255)
                                                                               aaCGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGGCGCGGCCCGCGGTGGGCGGGGCTACAA‑ACCCCGGGCAACGGGCGGCGGGGTGTGGCCCGCCCCGTGCAGGCCCCCCCAGGTCCGCGCAcc                                                                        >  1:1976653/1‑150 (MQ=255)
                                                                                 cGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCCCGGCCCCCGGTGGCCGCGGCGACAA‑CGCCGGGGCCAGGGCCGGCTGGGGGTCGCCCGCCCCCTGCCGGCCCCCCAACGGACGCGCAcccc                                                                      >  1:1228987/1‑149 (MQ=255)
                                                                                  ggCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAA‑GGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTg                                                                     >  2:3777112/1‑150 (MQ=255)
                                                                                   gCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAA‑GGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGc                                                                    <  1:561751/150‑1 (MQ=255)
                                                                                        ggCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCGCGGCCCGCGGTGGGCGCGGCGACAA‑GGCCGGGGCAAGGGCCCGCGGGGTGTCGCCGGCCCCAAGCAGGCCCCCCCACGTACGCTCACCCAGCCGGGc                                                               >  1:495877/1‑150 (MQ=255)
                                                                                          cAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGTTACAA‑GGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCCa                                                             >  1:3206226/1‑150 (MQ=255)
                                                                                            gCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCGGTGGCCGAGGCGACAA‑GGCCGGGTCAAGGGCCGGCGGGGTGTCGCCGGCCCCTTGCAGGCCCCGCCAGGTACGCTCACCCAGCAGGGCCAgc                                                           >  2:2369582/1‑150 (MQ=255)
                                                                                               tGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAA‑GGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCCAGCGCt                                                        <  2:5286229/150‑1 (MQ=255)
                                                                                                 cAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAA‑GGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCCAGCGCTTg                                                      <  1:4107858/150‑1 (MQ=255)
                                                                                                  aaCAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGGCGTGGGTAAAA‑TGCCGGGTGAAGGGGCGGCTGGTTTTCGCCGGGACCTTGCAGGCCCCGCAATGTAGGCTCACCCTGCAGGGCCAGCGCTTgg                                                     >  1:3638586/1‑150 (MQ=255)
                                                                                                  aaCAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAA‑GGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCCAGCGCTTgg                                                     <  2:2563797/150‑1 (MQ=255)
                                                                                                  aaCAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAA‑GGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCCAGCGCTTgg                                                     <  1:5407862/150‑1 (MQ=255)
                                                                                                   aCAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCGCGGCCAGCTGTGGCCGCGGCTACAA‑GGCCGGGTCAAGGGCCGGCTGGATTTCGCCGGCCCCTTGCAGGCCCCCCCCGGTACCCTCACCCAGCCGGGGCAGCGCTTGGc                                                    >  2:5168105/1‑150 (MQ=255)
                                                                                                   aCAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCGCGGCCAGCGGTGGGCCGGGGGACAACGGCCGGG‑GAAGGGGCCGCGGGATGTGGCCCGCCCCATGCAGGCCCCCCAACGGACGCGCACCCAGCAGGGCCCGCGCGTGGc                                                    >  1:2085992/1‑150 (MQ=255)
                                                                                                      gcgcCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAA‑GGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCCAGCGCTTGGCCCa                                                 <  1:2086606/150‑1 (MQ=255)
                                                                                                       cgcCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTCCAA‑GGCCGGGCCAAGGGCCGGCTGGCTGTCGCCGGCACCTTGCAGGCCCCGCCAGGTACGCGCACCCTGCAGGGCCCGCGCGTGGccccg                                                >  2:3287588/1‑148 (MQ=255)
                                                                                                        gcCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAA‑GGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCCAGCGCTTGGCCCAgg                                               <  1:576523/150‑1 (MQ=255)
                                                                                                        gcCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAA‑GGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCCAGCGCTTGGCCCAgg                                               <  1:5349983/150‑1 (MQ=255)
                                                                                                         ccTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAA‑GGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCCAGCGCTTGGCCCAGGc                                              <  2:1329476/150‑1 (MQ=255)
                                                                                                          ctggctggGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAA‑GGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCCAGCGCTTGGCCCAGGcc                                             <  2:2373125/150‑1 (MQ=255)
                                                                                                                gggCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAA‑GGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCCAGCGCTTGGCCCAGGCCGAAAAc                                       <  2:1241447/150‑1 (MQ=255)
                                                                                                                 ggCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAA‑GGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCCAGCGCTTGGCCCAGGCCGAAAAcc                                      <  2:3157598/150‑1 (MQ=255)
                                                                                                                  gCGTATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGGGGCCGAGGCTACAA‑GGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCCAGCGCTTGGCCCAGGCCGAAAACCg                                     <  2:3993841/150‑1 (MQ=255)
                                                                                                                     aaTCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAA‑GGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCCAGCGCTTGGCCCAGGCCGAAAACCGCCt                                  <  1:2716307/150‑1 (MQ=255)
                                                                                                                      aTCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAA‑GGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCCAGCGCTTGGCCCAGGCCGAAAACCGCCTg                                 <  1:1136201/150‑1 (MQ=255)
                                                                                                                       tCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAA‑GGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCCAGCGCtt                                                       >  2:2686323/1‑127 (MQ=255)
                                                                                                                       tCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAA‑GGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCCAGCGCtt                                                       <  1:2686323/127‑1 (MQ=255)
                                                                                                                       tCGACCTTGACTTGCCGCGCCGCGGCCCGCTGTGGCCGAGGCTACAA‑GGCCGGGTCAAGGGCCGGCGGGATGTCGCCGGCACCTTGCAGGCCCCCCAAGGTACGCTCACCCTGCAGGGCCAGCGCTTGGCCCCGGCCGAAAACCCCCGgc                                >  2:3213689/1‑150 (MQ=255)
                                                                                                                        cGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAA‑GGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCCAGCGCTTGGCCCAGGCCGAAAACCGCCTgca                               <  1:4783879/150‑1 (MQ=255)
                                                                                                                         gACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAA‑GGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCCAGCGCTTGGCCCAGGCCGAAAACCGCCTgcag                              <  1:5155624/150‑1 (MQ=255)
                                                                                                                           ccTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAA‑GGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCCAGCGCTTGGCCCAGGCCGAAAACCGCCTgcagca                            <  1:1530388/150‑1 (MQ=255)
                                                                                                                           ccTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAA‑GGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCCAGCGCTTGGCCCAGGCCGAAAACCGCCTgcagca                            <  1:1254726/150‑1 (MQ=255)
                                                                                                                           ccTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAA‑GGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCCAGCGCTTGGCCCAGGCCGAAAACCGCCTgcagca                            <  2:1716196/150‑1 (MQ=255)
                                                                                                                           ccTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAA‑GGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCCAGCGCTTGGCCCAGGCCGAAAACCGCCTgcagca                            <  2:2208327/150‑1 (MQ=255)
                                                                                                                             ttGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAA‑GGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCCAGCGCTTGGCCCAGGCCGAAAACCGCCTGCAGCAGt                          <  2:5228763/150‑1 (MQ=255)
                                                                                                                                 cTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAA‑GGCCGGGTCAAGGGCCGGCGGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCCAGCGCTTGGCCCAGGCCGAAAACCGCCTGCAGCAGTTggc                      <  1:3717815/150‑2 (MQ=255)
                                                                                                                                      cgcgCCTCGGCCAGCTGTGGCCGAGGCTACAA‑GGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCCAGCGCTTGGCCCAGGCCGAAAACCGCCTGCAGCAGTTGGACCTgg                 <  2:936372/150‑1 (MQ=255)
                                                                                                                                      cgcgCCTCGGCCAGCTGTGGCCGAGGCTACAA‑GGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCCAGCGCTTGGCCCAGGCCGAAAACCGCCTGCAGCAGTTGGACCTgg                 <  1:4330950/150‑1 (MQ=255)
                                                                                                                                      cgcgCCTCGGCCAGCTGTGGCCGAGGCTACAA‑GGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCCAGCGCTTGGCCCAGGCCGAAAACCGCCTGCAGCAGTTGGACCTgg                 <  1:3854210/150‑1 (MQ=255)
                                                                                                                                       gcgcCTCGGCCAGCTGTGGCCGAGGCTACAA‑GGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCCAGCGCTTGGCCCAGGCCGAAAACCGCCTGCAGCAGTTGGACCTGGa                <  1:502345/150‑1 (MQ=255)
                                                                                                                                       gcgcCTCGGCCAGCTGTGGCCGAGGCTACAA‑GGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCCAGCGCTTGGCCCAGGCCGAAAACCGCCTGCAGCAGTTGGACCTGGa                <  2:203943/150‑1 (MQ=255)
                                                                                                                                          ccTCGGCCAGCTGTGGCCGAGGCTACAA‑GGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCCAGCGCTTGGCCCAGGCCGAAAACCGCCTGCAGCAGTTGGACCTGGATGc             >  1:3885280/1‑150 (MQ=255)
                                                                                                                                               gCCAGCTGTGGCCGAGGCTACAA‑GGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCCAGCGCTTGGCCCAGGCCGAAAACCGCCTGCAGCAGTTGGACCTGGATGCCCGcc        >  2:386311/1‑150 (MQ=255)
                                                                                                                                                ccAGCTGTGGCCGAGGCTACAA‑GGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCCAGCGCTTGGCCCAGGCCGAAAACCGCCTGCAGCAGTTGGACCTGGATGCCCGCCt       <  2:4003491/150‑1 (MQ=255)
                                                                                                                                                ccAGCTGTGGCCGAGGCTACAA‑GGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTGGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCCAGCGCTTGGCCCAGGCCGAAAACCGCCTGCAGCAGTTGGACCTGGATGCCCGCCt       <  1:5238527/150‑1 (MQ=255)
                                                                                                                                                 cAGCTGTGGCCGAGGCTACAA‑GGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCCAGCGCTTGGCCCAGGCCGAAAACCGCCTGCAGCAGTTGGACCTGGATGCCCGCCTg      <  1:3275589/150‑1 (MQ=255)
                                                                                                                                                 cAGCTGTGGCCGAGGCTACAA‑GGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCCAGCGCTTGGCCCAGGCCGAAAACCGCCTGCAGCAGTTGGACCTGGATGCCCGCCTg      <  2:4791669/150‑1 (MQ=255)
                                                                                                                                                 cAGCTGTGGCCGAGGCTACAA‑GGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCCAGCGCTTGGCCCAGGCCGAAAACCGCCTGCAGCAGTTGGACCTGGATGCCCGCCTg      <  2:568114/150‑1 (MQ=255)
                                                                                                                                                     tgtgGCCGAGGCTACAA‑GGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCCAGCGCTTGGCCCAGGCCGAACACCGCCTGCAGCAGTTGGACCTGGATGCCCGCCTGGACa  >  2:3008998/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
GCTCAAGGCCGAGACCCAGGGGGCAGGGCAAAACTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAA‑GGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCCAGCGCTTGGCCCAGGCCGAAAACCGCCTGCAGCAGTTGGACCTGGATGCCCGCCTGGACA  >  NC_002947/3294139‑3294437

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑