Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 873,421 G→A 9.9% G18S (GGC→AGC)  PP_5457 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947873,4210GA 9.9% 304.7 / 19.0 111G18S (GGC→AGC) PP_5457hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (43/56);  new base A (4/7);  total (47/64)
Fisher's exact test for biased strand distribution p-value = 7.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

Only 100 of 121 total aligned reads displayed.

CCGAGGATTGCCAGGTGTGCTGTCAGCCTGTGGTGTGCATTCTGCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCTCGTCGTGGCTGATGAGCAGGCCGAATGCGCACGGCAGCTGATCGATGCGTACAA  >  NC_002947/873273‑873569
                                                                                                                                                    |                                                                                                                                                    
ccGAGGATTGCCAGGTGTGCTGTCAGCCTGTGGTGTGCATTCTGCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTggcagg                                                                                                                                                     >  2:1764370/1‑150 (MQ=255)
 cGAGGATTGCCAGGTGTGCTGTCAGCCTGTGGTGTGCATTCTGCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGAGCAAATGCTggcaggc                                                                                                                                                    >  1:704477/1‑150 (MQ=255)
   aGGATTGCCAGGTGTGCTGTCAGCCTGTGGTGTGCATTCTGCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAAGCAt                                                                                                                                                  >  2:1701576/1‑150 (MQ=255)
       ttGCCAGGTGTGCTGTCAGCCTGTGGTGTGCATTCTGCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCtt                                                                                                                                              <  2:260157/150‑1 (MQ=255)
         gCCAGGTGTGCTGTCAGCCTGTGGTGTGCATTCTGCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGc                                                                                                                                            <  2:180513/150‑1 (MQ=255)
          ccAGGTGTGCTGTCAGCCTGTGGTGTGCATTCTGCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGcc                                                                                                                                           >  1:2379145/1‑150 (MQ=255)
           cAGGTGTGCTGTCAGCCTGTGGTGTGCATTCTGCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAAGCATGCCTGCCa                                                                                                                                          <  2:857828/150‑1 (MQ=255)
           cAGGTGTGCTGTCAGCCTGTGGTGTGCATTCTGCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAAGCATGCCTGCCa                                                                                                                                          <  2:1456480/150‑1 (MQ=255)
               tgtgCTGTCAGCCTGTGGTGTGCATTCTGCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGa                                                                                                                                      <  1:1917354/150‑1 (MQ=255)
                gtgCTGTCAGCCTGTGGTGTGCATTCTGCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAAGCATGCCTGCCAGcatt                                                                                                                                     <  2:1450864/150‑4 (MQ=255)
                 tgCTGTCAGCCTGTGGTGTGCATTCTGCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAgg                                                                                                                                    >  2:583436/1‑150 (MQ=255)
                  gCTGTCAGCCTGTGGTGTGCATTCTGCAGGTGCATGGTGACGAATGGATGATTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAggg                                                                                                                                   <  1:121707/150‑1 (MQ=255)
                   cTGTCAGCCTGTGGTGTGCATTCTGCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAAGCATGCCTGCCAGcatttgc                                                                                                                                  <  1:1701576/150‑7 (MQ=255)
                    tGTCAGCCTGTGGTGTGCATTCTGCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGa                                                                                                                                 <  1:2139693/150‑1 (MQ=255)
                        aGCCTGTGGTGTGCATTCTGCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGa                                                                                                                             >  1:525463/1‑150 (MQ=255)
                         gCCTGTGGTGTGCATTCTGCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGAt                                                                                                                            <  2:2063710/150‑1 (MQ=255)
                             gtggtgTGCATTCTGCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTac                                                                                                                        >  1:1824825/1‑150 (MQ=255)
                               ggtgTGCATTCTGCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTacac                                                                                                                      <  1:583436/150‑1 (MQ=255)
                                gtgtGCATTCTGCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACg                                                                                                                     >  1:670397/1‑150 (MQ=255)
                                   tgCATTCTGCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCa                                                                                                                                          <  2:2026301/126‑1 (MQ=255)
                                   tgCATTCTGCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCa                                                                                                                                          >  1:2026301/1‑126 (MQ=255)
                                     cATTCTGCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCa                                                                                                                                          >  2:2056034/1‑124 (MQ=255)
                                     cATTCTGCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCa                                                                                                                                          <  1:2056034/124‑1 (MQ=255)
                                     cATTCTGCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTc                                                                                                                <  2:2214969/150‑1 (MQ=255)
                                       ttCTGCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCgg                                                                                                              <  1:158175/150‑1 (MQ=255)
                                       ttCTGCAGGTGCATGGGGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCgg                                                                                                              <  2:205536/150‑1 (MQ=255)
                                         cTGCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGcc                                                                                                            <  1:827418/150‑1 (MQ=255)
                                             aGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGa                                                                                                        >  1:2137202/1‑150 (MQ=255)
                                                 gCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCtt                                                                                                    >  1:1040690/1‑150 (MQ=255)
                                                   aTGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCAGTGTCCCTGt                                                                                                  <  2:1300588/150‑1 (MQ=255)
                                                      gTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTggcagg                                                                                                                                                     <  2:440564/96‑1 (MQ=255)
                                                      gTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTggcagg                                                                                                                                                     >  1:440564/1‑96 (MQ=255)
                                                       tGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAgg                                                                                                                                    >  1:1471637/1‑112 (MQ=255)
                                                       tGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAgg                                                                                                                                    <  2:1471637/112‑1 (MQ=255)
                                                         aCGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGcc                                                                                                            >  1:1782346/1‑134 (MQ=255)
                                                         aCGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGcc                                                                                                            <  2:1782346/134‑1 (MQ=255)
                                                         aCGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGt                                                                                              >  1:896870/1‑148 (MQ=255)
                                                          cGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTgcc                                                                                           >  2:513260/1‑150 (MQ=255)
                                                          cGAATGGATGCTTGAAGTCCGCCGCGAGGACGAAACCTGAGGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTgcc                                                                                           <  2:257082/150‑1 (MQ=255)
                                                                 aTGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGAt                                                                                                                            <  2:768111/110‑1 (MQ=255)
                                                                 aTGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGAt                                                                                                                            >  1:768111/1‑110 (MQ=255)
                                                                     ttGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGAt                                                                                                                            >  2:1602044/1‑106 (MQ=255)
                                                                     ttGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGAt                                                                                                                            <  1:1602044/106‑1 (MQ=255)
                                                                     ttGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCt                                                                                >  2:709972/1‑150 (MQ=255)
                                                                     ttGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAAGCATGCCTGCCa                                                                                                                                          >  1:451507/1‑92 (MQ=255)
                                                                     ttGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAAGCATGCCTGCCa                                                                                                                                          <  2:451507/92‑1 (MQ=255)
                                                                      tGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGa                                                                                                                             <  1:406562/104‑1 (MQ=255)
                                                                      tGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGa                                                                                                                             >  2:406562/1‑104 (MQ=255)
                                                                           tCCGCCGCGAGGACGACACCTGATGCATCGTATCTATGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTacac                                                                                                                      <  1:1020909/106‑1 (MQ=255)
                                                                           tCCGCCGCGAGGACGACACCTGATGCAGCGTATCTATGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTacac                                                                                                                      >  2:1020909/1‑106 (MQ=255)
                                                                                cgcgAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAgg                                                                     <  2:1040690/150‑1 (MQ=255)
                                                                                cgcgAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAgg                                                                     <  1:1057100/150‑1 (MQ=255)
                                                                                  tgAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGc                                                                   <  2:1097897/149‑1 (MQ=255)
                                                                                  cgAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGc                                                                   <  2:2379145/150‑1 (MQ=255)
                                                                                         gacACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTgg                                                            <  1:137802/150‑1 (MQ=255)
                                                                                           cacCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTgggg                                                          <  2:1569631/150‑1 (MQ=255)
                                                                                            acCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGAGCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGc                                                         >  2:1849896/1‑150 (MQ=255)
                                                                                              cTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCtc                                                       >  1:186172/1‑150 (MQ=255)
                                                                                               tGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCtcg                                                      >  2:1586288/1‑150 (MQ=255)
                                                                                                gATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCtcgt                                                     <  1:1586288/150‑1 (MQ=255)
                                                                                                 gtGCAGCGTATCTAAGAGCGGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCtcgtc                                                    <  2:621510/149‑1 (MQ=255)
                                                                                                  tGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAAGCATGCCTGCCAGcatt                                                                                                                                     <  1:696408/68‑4 (MQ=37)
                                                                                                  tGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAAGCATGCCTGCCAGcatt                                                                                                                                     >  2:696408/1‑65 (MQ=37)
                                                                                                   gCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCTcgtcgt                                                  <  2:2234957/150‑1 (MQ=255)
                                                                                                     aGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCTCGTCGTgg                                                <  1:1018812/150‑1 (MQ=255)
                                                                                                      gCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCGCGTCGTGGc                                               >  1:90232/1‑150 (MQ=255)
                                                                                                         tATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCTCGTCGTGGCtga                                            <  1:384355/150‑1 (MQ=255)
                                                                                                          aTCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAgg                                                                     <  1:569206/124‑1 (MQ=255)
                                                                                                          aTCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAgg                                                                     >  2:569206/1‑124 (MQ=255)
                                                                                                            cTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGa                                                                                                                                      >  2:1869257/1‑57 (MQ=255)
                                                                                                            cTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGa                                                                                                                                      <  1:1869257/57‑1 (MQ=255)
                                                                                                             tACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCTCGTCGTGGCTGATGAg                                        >  2:1889829/1‑150 (MQ=255)
                                                                                                             tACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCTCGTCGTGGCTGATGAg                                        >  2:2221064/1‑150 (MQ=255)
                                                                                                                 aaCCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCTCGTCGTGGCTGATGAGCAgg                                    <  2:1824482/150‑1 (MQ=255)
                                                                                                                  aCCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCTCGTCGTGGCTGATGAGCAGGc                                   <  2:1920255/150‑1 (MQ=255)
                                                                                                                     ggAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCTCGTCGTGGCTGATGAGCAGGCCGa                                >  2:452560/1‑150 (MQ=255)
                                                                                                                       aaaGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCTCGTCGTGGCTGATGAGCAGGCCGAAt                              <  1:1935649/150‑1 (MQ=255)
                                                                                                                        aaGCCTGCTTGAGGCTCAAATGCTGGTATGCATTCTTGCCAGCTATGTGATCGATGTACACGTGGTCGTCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGTGCCTGCTGGGGCTCGTCGTGGCTGATGATCAGTCCTAATg                             <  2:186172/150‑1 (MQ=255)
                                                                                                                           ccTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCTCGTCGTGGCTGATGAGCAGGc                                   >  1:2251784/1‑141 (MQ=255)
                                                                                                                           ccTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCTCGTCGTGGCTGATGAGCAGGc                                   <  2:2251784/141‑1 (MQ=255)
                                                                                                                               cTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCTCGTCGTGGCTGATGAGCAGGCCGAATGCGCACgg                      <  2:662590/150‑1 (MQ=255)
                                                                                                                                  gAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCTCGTCGTGGCTGATGAGCAGGCCGAATGCGCACGGCAg                   <  1:2221064/150‑1 (MQ=255)
                                                                                                                                    ggCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCTCGTCGTGGCTGATGAGCAGGCCGAATGCGCACGGCAGCt                 >  2:1984518/1‑150 (MQ=255)
                                                                                                                                     gcgcAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCTCGTCGTGGCTGATGAGCAGGCCGAATGCGCACGGCAGCTg                >  1:2131225/1‑150 (MQ=255)
                                                                                                                                       gcAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCTCGTCGTGGCTGATGAGCAGGCCGAATGCGCACGGc                     <  2:1633037/143‑1 (MQ=255)
                                                                                                                                       gcAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCTCGTCGTGGCTGATGAGCAGGCCGAATGCGCACGGc                     >  1:1633037/1‑143 (MQ=255)
                                                                                                                                       gcAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCTCGTCGTGGCTGATGAGCAGGCCGAATGCGCACGGCAGCTGAt              <  1:1526461/150‑1 (MQ=255)
                                                                                                                                        cAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCTCGTCGTGGCTGATGAGCAGGCCGAATGCGCACGGCAGCTGATc             <  1:2378797/150‑1 (MQ=255)
                                                                                                                                        cAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCTCGTCGTGGCTGATGAGCAGGCCGAATGCGCACGGCAGCTGATc             <  1:1889829/150‑1 (MQ=255)
                                                                                                                                              cTGGCAAGCATGCCTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTgggg                                                          <  2:308415/99‑1 (MQ=255)
                                                                                                                                              cTGGCAAGCATGCCTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTgggg                                                          >  1:308415/1‑99 (MQ=255)
                                                                                                                                               tGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCTCGTCGTg                                                 >  1:1491495/1‑107 (MQ=255)
                                                                                                                                               tGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCTCGTCGTg                                                 <  2:1491495/107‑1 (MQ=255)
                                                                                                                                               tGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCTCGTCGTGGCTGATGAGCAGGCCGAATGCGCACGGCAGCTGATc             >  2:2402102/1‑143 (MQ=255)
                                                                                                                                               tGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCTCGTCGTGGCTGATGAGCAGGCCGAATGCGCACGGCAGCTGATc             <  1:2402102/143‑1 (MQ=255)
                                                                                                                                                ggcaggcaTGCTTGCCAGCGAGGGGATCGATGTAAACGTAATTGGGCGACAGGTGTTCGGTGCCGCCGGCGACGACCCCTTGCAGGGCCTGCTGGGGCTCGTCGTGGCTGATGAGCAGGCCGAATGCGCACGGCAGCTGATCGATGCGTa     <  1:2341227/150‑1 (MQ=255)
                                                                                                                                                 gcaggcaTGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGc                                                                   <  1:2186404/87‑1 (MQ=255)
                                                                                                                                                 gcaggcaTGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCTCGTCGTGGCtgatg                                          <  2:639378/112‑1 (MQ=255)
                                                                                                                                                 gcaggcaTGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCTCGTCGTGGCtgatg                                          >  1:639378/1‑112 (MQ=255)
                                                                                                                                                   aggcaTGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCTCGTCGTGGCTGATGAGCAGGCCGAATGCGCACGGCAGCTGATCGATGCGTACaa  >  1:1391062/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CCGAGGATTGCCAGGTGTGCTGTCAGCCTGTGGTGTGCATTCTGCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCTCGTCGTGGCTGATGAGCAGGCCGAATGCGCACGGCAGCTGATCGATGCGTACAA  >  NC_002947/873273‑873569

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: