Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,761,360 T→A 6.9% L119Q (CTG→CAG)  gpmI → 2,3‑bisphosphoglycerate‑independent phosphoglycerate mutase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,761,3600TA 6.9% 256.3 / 9.0 87L119Q (CTG→CAG) gpmI2,3‑bisphosphoglycerate‑independent phosphoglycerate mutase
Reads supporting (aligned to +/- strand):  ref base T (40/40);  new base A (3/3);  total (43/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.10e-01

GAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTG‑‑GTGGCCATGGCCGAACTGGCCGCCCAGCGTGGCGCTGAAAAGATCTACCTGCACGCCTTCCTCGATGGCCGCGACACCCCGCCACGCAGCGCGCAGTCGTCCATCGAA  >  NC_002947/5761211‑5761505
                                                                                                                                                     |                                                                                                                                                   
gAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGActgct                                                                                                                                                     >  2:119443/1‑150 (MQ=255)
  aCCTCGGCGCAGGTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGt                                                                                                                                                   <  2:270682/150‑1 (MQ=255)
     tCGGCGCAGGTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCg                                                                                                                                                >  2:430549/1‑150 (MQ=255)
       ggCGCAGGTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGAt                                                                                                                                              <  1:1044717/150‑1 (MQ=255)
       ggCGCAGGTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGAt                                                                                                                                              <  1:1047139/150‑1 (MQ=255)
        gcgcAGGTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGAtg                                                                                                                                             <  1:1608595/150‑1 (MQ=255)
           cAGGTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGAtggtg                                                                                                                                          <  1:584370/150‑1 (MQ=255)
           cAGGTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGAtggtg                                                                                                                                          <  1:1035363/150‑1 (MQ=255)
           cAGGTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAACGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGAtggtg                                                                                                                                          <  1:346544/150‑1 (MQ=255)
             ggTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATggtggt                                                                                                                                        <  1:1215550/150‑1 (MQ=255)
               tCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGgtgt                                                                                                                                      >  1:45550/1‑150 (MQ=255)
                  cgtcgtGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTaca                                                                                                                                   <  1:181521/150‑1 (MQ=255)
                    tcgtGTATCCGGACTTCACGCGGGTGACCAAGGCGATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCGCGGACTGCTGTCCGATGGTGGTGTAcaca                                                                                                                                 >  2:1454986/1‑150 (MQ=255)
                    tcgtGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTAcaca                                                                                                                                 <  1:1080869/150‑1 (MQ=255)
                    tcgtGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTAcaca                                                                                                                                 <  2:45550/150‑1 (MQ=255)
                        gtATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGcca                                                                                                                             >  2:2220/1‑150 (MQ=255)
                           tCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGccacca                                                                                                                          >  1:564229/1‑150 (MQ=255)
                           tCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGccacca                                                                                                                          >  1:902624/1‑150 (MQ=255)
                           tCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGccacca                                                                                                                          >  1:1300977/1‑150 (MQ=255)
                           tCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCATCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGccacca                                                                                                                          >  2:740485/1‑150 (MQ=255)
                             aGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACAGCAGTCCGag                                                                                                                                              <  1:192537/128‑2 (MQ=255)
                             aGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACAGCAGTCCGag                                                                                                                                              >  2:192537/1‑127 (MQ=255)
                               gACTTCACGCGGGTGACCAAGGCCCTCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGAc                                                                                                                      <  2:714424/150‑1 (MQ=255)
                               gACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGccac                                                                                                                            >  1:220794/1‑144 (MQ=255)
                               gACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGccac                                                                                                                            <  2:220794/144‑1 (MQ=255)
                                  ttCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGccaacaattccac                                                                                                                   >  2:1572858/1‑143 (MQ=255)
                                  ttCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCAc                                                                                                                   >  1:690429/1‑150 (MQ=255)
                                    cACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACtt                                                                                                                 <  2:1249214/150‑1 (MQ=255)
                                     aCGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGCCCACTtg                                                                                                                >  2:44500/1‑150 (MQ=255)
                                     aCGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCANTtg                                                                                                                <  2:1194905/150‑1 (MQ=255)
                                     aCGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTtg                                                                                                                <  2:1854319/150‑1 (MQ=255)
                                       gcgGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTtg‑‑gt                                                                                                            >  1:860051/1‑150 (MQ=255)
                                        cgGGTGACCAAGGCCATCAGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGACAAGGCCGTGAACATCCTCGGACTGCGTTCCGATGGTGGTGTACACAGCCACCAGGACCACTtg‑‑gtg                                                                                                           <  2:1006217/150‑1 (MQ=255)
                                               ccAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTG‑‑GTGGCCATgg                                                                                                    >  2:83003/1‑150 (MQ=255)
                                               ccAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTG‑‑GTGGCCATgg                                                                                                    >  2:920169/1‑150 (MQ=255)
                                               ccAAGGCCATCCGCGACGGCGAGCTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCGGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTG‑‑GTGGCGATgg                                                                                                    >  2:945547/1‑150 (MQ=255)
                                                cAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTG‑‑GTGGCCATGGc                                                                                                   <  2:564229/150‑1 (MQ=255)
                                                     ccATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTG‑‑GTGGCCATGGCCGAAc                                                                                              >  1:1216051/1‑150 (MQ=255)
                                                      cATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACAGCAGTCCgcggat                                                                                                                                           <  2:1509213/106‑6 (MQ=255)
                                                      cATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACAGCAGTCCGaggat                                                                                                                                           >  1:1509213/1‑102 (MQ=255)
                                                       aTCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTG‑‑GTGGCCATGGCCGAACTg                                                                                            <  2:690429/150‑1 (MQ=255)
                                                       aTCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACAGCAGTCCGagg                                                                                                                                             <  2:1777453/103‑3 (MQ=255)
                                                       aTCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACAGCAGTCCGagg                                                                                                                                             >  1:1777453/1‑101 (MQ=255)
                                                        tCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTG‑‑GTGGCCATGGCCGAACTgg                                                                                           <  1:1803303/150‑1 (MQ=255)
                                                            cgACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTG‑‑GTGGCCATGGCCGAACTGgccgc                                                                                       <  1:753044/150‑1 (MQ=255)
                                                                 gCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAGGGCAGTCGGTGCTGGCAAGGCCGTACACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTG‑‑GTGGCCATGGCCGAACTGGCCGCCCAGc                                                                                  >  2:1658132/1‑150 (MQ=255)
                                                                   gAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTG‑‑GTGGCCATGGCCGAACTGGCCGCCCAGCGt                                                                                >  2:1594064/1‑150 (MQ=255)
                                                                      ttcttcGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGTAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTG‑‑GTGGCCATGGCCGAACTGGCCGCCCAGCGTGgc                                                                             >  2:423398/1‑150 (MQ=255)
                                                                      ttcttcGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGAtg                                                                                                                                             <  1:419235/88‑1 (MQ=255)
                                                                      ttcttcGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGAtg                                                                                                                                             >  2:419235/1‑88 (MQ=255)
                                                                      ttcttcGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTG‑‑GTGGCCATGGCCGAACTGGCCGCCCAGCGTGgc                                                                             >  2:1470190/1‑150 (MQ=255)
                                                                        cttcGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGAtggtg                                                                                                                                          >  2:302337/1‑89 (MQ=255)
                                                                        cttcGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGAtggtg                                                                                                                                          <  1:302337/89‑1 (MQ=255)
                                                                         ttcGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTG‑‑GTGGCCATGGCCGAACTGGCCGCCCAGCGTGGCGCt                                                                          >  2:1351689/1‑150 (MQ=255)
                                                                         ttcGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCGCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTG‑‑GTGGCCATGGCCGAACTGGCCGCCCAGCGTGGCGCt                                                                          >  2:1530583/1‑150 (MQ=255)
                                                                            gAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTG‑‑GTGGCCATGGCCGAACTGGCCGCCCAGCGTGGCGCTGaa                                                                       <  1:536252/150‑1 (MQ=255)
                                                                              aaaCCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTG‑‑GTGGCCATGGCCGAACTGGCCGCCCAGCGTGGCGCTGaaaa                                                                     <  1:288202/150‑1 (MQ=255)
                                                                                 cccGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTG‑‑GTGGCCATGGCCGAACTGGCCGCCCAGCGTGGCGCTGAAAAGAt                                                                  <  2:1216051/150‑1 (MQ=255)
                                                                                 cccGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTG‑‑GTGGCCATGGCCGAACTGGCCGCCCAGCGTGGCGCTGAAAAGAt                                                                  <  1:920169/150‑1 (MQ=255)
                                                                                  ccGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTG‑‑GTGGCCATGGCCGAACTGGCCGCCCAGCGTGGCGCt                                                                          <  2:1586463/141‑1 (MQ=255)
                                                                                  ccGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTG‑‑GTGGCCATGGCCGAACTGGCCGCCCAGCGTGGCGCt                                                                          >  1:1586463/1‑141 (MQ=255)
                                                                                     gTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTG‑‑GTGGCCATGGCCGAACTGGCCGCCCAGCGTGGCGCTGAAAAGATCTAc                                                              <  1:579778/150‑1 (MQ=255)
                                                                                     gTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTG‑‑GTGGCCATGGCCGAACTGGCCGCCCAGCGTGGCGCTGAAAAGATCTAc                                                              <  2:1300977/150‑1 (MQ=255)
                                                                                     gTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTG‑‑GTGGCCATGGCCGAACTGGCCGCCCAGCGTGGCGCTGAAAAGATCTAc                                                              <  1:119443/150‑1 (MQ=255)
                                                                                          gACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTG‑‑GTGGCCATGGCCGAACTGGCCGCCCAGCGTGGCGCTGAAAAGATCTACCTGCa                                                         >  2:926192/1‑150 (MQ=255)
                                                                                                   tgtgGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTG‑‑GTGGCCATGGCCGAACTGGCCGCCCAGCGTGGCGCTGAAAAGATCTACCTGCACGCCTTCCt                                                <  1:948990/150‑1 (MQ=255)
                                                                                                       gACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTG‑‑GTGGCCATGGCCGAACTGGCCGCCCAGCGTGGCGCTGAAAAGATCTACCTGCACGCCTTCCTCGAt                                            >  1:351176/1‑150 (MQ=255)
                                                                                                         cAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTG‑‑GTGGCCATGGCCGAACTGGCCGCCCAGCGTGGCGCTGAAAAGATCTACCTGCACGCCTTCCTCGATgg                                          <  2:860051/150‑1 (MQ=255)
                                                                                                          aaGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTG‑‑GTGGCCATGGCCGAACTGGCCGCCCAGCGTGGCGCTGAAAAGATCTACCTGCACGCCTTCCTCGATGGc                                         >  1:1914475/1‑150 (MQ=255)
                                                                                                            ggCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCAGCAGGACCACTTG‑‑GTGGCCATGGCCGAACTGGCCGCCCAGCGTGGCGCTGAAAAGATCTACCTGCACGCCTTCCTCGATGGCcg                                       >  2:1051214/1‑150 (MQ=255)
                                                                                                                         ggCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTG‑‑GTGGCCATGGCCGAACTGGCCGCCCAGCGTGGCGCTGAAAAGCTATGCCGGCACGCCTTCCTCGATGGCCGCGACACCCCGCCa                          >  1:1534647/1‑150 (MQ=255)
                                                                                                                          gCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTG‑‑GTGGCCATGGCCGAACTGGCCGCCCAGCGTGGCGCTGAAAAGATCTACCTGCACGCCTTCCTCGATGGCCGCGACACCCCGCCAc                         >  1:68477/1‑150 (MQ=255)
                                                                                                                          gCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTAGTGTACACAGCCACCAGGACCACTTG‑‑GTGGCCATGGCCGAACTGGCCGCCCAGCGTGGCGCTGAAAAGATCTACCTGCACGCCTTCCTCGATGGCCGCGACACCCCGCCAc                         >  2:1500530/1‑150 (MQ=255)
                                                                                                                            aaGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTG‑‑GTGGCCATGGCCGAACTGGCCGCCCAGCGTGGCGCTGAAAAGATCTACCTGCACGCCTTCCTCGATGGCCGCGACACCCCGCCACGc                       >  1:636490/1‑150 (MQ=255)
                                                                                                                               gCCGTGCACCTCCTGGGACTGCTGTCCGACGGAGGTGTACACAGCCACCAGGACCACTTGAAAAAATCAT‑GCTG‑GCTGGCCGCCCAGCGTGGCGCTGAAAAGATCTACCTGCACGCCTTCCTCGATGGCCGCGACACCCCGCCACGCAgc                    <  1:329256/150‑1 (MQ=255)
                                                                                                                                    gCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTG‑‑GTGGCCATGGCCGAACTGGCCGCCCAGCGTGGCGCTGAAAAGATCTACCTGCACGCCTTCCTCGATGGCCGCGACACCCCGCCACGCAGCGCGCa               >  1:1337704/1‑150 (MQ=255)
                                                                                                                                    gCACATCCTCGGACTGCTGTCCGATGGTGGTGTAAACAGCCACCAGGACCACTTG‑‑GTGGCCATGGCCGAACAAGCCGCCCAGCGTGGCGCTGAAAAGATCTACCTGCACGCCTTCCTCGATGGCCGCGACACCCCGCCACGCAGCGCGCa               >  1:1076031/1‑150 (MQ=255)
                                                                                                                                       caTCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTG‑‑GTGGCCATGGCCGAACTGGCCGCCCAGCGTGGCGCTGAAAAGATCTACCTGCACGCCTTCCTCGATGGCCGCGACACCCCGCCACGCAGCGCGCAgtc            <  1:83003/150‑1 (MQ=255)
                                                                                                                                        aTCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTG‑‑GTGGCCATGGCCGAACTGGCCGCCCAGCGTGGCGCTGAAAAGATCTACCTGCACGCCTTCCTCGATGTCCGCGACACCCCGCCACGCAGCGCGCa               <  1:1898497/146‑1 (MQ=255)
                                                                                                                                        aTCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTG‑‑GTGGCCATGGCCGAACTGGCCGCCCAGCGTGGCGCTGAAAAGATCTACCTGCACGCCTTCCTCGATGTCCGAGACACCCCGCCACGCAGCGCGCa               >  2:1898497/1‑146 (MQ=255)
                                                                                                                                        aTCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTG‑‑GTGGCCATGGCCGAACTGGCCGCCCAGCGTGGCGCTGAAAAGATCTACCTGCACGCCTTCCTCGATGGcc                                        >  1:1701975/1‑121 (MQ=255)
                                                                                                                                        aTCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTG‑‑GTGGCCATGGCCGAACTGGCCGCCCAGCGTGGCGCTGAAAAGATCTACCTGCACGCCTTCCTCGATGGcc                                        <  2:1701975/121‑1 (MQ=255)
                                                                                                                                        aTCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTG‑‑GTGGCCATGGCCGAACTGGCCGCCCAGCGTGGCGCTGAAAAGATCTACCTGCACGCCTTCCTCGATGGCCGCGACACCCCGCCACGCAGCGCGCAgtcg           >  2:513359/1‑150 (MQ=255)
                                                                                                                                          ccTCGGACTGCTGTCCGATGGTGGTGTACACAGCCAGCAGGACCACTTG‑‑GTGGCCATGGCCGAACTGGCCGCCCAGCGTGGCGCTGAAAAGATCTACCTGCACGCCTTCCTCGATGGCCGCGACACCCCGCCACGCAGCGCGCAgtcgtc         <  1:1051214/150‑1 (MQ=255)
                                                                                                                                           cTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTG‑‑GTGGCCATGGCCGAACTGGCCGCCCAGCGTGGCGCTGAAAAGATCTACCTGCACGCCTTCCTCGATGGCCGCGACACCCCGCCACGCAGCGCGCAGTCGTcc        >  1:1910097/1‑150 (MQ=255)
                                                                                                                                            tCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTG‑‑GTGGCCATGGCCGAACTGGCCGCCCAGCGTGGCGCTGAAAAGATCTACCTGCACGCCTTCCTCGATGGCCGCGACACCCCGCCACGCAGCGCGCAGTCGTCCa       <  1:1530583/150‑1 (MQ=255)
                                                                                                                                              ggACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTG‑‑GTGGCCATGGCCGAACTGGCCGCCCAGCGTGGCGCTGAAAAGATCTACCTGCACGCCTTCCTCGATGGCCGCGACACCCCGCCACGCAGCGCGCAGTCGTCCATc     <  2:636490/150‑1 (MQ=255)
                                                                                                                                              ggACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTG‑‑GTGGCCATGGCCGAACTGGCCGCCCAGCGTGGCGCTGAAAAGATCTACCTGCACGCCTTCCTCGATGGCCGCGACACCCCGCCACGCAGCGCGCAGTCGTCCATc     <  1:1658132/150‑1 (MQ=255)
                                                                                                                                                 ctgctgTCCGATGGTGGTGTACACAGCCACCAGGACCACTTG‑‑GTGGCCATGGCCGAACTGGCCGCCCAGCGTGGCGCTGAAAAGATCTACCTGCACGCCTTCCTCGATGGCCGCGACACCCCGCCACGCAGCGCGCAGTCGTCCATCGaa  <  1:1821181/150‑1 (MQ=255)
                                                                                                                                                    ctgTCCGATGGTGGTGTACACAGCCACCAGGACCACTTG‑‑GTGGCCATGGCCGAACTGGCCGCCCAGCGTGGCGCTGAAAAGATCTACCTGCACGCCTTCCTCGa                                             >  2:395010/1‑104 (MQ=255)
                                                                                                                                                    ctgTCCGATGGTGGTGTACACAGCCACCAGGACCACTTG‑‑GTGGCCATGGCCGAACTGGCCGCCCAGCGTGGCGCTGAAAAGATCTACCTGCACGCCTTCCTCGa                                             <  1:395010/104‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTG‑‑GTGGCCATGGCCGAACTGGCCGCCCAGCGTGGCGCTGAAAAGATCTACCTGCACGCCTTCCTCGATGGCCGCGACACCCCGCCACGCAGCGCGCAGTCGTCCATCGAA  >  NC_002947/5761211‑5761505

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