Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,916,151 A→G 6.5% intergenic (+65/+88) 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,1510AG 6.5% 356.8 / 8.5 126intergenic (+65/+88)spuI/argAglutamylpolyamine synthetase/amino acid N‑acetyltransferase
Reads supporting (aligned to +/- strand):  ref base A (56/54);  new base G (3/5);  total (64/62)
Fisher's exact test for biased strand distribution p-value = 7.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 131 total aligned reads displayed.

TCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGGCGGCCCAGGCTTACAGGGGCTTCTCGAAAATCTTCGAATTGCGCTGGTAGTTGTACAGCGAAGCCCGCGCCG  >  NC_002947/5916004‑5916300
                                                                                                                                                   |                                                                                                                                                      
tcgCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACG                                                                                                                                                      <  2:4147726‑M1/150‑1 (MQ=255)
 gcactgcgcactggttccgcgagcgcggctttgcgccgagtggcgtggagcggctgccggcggcgcgggcttcgctgtacaactaccagcgcaattcgaagattttcgagaagcccctgtaagcctgggcCGCCCTCTTCGCGGGCACGc                                                                                                                                                     <  1:3561638‑M1/20‑1 (MQ=255)
 gcactgcgcactggttccgcgagcgcggctttgcgccgagtggcgtggagcggctgccggcggcgcgggcttcgctgtacaactaccagcgcaattcgaagattttcgagaagcccctgtaagcctgggcCGCCCTCTTCGCGGGCACGc                                                                                                                                                     >  2:556111‑M1/131‑150 (MQ=255)
   cgTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCC                                                                                                                                                   <  1:1396186‑M1/150‑1 (MQ=255)
   cgTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCC                                                                                                                                                   <  1:2660046‑M1/150‑1 (MQ=255)
    ctgcgcactggttccgcgagcgcggctttgcgccgagtggcgtggagcggctgccggcggcgcgggcttcgctgtacaactaccagcgcaattcgaagattttcgagaagcccctgtaagcctgggcCGCCCTCTTCGCGGGCGCGCCCg                                                                                                                                                  <  1:5363774‑M1/23‑1 (MQ=255)
     tgcgcactggttccgcgagcgcggctttgcgccgagtggcgtggagcggctgccggcggcgcgggcttcgctgtacaactaccagcgcaattcgaagattttcgagaagcccctgtaagcctgggcCGCCCGCTTAGCGGGGccgcccgc                                                                                                                                                 <  1:3762855‑M1/24‑1 (MQ=255)
      gCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCT                                                                                                                                               <  1:2674811‑M1/150‑1 (MQ=255)
       cgcactggttccgcgagcgcggctttgcgccgagtggcgtggagcggctgccggcggcgcgggcttcgctgtacaactaccagcgcaattcgaagattttcgagaagcccctgtaagcctgggcCGCCCTCTTCGCGGGCACGCCCGCTc                                                                                                                                              >  2:873947‑M1/125‑150 (MQ=255)
       cgcactggttccgcgagcgcggctttgcgccgagtggcgtggagcggctgccggcggcgcgggcttcgctgtacaactaccagcgcaattcgaagattttcgagaagcccctgtaagcctgggcCGCCCTCTTCGCGGGCACGCCCGCTc                                                                                                                                              >  2:280767‑M1/125‑150 (MQ=255)
        aaGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCGTGCCCGCgaa                                                                                                                                              <  1:864372‑M1/150‑4 (MQ=255)
        aaGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCC                                                                                                                                             <  2:1321479‑M1/150‑1 (MQ=255)
         cactggttccgcgagcgcggctttgcgccgagtggcgtggagcggctgccggcggcgcgggcttcgctgtacaactaccagcgcaattcgaagattttcgagaagcccctgtaagcctgggcCGCCCTCTTCGCGGGCACGCCCGCTccc                                                                                                                                            >  1:4925342‑M1/123‑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)
                     gagcgcggctttgcgccgagtggcgtggagcggctgccggcggcgcgggcttcgctgtacaactaccagcgcaattcgaagattttcgagaagcccctgtaagcctgggcCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAAttt                                                                                                                                >  2:3225653‑M1/111‑150 (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:1622229/150‑1 (MQ=255)
                         agagTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGaa                                                                                                                            <  2:4458494/150‑1 (MQ=255)
                         agagTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGaa                                                                                                                            <  2:4338592/150‑1 (MQ=255)
                         agagTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGaa                                                                                                                            <  2:4053457/150‑1 (MQ=255)
                           agTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAgg                                                                                                                          <  1:689260/150‑1 (MQ=255)
                              tttgcgccgagtggcgtggagcggctgccggcggcgcgggcttcgctgtacaactaccagcgcaattcgaagattttcgagaagcccctgtaagcctgggcCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTg                                                                                                                       <  1:435617‑M1/49‑1 (MQ=255)
                                gAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGaa                                                                                                                            <  1:3747772/143‑1 (MQ=255)
                                gAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGaa                                                                                                                            >  2:3747772/1‑143 (MQ=255)
                                     cTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCtg                                                                                                                <  1:4952560/150‑1 (MQ=255)
                                      tCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTc                                                                                                                               >  2:558556/1‑134 (MQ=255)
                                      tCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTc                                                                                                                               <  1:558556/134‑1 (MQ=255)
                                      tCGATCTCCGACCTCGAGTACAACGGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCACGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGACGGCTGAAGCCtgt                                                                                                               <  2:4277790/150‑1 (MQ=255)
                                        gATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAATCGAAATCGCCGCCCCTCAAGGGCGGCGTGTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCCCGCCCGCTCCCGCAGGGTAATTTCGCAGGCTGAAGCCTGTGa                                                                                                             >  2:5235486/1‑150 (MQ=255)
                                           ctcCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACtt                                                                                                          <  2:1953376/150‑1 (MQ=255)
                                             gtggagcggctgccggcggcgcgggcttcgctgtacaactaccagcgcaattcgaagattttcgagaagcccctgtaagcctgggcCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTAc                                                                                                        <  2:4241382‑M1/64‑1 (MQ=255)
                                              cGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGGCCGCCCCCGCCCGCCCGGGTATTTTGGTTGCGGAAGCCAGTGGCTTAcc                                                                                                       >  1:2889119/1‑150 (MQ=255)
                                               ggagcggctgccggcggcgcgggcttcgctgtacaactaccagcgcaattcgaagattttcgagaagcccctgtaagcctgggcCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCt                                                                                                      >  2:3968193‑M1/85‑150 (MQ=255)
                                               gACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCt                                                                                                      <  1:1366919/150‑1 (MQ=255)
                                                aCCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCacaggtaagtcacaggcttcagccttcgaaattaccctg                                                                                                     <  1:2790650‑M1/150‑40 (MQ=255)
                                                 ccTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGATTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCtgt                                                                                                    <  2:3194385/150‑1 (MQ=255)
                                                  gcggctgccggcggcgcgggcttcgctgtacaactaccagcgcaattcgaagattttcgagaagcccctgtaagcctgggcCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCtgtg                                                                                                   >  1:4241382‑M1/82‑150 (MQ=255)
                                                  cTCGTGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTg                                                                                                   >  2:3060991‑M2/1‑149 (MQ=255)
                                                   tCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCCCGCCCGCGCCCCCACGGGAATTTCGAAGGGTGAAGCCTGTGGCTTAACAGTgg                                                                                                  <  1:2699769‑M2/150‑3 (MQ=255)
                                                   tCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGCCGCCCTCTTCGCGGGCACGCCCGCTCCCacaggtaagtcacaggcttcagccttcgaaattaccctgcgg                                                                                                  <  1:1764085‑M1/150‑43 (MQ=255)
                                                   cggctgccggcggcgcgggcttcgctgtacaactaccagcgcaattcgaagattttcgagaagcccctgtaagcctgggcCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTgg                                                                                                  >  2:1711223‑M1/81‑150 (MQ=255)
                                                     gcggccggcggcgcgggcttcgctgtacaactaccagcgcaattcgaagattttcgagaagccccggtaagcctgggcCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTACTTTCGCAGGCTGAAGCCTGTGACTTACCGGTGGGa                                                                                                >  2:5198825‑M1/79‑150 (MQ=255)
                                                        tACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTGGTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCACCCGCAGGGTAATTTCGAAGGCTGGAGCCGGTGACTTACCTGGgggcgcg                                                                                             <  2:3269329‑M2/150‑8 (MQ=255)
                                                         ccggcggcgcgggcttcgctgtacaactaccagcgcaattcgaagattttcgagaagcccctgtaagcctgggcCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCgg                                                                                            >  2:1920362‑M1/75‑150 (MQ=255)
                                                           aaCTGGTACCTGCATACCGTGTAAACGAACACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGCCGCCCTCTTCGCGGGCACGCCCGCTCCCacaggtaagtcacaggcttcagccttcgaaattaccctgcgggagcgggc                                                                                          <  1:1928202‑M1/150‑51 (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:3932980‑M1/1‑91 (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)
                                                                   ccTGCATACCGTGTAAACGAAACCGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCGCCCGCAGGGTAATCTCGACGGCTGAAGCCTGTGACTTCCCTGTgggagcgggcgtgcccgc                                                                                  >  2:4741769‑M2/1‑132 (MQ=255)
                                                                     tGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTgggagcgggcgtgcccgcga                                                                                >  2:4448852‑M2/1‑130 (MQ=255)
                                                                      cttcgctttacaactacccgcgcaattcgaagattttcgagaagcccctgtaagcctgggcCGCCCTCTTCGCGGGCACGCCCGCTCCCGCTGGTTAATTTCGAAGGCTGATGCCTGTGACTTACCTGTGGGAGCGGGCGTGCCCGCGaa                                                                               >  1:877401‑M1/62‑150 (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)
                                                                                 aaaCAAAAACGCCGCCCTTGAGGGGCGGCGTTTTCGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCCCGCCCGCCGGGTAATTTCGAAGGCTGAAGCCGGTGAATTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGGCGGccc                                                                    >  2:1545609/1‑150 (MQ=38)
                                                                                   aCGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGGCGGCCCAg                                                                  >  2:5363774/1‑150 (MQ=255)
                                                                                   aCGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCCCAGCCCCAGGGTAATATTGAATGCAGAAGACTTTTACATACCAGTGGGAGCGGGCGTGCCCGCGAAGAGGGCGGCCCAg                                                                  >  2:4189255/1‑150 (MQ=255)
                                                                                     gAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCGCGTCGCGGGCCCGCCCCCGCCCGCAGGGTAAATTCGGAGGCTGAAGCCAGGGGCATAACAGTGGGAGCGGGCGTGCCCGCGAAGAGGGCGGCCCAAGc                                                                >  1:1754926/1‑150 (MQ=255)
                                                                                       aaaCGCCGCCCTTGAGGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCGCCCGCAGGGTCATTTCGAAGGCTGAAGCCGGGGACGTACCGGTGGGAGCGGGCGTGCCCGCGAAGAGGGCGGCCCAGGCtt                                                              >  2:1339879/1‑150 (MQ=38)
                                                                                       aaaCGCCGCCCTTGAGGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCGCCCGCAGGGTAATTTCGAAGGCTGAAGCCGGGGACTTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGGCGGCCCAGGcgt                                                              >  2:1339935/1‑148 (MQ=38)
                                                                                        aaCGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCGCCCCCCCTCCCGCAGGGTTATTTCGAAGGCGGAAGCCTGTGTCATACCTGTGGGTGCGGGCGGGCCCCCGTAGAGGGCGGCCCAGGCGTa                                                             >  2:2051679/1‑150 (MQ=255)
                                                                                        aaCGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTgggagcgggcgt                                                                                        >  1:1201801‑M2/1‑111 (MQ=255)
                                                                                        aaCGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTgggagcgggcgt                                                                                        <  2:1201801‑M2/123‑13 (MQ=255)
                                                                                        aaCGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTgggagcggg                                                                                           >  1:5059419‑M2/1‑111 (MQ=255)
                                                                                        aaCGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTgggagcggg                                                                                           <  2:5059419‑M2/120‑10 (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)
                                                                                                   tCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCCCAGGGTAAATTCGAAGGCGGAAGCCGGTGACTTACCAGTGGGAGCGGGCGTGCCCCCGAAGAGGGCGGCCCAGGCATACACGGGCTtcg                                                  >  2:3379728/1‑149 (MQ=255)
                                                                                                   tCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCGTCGCCGGCACCCCCCCAACCCCCAGGGAATTTAGAAGGCGGAAGCCAGTGACAAAACATTGGGAGCGGGCGGGCCCCCGCATACGGCGGCCCACGCGTTCAGGGGCAtca                                                  >  2:4565691/1‑149 (MQ=255)
                                                                                                        ggcggcTTTTTGGTTTGAGCCGGTTGTCGCCCTCTTCGCGGGCACGCCCGCCCCGGCAGGGTATTTTCCATGGCTGAAGCGTTTGAGTTACCTGTGGGTGGGGGCGGGCCCGCGAAGAGGGCGGCCCAGGCTTACAGGGGCTTCTCGaaa                                             <  2:1021337/150‑1 (MQ=255)
                                                                                                                           aCCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCGGGCGGGCCCGCGAAGAGGGCGGCCCAGGCTTACAGGGGCTTCTCGAAAATCTTCGCATTGCGCTggg                          >  2:3581118/1‑149 (MQ=255)
                                                                                                                              ggTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGGCGGCCCAGGCTTACAGGGGCTTCTCGAAAATCTTCGAATTGCGCTGGTAGt                       >  2:3180174/1‑150 (MQ=255)
                                                                                                                                    GCCCTCTTCGCGGGCACGCCCGCGCCCGCAGGGTCATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGGCGGCCCCGGCGTACAGGGGCTTCTCGAAAAACTACGAATTGCGCTGGTAGTTGTACa                 <  2:2844936‑M1/150‑1 (MQ=255)
                                                                                                                                      CCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCAGTGGGAGCGGGCGTGCCCGCGAAGAGGGCGGCCCAGGCTTACAGGGGCTTCTCGAAAATCTTCGAATTGCGCTGGTAGTTGTACAGc               <  2:3044342‑M1/150‑1 (MQ=255)
                                                                                                                                        tctTCGCGGGCACGCCCGCTCCCacaggtaagtcacaggcttcagccttcgaaattaccctgcgggagcgggcgtgcccgcgaagagggcgacaccggtctcaaacaaaatcgccgcccttgaggggcggcgttttcgtttacacggtat             >  2:3125376‑M1/1‑23 (MQ=255)
                                                                                                                                        tctTCGCGGGCACGCCCGCGCCCccagggaagtcacaggcttcagccttcgaaattacccagcgggagcgggcgtgcccgcgaagagggcgacaccggtctcaaacaaaaacgccgcccttgaggggcggcgttttcgtttacacggtag             >  2:2339949‑M1/1‑23 (MQ=255)
                                                                                                                                        TCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGGCGGCCCAGGCTTACAGGGGCTTCTCGAAAATCTTCGAATTGCGCTGGTAGTTGTACAGCGa             <  2:3118812‑M1/150‑1 (MQ=255)
                                                                                                                                         cggCGCGGGCACGCCCGCCCCCGCAGGGTAATTTCGACGGCTGAAGCCTGTGAGTTCCGTGTGGGAGCGGGCGTGCCCGCGAAGAGGGCGGCCCAGGCTTACAGGGGCTTCTCGAAAATCTTCGAATTGCGCTGGTAGTTGTACAGCGaa            >  2:1546817‑M1/4‑150 (MQ=255)
                                                                                                                                          TTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGTGCGGGCGTGCCCGCGAAGAGGGCGGCCCAGGCTTACAGGGGCTTCTCGAAAATCTTCGAAGTGCGCTGGTAGTTGTACAGCGAAg           <  2:4778261‑M1/150‑1 (MQ=255)
                                                                                                                                           tCGCGGGCACGCCCGCTCCCacaggtaagtcacaggcttcagccttcgaaattaccctgcgggagcgggcgtgcccgcgaagagggcgacaccggtctcaaacaaaaacgccgcccttgaggggcggcgttttcgtttacacggtatgca          >  2:4448852‑M1/1‑20 (MQ=255)
                                                                                                                                             gcgGGCACGCCCGCTCCCacagggaagtcacaggcttcagccgtcgagattaccctgcgggcgcgggcgtgcccgcgaagagggcgacaccggtctcaaacaaaaacgccgcccttgaggggcggcggtttcgtttacacggtatgcagg        >  2:4741769‑M1/1‑18 (MQ=255)
                                                                                                                                              cgGGCACGCCCGGCCCCaccggtaaggcacgggctccacccttcgaaattaccctgggggtgcgggcgggcccgcgatttgggcgacaccggtctcaaacaaaaacgccgcccttgaggggcggcgttttcgtttacacggtatgcaggt       <  2:3335309‑M1/150‑134 (MQ=255)
                                                                                                                                                ggCACGCCCGCTCCCacaggtaagtcacaggcttcagccttcgaaattaccctgcgggagcgggcgtgcccgcgaagagggcgacaccggtctcaaacaaaaacgccgcccttgaggggcggcgttttcgtttacacggtatgcaggtac     >  1:2580550‑M1/1‑15 (MQ=255)
                                                                                                                                                 gCACGCCCGCTCCCacaggtaagtcacaggcttcagccttcgaaattaccctgcgggagcgggcgtgcccgcgaagagggcgacaccggtctcaaacaaaaacgccgcccttgaggggcggcgttttcgtttacacggtatgcaggtacc    >  2:4929133‑M1/1‑14 (MQ=255)
                                                                                                                                                   aCGCCCGCTCCCacaggtaagtcacaggcttcagccttcgaatttaccctgcgggagcgggcgtgcccgcgaagagggcgacaccggtctcaaacaaaaacgccgcccttgaggggcggcgttttcgtttacacggtatgcaggtaccag  <  2:3018463‑M1/150‑139 (MQ=255)
                                                                                                                                                   aCGCCCGCTCCCacaggtaagtcacaggcttcagccttcgaaattaccctgcgggagcgggcgtgcccgcgaagagggcgacaccggtctcaaacaaaaacgccgcccttgaggggcggcgtt                             <  2:1201801‑M1/123‑112 (MQ=255)
                                                                                                                                                   aCGCCCGCTCCCacaggtaagtcacaggcttcagccttcgaaattaccctgcgggagcgggcgtgcccgcgaagagggcgacaccggtctcaaacaaaaacgccgcccttgaggggcggcgtt                             >  1:1201801‑M1/1‑12 (MQ=255)
                                                                                                                                                   |                                                                                                                                                      
TCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGGCGGCCCAGGCTTACAGGGGCTTCTCGAAAATCTTCGAATTGCGCTGGTAGTTGTACAGCGAAGCCCGCGCCG  >  NC_002947/5916004‑5916300

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