Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,914,126 C→G 5.6% 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,914,1260CG 5.6% 266.7 / 3.5 84intergenic (+18/+135)spuI/argAglutamylpolyamine synthetase/amino acid N‑acetyltransferase
Reads supporting (aligned to +/- strand):  ref base C (40/38);  new base G (2/3);  total (42/42)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAGATCCTGAACAAGTACATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTC‑CGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGGCGGCCCAGGCTT  >  NC_002947/5913980‑5914262
                                                                                                                                                   |                                                                                                                                        
gagaTCCTGAACAAGTACATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTC‑CGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTc                                                                                                                                          >  2:291588/1‑147 (MQ=255)
 agaTCCTGAACAAGTACATCGATCCGAAGTACATCGTCATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTC‑CGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAg                                                                                                                                       >  1:61033/1‑149 (MQ=255)
  gaTCCTGAACAAGTACATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTC‑CGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAgg                                                                                                                                      >  1:126866/1‑149 (MQ=255)
      cTGAACAAGTACATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTC‑CGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGggcgg                                                                                                                                  >  1:52141/1‑149 (MQ=255)
       tGAACAAGTACATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTC‑CGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCTTGAGGggcgg                                                                                                                                  >  5:145964/1‑148 (MQ=255)
       tGAACAAGTACATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTC‑CGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGggcggc                                                                                                                                 >  2:309349/1‑149 (MQ=255)
       tGAACAAGTACATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTC‑CGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGggcggc                                                                                                                                 >  7:261528/1‑149 (MQ=255)
        gAACAAGTACATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTC‑CGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGgcggcg                                                                                                                                >  6:151248/1‑149 (MQ=255)
         aaCAAGTACATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTC‑CGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCTTGAGGggcgg                                                                                                                                  <  6:145964/146‑1 (MQ=255)
         aaCAAGTACATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTC‑CGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACAAAAACGCCGCCCTTGAGGGgcggcg                                                                                                                                <  5:228990/148‑1 (MQ=255)
          aCAAGTACATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTC‑CGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGtt                                                                                                                              <  5:151248/149‑1 (MQ=255)
               tACATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTC‑CGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTTAAGGGCGGCGTTtttgt                                                                                                                         <  2:10011/149‑1 (MQ=255)
               tACATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTC‑CGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTtttgt                                                                                                                         >  6:364812/1‑149 (MQ=255)
                aCATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTC‑CGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTtttgtt                                                                                                                        >  4:57043/1‑149 (MQ=255)
                aCATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTGCTAGCACTCGATCTC‑CGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTtttgt                                                                                                                         <  3:50219/148‑1 (MQ=255)
                 cATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTC‑CGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTtttgttt                                                                                                                       <  2:126866/149‑1 (MQ=255)
                 cATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTC‑CGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTtttgttt                                                                                                                       >  6:33538/1‑149 (MQ=255)
                 cATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTC‑CGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTtttgtt                                                                                                                        <  3:350198/148‑1 (MQ=255)
                  aTCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTC‑CGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTtttgttt                                                                                                                       >  7:310792/1‑148 (MQ=255)
                  aTCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTC‑CGACCTCGAGTAAAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTtttgttt                                                                                                                       >  7:197441/1‑148 (MQ=255)
                   tCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTC‑CGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTga                                                                                                                     <  6:155205/149‑1 (MQ=255)
                     gATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTC‑CGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTgaga                                                                                                                   >  8:144341/1‑149 (MQ=255)
                     gATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTC‑CGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTgaga                                                                                                                   <  5:364812/149‑1 (MQ=255)
                      aTCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTC‑CGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGAc                                                                                                                  >  8:171933/1‑149 (MQ=255)
                      aTCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAACCGAGCTGGAAGAGTTCGAGCACTCGATCTC‑CGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGAc                                                                                                                  <  1:309349/149‑1 (MQ=255)
                          gAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTC‑CGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGAc                                                                                                                  >  5:195638/1‑145 (MQ=255)
                          gAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTC‑CGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGAc                                                                                                                  <  6:195638/145‑1 (MQ=255)
                           aaGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCGAGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGgt                                                                                                              >  1:298689/1‑149 (MQ=255)
                                cATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTC‑CGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGcc                                                                                                        <  2:268097/149‑1 (MQ=255)
                                 aTCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTC‑CGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGccc                                                                                                       >  2:217891/1‑149 (MQ=255)
                                     aCATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTC‑CGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCtt                                                                                                   >  8:253105/1‑149 (MQ=255)
                                     aCATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTC‑CGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCtt                                                                                                   <  3:338639/149‑1 (MQ=255)
                                         cTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTC‑CGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCACTCAAGGGCGGCGTTTTTGTTTGAGACCGTTGTCGCCCTCTTcgcg                                                                                               >  6:54334/1‑149 (MQ=255)
                                          ttCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTC‑CGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCgg                                                                                              >  5:272552/1‑149 (MQ=255)
                                           tcgtcgCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTC‑CGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCgg                                                                                              >  1:243894/1‑148 (MQ=255)
                                           tcgtcgCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTC‑CGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCgg                                                                                              <  2:243894/148‑1 (MQ=255)
                                                gcgTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTC‑CGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGc                                                                                        <  3:57043/149‑1 (MQ=255)
                                                gcgTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTC‑CGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCATCGCGGGCACGc                                                                                        <  7:144341/149‑1 (MQ=255)
                                                        ggAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTC‑CGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCACTccc                                                                                >  1:106065/1‑149 (MQ=255)
                                                        ggAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTC‑CGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACACCCGCTccc                                                                                >  3:157847/1‑149 (MQ=255)
                                                            aGCGAGCTGGAAGAGTTCGAGCACTCGATCTC‑CGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAg                                                                            <  2:276359/149‑1 (MQ=255)
                                                             gCGAGCTGGAAGAGTTCGAGCACTCGATCTC‑CGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAgg                                                                           <  1:217891/149‑1 (MQ=255)
                                                             gCGAGCTGGAAGAGTTCGAGCACTCGATCTC‑CGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAgg                                                                           <  5:33538/149‑1 (MQ=255)
                                                                   tGGAAGAGTTCGAGCACTCGATCTC‑CGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTAGCGGGCACCCCCGCTCGCGCAGGGTAAtt                                                                     >  3:90640/1‑149 (MQ=255)
                                                                     gAAGAGTTCGAGCACTCGATCTC‑CGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTc                                                                   <  8:197441/149‑1 (MQ=255)
                                                                     gAAGAGTTCGAGCACTCGATCTC‑CGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTc                                                                   >  2:267344/1‑149 (MQ=255)
                                                                        gagTTCGAGCACTCGATCTC‑CGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGaa                                                                <  5:222146/149‑1 (MQ=255)
                                                                              gAGCACTCGATCTC‑CGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGGTg                                                           >  7:282353/1‑148 (MQ=255)
                                                                                gCACTCGATCTC‑CGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAg                                                        <  1:291588/149‑1 (MQ=255)
                                                                                  aCTCGATCTC‑CGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAgg                                                                           <  2:42569/128‑1 (MQ=255)
                                                                                   cTCGATCTC‑CGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAgg                                                              >  2:154196/1‑140 (MQ=255)
                                                                                   cTCGATCTC‑CGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAgg                                                              <  1:154196/140‑1 (MQ=255)
                                                                                      gATCTC‑CGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCtgt                                                   <  8:261528/148‑1 (MQ=255)
                                                                                             cGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCTCCCCTCAATGGCGTCGTTTTTGTTTGAGACCGGTGTCGACCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTAc                                            <  5:54334/149‑1 (MQ=255)
                                                                                             cGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGa                                                 <  7:253105/144‑1 (MQ=255)
                                                                                                 cTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGttttt                                                                                                                           <  6:363093/66‑1 (MQ=255)
                                                                                                 cTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGttttt                                                                                                                           >  5:363093/1‑66 (MQ=255)
                                                                                                 cTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCtgt                                        <  7:171933/149‑1 (MQ=255)
                                                                                                  tCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCtgtg                                       >  3:80465/1‑149 (MQ=255)
                                                                                                   cGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGa                                                 >  7:64364/1‑138 (MQ=255)
                                                                                                   cGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGa                                                 <  8:64364/138‑1 (MQ=255)
                                                                                                       tACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGc                                  <  2:61033/149‑1 (MQ=255)
                                                                                                         cAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTGGAAGGCTGAAGCCGGTGACTTACCGGTGGGAGCgg                                >  3:341141/1‑149 (MQ=255)
                                                                                                         cAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCgg                                <  1:267344/149‑1 (MQ=255)
                                                                                                         cAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACACCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCg                                 <  4:157847/148‑1 (MQ=255)
                                                                                                           aCTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCGGGc                              <  8:310792/149‑1 (MQ=255)
                                                                                                            cTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCAGAAGCCTGTGAGGTACCTGTGGGAGCGGGCg                             >  8:283816/1‑149 (MQ=255)
                                                                                                               gTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCGGGCGTGc                          <  2:298689/149‑1 (MQ=255)
                                                                                                               gTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACCCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCGGTGACGTACCTGTGGGAGCGGGCGTGc                          >  1:322285/1‑149 (MQ=255)
                                                                                                                   cTGCATACCGTGTAAACGAAAACGCCGCCCATCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGGAATTTCGAAGGCTGAAGCCTGTGACTTACCGGTGGGGGCGGGCGTGCCcgc                      >  1:32372/1‑149 (MQ=255)
                                                                                                                       aTACCGTGTAACCGAAATCCCCGCCCCTCACGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCCGGCACGCCCGCTCCCGCAGGGTAATTTCCAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCGGGCGTGCCCGCGAag                  <  7:283816/149‑1 (MQ=39)
                                                                                                                        tACCGTGTAAACGAAAACCCCGCCCCTCACGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCGGGCGTGCCCGCGAaga                 <  5:152478/149‑1 (MQ=255)
                                                                                                                         aCCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTgg                                      >  3:217176/1‑127 (MQ=255)
                                                                                                                         aCCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTgg                                      <  4:217176/127‑1 (MQ=255)
                                                                                                                          ccGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCGGGCGTGCCCGCGAagag                <  2:52141/148‑1 (MQ=255)
                                                                                                                              gtAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCACTTCGCGGGCCCGCCCCCGCCCCCAGGGGAATTTAGAAAGGGGAAGCCGGTGGCGTTACAGTGGGGGCGGGCGGGGCCGCGGAGACggggg           >  2:127076/1‑149 (MQ=255)
                                                                                                                                 aaCAAAAACGCCGCCCTTGAGGGGCGGCGTTTTTGATTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTAc                                            <  8:111370/113‑1 (MQ=38)
                                                                                                                                 aaCAAAAACGCCGCCCTTGAGGGGCGGCGTTTTTGATTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTAc                                            >  7:111370/1‑113 (MQ=38)
                                                                                                                                    gAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGGAATTTAGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGGCGGCCCAgg     >  2:66722/1‑149 (MQ=255)
                                                                                                                                       aaCGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGGAATTTCGAAGGCTGCAGCCGGTGACGTACCGGGGGGAGCGGGCGTGCCCGCGAAGAGGGCGGCCCCGGcg   >  1:295496/1‑147 (MQ=38)
                                                                                                                                       aaCGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCGCCCGCAGGGTTATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGGCGGCCCCGGcag  >  7:26698/1‑147 (MQ=255)
                                                                                                                                        aCGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGaa                                                         <  3:202803/93‑1 (MQ=255)
                                                                                                                                        aCGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGaa                                                         >  4:202803/1‑93 (MQ=255)
                                                                                                                                                 cTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGc                                  >  3:299416/1‑107 (MQ=255)
                                                                                                                                                 cTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGc                                  <  4:299416/107‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                        
GAGATCCTGAACAAGTACATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTC‑CGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGGCGGCCCAGGCTT  >  NC_002947/5913980‑5914262

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑