Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,274,242 G→A 5.5% G18D (GGT→GAT) ‡ metZ → O‑succinylhomoserine sulfhydrylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,274,2420GA 5.5% 869.7 / 24.3 271G18D (GGT→GAT) ‡metZO‑succinylhomoserine sulfhydrylase
Reads supporting (aligned to +/- strand):  ref base G (131/125);  new base A (4/11);  total (135/136)
Fisher's exact test for biased strand distribution p-value = 1.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 275 total aligned reads displayed.

CGCATCGAGCAGGCCCGTAACGACTTGGCCAAGGTCAAGAACCAGGCCGTCAGCGCGATCATCGACCTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTATGTGTTCCGCACGGCAGCCGATGCTGCTGCGCGCTTTGCCGGCGAAACGCCGGGCAACGTCTACT  >  NC_002947/2274093‑2274391
                                                                                                                                                     |                                                                                                                                                     
   aTCGAGCAGGCCCGTAACGACTTGGCCAAGGTCAAGAACCAGGCCGTCAGCGCGATCATCGACCTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGttt                                                                                                                                                    <  2:4310709/150‑1 (MQ=255)
   aTCGAGCAGGCCCGTAACGACTTGGCCAAGGTCAAGAACCAGGCCGTCAGCGCGATCATCGACCTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGttt                                                                                                                                                    >  1:702970/1‑150 (MQ=255)
    tCGAGCAGGCCCGTAACGACTTGGCCAAGGTCAAGAACCAGGCCGTCAGCGCGATCATCGACCTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTc                                                                                                                                                   <  1:739565/150‑1 (MQ=255)
    tCGAGCAGGCCCGTAACGACTTGGCCAAGGTCAAGAACCAGGCCGTCAGCGCGATCATCGACCTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTc                                                                                                                                                   <  2:4676852/150‑1 (MQ=255)
    acGAGCAGGCCCGTAACGACTTGGCCAAGGTCAAGAACCAGGCCGTCAGCGCGATCATCGACCTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTc                                                                                                                                                   <  2:3356019/149‑1 (MQ=255)
       aGCAGGCCCGTAACGACTTGGCCAAGGTCAAGAACCAGGCCGTCAGCGCGATCATCGACCTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGac                                                                                                                                                >  2:3593990/1‑150 (MQ=255)
        gCAGGCCCGTAACGACTTGGCCAAGGTCAAGAACCAGGCCGTCAGCGCGATCATCGACCTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGaca                                                                                                                                               <  2:403871/150‑1 (MQ=255)
            gCCCGTAACGACTTGGCCAAGGTCAAGAACCAGGCCGTCAGCGCGATCATCGACCTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGAAACCGACACCCt                                                                                                                                           <  2:2072501/150‑1 (MQ=255)
             cccGTAACGACTTGGCCAAGGTCAAGAACCAGGCCGTCAGCGCGATCATCGACCTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGAAACCGACACCCtc                                                                                                                                          <  2:1649541/150‑2 (MQ=255)
              ccGTAACGACTTGGCCAAGGTCAAGAACCAGGCCGTCAGCGCGATCATCGACCTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTgg                                                                                                                                         >  2:2934148/1‑150 (MQ=255)
              ccGTAACGACTTGGCCAAGGTCAAGAACCAGGCCGTCAGCGCGATCATCGACCTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTgg                                                                                                                                         >  1:4371433/1‑150 (MQ=255)
                gTAACGACTTGGCCAAGGTCAAGAACCAGGCCGTCAGCGCGATCATCGACCTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCGCCAGGGTGTCGAAACCGACACCCtcgag                                                                                                                                       <  2:551552/150‑5 (MQ=255)
                    cGACTTGGCCAAGGTCAAGAACCAGGCCTTCAGCGCGATAATCGACCTTTACAACAACTGATTGTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTAc                                                                                                                                   >  2:5312176/1‑150 (MQ=255)
                          ggCCAAGGTCAAGAACCAGGCCGTCAGCGCGATCATCGACCTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTg                                                                                                                             <  2:4835963/150‑1 (MQ=255)
                          ggCCAAGGTCAAGAACCAGGCCGTCAGCGCGATCATCGACCTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTg                                                                                                                             <  1:1947004/150‑1 (MQ=255)
                          ggCCAAGGTCAAGAACCAGGCCGTCAGCGCGATCATCGACCTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTg                                                                                                                             <  2:5135339/150‑1 (MQ=255)
                          ggCCAAGGTCAAGAACCAGGCCGTCAGCGCGATCATCGACCTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTg                                                                                                                             <  1:1259293/150‑1 (MQ=255)
                            ccAAGGTCAAGAACCAGGCCGTCAGCGCGATCATCGACCTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGc                                                                                                                           <  1:3328004/150‑1 (MQ=255)
                             cAAGGTCAAGAACCAGGCCGTCAGCGCGATCATCGACCTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGcc                                                                                                                          >  1:4185486/1‑150 (MQ=255)
                             cAAGGTCAAGAACCAGGCCGTCAGCGCGATCATCGACCTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGcc                                                                                                                          >  1:2755581/1‑150 (MQ=255)
                             cAAGGTCAAGAACCAGGCCGTCAGCGCGATCATCGACCTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGcc                                                                                                                          <  2:3558535/150‑1 (MQ=255)
                              gaGGTCAAGAACCAGGCCGTCAGCGCGATCATCGACCTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCa                                                                                                                         >  1:4861065/2‑150 (MQ=255)
                                 gTCAAGAACCAGGCCGTCAGCGCGATCATCGACCTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCaaaa                                                                                                                      <  2:2357912/150‑1 (MQ=255)
                                   cAAGAACCAGGCCGTCAGCGCGATCATCGACCTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAAcc                                                                                                                    <  1:1337944/150‑1 (MQ=255)
                                    aaGAACCAGGCCGTCAGCGCGATCATCGACCTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCg                                                                                                                   >  1:5281755/1‑150 (MQ=255)
                                    aaGAACCAGGCCGTCAGCGCGATCATCGACCTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCg                                                                                                                   >  2:2597904/1‑150 (MQ=255)
                                    aaGAACCAGGCCGTCAGCGCGATCATCGACCTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGGTTCGACACCCTGGCGGTACGCGCTGGCCACAACCg                                                                                                                   >  2:3507290/1‑150 (MQ=255)
                                      gAACCAGGCCGTCAGCGCGATCATCGACCTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTa                                                                                                                 >  2:5426626/1‑150 (MQ=255)
                                      gAACCAGGCCGTCAGCGCGATCATCGACCTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTa                                                                                                                 >  1:1274082/1‑150 (MQ=255)
                                       aaCCAGGCCGTCAGCGCGATCATCGACCTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTAc                                                                                                                <  2:275803/150‑1 (MQ=255)
                                         ccAGGCCGTCAGCGCGATCATCGACCTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTAccc                                                                                                              >  2:5001573/1‑150 (MQ=255)
                                           aGGCCGTCAGCGCGATCATCGACCTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCg                                                                                                            >  2:630988/1‑150 (MQ=255)
                                                   aGCGCGATCATCTACCTTTACAACAACTGATTTTGGGAGCGTCTGCATGACGGATCATTGGGTTGCCGGGCGACTGGTCAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCACGGAAGGCGa                                                                                                    <  2:2563144/150‑1 (MQ=255)
                                                   aGCGCGATCATCGACCTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGa                                                                                                    <  2:2753871/150‑1 (MQ=255)
                                                   aGCGCGATCATCGACCTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGa                                                                                                    <  1:5011782/150‑1 (MQ=255)
                                                    gcgcGATCATCGACCTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAg                                                                                                   <  1:858690/150‑1 (MQ=255)
                                                      gcgATCATCGACCTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGca                                                                                                 <  1:1208895/150‑1 (MQ=255)
                                                        gATCATCGACCTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGcaca                                                                                               >  2:3100605/1‑150 (MQ=255)
                                                        gATCATCGACCTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGcaca                                                                                               >  1:4217831/1‑150 (MQ=255)
                                                        gATCATCGACCTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGcaca                                                                                               <  2:1245637/150‑1 (MQ=255)
                                                             tcGACCTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGaa                                                                                          >  1:5451450/1‑150 (MQ=255)
                                                             tcGACCTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGaa                                                                                          >  2:3156673/1‑150 (MQ=255)
                                                              cGACCTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAg                                                                                         >  2:36144/1‑150 (MQ=255)
                                                                   tttACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTg                                                                                    <  1:171742/150‑1 (MQ=255)
                                                                      acaacaACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTc                                                                                 <  1:2668129/150‑1 (MQ=255)
                                                                      acaacaACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTc                                                                                 <  1:3917127/150‑1 (MQ=255)
                                                                        caccccTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACCGCGAAGCGCTGTTCCt                                                                               <  2:3808333/145‑1 (MQ=255)
                                                                        aacaacTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCt                                                                               <  2:2500003/150‑1 (MQ=255)
                                                                        aacaacTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCt                                                                               >  2:1667371/1‑150 (MQ=255)
                                                                           aacTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAGGGCGAGCACAGCGAAGCGCTGTTCCTGAc                                                                            >  1:1606671/1‑150 (MQ=255)
                                                                           aacTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGAc                                                                            >  1:4384180/1‑150 (MQ=255)
                                                                            acTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGAcc                                                                           >  1:2376255/1‑150 (MQ=255)
                                                                              tGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGt                                                                                   >  1:1073484/1‑140 (MQ=255)
                                                                              tGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGt                                                                                   <  2:1073484/140‑1 (MQ=255)
                                                                              tGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTc                                                                         >  1:2921032/1‑150 (MQ=255)
                                                                              tGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTc                                                                         >  2:4272871/1‑150 (MQ=255)
                                                                              tGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTc                                                                         >  1:3426930/1‑150 (MQ=255)
                                                                               gATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTcc                                                                        <  2:4247462/150‑1 (MQ=255)
                                                                                    tGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTa                                                                   <  2:278597/150‑1 (MQ=255)
                                                                                     gggAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTAt                                                                  <  2:1108826/150‑1 (MQ=255)
                                                                                       gAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTAtgt                                                                >  2:3664947/1‑150 (MQ=255)
                                                                                          cGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTATGTGtt                                                             >  2:391378/1‑150 (MQ=255)
                                                                                           gACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTATGTGTTc                                                            <  1:4182985/150‑1 (MQ=255)
                                                                                                cATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTATGTGTTCCGCAc                                                       <  2:3856236/150‑1 (MQ=255)
                                                                                                  tGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTAGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTATGTGTTCCGCACgg                                                     >  2:4584814/1‑150 (MQ=255)
                                                                                                   gACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTATGTGTTCCGCACGGc                                                    <  1:5113991/150‑1 (MQ=255)
                                                                                                    aCGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTATGTGTTCCGCACGGCa                                                   <  1:490456/150‑1 (MQ=255)
                                                                                                     cGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTATGTGTTCCGCACGGCAg                                                  <  1:3096666/150‑1 (MQ=255)
                                                                                                       gATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTATGTGTTCCGCACGGCAGcc                                                <  1:459861/150‑1 (MQ=255)
                                                                                                            aTGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTAt                                                                  >  1:908808/1‑127 (MQ=255)
                                                                                                            aTGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTATGTGTTCCGCACGGCAGCCGAtgc                                           >  2:422262/1‑150 (MQ=255)
                                                                                                            aTGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTATGTGTTCCGCACGGCAGCCGAtgc                                           <  2:794434/150‑1 (MQ=255)
                                                                                                              gggATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTATGTGTTCCGCACGGCAGCCGAtgctg                                         <  1:5568882/150‑1 (MQ=255)
                                                                                                                  tGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTATGTGTTCCGCACGGCAGCCGAtgctgctgc                                     <  1:5601433/150‑1 (MQ=255)
                                                                                                                   gCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTATGTGTTCCGCACGGCAGCCGATGCTGCTgcg                                    <  2:522217/150‑1 (MQ=255)
                                                                                                                      gggCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTATGTGTTCCGCACGGCAGCCGATGCTGCTgcgcgc                                 >  2:4106077/1‑150 (MQ=255)
                                                                                                                          gACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTATGTGTTCCGCACGGCAGCCGATGCTGCTGCGCGCTTTg                             >  2:1233189/1‑150 (MQ=255)
                                                                                                                          gACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTATGAGTTCCGCACGGCAGCCGATGCTGCTGCGCGCTTTg                             >  1:3532076/1‑150 (MQ=255)
                                                                                                                           aCTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTATGTGTTCCGCACGGCAGCCGATGCTGCTGCGCGCTTTGc                            <  2:1282236/150‑1 (MQ=255)
                                                                                                                               gACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTATGTGTTCCGCACGGCAGCCGATGCTGCTGCGCGCTTTGCCGGc                        >  1:139747/1‑150 (MQ=255)
                                                                                                                                 cAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTATGTGTTCCGCACGGCAGCCGATGCTGCTGCGCGCTTTGCCGGCGa                      <  1:370411/150‑1 (MQ=255)
                                                                                                                                        cgccaGGGTGTCGAAACCGACACCCTCGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTATGTGTTCCGCACGGCAGCCGATGCTGCTGCGCGCTTTGCCGGCGAAACGCCg               >  1:25889/5‑150 (MQ=255)
                                                                                                                                        cTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTATGTGTTCCGCACGGCAGCCGATGCTGCTGCGCGCTTTGCCGGCGAAACGCCg               <  1:3338236/150‑1 (MQ=255)
                                                                                                                                        cTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTATGTGTTCCGCACGGCAGCCGATGCTGCTGCGCGCTTTGCCGGCGAAACGCCg               <  1:5316771/150‑1 (MQ=255)
                                                                                                                                        cTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTATGTGTTCCGCACGGCAGCCGATGCTGCTGCGCGCTTTGCCGGCGAAACGCCg               <  2:489958/150‑1 (MQ=255)
                                                                                                                                         tCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTATGTGTTCCGCACGGCAGCCGATGCTGCTGCGCGCTTTGCCGGCGAAACGCCgg              <  1:1113736/150‑1 (MQ=255)
                                                                                                                                         tCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTATGTGTTCCGCACGGCAGCCGATGCTGCTGCGCGCTTTGCCGGCGAAACGCCgg              <  2:1851223/150‑1 (MQ=255)
                                                                                                                                         tCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTATGTGTTCCGCACGGCAGCCGATGCTGCTGCGCGCTTTGCCGGCGAAACGCCgg              <  1:2493716/150‑1 (MQ=255)
                                                                                                                                         gccaGGGTGTCGAAACCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTATGTGTTCCGCACGGCAGCCGATGCTGCTGCGCGCTTTGCCGGCGAAACGCCgg              >  2:3763324/4‑150 (MQ=255)
                                                                                                                                         gccaGGGTGTCGAAACCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTATGTGTTCCGCACGGCAGCCGATGCTGCTGCGCGCTTTGCCGGCGAAACGCCgg              >  1:4006283/4‑150 (MQ=255)
                                                                                                                                          cGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTATGTGTTCCGCACGGCAGCCGATGCTGCTGCGCGCTTTGCCGGCGAAACGCCggg             <  2:903318/150‑1 (MQ=255)
                                                                                                                                           gAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTATGTGTTCCGCACGGCAGCCGATGCTGCTGCGCGCTTTGCCGGCGAAACGCCGGGc            >  1:1756347/1‑150 (MQ=255)
                                                                                                                                            aGGGTGTCGAAACCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTATGTGTTCCGCACGGCAGCCGATGCTGCTGCGCGCTTTGCCGGCGAAACGCCGGGCa           >  1:3053819/1‑150 (MQ=255)
                                                                                                                                             gggTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTATGTGTTCCGCACGGCAGCCGATGCTGCTGCGCGCTTTGCCGGCGTTTCGCCggcaaa          <  2:1513353/150‑5 (MQ=255)
                                                                                                                                              ggTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTATGTGTTCCGCACGGCAGCCGATGCTGCTGCGCGCTTTGCCGGCGAAACGCCGGGCAAc         >  1:5132533/1‑150 (MQ=255)
                                                                                                                                              ggTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTATGTGTTCCGCACGGCAGCCGATGCTGCTGCGCGCTTTGCCGGCGAAACGCCGGGCAAc         <  1:2887902/150‑1 (MQ=255)
                                                                                                                                              ggTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTATGTGTTCCGCACGGCAGCCGATGCTGCTGCGCGCTTTACCGGCGAAACGCCGGGCAAc         <  1:3655751/150‑1 (MQ=255)
                                                                                                                                                tgtCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTATGTGTTCCGCACGGCAGCCGATGCTGCTGCGCGCTTTGCCGGCGAAACGCCGGGCAACGt       <  2:5550118/150‑1 (MQ=255)
                                                                                                                                                 gtCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTATGTGTTCCGCACGGCAGCCGATGCTGCTGCGCGCTTTGCCGGCGAAACGCCGGGCAACGTc      >  1:381479/1‑150 (MQ=255)
                                                                                                                                                    ggTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTATGTGTTCCGCACGGCAGCCGATGCTGCTGCGCGCTTTGCCGGCGAAACGCCGGGCAACGTCTAc   <  1:3777186/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGCATCGAGCAGGCCCGTAACGACTTGGCCAAGGTCAAGAACCAGGCCGTCAGCGCGATCATCGACCTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTATGTGTTCCGCACGGCAGCCGATGCTGCTGCGCGCTTTGCCGGCGAAACGCCGGGCAACGTCTACT  >  NC_002947/2274093‑2274391

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: