Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 953,241 G→A 5.0% G16E (GGG→GAG)  cyoC → cytochrome bo terminal oxidase subunit III

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947953,2410GA 5.0% 306.8 / 2.9 99G16E (GGG→GAG) cyoCcytochrome bo terminal oxidase subunit III
Reads supporting (aligned to +/- strand):  ref base G (48/46);  new base A (3/2);  total (51/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

AAGAAGTGGCGCGCATCGAAGGCGAGCGTATGAAAGCGCTGGCCAAAGCAGGTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGTTTGCGTCGCTCTTCGCCACCTACGCGGTGCTGTCCGGCAGTTTTGCCGGCGGCCCGTCGGGTCATGA  >  NC_002947/953094‑953388
                                                                                                                                                   |                                                                                                                                                   
aagaagTGGCGCGCATCGAAGGCGAGCGTATGAAAGCGCTGGCCAAAGCAGGTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGc                                                                                                                                                   >  1:1918897/1‑150 (MQ=255)
   aagTGGCGCGCATCGAAGGCGAGCGTATGAAAGCGCTGGCCAAAGCAGGTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATg                                                                                                                                                >  1:928094/1‑150 (MQ=255)
    agTGGCGCGCATCGAAGGCGAGCGTATGAAAGCGCTGGCCAAAGCAGGTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATga                                                                                                                                               <  1:299455/150‑1 (MQ=255)
    agTGGCGCGCATCGAAGGCGAGCGTATGAAAGCGCTGGCCAAAGCAGGTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATga                                                                                                                                               <  1:1254640/150‑1 (MQ=255)
    agTGGCGCGCATCGAAGGCGAGCGTATGAAAGCGCTGGCCAAAGCAGGTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCCCATGGTCGTGACCATGAGCagcg                                                                                                                                               <  2:392372/150‑4 (MQ=255)
          gcgcATCGAAGGCGAGCGTATGAAAGCGCTGGCCAAAGCAGGTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTTTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGacca                                                                                                                                         >  1:119793/1‑150 (MQ=255)
          gcgcATCGAAGGCGAGCGTATGAAAGCGCTGGCCAAAGCAGGTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGacca                                                                                                                                         >  1:948420/1‑150 (MQ=255)
           cgcATCGAAGGCGAGCGTATGAAAGCGCTGGCCAAAGCAGGTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGaccac                                                                                                                                        <  1:1624551/150‑1 (MQ=255)
              aTCGAAGGCGAGCGTATGAAAGCGCTGGCCAAAGCAGGTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGaccaccac                                                                                                                                     <  2:120547/150‑1 (MQ=255)
                cGAAGGCGAGCGTATGAAAGCGCTGGCCAAAGCAGGTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAGTGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGAccaccacca                                                                                                                                   <  1:175955/150‑1 (MQ=255)
                     gCGAGCGTATGAAAGCGCTGGCCAAAGCAGGTGCTCTGCCTGCCTGCGCACGTGTCGAATCGTTTGAACGGGGGTAATCACTGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGAAGACCACCACCACGACt                                                                                                                              >  1:1571605/1‑150 (MQ=255)
                       gAGCGTATGAAAGCGCTGGCCAAAGCAGGTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCg                                                                                                                            >  2:226441/1‑150 (MQ=255)
                        aGCGTATGAAAGCGCTGGCCAAAGCAGGTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCgg                                                                                                                           <  1:226441/150‑1 (MQ=255)
                        aGCGTATGAAAGCGCTGGCCAAAGCAGGTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCgg                                                                                                                           <  2:624789/150‑1 (MQ=255)
                           gTATGAAAGCGCTGGCCAAAGCAGGTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGcc                                                                                                                        <  2:1403340/150‑1 (MQ=255)
                              tGAAAGCGCTGGCCAAAGCAGGTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGa                                                                                                                     <  1:150554/150‑1 (MQ=255)
                               gAAAGCGCTGGCCAAAGCAGGTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGATGCGCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGAt                                                                                                                    <  1:526383/150‑1 (MQ=255)
                                       tGGCCAAAGCAGGTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTAc                                                                                                            >  2:1732548/1‑150 (MQ=255)
                                       tGGCCAAAGCAGGTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTAc                                                                                                            >  1:1285297/1‑150 (MQ=255)
                                         gCCAAAGCAGGTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTg                                                                                                          <  1:149114/150‑1 (MQ=255)
                                           cAAAGCAGGTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTggg                                                                                                        >  1:308480/1‑150 (MQ=255)
                                           cAAAGCAGGTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATCGTCACGACCATGGGAATAACCTCCTCCACCACGACTCGGGCCAGATGACCGTACTggg                                                                                                        >  2:912690/1‑150 (MQ=255)
                                            aaaGCAGGTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGc                                                                                                       <  1:378235/150‑1 (MQ=255)
                                              aGCAGGGGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCtt                                                                                                     <  2:811167/150‑1 (MQ=255)
                                                  ggTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTgg                                                                                                 <  2:1560202/150‑1 (MQ=255)
                                                  ggTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTgg                                                                                                 <  1:1721694/150‑1 (MQ=255)
                                                  ggTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTgg                                                                                                 >  2:282714/1‑150 (MQ=255)
                                                   gTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGGCGGCCGACACGACCATGGGCATGACGACCACCACCACGACTCGGGATAGATGACCGTACTGGGCTTCTGGc                                                                                                >  1:118071/1‑150 (MQ=255)
                                                   gTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGc                                                                                                <  1:755976/150‑1 (MQ=255)
                                                     gCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTg                                                                                              >  1:1199541/1‑150 (MQ=255)
                                                       tctGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTa                                                                                            <  1:252943/150‑1 (MQ=255)
                                                       tctGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACAACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGtt                                                                                            >  1:1199667/1‑149 (MQ=255)
                                                        ctgcctgcCGGCGCACGTGTCCAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTgaat                                                                                           <  1:1809468/150‑4 (MQ=255)
                                                           cctgccGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCtg                                                                                        <  1:795226/150‑1 (MQ=255)
                                                             tgccGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCtgat                                                                                      <  1:1417661/150‑1 (MQ=255)
                                                                cGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTCGGCTTCTGGCTGTACCTGATGAc                                                                                   <  2:948420/150‑1 (MQ=255)
                                                                  gcgcACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCCCATGGTCGTGACCATGAGCa                                                                                                                                                  <  1:1239231/85‑1 (MQ=38)
                                                                  gcgcACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCCCATGGTCGTGACCATGAGCa                                                                                                                                                  >  2:1239231/1‑85 (MQ=38)
                                                                    gcACGTGTCGAAGCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGAc                                                                               <  2:193279/150‑1 (MQ=255)
                                                                        gtgtCGAATCGTTTGAGCGGGAATAATCAAGGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCa                                                                           >  2:1658458/1‑150 (MQ=255)
                                                                           tCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATcc                                                                        >  1:1105176/1‑150 (MQ=255)
                                                                             gAANCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTg                                                                      >  2:1200601/1‑150 (MQ=255)
                                                                              aaTCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACCAGCGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGt                                                                     >  1:497312/1‑150 (MQ=255)
                                                                                tCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGAc                                                                                   <  2:1694143/134‑1 (MQ=255)
                                                                                tCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGAc                                                                                   >  1:1694143/1‑134 (MQ=255)
                                                                                tCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGttt                                                                   <  2:1442368/150‑1 (MQ=255)
                                                                                 cGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAAGGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGTTTg                                                                  <  1:1277157/150‑1 (MQ=255)
                                                                                   tttGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGTTTGCg                                                                >  2:1308096/1‑150 (MQ=255)
                                                                                   tttGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGTTTGCg                                                                >  1:78395/1‑150 (MQ=255)
                                                                                    ttGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGTTTGCGt                                                               <  1:1308096/150‑1 (MQ=255)
                                                                                     tGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGTGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGTTTGCGTc                                                              <  2:1417156/150‑1 (MQ=255)
                                                                                       aGCGGGTGTAATCAATGTCCAGTGAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGTTTGCGTCGc                                                            >  2:942985/1‑150 (MQ=255)
                                                                                          gggTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCAt                                                                                                                                                 >  2:896726/1‑62 (MQ=255)
                                                                                          gggTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCAt                                                                                                                                                 <  1:896726/62‑1 (MQ=255)
                                                                                            gtgtAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCt                                                                                               >  1:837368/1‑110 (MQ=255)
                                                                                            gtgtAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCt                                                                                               <  2:837368/110‑1 (MQ=255)
                                                                                            gtgtAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCGCATGGTCACGACCATGGGCATGACGACCACCCCCAACGCTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGTTTGCGTCGCTCTTc                                                       <  1:1869230/150‑1 (MQ=255)
                                                                                             tgtAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGTTTGCGTCGCTCTTCg                                                      >  2:293098/1‑150 (MQ=255)
                                                                                             cgtAATTAATCTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGTTTGCGTCGCTCTTCg                                                      >  2:471955/2‑150 (MQ=255)
                                                                                              gtAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGTTTGCGTCGCTCTTCGc                                                     <  1:1892355/150‑1 (MQ=255)
                                                                                               tAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGTTTGCGTCGCTCTTCGcc                                                    >  2:1423885/1‑150 (MQ=255)
                                                                                                   cAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGTTTGCGTCGCTCTTCGCCACCt                                                >  2:1375176/1‑150 (MQ=255)
                                                                                                    aaTGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGa                                                                                                                     <  2:1384774/80‑1 (MQ=255)
                                                                                                    aaTGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGa                                                                                                                     >  1:1384774/1‑80 (MQ=255)
                                                                                                        tCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGTTTGCGTCGCTCTTCGCCACCTAcgcg                                           >  2:52431/1‑150 (MQ=255)
                                                                                                          cAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCAACCATTAATGACCGACTGCATCCTGTTTGCGTCGCTCTTCGCCACCTACGCGGt                                         >  1:593496/1‑150 (MQ=255)
                                                                                                           aGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGTTTGCGTCGCTCTTCGCCACCTACGCGGTg                                        >  2:1644299/1‑150 (MQ=255)
                                                                                                           aGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGTTTGCGTCGCTCTTCGCCACCTACGCGGTg                                        <  1:139846/150‑1 (MQ=255)
                                                                                                             tCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGTTTGCGTCGCTCTTCGCCACCTACGCGGTGCt                                      >  1:1220110/1‑150 (MQ=255)
                                                                                                             tCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGTTTGCGTCGCTCTTCGCCACCTACGCGGTGCt                                      >  1:345257/1‑150 (MQ=255)
                                                                                                               aaGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTAc                                                                                           <  1:873/95‑1 (MQ=255)
                                                                                                               aaGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTAc                                                                                           >  2:873/1‑95 (MQ=255)
                                                                                                                aGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGCCGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGTTTGCGTCGCTCTTCGCCACCTACGCGGTGCTGTc                                   <  1:107765/150‑1 (MQ=255)
                                                                                                                aGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGTTTGCGTCGCTCTTCGCCACCTACGCGGTGCTGTc                                   >  1:828969/1‑150 (MQ=255)
                                                                                                                aGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGTTTGCGTCGCTCTTCGCCACCTACGCGGTGCTGTc                                   <  1:282714/150‑1 (MQ=255)
                                                                                                                  tAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTAGCTGATGACCGAGTGGATCCTGGTTGTGTTGTTCTTCGCTATGTACATGGTGCTggtcg                                 >  2:10946/1‑146 (MQ=255)
                                                                                                                  tAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGTTTGCGTCGCTCTTCGCCACCTACGCGGTGCTGTCCg                                 <  2:1918897/150‑1 (MQ=255)
                                                                                                                  tAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGTTTGCGTCGCTCTTCGCCACCTACGCGGTGCTGTCCg                                 >  1:1422255/1‑150 (MQ=255)
                                                                                                                   aaTGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGTTTGCGTCGCTCTTc                                                       >  2:1239712/1‑127 (MQ=255)
                                                                                                                   aaTGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGTTTGCGTCGCTCTTc                                                       <  1:1239712/127‑1 (MQ=255)
                                                                                                                      gCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGTTTGCGTCGCTCTTCGCCACCTACGCGGTGCTGTCCGGCAg                             <  2:1571605/150‑1 (MQ=255)
                                                                                                                       cACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGTTTGCGTCGCTCTTCGCCACCTACGCGGTGCTGTCCgg                                >  1:1456139/1‑146 (MQ=255)
                                                                                                                       cACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGTTTGCGTCGCTCTTCGCCACCTACGCGGTGCTGTCCgg                                <  2:1456139/146‑1 (MQ=255)
                                                                                                                          ggCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGTTTGCGTCGCTCTTCGCCACCTACGCGGTGCTGTCCGGCAGtttt                         <  1:1423885/150‑1 (MQ=255)
                                                                                                                          ggCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGTTTGCGTCGCTCTTCGCCACCTACGCGGTGCTGGCCGGCAGtttt                         >  1:1340094/1‑150 (MQ=255)
                                                                                                                             gctgctCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGTTTGCGTCGCTCTTCGCCACCTACGCGGTGCTGTCCGGCAGTTTTGcc                      >  1:581374/1‑150 (MQ=255)
                                                                                                                             gctgctCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGTTTGCGTCGCTCTTCGCCACCTACGCGGTGCTGTCCGGCAGTTTTGcc                      >  1:454297/1‑150 (MQ=255)
                                                                                                                               tgctCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGTTTGCGTCGCTCTTCGCCACCTACGCGGTGCTGTCCGGCAGTTTTGCcgg                    <  2:78395/150‑1 (MQ=255)
                                                                                                                                    aTGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGTTTGCGTCGCTCTTCGCCACCTACGCGGTg                                        <  2:466033/125‑1 (MQ=255)
                                                                                                                                    aTGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGTTTGCGTCGCTCTTCGCCACCTACGCGGTg                                        >  1:466033/1‑125 (MQ=255)
                                                                                                                                    aTGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGTTTGCGTCGCTCTTCGCCACCTACGCGGTGCTGTCCGGCAGTTTTGCCGGCGGcc               >  2:225690/1‑150 (MQ=255)
                                                                                                                                    aTGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGTTTGCGTCGCTCTTCGCCACCTACGCGGTGCTGTCCGGCAGTTTTGCCGGCGGcc               >  2:805995/1‑150 (MQ=255)
                                                                                                                                        tCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGTTTGCGTCGCTCTTCGCCACCTACGCGGTGCTGTCCGGCAGTTTTGCCGGCGGCCCGTc           <  1:200090/150‑1 (MQ=255)
                                                                                                                                          gtgACCATGAGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGTTTGCGTCGCTCTTCGCCACCTACGCGGTGCTGTCCGGCAGTTTTGCCGGCGGCCCGTCgg         >  1:859856/3‑150 (MQ=255)
                                                                                                                                          gtgACCATGAGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGGACCTGATGACCGACTGCATCCTGTTTGCGTCGCTCTTCGCCACCTACGCGGTGCTGTCCGGCAGTTTTGCCGGCGGCCCGTCgg         >  1:842502/3‑150 (MQ=255)
                                                                                                                                          aCGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGTTTGCGTCGCTCTTCGCCACCTACGCGGTGCTGTCCGGCAGTTTTGCCGGCGGCCCGTCgg         <  1:1732548/150‑1 (MQ=255)
                                                                                                                                           cGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGTTTGCGTCGCTCTTCGCCACCTACGCGGTGCTGTCCGGCAGTTTTGCCGGCGGCCCGTCggg        >  1:1151125/1‑150 (MQ=255)
                                                                                                                                               cATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGTTTGCGTCGCTCTTCGCCACCTACGCGGTGCTGTCCGGCAGTTTTGCCGGCGGCCCGTCGGGTCAt    >  2:702688/1‑150 (MQ=255)
                                                                                                                                                 tGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGTTTGCGACAGCCATCACCGCGTACGCGGTGCTGTCCGGCAGTTTTGCCGGCGGCCCGTCGGGTCATGa  <  2:928094/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
AAGAAGTGGCGCGCATCGAAGGCGAGCGTATGAAAGCGCTGGCCAAAGCAGGTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGTTTGCGTCGCTCTTCGCCACCTACGCGGTGCTGTCCGGCAGTTTTGCCGGCGGCCCGTCGGGTCATGA  >  NC_002947/953094‑953388

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: