Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,916,102 C→T 8.0% intergenic (+16/+137) spuI → / ← argA glutamylpolyamine synthetase/amino acid N‑acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,916,1020CT 8.0% 406.9 / 17.3 137intergenic (+16/+137)spuI/argAglutamylpolyamine synthetase/amino acid N‑acetyltransferase
Reads supporting (aligned to +/- strand):  ref base C (47/79);  new base T (4/7);  total (51/86)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 140 total aligned reads displayed.

CAGCGAGATCCTGAACAAGTACATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGGCGGCCCAGGCTTACAGGGGCTTC  >  NC_002947/5915954‑5916251
                                                                                                                                                    |                                                                                                                                                     
cAGCGAGATCCTGAACAAGTACATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCt                                                                                                                                                      >  1:4197256/1‑150 (MQ=255)
 aGCGAGATCCTGAACAAGTACATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTc                                                                                                                                                     >  1:1021337/1‑150 (MQ=255)
   cGAGATCCTGAACAAGTACATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCaa                                                                                                                                                   <  2:5316985/150‑1 (MQ=255)
    gagaTCCTGAACAAGTACATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAg                                                                                                                                                  <  2:1448289/150‑1 (MQ=255)
    gagaTCCTGAACAAGTACATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAg                                                                                                                                                  >  1:4929133/1‑150 (MQ=255)
    gagaTCCTGAACAAGTACATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAg                                                                                                                                                  <  2:1630702/150‑1 (MQ=255)
    gagaTCCTGAACAAGTACATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAg                                                                                                                                                  >  1:2817926/1‑150 (MQ=255)
    gagaTCCTGAACAAGTACATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAg                                                                                                                                                  <  1:2796702/150‑1 (MQ=255)
    gagaTCCTGAACAAGTACATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAg                                                                                                                                                  <  1:2796664/150‑1 (MQ=255)
     agaTCCTGAACAAGTACATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAgg                                                                                                                                                 >  1:4436936/1‑150 (MQ=255)
        tCCTGAACAAGTACATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGggcg                                                                                                                                              <  2:1559771/150‑1 (MQ=255)
         ccTGAACAAGTACATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGggcgg                                                                                                                                             >  2:1366919/1‑150 (MQ=255)
          cTGAACAAGTACATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGggcggc                                                                                                                                            <  2:4244967/150‑1 (MQ=255)
          cTGAACAAGTACATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCGGCATACCGTGTAAACAAAAACGCCGCCCTTGAGGggcggc                                                                                                                                            <  2:1124963/150‑1 (MQ=255)
           tGAACAAGTACATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGgcggcg                                                                                                                                           >  1:53314/1‑150 (MQ=255)
            gAACAAGTACATCGATCCGAAGTACATCGACATCTGCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACAAAAACGCCGCCCTTGAGGGGCGGCGt                                                                                                                                          <  1:3348749/150‑1 (MQ=38)
             aaCAAGTACATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGtt                                                                                                                                         <  1:3864037/150‑1 (MQ=255)
             aaCAAGTACATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGtt                                                                                                                                         >  2:2580550/1‑150 (MQ=255)
                aaGTACATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACAAAAACGCCGCCCTTGAGGGGCGGCGttttc                                                                                                                                      <  2:1303327/150‑2 (MQ=255)
                 aGTACATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTtttg                                                                                                                                     <  2:4436936/150‑1 (MQ=255)
                 aGTACATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACAAAAACGCCGCCCTTGAGGGGCGGCGttttcg                                                                                                                                     <  2:2244573/150‑3 (MQ=255)
                    aCATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTtttgttt                                                                                                                                  <  2:3142105/150‑1 (MQ=255)
                    aCATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTtttgttt                                                                                                                                  <  2:2984420/150‑1 (MQ=255)
                      aTCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTtttgtggga                                                                                                                                >  2:168045/1‑146 (MQ=255)
                      aTCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTga                                                                                                                                <  1:26427/150‑1 (MQ=255)
                      aTCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTga                                                                                                                                >  1:1546817/1‑150 (MQ=255)
                       tCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACAAAAACGCCGCCCTTGAGGGGCGGCGTTttcgtttaca                                                                                                                               <  1:1693064/150‑4 (MQ=38)
                        cGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTgaga                                                                                                                              <  1:2335505/150‑1 (MQ=255)
                           tCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCg                                                                                                                           >  1:2339949/1‑150 (MQ=255)
                           tCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGGACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTATTTGTTTGAGACCg                                                                                                                           <  2:3949835/150‑1 (MQ=255)
                              gAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGttttt                                                                                                                                      <  2:3814066/136‑1 (MQ=255)
                                aGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTC                                                                                                                      <  1:186244‑M1/150‑1 (MQ=255)
                                 gTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCG                                                                                                                     >  2:1956672‑M1/1‑150 (MQ=255)
                                      tCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTtttgtttgca                                                                                                                               >  2:5345786/1‑133 (MQ=255)
                                      tCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCac                                                                                                                >  1:2741960‑M1/1‑148 (MQ=255)
                                       cGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCT                                                                                                               <  1:2694587‑M1/150‑1 (MQ=255)
                                       cGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCT                                                                                                               <  1:2694583‑M1/150‑1 (MQ=255)
                                        gACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTT                                                                                                              <  1:1435991‑M1/150‑1 (MQ=255)
                                             cTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGG                                                                                                         <  2:1967827‑M1/150‑1 (MQ=255)
                                                  tcgCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACG                                                                                                    <  2:4147726‑M1/150‑1 (MQ=255)
                                                     cgTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCC                                                                                                 <  1:1396186‑M1/150‑1 (MQ=255)
                                                        gCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCT                                                                                              <  1:2674811‑M1/150‑1 (MQ=255)
                                                          aaGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCC                                                                                            <  2:1321479‑M1/150‑1 (MQ=255)
                                                          aaGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGC                                                                                                       <  1:1222367‑M1/139‑1 (MQ=255)
                                                          aaGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGC                                                                                                       >  2:1222367‑M1/1‑139 (MQ=255)
                                                            ggAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGttt                                                                                                                                        <  1:1788331/104‑1 (MQ=255)
                                                            ggAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGttt                                                                                                                                        >  2:1788331/1‑104 (MQ=255)
                                                            ggAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTACATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCGTGCCCGCgaaga                                                                                          <  2:4453421‑M1/150‑6 (MQ=255)
                                                               aaGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCacag                                                                                       <  1:2282838‑M1/150‑5 (MQ=255)
                                                                aGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGATTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAgg                                                                                      <  1:4061372/150‑1 (MQ=255)
                                                                  cGAGCTGGAAGAGTGCGAGCACTCGATCTCCGCCCTCGAGTACAACTGGTACCTGCAGACCGTGTAAACGAAAACGCCGCCCCTCACGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGt                                                                                    <  1:3891431/150‑1 (MQ=255)
                                                                       tGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAAttt                                                                               <  1:347878/150‑1 (MQ=255)
                                                                        ggAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCCTACCGTGTAAACAAAAACGCCGCCCCTCCAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTc                                                                              <  2:4171855/150‑1 (MQ=255)
                                                                        ggAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCGCTTCGCGGGGACGCCCCCCCCCGCCGGGGAAtgtg                                                                              >  1:4361576/1‑147 (MQ=255)
                                                                          aaGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGa                                                                            <  1:3250135/150‑1 (MQ=255)
                                                                           agagTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGaa                                                                           <  2:4053457/150‑1 (MQ=255)
                                                                           agagTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGaa                                                                           <  2:4338592/150‑1 (MQ=255)
                                                                           agagTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGaa                                                                           <  2:4458494/150‑1 (MQ=255)
                                                                             agTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAgg                                                                         <  1:689260/150‑1 (MQ=255)
                                                                                  gAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGaa                                                                           >  2:3747772/1‑143 (MQ=255)
                                                                                  gAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGaa                                                                           <  1:3747772/143‑1 (MQ=255)
                                                                                    gCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGc                                                                  <  2:5364261/150‑1 (MQ=255)
                                                                                       cTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCtg                                                               <  1:4952560/150‑1 (MQ=255)
                                                                                        tCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTc                                                                              >  2:558556/1‑134 (MQ=255)
                                                                                        tCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTc                                                                              <  1:558556/134‑1 (MQ=255)
                                                                                        tCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCtgt                                                              <  1:3846467/150‑1 (MQ=255)
                                                                                        tCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCtgt                                                              <  2:4197256/150‑1 (MQ=255)
                                                                                        tCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCtgt                                                              >  2:3007553/1‑150 (MQ=255)
                                                                                          gATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAATCGAAATCGCCGCCCCTCAAGGGCGGCGTGTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCCCGCCCGCTCCCGCAGGGTAATTTCGCAGGCTGAAGCCTGTGa                                                            >  2:5235486/1‑150 (MQ=255)
                                                                                          gATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCT                                                                                                                 >  2:322761‑M1/1‑97 (MQ=255)
                                                                                          gATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGGGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCT                                                                                                                 <  1:322761‑M1/97‑1 (MQ=255)
                                                                                                cGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGGCCGCCCCCGCCCGCCCGGGTATTTTGGTTGCGGAAGCCAGTGGCTTAcc                                                      >  1:2889119/1‑150 (MQ=255)
                                                                                                  aCCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCacaggtaagtcacaggcttcagccttcgaaattaccctg                                                    <  1:2790650‑M1/150‑40 (MQ=255)
                                                                                                   ccTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGATTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCtgt                                                   <  2:3194385/150‑1 (MQ=255)
                                                                                                    cTCGTGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTg                                                  >  2:3060991‑M2/1‑149 (MQ=255)
                                                                                                          tACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTGGTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCACCCGCAGGGTAATTTCGAAGGCTGGAGCCGGTGACTTACCTGGgggcgcg                                            <  2:3269329‑M2/150‑8 (MQ=255)
                                                                                                           aCAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTgggagcgg                                           >  2:4957686‑M2/1‑142 (MQ=255)
                                                                                                             aaCTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCacaggtaagtcacaggcttcagccttcgaaattaccctgcgggagcgggc                                         <  1:5157354‑M1/150‑51 (MQ=255)
                                                                                                                tGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCC                                                                                                 >  1:5291000‑M1/1‑91 (MQ=255)
                                                                                                                 ggTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTgggagcgggcgtgc                                     >  2:4929133‑M2/1‑136 (MQ=255)
                                                                                                                  gTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCGAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCGTGCCCGCgaa                                                                                            <  1:3932980‑M1/94‑4 (MQ=255)
                                                                                                                  gTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCGAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCGTGCCCGCgaa                                                                                            >  2:3932935‑M1/1‑91 (MQ=255)
                                                                                                                     ccTGCATACCGTGTAAACGAAACCGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCGCCCGCAGGGTAATCTCGACGGCTGAAGCCTGTGACTTCCCTGTgggagcgggcgtgcccgc                                 >  2:4741769‑M2/1‑132 (MQ=255)
                                                                                                                          ctaCCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCTGGGTAATTTCGAAGGCTGAAGCCTGTGACTTCCCTGGgggcgcgggcgtgcccgcgaaga                            >  2:2339949‑M2/2‑127 (MQ=255)
                                                                                                                          aTACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGATTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTgggagcgggcgtgcccgcgaaga                            >  2:3125376‑M2/1‑127 (MQ=255)
                                                                                                                               gtgtAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCacaggtaagtcacaggcttcagccttcgaaattaccctgcgggagcgg                                           <  2:2633628‑M1/130‑49 (MQ=255)
                                                                                                                               gtgtAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCacaggtaagtcacaggcttcagccttcgaaattaccctgcgggagcgg                                           >  1:2633628‑M1/1‑82 (MQ=255)
                                                                                                                                   aaaCAAAAACGCCGCCCTTGAGGGGCGGCGTTTTCGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCCCGCCCGCCGGGTAATTTCGAAGGCTGAAGCCGGTGAATTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGGCGGccc                   >  2:1545609/1‑150 (MQ=38)
                                                                                                                                   aaaCAAAAACGCCGCCCTTGAGGGGCGGCGTTTTCGTTTGAGACCGGTGTCGCCCTCATAGCGGGCCCGCCCGCTCCCCCAGGGTAAATTCGAAGGCGGGACCCAGTGACGTACCAGTGGGAGCGGGCGTGCCCGCGCAGAGGGCGGccc                   >  1:2723062/1‑150 (MQ=25)
                                                                                                                                    aaCGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTAGCGGGCCACCCCGCACCCGCACGGTAATTTCGAAGGCTGGAGCCGGTGGCTTACCGGTTGGAGCCGGGGGGCCCCCGGAGGGGGCGGcccc                  >  2:52633/1‑149 (MQ=255)
                                                                                                                                     aCGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCCCAGCCCCAGGGTAATATTGAATGCAGAAGACTTTTACATACCAGTGGGAGCGGGCGTGCCCGCGAAGAGGGCGGCCCAg                 >  2:4189255/1‑150 (MQ=255)
                                                                                                                                       gAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCGCGTCGCGGGCCCGCCCCCGCCCGCAGGGTAAATTCGGAGGCTGAAGCCAGGGGCATAACAGTGGGAGCGGGCGTGCCCGCGAAGAGGGCGGCCCAAGc               >  1:1754926/1‑150 (MQ=255)
                                                                                                                                         aaaCGCCGCCCTTGAGGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCGCCCGCAGGGTAATTTCGAAGGCTGAAGCCGGGGACTTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGGCGGCCCAGGcgt             >  2:1339935/1‑148 (MQ=38)
                                                                                                                                          aaCGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCGCCCCCCCTCCCGCAGGGTTATTTCGAAGGCGGAAGCCTGTGTCATACCTGTGGGTGCGGGCGGGCCCCCGTAGAGGGCGGCCCAGGCGTa            >  2:2051679/1‑150 (MQ=255)
                                                                                                                                          aaCGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTgggagcgggcgt                                       <  2:1201801‑M2/123‑13 (MQ=255)
                                                                                                                                          aaCGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTgggagcgggcgt                                       >  1:1201801‑M2/1‑111 (MQ=255)
                                                                                                                                            cgccgcCCCTCCCGGGCGGCGTTTTTTTTTGCGACCGGTTTCGCCCTCTGCGGGGGCCCGCCCGCTCCCGCTGGTTATTTTCGAGGGCTCACGCCTTTGAGTTACCTGTGGGCGGGGGGGTCCCCGCGAAGAGGGCGGCCCAGGCTTACa          <  2:2817926/150‑1 (MQ=255)
                                                                                                                                            cgccgcCCCTACCGGGCGGCTTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGAAGGGTAATTTCGAAGGCTGAAGCCCGTGACTTACCTGTGGTAGGGGGCGTGCCCGCGAAGAGGGCGGCCCAGGCTTACa          <  2:4361576/150‑1 (MQ=255)
                                                                                                                                                    cTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCAAGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCGGGGACGTAACTGTGGGAGCGGGCGGGCCCGCGAAGAGGGCGGCCCAGGCGTACAGGGGCGtc  >  2:327890/1‑150 (MQ=255)
                                                                                                                                                    cTCAAGGGCGGCGGTTTTGTTAGAGACCGGAGTCGCCCTCATCTCGGGCACGCCCCCTCCCGCAGGGTAATTTCGAAGGCGGAACCCTGTGACTTAACGGGTTGAGCCGCCGTGGCCCCGAAGGGGGCGGCCCAGCCGTACGCGGGCttg  >  2:5315398/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CAGCGAGATCCTGAACAAGTACATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGGCGGCCCAGGCTTACAGGGGCTTC  >  NC_002947/5915954‑5916251

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