Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,579,942 G→C 5.2% S71R (AGC→AGG glgB ← 1,4‑alpha‑glucan branching enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,579,9420GC 5.2% 303.3 / 6.8 97S71R (AGC→AGGglgB1,4‑alpha‑glucan branching enzyme
Reads supporting (aligned to +/- strand):  ref base G (49/42);  new base C (2/3);  total (51/46)
Fisher's exact test for biased strand distribution p-value = 6.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.75e-01

TGCCTTCCGCGAACAAGTGAAGGTCCATGTCTCCCAATTGCGGGCCGAAGCTGTAGGGGTCTTCGGTAACTTGCTCGCCACCGGCCCAGACGATGTGAAGCTGATAGGGCTGTGCCTTATCGAGGTGCGCGCTGAACAACCCGGGCAGGCTGCCCTGCTCCATTTCGGCAAGCACCCGCCCGTCATCCCGGGCCTGAATGCGCACTTTCAGGGCATTGGGCAGAAAGGCGCGGACCACCTGCCCGCCTGCGCCGTCGCCGTGGGGGCCGAGTACCGCAAAAGGATCGGCGTGTTCG  >  NC_002947/4579794‑4580089
                                                                                                                                                    |                                                                                                                                                   
tGCCTTCCGCGAACAAGTGAAGGTCCATGTCTCCCAATTGCGGGCCGAAGCTGTAGGGGTCTTCGGTAACTTGCTCGCCACCGGCCCAGACGATGTGAAGCTGATAGGGCTGTGCCTTATCGAGGTGCGCGCTGAACAACCCGGGCAGGc                                                                                                                                                    >  1:601883/1‑150 (MQ=255)
        gcgAACAAGTGAAGGTCCATGTCTCCCAATTGCGGGCCGAAGCTGTAGGGGTCTTCGGTAACTTGCTCGCCACCGGCCCAGACGATGTGAAGCTGATAGGGCTGTGCCTTATCGAGGTGCGCGCTGAACAACCCGGGCAGGCTGCCCTGc                                                                                                                                            <  1:2212917/150‑1 (MQ=255)
        gcgAACAAGTGAAGGTCCATGTCTCCCAATTGCGGGCCGAAGCTGTAGGGGTCTTCGGTAACTTGCTCGCCACCGGCCCAGACGATGTGAAGCTGATAGGGCTGTGCCTTATCGAGGTGCGCGCTGAACAACCCGGGCAGGCTGCCCTGc                                                                                                                                            <  2:1943246/150‑1 (MQ=255)
        gcgAACAAGAGTAGGTCCATGTCTCCCAATTGCGGGCCGAAGCTGTAGGGGTCTTCGGTACGTAGCTCACGACCGGCCATAACGATGTGAAGCTGCTAGGGCTGTGCCTTATCGAGGTGCGCGCTGAACAACCCGGGCAGGCTGCCCTGc                                                                                                                                            <  1:1355597/150‑1 (MQ=255)
         cgAACAAGTGAAGGTCCATGTCTCCCAATTGCGGGCCGAAGCTGTAGGGGTCTTCGGTAACTTGCTCGCCACCGGCCCAGACGATGTGAAGCTGATAGGGCTGTGCCTTATCGAGGTGCGCGCTGAACAACCCGGGCAGCCTGcccgggt                                                                                                                                           <  2:212107/150‑5 (MQ=255)
           aaCAAGTGAAGGTCCATGTCTCCCACTTGCGGGCCGAAGCTGTAGGGGTCTTCGGTAACTTGCTCGCCACCGGCCCAGACGATGTGAAGCTGATAGGGCTGTGCCTTATCGAGGTGCGCGCTGAACAACCCGGGCAGGCTGCCCTGCTcc                                                                                                                                         <  2:555975/150‑1 (MQ=255)
              aaGTGAAGGTCCATGTCTCCCAATTGCGGGCCGAAGCTGTAGGGGTCTTCGGTAACTTGCTCGCCACCGGCCCAGACGATGTGAAGCTGATAGGGCTGTGCCTTATCGAGGTGCGCGCTGAACAACCCGGGCAGGCTGCCCTGCTCCAtt                                                                                                                                      >  2:721753/1‑150 (MQ=255)
                gTGAAGGTCCATGTCTCCCAATTGCGGGCCGAAGCTGTAGGGGTCTTCGGTAACTTGCTCGCCACCGGCCCAGACGATGTGAAGCTGATAGGGCTGTGCCTTATCGAGGTGCGCGCTGAACAACCCGGGCAGGCTGCCCTGCTCCATTTc                                                                                                                                    >  1:2319766/1‑150 (MQ=255)
                gTGAAGGTCCATGTCTCCCAATTGCGGGCCGAAGCTGTAGGGGTCTTCGGTAACTTGCTCGCCACCGGCCCAGACGATGTGAAGCTGATAGGGCTGTGCCTTATCGAGGTGCGCGCTGAACAACCCGGGCAGGCTGCCCTGCTCCATTTc                                                                                                                                    <  2:1865730/150‑1 (MQ=255)
                     ggTCCATGTCTCCCAATTGCGGGCCGAAGCTGTAGGGGTCTTCGGTAACTTGCTCGCCACCGGCCCAGACGATGTGAAGCTGATAGGGCTGTGCCTTATCGAGGTGCGCGCTGAACAACCCGGGCAGGCTGCCCTGCTCCATTTCGGCaa                                                                                                                               <  1:1721761/150‑1 (MQ=255)
                      gTCCATGTCTCCCAATTGCGGGCCGAAGCTGTAGGGGTCTTCGGTAACTTGCTCGCCACCGGCCCAGACGATGTGAAGCTGATAGGGCTGTGCCTTATCGAGGTGCGCGCTGAACAACCCGGGCAGGCTGCCCTGCTCCATTTCGGCAAg                                                                                                                              >  1:1169774/1‑150 (MQ=255)
                       tCCATGTCTCCCAATTGCGGGCCGAAGCTGTAGGGGTCTTCGGTAACTTGCTCGCCACCGGCCCAGACGATGTGAAGCTGATAGGGCTGTGCCTTATCGAGGTGCGCGCTGAACAACCCGGGCAGGCTGCCCTGCTCCATTTCGGCAAGc                                                                                                                             >  2:807973/1‑150 (MQ=255)
                            gTCTCCCAATTGCGGGCCGAAGCTGTAGGGGTCTTCGGTAACTTGCTCGCCACCGGCCCAGACGATGTGGATCTGATAGGGCCGTCCCTTATCGAGGTGCGCGCTGAACAACCCGGGCAGGCTGCCCTGCTCCATTTCGGCAAGCAcccg                                                                                                                        >  1:166712/1‑150 (MQ=255)
                               tcCCAATTGCGGGCCGAAGCTGTAGGGGTCTTCGGTAACTTGCTCGCCACCGGCCCAGACGATGTGAAGCTGATAGGGCTGTGCCTTATCGAGGTGCGCGCTGAACAACCCGGGCAGGCTg                                                                                                                                                  >  1:1233051/1‑121 (MQ=255)
                               tcCCAATTGCGGGCCGAAGCTGTAGGGGTCTTCGGTAACTTGCTCGCCACCGGCCCAGACGATGTGAAGCTGATAGGGCTGTGCCTTATCGAGGTGCGCGCTGAACAACCCGGGCAGGCTGCCCTGCTCCATTTCGGCAAGCAcccgccc                                                                                                                     >  1:1278288/1‑150 (MQ=255)
                               tcCCAATTGCGGGCCGAAGCTGTAGGGGTCTTCGGTAACTTGCTCGCCACCGGCCCAGACGATGTGAAGCTGATAGGGCTGTGCCTTATCGAGGTGCGCGCTGAACAACCCCGGCAGGCTGCCCTGCTCCATTTCGGCAAGCAcccgccc                                                                                                                     >  1:2060150/1‑150 (MQ=255)
                               tcCCAATTGCGGGCCGAAGCTGTAGGGGTCTTCGGTAACTTGCTCGCCACCGGCCCAGACGATGTGAAGCTGACAGGGCTGTGCCTTATCGAGGTGCGCGCTGAACAACCCGGGCAGGCTg                                                                                                                                                  <  2:1233051/121‑1 (MQ=255)
                                   aaTTGCGGGCCGAAGCTGTAGGGGTCTTCGGTAACTTGCTCGCCACCGGCCCAGACGATGTGAAGCTGATAGGGCTGTGCCTTATCGAGGTGCGCGCTGAACAACCCGGGCAGGCTGCCCTGCTCCATTTCGGCAAGCACCCGCCCGTCa                                                                                                                 >  1:1783785/1‑150 (MQ=255)
                                   aaTTGCGGGCCGAAGCTGTAGGGGTCTTCGGTAACTTGCTCGCCACCGGCCCAGACGATGTGAAGCTGATAGGGCTGTGCCTTATCGAGGTGCGCGCTGAACAACCCGGGCAGGCTGCCCTGCTCCATTTCGGCAAGCACCCGCCCGTCa                                                                                                                 >  1:1783775/1‑150 (MQ=255)
                                   aaTTGCGGGCCGAAGCTGTAGGGGTCTTCGGTAACTTGCTCGCCACCGGCCCAGACGATGTGAAGCTGATAGGGCTGTGCCTTATCGAGGGGCGCGCTGAACAACCCGGGCAGGCTGCCCTGCTCCATTTCGGCAAGCACCCGCCCGTCa                                                                                                                 >  1:616770/1‑150 (MQ=255)
                                          ggCCGAAGCGGTAGGGGTCTTCGGTAACTTGCTCGCCACCGGCCCTGACGATGTGAAGCTGATAGGGCTGTGCCTTATCGAGGTGCGCGCTGAACAACCCGGTCAGGCTGCCCTGCTCCATTTCGGCAAGCACCCGCCCGTCATCCCggg                                                                                                          <  2:321251/150‑1 (MQ=255)
                                           gCCGAAGCTGTAGGGGTCTTCGGTAACTTGCTCGCCACCGGCCCAGACGATGTGAAGCTGATAGGGCTGTGCCTTATCGAGGTGCGCGCTGAACAACCCGGGNAGGCTGNCCTGCTCCATTTCGGCAAGCACCCGCCCGTCATGCCGNGc                                                                                                         <  1:2184680/150‑1 (MQ=255)
                                           gCCGAAGCTGTAGGGGTCTTCGGTAACTTGCTCGCCACCGGCCCAGACGATGTGAAGCTGATAGGGCTGTGCCTTATCGAGGTGCGCGCTGAACAACCCGGGCAGGCTGCCCTGCTCCATTTCGGCAAGCACCCGCCCGTCATCCCGGGc                                                                                                         >  2:1413212/1‑150 (MQ=255)
                                           gCCGAAGCTGTAGCGGTCTTCGGTAACTTGCTCGCCACCGGCCCAGACGATGTGAAGCTGATAGGGCCGTGCCTTATCGAGGTGCGCGCTGAACAACCCGGGCAGGCTGCCCTGCTCCATTTCGGCAAGCACCCGCCCGTCATCCCGGGc                                                                                                         <  2:1295824/150‑1 (MQ=255)
                                                aGCTGTAGGGGTCTTCGGTAACTTGCTCGCCACCGGCCCAGACGATGTGAAGCTGATAGGGCTGTGCCTTATCGAGGTGCGCGCTGAACAACCCGGGCAGGCTGCCCTGCTCCATTTCGGCAAGCACCCGCCCGTCATCCCGGGCCTGaa                                                                                                    <  1:1754984/150‑1 (MQ=255)
                                                   tGTAGGGGTCTTCGGTAACTTGCTCGCCACCGGGCCAGACGATGTGAAGCTGATAGGGCTGTGCCTTATCGAGGTGCGCGCTGAACAACCCGGGCAGGCTGCCCTGCTCCATTTCGGCAAGCACCCGCCCGTCATCCCGGGCCTGAATgc                                                                                                 >  2:79144/1‑150 (MQ=255)
                                                    gTAGGGGTCTTCGGTAACTTGCTCGCCACCGGCCCAGACGATGTGAAGCTGATAGGGCTGTGCCTTATCGAGGTGCGCGCTGAACAACCCGGGCAGGCTGCCCTGCTCCATTTCGGCAAGCACCCGCCCGTCATCCCGGGCCTGAATgcg                                                                                                <  2:1250687/150‑1 (MQ=255)
                                                    gTAGGGGTCTTCGGTAACTTGCTCGCCACCGGCCCAGACGATGTGAAGCTGATAGGGCTGTGCCTTATCGAGGTGCGCGCTGAACAACCCGGGCAGGCTGCCCTGCTACATTTCGGCAAGCACCCGCCCGTCATCCCGGGCCTGAATgcg                                                                                                >  1:1690398/1‑150 (MQ=255)
                                                     tAGGGGTCTTCGGTAACTTGCTCGCCACCGGCCCAGACGATGTGAAGCTGATAGGGCTGTGCCTTATCGAGGTGCGCGCTGAACAACCCGGGCAGGCTGCCCTGCTCCATTTCGGCAAGCACCCGCCCGTCATCCCGGGCCTGAATgcgc                                                                                               <  2:1278288/150‑1 (MQ=255)
                                                     tAGGGGTCTTCGGTAACTTGCTCGCCACCGGCCCAGACGATGTGAAGCTGATAGGGCTGTGCCTTATCGAGGTGCGCGCTGAACAACCCGGGCAGGCTGCCCTGCTCCATTTCGGCAAGCACCCGCCCGTCATCCCGGGCCTGAATgcgc                                                                                               >  2:288253/1‑150 (MQ=255)
                                                        gggTCTTCGGTAACTTGCTCGCCACCGGCCCAGACGATGTGAAGCTGATAGGGCTGTGCCTTATCGAGGTGCGCGCTGAACAACCCGGGCAGGCTGCCCTGCTCCATTTCGGCAAGCACCCGCCCGTCATCCCGGGCCTGAATGCGCACt                                                                                            <  2:2303042/150‑1 (MQ=255)
                                                        gggTCTTCGGTAACTTGCTCGCCACCGGCCCAGACGATGTGAAGCTGATAGGGCTGTGCCTTATCGAGGTGCGCGCTGAACAACCCGGGCAGGCTGCCCTGCTCCATTTCGGCAAGCACCCGCCCGTCATCCCGGGCCTGAATGCGCACt                                                                                            <  1:1196815/150‑1 (MQ=255)
                                                        gggTCTTCGGTAACTTGCTCGCCACCGGCCCAGACGATGTGAAGCTGATAGGGCTGTGCCTTATCGAGGTGCGCGCTGAACAACCCGGGCAGGCTGCCCTGCTCCATTTCGGCAAGCACCCGCCCGTCATCCCGGGCCTGAATGCGCACt                                                                                            >  1:434445/1‑150 (MQ=255)
                                                        gggTCTTCGGTAACTTGCTCGCCACCGGCCCAGACGATGTGAAGCTGATAGGGCTGTGCCTTATCGAGGTGCGCGCTGAACAACCCGGGCAGCCTGcccggg                                                                                                                                            <  1:2116958/102‑4 (MQ=255)
                                                        gggTCTTCGGTAACTTGCTCGCCACCGGCCCAGACGATGTGAAGCTGATAGGGCTGTGCCTTATCGAGGTGCGCGCTGAACAACCCGGGCAGCCTGcccggg                                                                                                                                            >  2:2116958/1‑99 (MQ=255)
                                                          gTCTTCGGTAACTTGCTCGCCACCGGCCCAGACGATGTGAAGCTGATAGGGCTGTGCCTTATCGAGGTGCGCGCTGAACAACCCGGGCAGGCTGCCCTGCTCCATTTCGGCAAGCACCCGCCCGTCATCCCGGGCCTGAATGCGCACttt                                                                                          >  1:949193/1‑150 (MQ=255)
                                                          gTCTTCGGTAACTTGCTCGCCACCGGCCCAGACGATGTGAAGCTGATAGGGCTGTGCCTTATCGAGGTGCGCGCTGAACAACCCGGGCAGCCTGcccgggt                                                                                                                                           <  1:1823155/101‑5 (MQ=255)
                                                          gTCTTCGGTAACTTGCTCGCCACCGGCCCAGACGATGTGAAGCTGATAGGGCTGTGCCTTATCGAGGTGCGCGCTGAACAACCCGGGCAGCCTGcccgggt                                                                                                                                           >  2:1823155/1‑97 (MQ=255)
                                                            cTTCGGTAACTTGCTCGCCACCGGCCCAGACGATGTGAAGCTGATAGGGCTGTGCCTTATCGAGGTGCGCGCTGAACCAACCGGGCAGGCTGCCCTGCTCCATTTCGGCAAGCACCCGCCCGTCATCCCGGGCCTGAATGCGCACTTTCa                                                                                        <  1:1353444/150‑1 (MQ=255)
                                                            cTTCGGTAACTTGCTCGCCACCGGCCCAGACGATGTGAAGCTGATAGGGCTGTGCCTTATCGAGGTGCGCGCTGAACAACCCGGGCAGGCTGCCCTGCTCCATTTCGGCAAGCACCCGCCCGTCATCCCGGGCCTGAATGCGCACTTTCa                                                                                        <  2:2337895/150‑1 (MQ=255)
                                                            cTTCGGTAACTTGCTCGCCACCGGCCCAGACGATGTGAAGCTGATAGGGCTGTGCCTTATCGAGGTGCGCGCTGAACAACCCGGGCAGGCTGCCCTGCTCCATTTCGGCAAGCACCCGCCCGTCATCCCGGGCCTGAATGCGCACTTTCa                                                                                        >  1:548334/1‑150 (MQ=255)
                                                               cGGTAACTTGCTCGCCACCGGCCCAGACGATGTGAAGCTGATAGGGCTGTGCCTTATCGAGGTGCGCGCTGAACAACCCGGGCAGGCTGCCCTGCTTCTTGTCGGCAAGCACCCGCCCGTCATCCCGGGCCTGAATGCGCACTTTCAggg                                                                                     <  1:1995634/150‑1 (MQ=255)
                                                               cGGTAACTTGCTCGCCACCGGCCCAGACGATGTGAAGCTGATAGGGCTGTGCCTTATCGAGGTGCGCGCTGAACAACCCGGGCAGGCTGCCCTGCTCCATTTCGGCAAGCACCCGCCCGTCATCCCGGGCCTGAATGCGCACTTTCAggg                                                                                     >  1:569564/1‑150 (MQ=255)
                                                                    aCTTGCTCGCCACCGGCCCAGACGATGTGAAGCTGATAGGGCTGTGCCTTATCGAGGTGCGCGCTGAACAACCCGGGCAGGCTGCCCTGCTCCATTTCGGCAAGCACCCGCCCGTCATCCCGGGCCTGAATGCGCACTTTCAGGGCATTg                                                                                >  1:1363974/1‑150 (MQ=255)
                                                                      ttGCTCGCCACCGGCCCAGACGATGTGAAGCTGATAGGGCTGTGCCTTATCGAGGTGCGCGCTGAACAACCCGGGCAGGCTGCCCTGCTCCATTTCGGCAAGCACCCGCCCGTCATCCCGGGCCTGAATGCGCACTTTCAg                                                                                       >  2:2316092/1‑141 (MQ=255)
                                                                      ttGCTCGCCACCGGCCCAGACGATGTGAAGCTGATAGGGCTGTGCCTTATCGAGGTGCGCGCTGAACAACCCGGGCAGGCTGCCCTGCTCCATTTCGGCAAGCACCCGCCCGTCATCCCGGGCCTGAATGCGCACTTTCAg                                                                                       <  1:2316092/141‑1 (MQ=255)
                                                                      ttGCTCGCCACCGGCCCAGACGATGTGAAGCTGATAGGGCTGTGCCTTATCGAGGTGCGCGCTGAACAACCCGGGCAGGCTGCCCTGCTCCATTTCGGCAAGCACCCGCCCGTCATCCCGGGCCTGAATGCGCACTTTCAGGGCATTggg                                                                              >  2:1091236/1‑150 (MQ=255)
                                                                          tCGCCACCGGCCCAGACGATGTGAAGCTGATAGGGCTGTGCCTTATCGAGGTGCGCGCTGAACAACCCGGGCAGGCTGCCCTGCTCCATTTCGGCAAGCACCCGCCCGTCATCCCGGGCCTGAATGCGCACTTTCAGGGCATTGGGCAgg                                                                          >  1:1472917/1‑149 (MQ=255)
                                                                          tCGCCACCGGCCCAGACGATGTGAAGCTGATAGGGCTGTGCCTTATCGAGGTGCGCGCTGAACAACCCGGGCAGGCTGCCCTGCTCCATTTCGGCAAGCACCCGCCCGTCATCCCGGGCCTGAATGCGCACTTTCAGGGCATTGGGCAGa                                                                          >  1:1129674/1‑150 (MQ=255)
                                                                          tCGCCACCGGCCCAGACGATGTGAAGCTGATAGGGCTGTGCCTTATCGAGGTGCGCGCTGAACAACCCGGGCAGGCTGCCCTGCTCCATTTCGGCAAGCACCCGCCCGTCATCCCGGGCCTGAATGCGCACTTTCAGGGCATTGGGCAGa                                                                          >  1:1473428/1‑150 (MQ=255)
                                                                            gCCACCGGCCCAGACGATGTGAAGCTGATAGGGCTGTGCCTTATCGAGGTGCGCGCTGAACAACCCGGGCAGGCTGCCCTGCTCCATTTCGGCaa                                                                                                                               <  2:2348303/95‑1 (MQ=255)
                                                                            gCCACCGGCCCAGACGATGTGAAGCTGATAGGGCTGTGCCTTATCGAGGTGCGCGCTGAACAACCCGGGCAGGCTGCCCTGCTCCATTTCGGCaa                                                                                                                               >  1:2348303/1‑95 (MQ=255)
                                                                               aCCGGCCCAGACGATGTGAAGCTGATAGGGCTGTGCCTTATCGAGGTGCGCGCTGAACAACCCGGGCAGGCTGCCCTGCTCCATTTCGGCAAGCACCCGCCCGTCATCCCGGGCCTGAATGCGCACTTTCAGGGCATTGGGCAGAAAGgc                                                                     >  2:1581379/1‑150 (MQ=255)
                                                                                 cGGCCCAGACGATGTGAAGCTGATAGGGCTGTGCCTTATCGAGGTGCGCGCTGAACAACCCGGGCAGGCTGCCCTGCTCCATTTCGGCAAGCACCCGCCCGTCATCCCGGGCCTGAATGCGCACTTTCAGGGCATTGGGCAGAAAGgcgc                                                                   >  1:783863/1‑150 (MQ=255)
                                                                                  ggCCCAGACGATGTGAAGCTGATAGGGCTGTGCCTTATCGAGGTGCGCGCTGAACAACCCGGGCAGGCTGCCCTGCTCAATTTCGGCAAGCACCCGCCCGTCATCCCGGGCCTGAATGCGCACTTTCAGGGCATTGGGCAGAAAGGcgcg                                                                  <  1:1993622/150‑1 (MQ=255)
                                                                                  ggCCCACACGATGGGAAGCTGATAGGGCTGTGCCTTATCGAGGTGCGCGCTGAACAACCCGGGCAGGCTGCCCTGCTCCATTTCGGCAAGCACCCGCCCGTCATCCCGGGCCTGAATGCGCACTTTCAGGGCATTGGGCAGAAAGGcgcg                                                                  <  2:60003/150‑1 (MQ=255)
                                                                                   gCCCAGACGATGTGAAGCTGATAGGGCTGTGCCTTATCGAGGTGCGCGCTGAACAACCCGGGCAGGCTGCCCTGCTCCATTTCGGCAAGCACCCGCCCGTCATCCCGGGCCTGAATGCGCACTTTCAGGGCATTGGGCAGAAAGGCGCgg                                                                 >  1:662909/1‑150 (MQ=255)
                                                                                    cccAGACGATGTGAAGCTGATAGGGCTGTGCCTTATCGAGGTGCGCGCTGAACAACCCGGGCAGGCTGCCCTGCTCCATTTCGGCAAGCACCCGCCCGTCATCCCGGGCCTGAATGCGCACTTTCAGGGCATTGGGCAGAAAGGCGCGGa                                                                <  2:1169774/150‑1 (MQ=255)
                                                                                     ccAGACGATGTGAAGCTGATAGGGCTGTGCCTTATCGAGGTGCGCGCTGAACAACCCGGGCAGGCTGCCCTGCTCCATTTCGGCAAGCACCCGCCCGTCATCCCGGGCCTGAAt                                                                                                   >  1:1000999/1‑114 (MQ=255)
                                                                                     ccAGACGATGTGAAGCTGATAGGGCTGTGCCTTATCGAGGTGCGCGCTGAACAACCCGGGCAGGCTGCCCTGCTCCATTTCGGCAAGCACCCGCCCGTCATCCCGGGCCTGAAt                                                                                                   <  2:1000999/114‑1 (MQ=255)
                                                                                      cAGACGATGTGAAGCTGATAGGGCTGTGCCTTATCGAGGTGCGCGCTGAACAACCCGGGCAGGCTGCCCTGCTCCATTTCGGCAAGCACCCGCCCGTCATCCCGGGCCTGAATGCGCACTTTCAGGGCATTGGGCAGAAAGGCGCGGacc                                                              >  2:1585927/1‑150 (MQ=255)
                                                                                        gACGATGTGAAGCTGATAGGGCTGTGCCTTATCGAGGTGCGCGCTGAACAACCCGGGCAGGCTGCCCTGCTCCATTTCGGCAAGCACCCGCCCGTCATCCCGGGCCTGAATGCGCACTTTCAGGGCATTGGGCAGAAAGGCGCGGaccac                                                            >  2:772981/1‑150 (MQ=255)
                                                                                        gACGATGTGAAGCTGATAGGGCTGTGCCTTATCGAGGTGCGCGCTGAACAACCCGGGCAGGCTGCCCTGCTCCATTTCGGCAAGCACCCGCCCGTCATCCCGGGCCTGAATGCGCACTTTCAGGGCATTGGGCAGAAAGGCGCGGaccac                                                            <  2:783863/150‑1 (MQ=255)
                                                                                              gtgAAGCTGATAGGGCTGTGCCTTATCGAGGTGCGCGCTGAACAACCCGGGCCGGGCGCTGTCGTCCGATTGTTTAGGCGCCCGCCCGTCATCCCGGGCCTGAATGCGCACTTTCAGGGCATTGGGCAGAAAGGCGCGGACCACCTGccc                                                      <  1:1750578/150‑1 (MQ=255)
                                                                                              gtgAAGCTGATAGGGCTGTGCCTTATCGAGGTGCGCGCTGAACAACCCGGGAAGAATCGCGTTCGATGTGACTATGATCGCCCGCCCTTCATCCCGGGCCTGAATGCGCACTTTCAGGGCATTGGGCAGAAAGGCGCGGACCACCTGccc                                                      <  1:1750624/150‑1 (MQ=255)
                                                                                                gAAGCTGATAGGGCTGTGCCTTATCGAGGTGCGCGCTGAACAACCCGGGCAGGCTGCCCTGCTCCATTTCGGCAAGCACCCGCCCGTCATCCCGGGCCTGAATGCGCACTTTCAGGGCATTGGGCAGAAAGGCGCGGACCACCTGCCCGc                                                    <  1:1413212/150‑1 (MQ=255)
                                                                                                   gcttagaGGGCTGTGCCTTATCGAGGTGCGCGCTGAACAACCCGGGCAGGCTGCCCTGCTCCATTTCGGCAAGCACCCGCCCGTCATCCCGGGCCTGAATGCGCACTTTCAGGGCATTGGGCAGAAAGGCGCGGACCACCTGCCCGCCTg                                                 >  2:595551/7‑150 (MQ=255)
                                                                                                   gCTGATAGGGCTGTGCCTTATCGAGGTGCGCGCTGAACAACCCGGGCAGGCTGCCCTGCTCCATTTCGGCAAGCACCCGCCCGTCATCCCGGGCCTGAATGCGCACTTTCAGGGCATTGGGCAGAAAGGCGCGGACCACCTGCCCGCCtt                                                 >  2:516407/1‑149 (MQ=255)
                                                                                                   gCTGATAGGGCTGTGCCTTATCGAGGTGCGCGCTGAACAACCCGGGCAGGCTGCCCTGCTCCATTTCGGCAAGCACCCGCCCGTCATCCCGGGCCTGAATGCGCACTTTCAGGGCATTGGGCAGAAAGGCGCGGACCACCTGCCCGCCTg                                                 >  1:1691612/1‑150 (MQ=255)
                                                                                                     tGATAGGGCTGTGCCTTATCGAGGTGCGCGCTGAACAACCCGGGCAGGCTGCCCTGCTCCATTTCGGCAAGCACCCGCCCGTCATCCCGGGCCTGAATGCGCACTTTCAGGGCATTGGGCAGAAAGGCGCGGACCACCTGCCCGCCTgcg                                               >  1:1674712/1‑150 (MQ=255)
                                                                                                        tAGGGCTGTGCCTTATCGAGGTGCGCGCTGAACAACCCGGGCAGGCTGCCCTGCTCCATTTCGGCAAGCACCCGCCCGTCATCCCGGGCCTGAATgcgc                                                                                               >  1:1374707/1‑99 (MQ=255)
                                                                                                        tAGGGCTGTGCCTTATCGAGGTGCGCGCTGAACAACCCGGGCAGGCTGCCCTGCTCCATTTCGGCAAGCACCCGCCCGTCATCCCGGGCCTGAATgcgc                                                                                               <  2:1374707/99‑1 (MQ=255)
                                                                                                           ggCTGTGCCTTATCGAGGTGCGCGCTGAACAACCCGGGCAGGCTGCCCTGCTCCATTTCGGCAAGCACCCGCCCGTCATCCCGGGCCTGAATGCGCACTTTCAGGGCATTTGGCAGAAAGGCGCGGACCACCTGCCCGCCTGcgccgtcg                                         >  1:301251/1‑150 (MQ=255)
                                                                                                                tgCCTTATCGAGGTGCGCGCTGAACAACCCGGGCAGGCTGCCCTGCTCCATTTCGGCAAGCACCCGCCCGTCATCCCGGGCCTGAATGCGCACTTTCAGGGCATTGGGCAGAAAGGCGCGGACCACCTGCCCGCCTGCGCCGTCGCCGTg                                    >  2:52741/1‑150 (MQ=255)
                                                                                                                tgCCTTATCGAGGTGCGCGCTGAACAACCCGGGCAGGCTGCCCTGCTCCATTTCGGCAAGAAACCGCCCGTCAATTGGGGCCTTAGGGCGCACTTTCAGGGCATTGGGCAGAAAGGCGCGGACCACCTGCCCGCCTGCGCCGTCGCCGTg                                    >  1:1484562/1‑150 (MQ=255)
                                                                                                                       tCGAGGTGCGCGCTGAACAACCCGGGCAGGCTGCCCTGCTCCATTTCGGCAAGCACCCGCCCGTCATCCCGGGCCTGAATgc                                                                                                 <  1:1763228/82‑1 (MQ=255)
                                                                                                                       tCGAGGTGCGCGCTGAACAACCCGGGCAGGCTGCCCTGCTCCATTTCGGCAAGCACCCGCCCGTCATCCCGGGCCTGAATgc                                                                                                 >  2:1763228/1‑82 (MQ=255)
                                                                                                                        cGAGGTGCGCGCTGAACAACCCGGGCAGGCTGCCCTGCTCCATTTCGGCAAGCACCCGCCCGTCATCCCGGGCCTGAATGCGCACTTTCAGGGCATTGGGCAGAAAGGCGCGGACCACCTGCCCGCCTGCGCCGTCGCCGTGGGGGCCGa                            <  2:1783775/150‑1 (MQ=255)
                                                                                                                        cGAGGTGCGCGCTGAACAACCCGGGCAGGCTGCCCTGCTCCATTTCGGCAAGCACCCGCCCGTCATCCCGGGCCTGAATGCGCACTTTCAGGGCATTGGGCAGAAAGGCGCGGACCACCTGCCCGCCTGCGCCGTCGCCGTGGGGGCCGa                            <  1:2429156/150‑1 (MQ=255)
                                                                                                                        cGAGGTGCGCGCTGAACAACCCGGGCAGGCTGCCCTGCTCCATTTCGGCAAGCACCCGCCCGTCATCCCGGGCCTGAATGCGCACTTTCAGGGCATTGGGCAGAAAGGCGCGGACCACCTGCCCGCCTGCGCCGTCGCCGTGGGGGCCGa                            <  2:1783785/150‑1 (MQ=255)
                                                                                                                         gAGGTGCGCGCTGAACAACCCGGGCAGGCTGCCCTGCTCCATTTCGGCAAGCACCCGCCCGTCATCCCGGGCCTGAATGCGCACTTTCAGGGCATTGGGCAGAAAGGCGCGGACCACCTGCCCGCTGGCGCCGTCGCCGTGGGGGCCGAg                           <  1:2156180/150‑1 (MQ=255)
                                                                                                                            gTGCGCGCTGAACAACCCGGGCAGGCTGCCCTGCTCCATTTCGGCAAGCACCCGCCCGTCATCCCGGGCCTGAATGCGCACTTTCAGGGCATTg                                                                                >  2:971340/1‑94 (MQ=255)
                                                                                                                            gTGCGCGCTGAACAACCCGGGCAGGCTGCCCTGCTCCATTTCGGCAAGCACCCGCCCGTCATCCCGGGCCTGAATGCGCACTTTCAGGGCATTg                                                                                <  1:971340/94‑1 (MQ=255)
                                                                                                                             tGCGCGCTGAACAACCCGGGCAGGCTGCCCTGCTCCATTTCGGCAAGCACCCGCCCGTCATCCCGGGCCTGAATGCGCACTTTCAGGGCATTGGGCAGAAAGGCGCGGACCACCTGCCCGCCTGCGCCGTCGCCGTGGGGGCCGAGTAcc                       <  2:569564/150‑1 (MQ=255)
                                                                                                                                  gcTGAACAACCCGGGCAGGCTGCCCTGCTCCATTTCGGCAAGCACCCGCCCGTCATCCCGGGCCTGAATGCGCACTTTCAGGGCATTGGGCAGAAAGGCGCGGACCACCTGCCCGCCTGCGCCGTCGCCGTGGGGGCCGAGTACCGCaaa                  >  2:1648864/1‑150 (MQ=255)
                                                                                                                                    tGAACAACCCGGGCAGGCTGCCCTGCTCCATTTCGGCAAGCACCCGCCCGTCATCCCGGGCCTGAATGCGCACTTTCAg                                                                                       <  1:972899/79‑1 (MQ=255)
                                                                                                                                    tGAACAACCCGGGCAGGCTGCCCTGCTCCATTTCGGCAAGCACCCGCCCGTCATCCCGGGCCTGAATGCGCACTTTCAg                                                                                       >  2:972899/1‑79 (MQ=255)
                                                                                                                                    tGAACAACCCGGGCAGGCTGCCCTGCTCCATTTCGGCAAGCACCCGCCCGTCATCCCGGGCCTGAATGCGCACTTTCAGGGCATTGGGCAGAAAGGCGCGGACCACCTGCCCGCCTGCGCCGTCGCCGTGGGGGCCGAGTACCGCAAAAg                <  2:1363974/150‑1 (MQ=255)
                                                                                                                                      aacaacCCGGGCAGGATGCCCTGCGTCATTTCGGCAAGCACCCGCCCGTCATCCCGGGCCTGAATGCGCACTTTCAGGGCATTGGGCAGAAAGGCGCGGACCACCTGCCCGCCTGCGCCGTCGCCGTGGGGGCCGAGTACCGCAAAAGGa              >  2:1971203/1‑150 (MQ=255)
                                                                                                                                       acaacGCGGGCGTGCTTCCCTGATCCATTTCGGCAAGCACCCGCCCGTCATCCCGGTCCTGAACGCGCACTTTCAGGGCATTGGGCAGAAAGGCGCGGACCACCTGCCCGCCTGCGCCGTCGCCGTGGGGGCCGAGTACCGCAAAAGGAt             <  2:301251/150‑1 (MQ=255)
                                                                                                                                       acaacCCGGGCAGGCTGCCCTGCTCCATTTCGGCAAGCACCCGCCCGTCATCCCGGGCCTGAATGCGCACTTTCAGGGCATTGGGCAGAAAGGCGCGGACCACCTGCCCGCCTGCGCCGTCGCCGTGGGGGCCGAGTACCGCAAAAGGAt             >  2:816784/1‑150 (MQ=255)
                                                                                                                                        caacCCGGGCAGGCTGCCCTGCTCCATTTCGGCAAGCACCCGCCCGTCATCCCGGGCCTGAATGCGCACTTTCAGGGCATTGGGCAGAAAGGCGCGGACCACCTGCCCGCCTGCGCCGTCGCCGTGGGGGCCGAGTACCGCAAAAGGATc            <  1:1416179/150‑1 (MQ=255)
                                                                                                                                           cccGGGCAGGCTGCCCTGCTCCATTTCGGCAAGCACCCGCCCGTCATCCCGGGCCTGAATGCGCACTTTCAGGGCATTGGGCAGAAAGGCGCGGACCACCTGCCCGCCTGCGCCGTCGCCGTGGGGGTCGAGTACCGCAAAAGGATCGGc         <  1:52741/150‑1 (MQ=255)
                                                                                                                                           cccGGGCAGGCTGCCCTGCTCCATTTCGGCAAGCACCCGCCCGTCATCCCGGGCCTGAATGCGCACTTTCAGGGCATTGGGCAGAAAGGCGCGGACCACCTGCCCGCCTGCGCCGTCGCCGTGGGGGCCGAGTACCGCAAAAGGATCGGc         >  2:2233756/1‑150 (MQ=255)
                                                                                                                                           cccGGGCAGGCTGCCCTGCTCCATTTCGGCAAGCACCCGCCCGTCATCCCGGGCCTGAATGCGCACTTTCAGGGCATTGGGCAGAAAGGCGCGGACCACCTGCCCGCCTGCGCCGTCGCCGTGGGGGCCGAGTACCGCAAAAGGATCGGc         <  2:166712/150‑1 (MQ=255)
                                                                                                                                           cccGCGCAGGCTGCCCTGCTCCATTGCGGCAAGCACCCGCCCGTCATCCCGGGCCTGAATGCGCACTTTCAGGGCACTGGGCAGAAAGGCGCGGACCACCTGCCCGCCTGCGCCGTCGCCGTGGGGCCCGAGTACCGCAAAAGGATCGGc         <  2:1190570/150‑1 (MQ=255)
                                                                                                                                                  aGGCTGCCCTGCTCCATTTCGGCAAGCACCCGCCCGTCATCCCGGGCCTGAATGCGCACTTTCAGGGCATTGGGCAGAAAGGCGCGGACCACCTGCCCGCCTGCGCCGTCGCCGTGGGGGCCGAGTACCGCAAAAGGATCGGCGTGTTCg  >  2:1926184/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TGCCTTCCGCGAACAAGTGAAGGTCCATGTCTCCCAATTGCGGGCCGAAGCTGTAGGGGTCTTCGGTAACTTGCTCGCCACCGGCCCAGACGATGTGAAGCTGATAGGGCTGTGCCTTATCGAGGTGCGCGCTGAACAACCCGGGCAGGCTGCCCTGCTCCATTTCGGCAAGCACCCGCCCGTCATCCCGGGCCTGAATGCGCACTTTCAGGGCATTGGGCAGAAAGGCGCGGACCACCTGCCCGCCTGCGCCGTCGCCGTGGGGGCCGAGTACCGCAAAAGGATCGGCGTGTTCG  >  NC_002947/4579794‑4580089

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: