Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 953,232 A→G 5.3% H13R (CAC→CGC) ‡ 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,2320AG 5.3% 2369.4 / 70.8 731H13R (CAC→CGC) ‡cyoCcytochrome bo terminal oxidase subunit III
Reads supporting (aligned to +/- strand):  ref base A (354/336);  new base G (19/20);  total (374/357)
Fisher's exact test for biased strand distribution p-value = 8.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.96e-01

Only 100 of 752 total aligned reads displayed.

GGTTCCGGCCGAAGAAGTGGCGCGCATCGAAGGCGAGCGTATGAAAGCGCTGGCCAAAGCAGGTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGTTTGCGTCGCTCTTCGCCACCTACGCGGTGCTGTCCGGCAGTTTTGCCGGCGGCCCGTCGG  >  NC_002947/953083‑953381
                                                                                                                                                     |                                                                                                                                                     
 gTTCCGGCCGAAGAAGTGGCGCGCATCGAAGGCGAGCGTATGAAAGCGCTGGCCAAAGCAGGTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCAc                                                                                                                                                      <  2:717442/150‑1 (MQ=255)
   tCCGGCCGAAGAAGTGGCGCGCATCGAAGGCGAGCGTATGAAAGCGCTGGCCAAAGCAGGTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGa                                                                                                                                                    <  1:6803778/150‑1 (MQ=255)
       gCCGAAGAAGTGGCGCGCATCGAAGGCGAGCGTATGAAAGCGCTGGCCAAAGCAGGTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCAt                                                                                                                                                <  1:5765281/150‑1 (MQ=255)
        ccGAAGAAGTGGCGCGCATCGAAGGCGAGCGTATGAAAGCGCTGGCCAAAGCAGGTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATg                                                                                                                                               <  1:1168934/150‑1 (MQ=255)
         cGAAGAAGTGGCGCGCATCGAAGGCGAGCGTATGAAAGCGCTGGCCAAAGCAGGTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATgg                                                                                                                                              <  1:5356590/150‑1 (MQ=255)
         cGAAGAAGTGGCGCGCATCGAAGGCGAGCGTATGAAAGCGCTGGCCAAAGCAGGTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATgg                                                                                                                                              >  2:549255/1‑150 (MQ=255)
           aagaagTGGCGCGCATCGAAGGCGAGCGTATGAAAGCGCTGGCCAAAGCAGGAGCTCAGCCAGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCAAGGCGCTGCTCATGGTCACGACCAAGGGc                                                                                                                                            >  1:18310/1‑150 (MQ=255)
            agaagTGGCGCGCATCGAAGGCGAGCGTATGAAAGCGCTGGCCAAAGCAGGTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCa                                                                                                                                           <  1:671725/150‑1 (MQ=255)
               agTGGCGCGCATCGAAGGCGAGCGTATGAAAGCGCTGGCCAAAGCAGGTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATga                                                                                                                                        <  1:12829/150‑1 (MQ=255)
                 tGGCGCGCATCGAAGGCGAGCGTATGAAAGCGCTGGCCAAAGCAGGTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATgacg                                                                                                                                      >  1:2844280/1‑150 (MQ=255)
                  ggCGCGCATCGAAGGCGAGCGTATGAAAGCGCTGGCCAAAGCAGGTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATgacga                                                                                                                                     <  2:283267/150‑1 (MQ=255)
                   gcgcgcATCGAAGGCGAGCGTATGAAAGCGCTGGCCAAAGCAGGTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATgacgac                                                                                                                                    >  2:2685095/1‑150 (MQ=255)
                    cgcgcATCGAAGGCGAGCGTATGAAAGCGCTGGCCAAAGCAGGTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGacc                                                                                                                                   <  1:786140/150‑1 (MQ=255)
                     gcgcATCGAAGGCGAGCGTATGAAAGCGCTGGCCAAAGCAGGTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCCCATGGTCGTGACCATgag                                                                                                                                             <  2:1501778/139‑3 (MQ=255)
                     gcgcATCGAAGGCGAGCGTATGAAAGCGCTGGCCAAAGCAGGTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCCCATGGTCGTGACCATgag                                                                                                                                             >  1:1501778/1‑137 (MQ=255)
                      cgcATCGAAGGCGAGCGTATGAAAGCGCTGGCCAAAGCAGGTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGaccac                                                                                                                                 <  2:2334860/150‑1 (MQ=255)
                      cgcATCGAAGGCGAGCGTATGAAAGCGCTGGCCAAAGCAGGTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGaccac                                                                                                                                 <  1:4785351/150‑1 (MQ=255)
                         aTCGAAGGCGAGCGTATGAAAGCGCTGGCCAAAGCAGGTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGaccaccac                                                                                                                              >  1:1573141/1‑150 (MQ=255)
                         aTCGAAGGCGAGCGTATGAAAGCGCTGGCCAAAGCAGGTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGaccaccac                                                                                                                              >  2:4929700/1‑150 (MQ=255)
                           cGAAGGCGAGCGTATGAAAGCGCTGGCCAAAGCAGGTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGAccaccacca                                                                                                                            >  1:3342766/1‑150 (MQ=255)
                           cGAAGGCGAGCGTATGAAAGCGCTGGCCAAAGCAGGTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGAccaccacca                                                                                                                            <  1:5986700/150‑1 (MQ=255)
                            gAAGGCGAGCGTATGAAAGCGCTGGCCAAAGCAGGTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACcaccaccac                                                                                                                           >  2:1728401/1‑150 (MQ=255)
                                   aGCGTATGAAAGCGCTGGCCAAAGCAGGTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCgg                                                                                                                    <  2:3607852/150‑1 (MQ=255)
                                   aGCGTATGAAAGCGCTGGCCAAAGCAGGTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCgg                                                                                                                    <  2:8453276/150‑1 (MQ=255)
                                   aGCGTATGAAAGCGCTGGCCAAAGCAGGTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCgg                                                                                                                    <  2:6418780/150‑1 (MQ=255)
                                    gCGTATGAAAGCGCTGGCCAAAGCAGGTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCggg                                                                                                                   >  1:1038914/1‑150 (MQ=255)
                                         tGAAAGCGCTGGCCAAAGCAGGTGCTCTGCCTGCCGGCGCAAGAGACGAAAAGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGAACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGa                                                                                                              >  1:166395/1‑150 (MQ=255)
                                           aaaGCGCTGGCCAAAGCAGGTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATg                                                                                                            <  1:1728401/150‑1 (MQ=255)
                                            aaGCGCTGGCCAAAGCAGGTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTTTAATCAATGTCCAGTCTAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGTCATGACGACCACCACCACGACTCGGGCCAGATGa                                                                                                           <  1:2635736/150‑1 (MQ=255)
                                            aaGCGCTGGCCAAAGCAGGTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGa                                                                                                           >  2:6550411/1‑150 (MQ=255)
                                               cgcTGGCCAAAGCAGGTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCg                                                                                                        <  2:6982494/150‑1 (MQ=255)
                                                gcTGGCCAAAGCAGGTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGt                                                                                                       >  2:3133594/1‑150 (MQ=255)
                                                  tGGCCAAAGCAGGTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTAc                                                                                                     >  2:477643/1‑150 (MQ=255)
                                                   ggCCAACGCAGGTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACt                                                                                                    <  1:4414613/150‑1 (MQ=255)
                                                    gCCAAAGCAGGTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCCCGACCATGGGCAGGACGACCACCACCACGACTCGGGCCAGATGACCGTAcgg                                                                                                   >  1:5542674/1‑148 (MQ=255)
                                                    gCCAAAGCAGGTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAAACAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTg                                                                                                   >  1:6532943/1‑150 (MQ=255)
                                                     ccAAAGCAGGTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTgg                                                                                                  <  2:6638357/150‑1 (MQ=255)
                                                       aaaGCAGGTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGc                                                                                                >  1:2072729/1‑150 (MQ=255)
                                                       aaaGCAGGTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGc                                                                                                >  2:4096783/1‑150 (MQ=255)
                                                       aaaGCAGGTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGc                                                                                                >  2:2191364/1‑150 (MQ=255)
                                                          gCAGGTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTc                                                                                             >  1:2897117/1‑150 (MQ=255)
                                                            aGGTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTg                                                                                           <  1:3816120/150‑1 (MQ=255)
                                                            aGGTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTg                                                                                           <  2:8310615/150‑1 (MQ=255)
                                                             ggTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTgg                                                                                          <  2:4226627/150‑1 (MQ=255)
                                                             ggTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTgg                                                                                          >  2:5271782/1‑150 (MQ=255)
                                                             ggTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTgg                                                                                          <  2:755367/150‑1 (MQ=255)
                                                              gTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGc                                                                                         <  2:826867/150‑1 (MQ=255)
                                                                        tgccGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCtgat                                                                               <  1:4212319/150‑1 (MQ=255)
                                                                        tgccGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCtgat                                                                               >  1:5821924/1‑150 (MQ=255)
                                                                        tgccGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAAGGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCtgat                                                                               <  2:7872755/150‑1 (MQ=255)
                                                                                   gtgtCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCa                                                                    <  1:4694215/150‑1 (MQ=255)
                                                                                      tCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATcc                                                                 >  1:8461803/1‑150 (MQ=255)
                                                                                      tCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATcc                                                                 >  2:3394892/1‑150 (MQ=255)
                                                                                        gAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTg                                                               >  1:8601686/1‑150 (MQ=255)
                                                                                         aaTCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGt                                                              >  2:2011096/1‑150 (MQ=255)
                                                                                         aaTCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGt                                                              >  2:611861/1‑150 (MQ=255)
                                                                                          aTCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATTACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGtt                                                             <  2:8106648/150‑1 (MQ=255)
                                                                                           tCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGttt                                                            <  2:7742735/150‑1 (MQ=255)
                                                                                              tttGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGTTTGCg                                                         >  2:3817810/1‑150 (MQ=255)
                                                                                                   gCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGTTTGCGTCGct                                                    >  1:993707/1‑150 (MQ=255)
                                                                                                     gggTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGTTTGCGTCGctct                                                  >  2:4544813/1‑150 (MQ=255)
                                                                                                      ggTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGTTTGCGTCGCTCtt                                                 <  1:8079767/150‑1 (MQ=255)
                                                                                                       gtgtAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGTTTGCGTCGCTCTTc                                                <  1:635035/150‑1 (MQ=255)
                                                                                                       gtgtAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGTTTGCGTCGCTCTTc                                                <  2:5701776/150‑1 (MQ=255)
                                                                                                          tAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGTTTGCGTCGCTCTTCGcc                                             >  2:5368034/1‑150 (MQ=255)
                                                                                                           aaTCAAGGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGTTTGCGTCGCTCTTCGCCa                                            <  2:6004605/150‑1 (MQ=255)
                                                                                                            aTCAATGTCCAGTCTAGTAATGCACGGCGGTGCTCATGGTCTCGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGTCCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGTTTGCGTCGCTCTTCGCCAc                                           <  2:7308243/150‑1 (MQ=255)
                                                                                                            aTCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGTTTGCGTCGCTCTTCGCCAc                                           >  1:701955/1‑150 (MQ=255)
                                                                                                             tCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGTTTGCGTCGCTCTTCGCCAcc                                          >  2:8584154/1‑150 (MQ=255)
                                                                                                             tCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGTTTGCGTCGCTCTTCGCCAcc                                          >  2:247934/1‑150 (MQ=255)
                                                                                                               aaTGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGTTTGCGTCGCTCTTCGCCACCTa                                        >  1:7671119/1‑150 (MQ=255)
                                                                                                                aTGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGTTTGCGTCGCTCTTCGCCACCTAc                                       <  2:6306207/150‑1 (MQ=255)
                                                                                                                aTGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGTTTGCGTCGCTCTTCGCCACCTAc                                       <  2:1175204/150‑1 (MQ=255)
                                                                                                                   tCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGTTTGCGTCGCTCTTCGCCACCTAcgcg                                    >  1:4620839/1‑150 (MQ=255)
                                                                                                                      aGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGTTTGCGTCGCTCTTCGCCACCTACGCGGTg                                 <  1:6147387/150‑1 (MQ=255)
                                                                                                                      aGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGTTTGCGTCGCTCTTCGCCACCTACGCGGTg                                 <  1:3099883/150‑1 (MQ=255)
                                                                                                                      aGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGTTTGCGTCGCTCTTCGCCACCTACGCGGTg                                 >  1:8385531/1‑150 (MQ=255)
                                                                                                                        tCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGTTTGCGTCGCTCTTCGCCACCTACGCGGTGCt                               >  1:3712294/1‑150 (MQ=255)
                                                                                                                           aGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGTTTGCGTCGCTCTTCGCCACCTACGCGGTGCTGTc                            >  1:402568/1‑150 (MQ=255)
                                                                                                                            gTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGTTTGCGTCGCTCTTCGCCACCTACGCGGTGCTGTcc                           >  1:7565729/1‑150 (MQ=255)
                                                                                                                             tAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGTTTGCGTCGCTCTTCGCCACCTACGCGGTGCTGTCCg                          >  2:4502999/1‑150 (MQ=255)
                                                                                                                             tAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGTTTGCGTCGCTCTTCGCCACCTACGCGGTGCTGTCCg                          >  1:867628/1‑150 (MQ=255)
                                                                                                                                tGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGTTTGCGTCGCTCTTCGCCACCTACGCGGTGCTGTCCGGCa                       <  2:5160574/150‑1 (MQ=255)
                                                                                                                                    cGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGTTTGCGTCGCTCTTCGCCACCTACGCGGTGCTGTCCGGCAGttt                   <  1:7286387/150‑1 (MQ=255)
                                                                                                                                     gtagctgctCATGGTGAAGACCTTTGGCATGACGAACTCCACGACGACTCGTTCGTGATGACGGTACTTGGCTTCTGGCTGTACCTGATGACCGTCTTCATCCTGTTTGCGTCTCTCTTCGCCACCTACGCTGTGCTGTCCGTCAGtttt                  <  2:3273090/147‑1 (MQ=255)
                                                                                                                                        gctgctCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGTTTGCGTCGCTCTTCGCCACCTACGCGGTGCTGTCCGGCAGTTTTGcc               <  2:6783078/150‑1 (MQ=255)
                                                                                                                                        gctgctCATGGTCACGACCATGGGCATGACGACCAACACCACGACTCGGGCCAGATGACCGTACTGGTCTTCTGGCTGTACCTGATGACCGACTGCATCCTGTTTGCGTCGCTCTTCGCCACCTACGCGGTGCTGTCCGGCAGTTTTGcc               <  1:4528685/150‑1 (MQ=255)
                                                                                                                                        catgccCATGGTCGTGACCATGAGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGTTTGCGTCGCTCTTCGCCACCTACGCGGTGCTGTCCGGCAGTTTTGcc               >  1:5042873/3‑150 (MQ=255)
                                                                                                                                        catgccCATGGTCGTGACCATGAGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGTTTGCGTCGCTCTTCGCCACCTACGCGGTGCTGTCCGGCAGTTTTGcc               >  1:648887/3‑150 (MQ=255)
                                                                                                                                          tgccCATGGTCGTGACCATGAGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGTTTGCGTCGCTCTTCGCCACCTACGCGGTGCTGTCCGGCAGTTTTGCcgg             >  1:2738352/1‑150 (MQ=255)
                                                                                                                                            ctCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGTTTGCGTCGCTCTTCGCCACCTACGCGGTGCTGTCCGGCAGTTTTGCcggcg           >  1:3555954/1‑150 (MQ=255)
                                                                                                                                             tCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGTTTGCGTCGCTCTTCGCCACCTACGCGGTGCTGTCCGGCAGTTTTGCcggcgg          <  1:5151826/150‑1 (MQ=255)
                                                                                                                                               aTGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGTTTGCGTCGCTCTTCGCCACCTACGCGGTGCTGTCCGGCAGTTTTGCCGGCGGcc        >  2:3990051/1‑150 (MQ=255)
                                                                                                                                                tGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGTTTGCGTCGCTCTTCGCCACCTACGCGGTGCTGTCCGGCAGTTTTGCCGGCGGccc       >  1:3220426/1‑150 (MQ=255)
                                                                                                                                                 ggTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGTTTGCGTCGCTCTTCGCCACCTACGCGGTGCTGTCCGGCAGTTTTGCCGGCGGCCCg      >  2:1426079/1‑150 (MQ=255)
                                                                                                                                                  gTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGAc                                                                        <  1:4348979/83‑1 (MQ=255)
                                                                                                                                                   tCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGTTTGCGTCGCTCTTCGCCACCTACGCGGTGCTGTCCGGCAGTTTTGCCGGCGGCCCGTc    >  1:2365486/1‑150 (MQ=255)
                                                                                                                                                     aCGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGTTTGCGTCGCTCTTCGCCACCTACGCGGTGCTGTCCGGCAGTTTTGCCGGCGGCCCGTCgg  <  1:8245497/150‑1 (MQ=255)
                                                                                                                                                     aCGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGTTTGCGTCGCTCTTCGCCACCTACGCGGTGCTGTCCGGCAGTTTTGCCGGCGGCCCGTCgg  >  2:1351101/1‑150 (MQ=255)
                                                                                                                                                     aCGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGTTTGCGTCGCTCTTCGCCACCTACGCGGTGCTGTCCGGCAGTTTTGCCGGCGGCCCGTCgg  <  2:8739017/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGTTCCGGCCGAAGAAGTGGCGCGCATCGAAGGCGAGCGTATGAAAGCGCTGGCCAAAGCAGGTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGTTTGCGTCGCTCTTCGCCACCTACGCGGTGCTGTCCGGCAGTTTTGCCGGCGGCCCGTCGG  >  NC_002947/953083‑953381

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