Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 99,268 A→G 5.8% R232G (AGG→GGG)  PP_0093 → zinc protease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00294799,2680AG 5.8% 213.0 / 2.2 69R232G (AGG→GGG) PP_0093zinc protease
Reads supporting (aligned to +/- strand):  ref base A (37/28);  new base G (2/2);  total (39/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.33e-01

CTACGTTGATCTTTCATGAACTGGCCCACCAGCGCTTCTATGTGCAGGACGATACCGAGTTCAATGAATCGTTTGCCTCGTTCGTCGAGCAGGAAGGCACGCGGCAATGGCGGGCGGCGCGTGGGCTGGCGGCGGTTGGCGCAGGGCAGAGGCAGCAGCAGCGTGACCAGTTCATCCGCCTGGTGCTGGCCTGTCGTGAGCGGTTGCAGGCCATTTATGCCGGGCCTTTGGATGATGCGCACAAGCGCCGGGCCAAGCAGGCCGAGTTTGAGCGTTTGAGGCGTGAATACC  >  NC_002947/99119‑99409
                                                                                                                                                     |                                                                                                                                             
cTACGTTGATCTTTCATGAACTGGCCCACCAGCGCTTCTATGTGCAGGACGATACCGAGTTCAATGAATCGTTTGCCTCGTTCGTCGAGCAGGAAGGCACGCGGCAATGGCGGGCGGCGCGTGGGCTGGCGGCGGTTGGCGCAGGGCaga                                                                                                                                               <  2:1757641/150‑1 (MQ=255)
           tttCATGAACTGGCCCACCAGCGCTTCTATGTGCAGGACGATACCGAGTTCAATGAATCGTTTGCCTCGTTCGTCGAGCAGGAAGGCACGCGGCAATGGCGGGCGGCGCGTGGGCTGGCGGCGGTTGGCGCAGGGCAGAGGcagcagcag                                                                                                                                    >  1:203378/1‑150 (MQ=255)
            ttCATGAACTGGCCCACCAGCGCTTCTATGTGCAGGACGATACCGAGTTCAATGAATCGTTTGCCTCGTTCGTCGAGCAGGAAGGCACGCGGCAATGGCGGGCGGCGCGTGGGCTGGCGGCGGTTGGCGCAGGGCAGAGGCagcagcagc                                                                                                                                   <  1:2240740/150‑1 (MQ=255)
            ttCATGAACTGGCCCACCAGCGCTTCTATGTGCAGGACGATACCGAGTTCAATGAATCGTTTGCCTCGTTCGTCGAGCAGGAAGGCACGCGGCAATGGCGGGCGGCGCGTGGGCTGGCGGCGGTTGGCGCAGGGCAGAGGCagcagcagc                                                                                                                                   <  2:1779800/150‑1 (MQ=255)
               aTGAACTGGCCCACCAGCGCTTCTATGTGCAGGACGATACCGAGTTCAATGAATCGTTTGCCTCGTTCGTCGAGCAGGAAGGCACGCGGCAATGGCGGGCGGCGCGTGGGCTGGCGGCGGTTGGCGCAGGGCAGAGGCAGCAGCAGCGTg                                                                                                                                >  2:1065279/1‑150 (MQ=255)
                 gAACTGGCCCACCAGCGCTTCTATGTGCAGGACGATACCGAGTTCAATGAATCGTTTGCCTCGTTCGTCGAGCAGGAAGGCACGCGGCAATGGCGGGCGGCGCGTGGGCTGGCGGCGGTTGGCGCAGGGCAGAGGCAGCAGCAGCGTGAc                                                                                                                              <  2:2187089/150‑1 (MQ=255)
                        cccACCAGCGCTTCTATGTGCAGGACGATACCGAGTTCAATGAATCGTTTGCCTCGTTCGTCGAGCAGGAAGGCACGCGGCAATGGCGGGCGGCGCGTGGGGTGGAGGCGGTTGGAGCAGGGCAGAGGCAGCAGCCGCGTGACCAGGTCa                                                                                                                       >  1:417203/1‑150 (MQ=255)
                         ccaccaGCGCTTCTATGTGCAGGACGATACCGAGTTCAATGAATCGTTTGCCTCGTTCGTCGAGCAGGAAGGCACGCGGCAATGGCGGGCGGCGCGTGGGCTGGCGGCGGTTGGCGCAGGGCAGAGGCAGCAGCAGCGTGACCAGTTCAt                                                                                                                      >  1:935677/1‑150 (MQ=255)
                         ccaccaGCGCTTCTATGTGCAGGACGATACCGAGTTCAATGAATCGTTTGCCTCGTTCGTCGAGCAGGAAGGCACGCGGCAATGGCGGGCGGCGCGTGGGCTGGCGGCGGTTGGCGCAGGGCAGAGGCAGCAGCAGCGTGACCAGTTCAt                                                                                                                      >  2:855683/1‑150 (MQ=255)
                           accaGCGCTTCTATGTGCAGGACGATACCGAGTTCAATGAATCGTTTGCCTCGTTCGTCGAGCAGGAAGGCACGCGGCAATGGCGGGCGGCGCGTGGGCTGGCGGCGGTTGGCGCAGGGCAGAGGCAGCAGCAGCGTGACCAGTTCATcc                                                                                                                    <  1:1453483/150‑1 (MQ=255)
                               gcgcTTCTATGTGCAGGACGATACCGAGTTCAATGAATCGTTTGCCTCGTTCGTCGAGCAGGAAGGCACGCGGCAATGGCGGGCGGCGCGTGGGCTGGCGGCGGTTGGCGCAGGGCAGAGgcagc                                                                                                                                         <  1:1404208/125‑1 (MQ=255)
                               gcgcTTCTATGTGCAGGACGATACCGAGTTCAATGAATCGTTTGCCTCGTTCGTCGAGCAGGAAGGCACGCGGCAATGGCGGGCGGCGCGTGGGCTGGCGGCGGTTGGCGCAGGGCAGAGgcagc                                                                                                                                         >  2:1404208/1‑125 (MQ=255)
                               gcgcTTCTATGTGCAGGACGATACCGAGTTCAATGAATCGTTTGCCTCGTTCGTCGAGCAGGAAGGCACGCGGCAATGGCGGGCGGCGCGTGGGCTGGCGGCGGTTGGCGCAGGGCAGAGGCAGCAGCAGCGTGACCAGTTCATCCGCCt                                                                                                                >  1:958452/1‑150 (MQ=255)
                                  cTTCTATGTGCAGGACGATACCGAGTTCAATGAATCGTTTGCCTCGTTCGTCGAGCAGGAAGGCACGCGGCAATGGCGGGCGGCGCGTGGGCTGGCGGCGGTTGGCGCAGGGCAGAGGCAGCAGCAGCGTGACCAGTTCATCCGCCTGGt                                                                                                             <  1:855683/150‑1 (MQ=255)
                                     cTATGTGCAGGACGATACCGAGTTCAATGAATCGTTTGCCTCGTTCGTCGAGCAGGAAGGCACGCGGCAATGGCGGGCGGCGCGTGGGCTGGCGGCGGTTGGCGCAGGGCAGAGGCAGCAGCAGCGTGACCAGTTCATCCGCCTGGTGCt                                                                                                          >  2:78213/1‑150 (MQ=255)
                                          tgCAGGACGATACCGAGTTCAATGAATCGTTTGCCTCGTTCGTCGAGCAGGACGGCAAGCGGCAATGGCGGGCGGCGCGTGGGCTGGCGGCGGTTGGCGCAGGGCAGAGGCAGCAGCAGCGTGACCAGTTCCTCCGCCTGGTGCTGGCCt                                                                                                     >  2:467506/1‑150 (MQ=255)
                                          tgCAGGACGATACCGAGTTCAATGAATCGTTTGCCTCGTTCGTCGAGCAGGAAGGCACGCGGCAATGGCGGGCGGCGCGTGGGCTGGCGGCGGTTGGCGCAGGGCAGGGGCAGCAGCAGCGTGACCAGTTCATCCGCCTGGTGCTGGCCt                                                                                                     >  1:2010100/1‑150 (MQ=255)
                                             aGGACGATACCGAGTTCAATGAATCGTTTGCCTCGTTCGTCGAGCAGGAAGGCACGCGGCAATGGCGGGCGGCGCGTGGGCTGGCGGCGGTTGGCGCAGGGCAGAGGCAGCAGCAGCGt                                                                                                                                 <  1:1294798/119‑1 (MQ=255)
                                             aGGACGATACCGAGTTCAATGAATCGTTTGCCTCGTTCGTCGAGCAGGAAGGCACGCGGCAATGGCGGGCGGCGCGTGGGCTGGCGGCGGTTGGCGCAGGGCAGAGGCAGCAGCAGCGt                                                                                                                                 >  2:1294798/1‑119 (MQ=255)
                                                   aTACCGAGTTCAATGAATCGTTTGCCTCGTTCGTCGAGCAGGAAGGCACGCGGCAATGGCGGGCGGCGCGTGGGCTGGCGGCGGTTGGCGCAGGGCAGAGGCAGCAGCAGCGTGACCAGTTCATCCGCCTGGTGCTGGCCTGTCGTGAGc                                                                                            >  2:1241259/1‑150 (MQ=255)
                                                       cGAGTTCAATGAATCGTTTGCCTCGTTCGTCGAGCAGGAAGGCACGCGGCAATGGCGGGCGGCGCGTGGGCTGGCGGGGGTTGGCGCAGGGCAGAGGCAGCAGCAGCGTGACCAGTTCATCCGCCTGGTGCTGGCGTGTCGTGCGCGGtt                                                                                        >  1:410409/1‑150 (MQ=255)
                                                        gAGTTCAATGAATCGTTTGCCTCGTTCGTCGAGCAGGAAGGCACGCGGCAATGGCGGGCGGCGCGTGGGCTGGCGGCGGTTGGCGCAGGGCAGAGGCAGCAGCAGCGTGACCAGTTCATCCGCCTGGTg                                                                                                            >  2:1240130/1‑129 (MQ=255)
                                                        gAGTTCAATGAATCGTTTGCCTCGTTCGTCGAGCAGGAAGGCACGCGGCAATGGCGGGCGGCGCGTGGGCTGGCGGCGGTTGGCGCAGGGCAGAGGCAGCAGCAGCGTGACCAGTTCATCCGCCTGGTg                                                                                                            <  1:1240130/129‑1 (MQ=255)
                                                              tgcctgtCGTTTGCCTCGTTCGTCGAGCAGGAAGGCACGCGGCAATGGCGGGCGGCGCGTGGGCTGGCGGCGGTTGGCGCAGGGCAGAGGCAGCAGCAGCGTGACCAGTTCATCCGCCTGGTGCTGGCCTGTCGTGAGCGGTTGCAGGcc                                                                                 >  1:83988/7‑150 (MQ=255)
                                                              aaTGAATCGTTTGCCTCGTTCGTCGAGCAGGAAGGCACGCGGCAATGGCGGGCGGCGCGTGGGCTGGCGGCGGTTGGCGCAGGGCAGAGGCAGCAGCAGCGTGACCAGTTCATCCGCCTGGTGCTGGCCTGTCGTGAGCGGTTGCAGGcc                                                                                 >  1:300577/1‑150 (MQ=255)
                                                               aTGAATCGTTTGCCTCGTTCGTCGAGCAGGAAGGCACGCGGCAATGGCGGGCGGCGCGTGGGCTGGCGGCGGTTGGCGCAGGGCAGAGGCAGCAGCAGCGTGACCAGTTCATCCGCCTGGTGCTGGCCTGTCGTGAg                                                                                             >  1:69932/1‑137 (MQ=255)
                                                               aTGAATCGTTTGCCTCGTTCGTCGAGCAGGAAGGCACGCGGCAATGGCGGGCGGCGCGTGGGCTGGCGGCGGTTGGCGCAGGGCAGAGGCAGCAGCAGCGTGACCAGTTCATCCGCCTGGTGCTGGCCTGTCGTGAg                                                                                             <  2:69932/137‑1 (MQ=255)
                                                               aTGAATCGTTTGCCTCGTTCGTCGAGCAGGAAGGCACGCGGCAATGGCGGGCGGCGCGTGGGCTGGCGGCGGTTGGCGCAGGGCAGAGGCAGCAGCAGCGTGACCAGTTCATCCGCCTGGTGCTGGCCTGTCGTGAGCGGTTGCAGGCCa                                                                                <  2:300577/150‑1 (MQ=255)
                                                                  aaTCGTTTGCCTCGTTCGTCGAGCAGGAAGGCACGCGGCAATGGCGGGCGGCGCGTGGGCTGGCGGCGGTTGGCGCAGGGCAGAGGCAGCAGCAGCGTGACCAGTTCATCCGCCTGGTGCTGGCCTGTCGTGAGCGGTTGCAGGCCAttt                                                                             <  1:1241259/150‑1 (MQ=255)
                                                                    tCGTTTGCCTCGTTCGTCGAGCAGGAAGGCACGCGGCAATGGCGGGCGGCGCGTGGGCTGGCGGCGGTTGGCGCAGGGCAGGGGCAGCAGCAGCGTGACCAGTTCATCCGCCTGGTGCTGGCCt                                                                                                     >  1:765933/1‑124 (MQ=255)
                                                                    tCGTTTGCCTCGTTCGTCGAGCAGGAAGGCACGCGGCAATGGCGGGCGGCGCGTGGGCTGGCGGCGGTTGGCGCAGGGCAGGGGCAGCAGCAGCGTGACCAGTTCATCCGCCTGGTGCTGGCCt                                                                                                     <  2:765933/124‑1 (MQ=255)
                                                                       tttGCCTCGTTCGTCGAGCAGGAAGGCACGCGGCAATGGCGGGCGGCGCGTGGGCTGGCGGCGGTTGGCGCAGGGCAGAGGCAGCAGCAGCGTGACCAGTTCATCCGCCTGGTGCTGGCCTGTCGTGAGCGGTTGCAGGCCATTTATGcc                                                                        >  2:1973130/1‑150 (MQ=255)
                                                                        ttGCCTCGTTCGTCGAGCAGGAAGGCACGCGGCAATGGCGGGCGGCGCGTGGGCTGGCGGCGGTTGGCGCAGGGCAGAGGCAGCAGCAGCGTGACCAGTTCATCCGCCTGGTGCTGGCCTGTCGTGAGCGGTTGCAGGCCATTTATGCCg                                                                       <  2:958452/150‑1 (MQ=255)
                                                                            cTCGTTCGTCGAGCAGGAAGGCACGCGGCAATGGCGGGCGGCGCGTGGGCTGGCGGCGGTTGGCGCAGGGGAGAGgca                                                                                                                                           <  2:1787846/78‑1 (MQ=255)
                                                                            cTCGTTCGTCGAGCAGGAAGGCACGCGGCAATGGCGGGCGGCGCGTGGGCTGGCGGCGGTTGGCGCAGGGGAGAGgca                                                                                                                                           >  1:1787846/1‑78 (MQ=255)
                                                                            cTCGTTCGTCGAGCAGGAAGGCACGCGGCAATGGCGGGCGGCGCGTGGGCTGGCGGCGGTTGGCGCAGGGCAGAGGCAGCAGCAGCGTGACCAGTTCATCCGCCTGGTGCTGGCCTGTCGTGAGCGGTTGCAGGCCATTTATGCCGGGcc                                                                   >  1:2191111/1‑150 (MQ=255)
                                                                            cTCGTTCGTCGAGCAGGAAGGCACGCGGCAATGGCGGGCGGCGCGTGGGCTGGCGGCGGTTGGCGCAGGGCAGAGGCAGCAGCAGCGTGACCAGTTCATCCGCCTGGTGCTGGCCTGTCGTGAGCGGTTGCAGGCCATTTATGCCGGGcc                                                                   >  1:515242/1‑150 (MQ=255)
                                                                             tcgttcgtCGAGCAGGAAGGCACGCGGCAATGGCGGGCGGCGCGTGGGCTGGCGGCGGTTGGCGCAGGGCAGAGGCAGCAGCAGCGTGACCAGTTCATCCGCCTGGTGCTGGCCTGTCGTGAGCGGTTGCAGGCCATTTATGCCGGGCCt                                                                  >  2:1862540/1‑150 (MQ=255)
                                                                             tcgttcgtCGAGCAGGAAGGCACGCGGCAATGGCGGGCGGCGCGTGGGCTGGCGGCGGTTGGCGCAGGGCAGAGGCAGCAGCAGCGTGACCAGTTCATCCGCCTGGTGCTGGCCTGTCGTGAGCGGTTGCAGGCCATTTATGCCGGGCCt                                                                  >  1:2158618/1‑150 (MQ=255)
                                                                               gttcgtCGAGCAGGAAGGCACGCGGCAATGGCGGGCGGCGCGTGGGCTGGCGGCGGTTGGCGCAGGGCAGAGGCAGCAGCAGCGTGACCAGTTCATCCGCCTGGTGCTGGCCTGTCGTGAGCGGTTGCAGGCCATTTATGCCGGGCCttt                                                                >  2:2255237/1‑150 (MQ=255)
                                                                                 tcgtcgAGCAGGAAGGCACGCGGCAATGGCGGGCGGCGCGTGGGCTGGCGGCGGTTGGCGCAGGGCAGAGGCAGCAGCAGCGTGACCAGTTCATCCGCCTGGTGCTGGCCTGTCGTGAGCGGTTGCAGGCCATTTATGCCGGGCCTTTgg                                                              <  2:1379465/150‑1 (MQ=255)
                                                                                   gtcgAGCAGGAAGGCACGCGGCAATGGCGGCCGGCCCACGCGGCGCAGGATAGTCCCTCTGGACAGACTCATCAGTAGTGTGACCAGTTCATCCGACTGGTGCTGGCGTGTCGTGAGCGGTTGCAGTCCATTTATGCCGGGCCTTTGgat                                                            >  1:656725/1‑150 (MQ=255)
                                                                                    tcgAGCAGGAAGGCACGCGGCAATGGCGGGCGGCGCGTGGGCTGGCGGCGGTTGGCGCAGGGCAGAGGCAGCAGCAGCGTGACCAGTTCATCCGCCTGGTGCTGGCCTGTc                                                                                                  >  2:1116366/1‑111 (MQ=255)
                                                                                    tcgAGCAGGAAGGCACGCGGCAATGGCGGGCGGCGCGTGGGCTGGCGGCGGTTGGCGCAGGGCAGAGGCAGCAGCAGCGTGACCAGTTCATCCGCCTGGTGCTGGCCTGTc                                                                                                  <  1:1116366/111‑1 (MQ=255)
                                                                                        gCAGGAAGGCACGCGGCAATGGCGGGCGGCGCGTGGGCTGGCGGCGGTTGGCGCAGGGCAGAGGCAGCAGCAGCGTGACCAGTTCATCCGCCTGGTGCTGGCCTGTCGTGAGCGGTTGCAGGCCATTTATGCCGGGCCTTTGGATGATgc                                                       >  2:1423042/1‑150 (MQ=255)
                                                                                          aGGAAGGCACGCGGCAATGGCGGGCGGCGCGTGGGCTGGCGGCGGTTGGCGCAGGGCAGAGGCAGCAGCAGCGTGACCAGtt                                                                                                                         >  2:1705220/1‑82 (MQ=255)
                                                                                          aGGAAGGCACGCGGCAATGGCGGGCGGCGCGTGGGCTGGCGGCGGTTGGCGCAGGGCAGAGGCAGCAGCAGCGTGACCAGtt                                                                                                                         <  1:1705220/82‑1 (MQ=255)
                                                                                           ggAAGGCACGCGGCAATGGCGGGCGGCGCGTGGGCTGGCGGCGGTTGGCGCAGGGCAGAGGCAGCAGCAGCGTGACCAGTTCATCCGCCTGGTGCTGGCCTGTCGTGAGCGGTTGCAGGCCATTTATGCCGGGCCTTTGGATGATGCGca                                                    <  1:1423042/150‑1 (MQ=255)
                                                                                            gAAGGCACGCGGCAATGGCGGGCGGCGCGTGGGCTGGCGGCGGTTGGCGCAGGGCAGAGGCAGCAGCAGCGTGACCAGTTCATCCGCCTGGTGCTGGCCTGTCGTGAGCGGTTGCAGGCCATTTATGCCGGGCCTTTGGATGATGCGcac                                                   >  2:722428/1‑150 (MQ=255)
                                                                                               ggCACGCGGCAATGGCGGGCGGCGCGTGGGCTGGCGGCGGTTGGCGCAGGGCAGAGGCAGCAGCAGCGTGACCAGTTCATCCGCCTGGTGCTGGCCTGTCGTGAGCGGTTGCAGGCCATTTATGCCGGGCCTTTGGATGATGCGCACAAg                                                >  1:1616987/1‑150 (MQ=255)
                                                                                                 cACGCGGCAATGGCGGGCGGCGCGTGGGCTGGCGGCGGTTGGCGCAGGGCAGAGGCAGCAGCAGCGTGACCAGTTCATCCGCCTGGTGCTGGCCTGTCGTGAGCGGTTGCAGGCCATTTATGCCGGGCCTTTGGATGATGCGCACAAgcg                                              >  2:19760/1‑150 (MQ=255)
                                                                                                    gcgGCAATGGCGGGCGGCGCGTGGGCTGGCGGCGGTTGGCGCAGGGCAGAGGCAGCAGCAGCGt                                                                                                                                 <  1:272314/64‑1 (MQ=255)
                                                                                                    gcgGCAATGGCGGGCGGCGCGTGGGCTGGCGGCGGTTGGCGCAGGGCAGAGGCAGCAGCAGCGt                                                                                                                                 >  2:272314/1‑64 (MQ=255)
                                                                                                      ggCAATGGCGGGCGGCGCGTGGGCTGGCGGCGGTTGGCGCAGGGAAGAGGCAGCAGCAGCGTGACCAGTTCATCCGCCTGGTGCTGGCCTGTCGTGAGCGGTTGCAGGCCATTTATGCCGGGCCTTTGGATGATGCGCACAAGCGCCggg                                         <  2:271012/150‑1 (MQ=255)
                                                                                                         aaTGGCGGGCGGCGCGTGGGCTGGCGGCGGTTGGCGCAGGGCAGAGGCAGCAGCAGCGTGACCAGTTCATCCGCCTGGTGCTGGCCTGTCGTGa                                                                                              >  1:1321666/1‑94 (MQ=255)
                                                                                                         aaTGGCGGGCGGCGCGTGGGCTGGCGGCGGTTGGCGCAGGGCAGAGGCAGCAGCAGCGTGACCAGTTCATCCGCCTGGTGCTGGCCTGTCGTGa                                                                                              <  2:1321666/94‑1 (MQ=255)
                                                                                                                   ggcgCGTGGGCTGGCGGCGGTTGGCGCAGGGCAGAGGCAGCAGGAGCGTGACCAGTTCATCCGCCTGGTGCTGGCCTGTCGTGAGCGGTTGCAGGCCATTTATGCCGGGCCTTTGGATGATGCGCACAAGCGCCGGGCCAAGCAGGCCGa                            <  2:935677/150‑1 (MQ=255)
                                                                                                                   ggcgCGTGGGCTGGCGGCGGTTGGCGCAGGGCAGAGGCAGCAGCAGCGTGACCAGTTCATCCGCCTGGTGCTGGCCTGTCGTGAGCGGTTGCAGGCCATTTATGCCgg                                                                      >  1:880899/1‑108 (MQ=255)
                                                                                                                   ggcgCGTGGGCTGGCGGCGGTTGGCGCAGGGCAGAGGCAGCAGCAGCGTGACCAGTTCATCCGCCTGGTGCTGGCCTGTCGTGAGCGGTTGCAGGCCATTTATGCCgg                                                                      <  2:880899/108‑1 (MQ=255)
                                                                                                                      gcgTGGGCTGGCGGCGGTTGGCGCAGGGCAGAGGCAGCAGCAGCGTGACCAGTTCATCCGCCTGGTGCTGGCCTGTCGTGAGCGGTTGCAGGCCATTTATGCCGGGCCTTTGGATGATGCGCACAAGCGCCGGGCCAAGCAGGCCGAGtt                         >  1:1115447/1‑150 (MQ=255)
                                                                                                                         tGGGCTGGCGGCGGTTGGCGCAGGGCAGGGGCAGCAGCAGCGTGACCAGTTCATCCGCCTGGTGCTGGCCTGTCGTGAGCGGTTGCAGGCCATTTATGCCGGGCCTTTGGATGATGCGCACAAGCGCCGGGCCAAGCAGGCCGAGTTTGa                      <  2:2010100/150‑1 (MQ=255)
                                                                                                                         tGGGCTGGCGGCGGTTGGCGCAGGGCAGAGGCAGCAGCAGCGTGACCAGTTCATCCGCCTGGTGCTGGCCTGTCGTGAGCGGTTGCAGGCCATTTATGCCGGGCCTTTGGATGATGCGCACAAGCGCCGGGCCAAGCAGGCCGAGTTTGa                      <  1:1808315/150‑1 (MQ=255)
                                                                                                                             cTGGCGGCGGTTGGCGCAGGGCAGAGGCAGCAGCAGCGTGACCAGTTCATCCGCCTGGTGCTGGCCTGTCGTGAGCGGTTGCAGGCCATTTATGCCGGGCCTTTGGATGATGCGCACAAGCGCCGGGCCAAGCAGGCCGAGTTTGAGCGt                  <  2:2191111/150‑1 (MQ=255)
                                                                                                                               ggcggcGGTTGGCGCAGGGCAGAGGCAGCAGCAGCGTGACCAGTTCATCCGCCTGGTGCTGGCCTGTCGTGAGCGGTTGCAGGCCATTTATGCCGGGCCTTTGGATGATGCGCACACGCGCCGGGCCAAGCAGGCCGAGTTTGAGCGttt                <  1:1862540/150‑1 (MQ=255)
                                                                                                                                 cggcggTTGGCGCAGGGCAGAGGCAGCAGCAGCGTGACCAGTTCATCCGCCTGGTGCTGGCCTGTCGTGAGCGGTTGCAGGCCATTTATGCCGGGCCTTTGGATGATGCGCACAAGCGCCGGGCCAAGCAGGCCGAGTTTGAGCGTTTGa              <  2:2158618/150‑1 (MQ=255)
                                                                                                                                  ggcggTTGGCGCAGGGCAGAGGCAGCAGCAGCGTGACCAGTTCATCCGCCTGGTGCTGGCCTGTCGTGAGCGGTTGCAGGCCATTTATGCCGGGCCTTTGGATGATGCGCACAAGCGCCGGGCCAAGCAGGCCGAGTTTGAGCGTTTGAg             >  2:2388379/1‑150 (MQ=255)
                                                                                                                                      gTTGGCGCAGGGCAGAGGCAGCAGCAGCGTGACCAGTTCATCCGCCTGGTGCTGGTCTGTCGTGAGCGGTTGCAGGCCATTTATGCCGGGCCTTTGGATGATGCGCACAAGCGCCGGGCCAAGCAGGCCGAGTTTGAGCGTTTGAGGCGt         <  2:1959324/150‑1 (MQ=255)
                                                                                                                                           cgcAGGGCAGAGGCAGCAGCAGCGTGACCAGTTCATCCGCCTGGTGCTGGCCTGTCGTGAGCGGTTGCAGGCCATTTATGCCGGGCCTTTGGATGATGCGCACAAGCGCCGGGCCAAGCAGGCCGAGTTTGAGCGTTTGAGGCGTGAATa    >  1:1775019/1‑150 (MQ=255)
                                                                                                                                           cgcAGGGCAGAGGCAGCAGCAGCGTGACCAGTTCATCCGCCTGGTGCTGGCCTGTCGTGAGCGGTTGCAGGCCATTTATGCCGGGCCTTTGGATGATGCGCACAAGCGCCGGGCCAAGCAGGCCGAGTTTGAGCGTTTGAGGCGTGAATa    >  1:317013/1‑150 (MQ=255)
                                                                                                                                             cAGGGCAGAGGCAGCAGCAGCGTGACCAGTTCATCCGCCTGGTGCTGGCCTGTCGTGAGCGGTTGCAGGCCATTTATGCCGGGCCTTTGGATGATGCGCACAAGCGCCGGGCCAAGCAGGCCGAGTTTGAGCGTTTGAGGCGTGAATAcc  <  2:1616987/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                             
CTACGTTGATCTTTCATGAACTGGCCCACCAGCGCTTCTATGTGCAGGACGATACCGAGTTCAATGAATCGTTTGCCTCGTTCGTCGAGCAGGAAGGCACGCGGCAATGGCGGGCGGCGCGTGGGCTGGCGGCGGTTGGCGCAGGGCAGAGGCAGCAGCAGCGTGACCAGTTCATCCGCCTGGTGCTGGCCTGTCGTGAGCGGTTGCAGGCCATTTATGCCGGGCCTTTGGATGATGCGCACAAGCGCCGGGCCAAGCAGGCCGAGTTTGAGCGTTTGAGGCGTGAATACC  >  NC_002947/99119‑99409

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