Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,916,104 C→G 7.9% intergenic (+18/+135) 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,1040CG 7.9% 450.9 / 22.6 140intergenic (+18/+135)spuI/argAglutamylpolyamine synthetase/amino acid N‑acetyltransferase
Reads supporting (aligned to +/- strand):  ref base C (49/79);  new base G (4/7);  total (53/87)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

Only 100 of 142 total aligned reads displayed.

AGCGAGATCCTGAACAAGTACATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGGCGGCCCAGGCTTACAGGGGCTTCTC  >  NC_002947/5915955‑5916253
                                                                                                                                                     |                                                                                                                                                      
aGCGAGATCCTGAACAAGTACATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTc                                                                                                                                                        >  1:1021337/1‑150 (MQ=255)
  cGAGATCCTGAACAAGTACATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCaa                                                                                                                                                      <  2:2422862/150‑1 (MQ=255)
  cGAGATCCTGAACAAGTACATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCaa                                                                                                                                                      <  2:5316985/150‑1 (MQ=255)
   gagaTCCTGAACAAGTACATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAg                                                                                                                                                     <  2:1630702/150‑1 (MQ=255)
   gagaTCCTGAACAAGTACATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAg                                                                                                                                                     <  2:1448289/150‑1 (MQ=255)
   gagaTCCTGAACAAGTACATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAg                                                                                                                                                     >  1:2817926/1‑150 (MQ=255)
   gagaTCCTGAACAAGTACATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAg                                                                                                                                                     <  1:2796664/150‑1 (MQ=255)
   gagaTCCTGAACAAGTACATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAg                                                                                                                                                     <  2:4194969/150‑1 (MQ=255)
   gagaTCCTGAACAAGTACATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAg                                                                                                                                                     >  1:4929133/1‑150 (MQ=255)
         cTGAACAAGTACATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGggcggc                                                                                                                                               <  2:4244967/150‑1 (MQ=255)
         cTGAACAAGTACATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACAAAAACGCCGCCCTTGAGGggcggc                                                                                                                                               <  1:2809934/150‑1 (MQ=255)
         cTGAACAAGTACATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCGGCATACCGTGTAAACAAAAACGCCGCCCTTGAGGggcggc                                                                                                                                               <  2:1124963/150‑1 (MQ=255)
            aaCAAGTACATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGtt                                                                                                                                            <  1:41455/150‑1 (MQ=255)
               aaGTACATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACAAAAACGCCGCCCTTGAGGGGCGGCGttttc                                                                                                                                         <  2:1303327/150‑2 (MQ=255)
                  tACATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGtttttttt                                                                                                                                      >  1:3000112/1‑150 (MQ=255)
                   aCATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTtttgttt                                                                                                                                     <  2:3142105/150‑1 (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)
                          tCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGGACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTATTTGTTTGAGACCg                                                                                                                              <  2:3949835/150‑1 (MQ=255)
                             gAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGttttt                                                                                                                                         <  2:3814066/136‑1 (MQ=255)
                             gAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGttttt                                                                                                                                         >  1:3814066/1‑136 (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)
                                     tCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTtttgtttgca                                                                                                                                  <  1:5345786/135‑3 (MQ=255)
                                     tCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCac                                                                                                                   >  1:2741960‑M1/1‑148 (MQ=255)
                                      cGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCT                                                                                                                  <  1:2694583‑M1/150‑1 (MQ=255)
                                       gACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTT                                                                                                                 <  1:1435991‑M1/150‑1 (MQ=255)
                                            cTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCACGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGG                                                                                                            <  2:1968080‑M1/150‑1 (MQ=255)
                                             ttCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGG                                                                                                           <  1:4751018‑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)
                                                         aaGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCGTGCCCGCgaa                                                                                               <  1:864372‑M1/150‑4 (MQ=255)
                                                         aaGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGC                                                                                                          <  1:1222367‑M1/139‑1 (MQ=255)
                                                           ggAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGttt                                                                                                                                           <  1:1788331/104‑1 (MQ=255)
                                                           ggAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGttt                                                                                                                                           >  2:1788331/1‑104 (MQ=255)
                                                           ggAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTACATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCGTGCCCGCgaaga                                                                                             <  2:4453421‑M1/150‑6 (MQ=255)
                                                            gAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGc                                                                                           >  2:3250135/1‑150 (MQ=255)
                                                             aaaGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCa                                                                                          <  1:343209/150‑1 (MQ=255)
                                                              aaGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCacag                                                                                         <  1:2282838‑M1/150‑5 (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:1622229/150‑1 (MQ=255)
                                                                          agagTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGaa                                                                             <  2:4458494/150‑1 (MQ=255)
                                                                          agagTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGaa                                                                             <  2:4338592/150‑1 (MQ=255)
                                                                           gagTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAg                                                                            <  2:1530675/150‑1 (MQ=255)
                                                                                 gAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGaa                                                                             >  2:3747772/1‑143 (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                                                                <  2:4197256/150‑1 (MQ=255)
                                                                                       tCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCtgt                                                                <  1:3846467/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)
                                                                                            ctcCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACtt                                                           <  2:1953376/150‑1 (MQ=255)
                                                                                                 aCCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCacaggtaagtcacaggcttcagccttcgaaattaccctg                                                      <  1:2790650‑M1/150‑40 (MQ=255)
                                                                                                    tCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCCCGCCCGCGCCCCCACGGGAATTTCGAAGGGTGAAGCCTGTGGCTTAACAGTgg                                                   <  1:2699769‑M2/150‑3 (MQ=255)
                                                                                                    tCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCGCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACGTACCGGTgg                                                   <  2:313428‑M2/150‑3 (MQ=255)
                                                                                                          aCAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTgggagcgg                                             >  2:4957686‑M2/1‑142 (MQ=255)
                                                                                                            aaCTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCacaggtaagtcacaggcttcagccttcgaaattaccctgcgggagcgggc                                           <  1:5157354‑M1/150‑51 (MQ=255)
                                                                                                              cTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAAATTCGAAGGCTGAAGCCTGTGACTTACCTGTgggagcgggcgt                                         <  2:3018463‑M2/150‑13 (MQ=255)
                                                                                                               tGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCC                                                                                                    <  2:5291000‑M1/91‑1 (MQ=255)
                                                                                                               tGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCC                                                                                                    >  1:5291000‑M1/1‑91 (MQ=255)
                                                                                                                ggTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTgggagcgggcgtgc                                       >  2:4929133‑M2/1‑136 (MQ=255)
                                                                                                                 gTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCGAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCGTGCCCGCgaa                                                                                               >  2:3932935‑M1/1‑91 (MQ=255)
                                                                                                                 gTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCGAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCGTGCCCGCgaa                                                                                               <  1:3932980‑M1/94‑4 (MQ=255)
                                                                                                                 gTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCGAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCGTGCCCGCgaa                                                                                               <  1:3932935‑M1/94‑4 (MQ=255)
                                                                                                                 gTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTgggagcgggcgtgcc                                      >  1:2580550‑M2/1‑135 (MQ=255)
                                                                                                                   aCCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCAAATCGCGGGCCCGCCCGCACCCCCAGGGTAATTTCGAAGGGTGGAGCCCGTGCCTTACCGGTggggccgggcgtgcccg                                    <  2:3335309‑M2/150‑18 (MQ=255)
                                                                                                                    ccTGCATACCGTGTAAACGAAACCGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCGCCCGCAGGGTAATCTCGACGGCTGAAGCCTGTGACTTCCCTGTgggagcgggcgtgcccgc                                   >  2:4741769‑M2/1‑132 (MQ=255)
                                                                                                                      tGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTgggagcgggcgtgcccgcga                                 >  2:4448852‑M2/1‑130 (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)
                                                                                                                                      gAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCGCGTCGCGGGCCCGCCCCCGCCCGCAGGGTAAATTCGGAGGCTGAAGCCAGGGGCATAACAGTGGGAGCGGGCGTGCCCGCGAAGAGGGCGGCCCAAGc                 >  1:1754926/1‑150 (MQ=255)
                                                                                                                                        aaaCGCCGCCCTTGAGGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCGCCCGCAGGGTCATTTCGAAGGCTGAAGCCGGGGACGTACCGGTGGGAGCGGGCGTGCCCGCGAAGAGGGCGGCCCAGGCtt               >  2:1339879/1‑150 (MQ=38)
                                                                                                                                        aaaCGCCGCCCTTGAGGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCGCCCGCAGGGTAATTTCGAAGGCTGAAGCCGGGGACTTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGGCGGCCCAGGcgt               >  2:1339935/1‑148 (MQ=38)
                                                                                                                                        aaaCGCCGCCCCTCAAGGGCGGCGTTTTTGGTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCCCACCCCCAGGGTAAGTTCGAAGGCGGTAGCCGGTTACTTACCAGGGGGGGCGGGCGTGCCCGCGGAGTGGGCGGCCCAGGCtt               >  2:125969/1‑150 (MQ=255)
                                                                                                                                         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)
                                                                                                                                                   cTCATGGGCGGCTTTTTTTTTTGTGACCGGTGTCGCCCTATTCGCGGGCACGCCCGCTCCCGCAGGTTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGGCGGCCCAGGCTTACAGGGGCTtc    <  2:2683839/150‑1 (MQ=38)
                                                                                                                                                   cTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCAAGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCGGGGACGTAACTGTGGGAGCGGGCGGGCCCGCGAAGAGGGCGGCCCAGGCGTACAGGGGCGtc    >  2:327890/1‑150 (MQ=255)
                                                                                                                                                    tCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCGTCGCCGGCACCCCCCCAACCCCCAGGGAATTTAGAAGGCGGAAGCCAGTGACAAAACATTGGGAGCGGGCGGGCCCCCGCATACGGCGGCCCACGCGTTCAGGGGCAtc  >  2:4565691/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
AGCGAGATCCTGAACAAGTACATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGGCGGCCCAGGCTTACAGGGGCTTCTC  >  NC_002947/5915955‑5916253

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: