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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,914,1240CT 8.9% 219.0 / 7.1 78intergenic (+16/+137)spuI/argAglutamylpolyamine synthetase/amino acid N‑acetyltransferase
Reads supporting (aligned to +/- strand):  ref base C (37/34);  new base T (4/3);  total (41/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCGAGATCCTGAACAAGTACATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTA‑CAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGGCGGCCCAGGCTTA  >  NC_002947/5913977‑5914263
                                                                                                                                                    |                                                                                                                                           
aGCGAGATCCTGAACAAGTACATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTA‑CAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCt                                                                                                                                            <  6:16623/149‑1 (MQ=255)
aGCGAGATCCTGAACAAGTACATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTA‑CAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCt                                                                                                                                            >  3:316255/1‑149 (MQ=255)
   gagaTCCTGAACAAGTACATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTA‑CAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCaa                                                                                                                                         >  4:209869/1‑149 (MQ=255)
       tCCTGAACAAGTACATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTA‑CAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGggc                                                                                                                                     <  5:260705/149‑1 (MQ=255)
       tCCTGAACAAGTACATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTA‑CAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGggc                                                                                                                                     <  3:209869/149‑1 (MQ=255)
           gAACAAGTACATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTA‑CAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGgcggcg                                                                                                                                 >  4:73983/1‑149 (MQ=255)
             aCAAGTACATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTA‑CAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGtt                                                                                                                               <  8:227079/149‑1 (MQ=255)
              cAAGTACATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTA‑CAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGttt                                                                                                                              >  7:197374/1‑149 (MQ=255)
                aGTACATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTA‑CAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGttttt                                                                                                                            <  8:197374/149‑1 (MQ=255)
                    cATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTA‑CAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAgg                                                                                                                                       >  3:120906/1‑134 (MQ=255)
                    cATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTA‑CAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAgg                                                                                                                                       <  4:120906/134‑1 (MQ=255)
                    cATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTA‑CAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTtttgttt                                                                                                                        <  5:80155/149‑1 (MQ=255)
                    cATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTA‑CAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTtttgttt                                                                                                                        >  5:153606/1‑149 (MQ=255)
                    cATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTA‑CAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTtttgtt                                                                                                                         <  7:188933/148‑1 (MQ=255)
                       cGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTA‑CAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTtttgttt                                                                                                                        >  3:289443/1‑146 (MQ=255)
                        gATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTA‑CAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCTTGAGGGGCGGCGTTttcgtttaca                                                                                                                     <  3:29135/148‑4 (MQ=39)
                        gATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTA‑CAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTgaga                                                                                                                    <  1:157957/149‑1 (MQ=255)
                          tCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTA‑CAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGAcc                                                                                                                  >  6:144072/1‑149 (MQ=255)
                            cGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTA‑CAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCg                                                                                                                 <  5:1483/148‑1 (MQ=255)
                                gTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTA‑CAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTtttgttt                                                                                                                        >  5:190124/1‑137 (MQ=255)
                                gTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTA‑CAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTtttgttt                                                                                                                        <  6:190124/137‑1 (MQ=255)
                                     tCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTA‑CAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCct                                                                                                       <  5:144072/149‑1 (MQ=255)
                                        aCATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTAGCAACTGGTACATGAATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCTTTTTTGTTTGAGACCGGTGTCGCCctct                                                                                                     >  1:275325/1‑149 (MQ=255)
                                                          aGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTA‑CAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTcc                                                                                  >  2:188387/1‑149 (MQ=255)
                                                          aGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTA‑CAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTcc                                                                                  >  4:274216/1‑149 (MQ=255)
                                                           ggAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTA‑CAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTccc                                                                                 >  7:131724/1‑149 (MQ=255)
                                                           ggAAAGCGAGCTGGAAGAGTTCGAGCACTAGATCTCCGAACTCGAGTA‑CAACTGGTACATGCATACCGAGTAAAAGAAAACGCCGCCCCTCAAGGCCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCTCCCACGCCCGCAcc                                                                                  >  8:200147/1‑148 (MQ=255)
                                                             aaaTCGATCTGGAAGTGTTCGATCACTCGATATCCGACCTCGAGTA‑CAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGAGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCg                                                                                <  7:200147/148‑1 (MQ=255)
                                                              aaGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTA‑CAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTcgc                                                                                                 <  8:207381/130‑1 (MQ=255)
                                                              aaGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTA‑CAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTcgc                                                                                                 >  7:207381/1‑130 (MQ=255)
                                                                gCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTA‑CAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAgg                                                                            >  5:149639/1‑149 (MQ=255)
                                                                gCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTA‑CAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAgg                                                                            <  6:197707/149‑1 (MQ=255)
                                                                       ggAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTA‑CAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAAttt                                                                     >  7:198001/1‑149 (MQ=255)
                                                                        gAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTA‑CAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGccc                                                                                       <  6:54526/130‑1 (MQ=255)
                                                                        gAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTA‑CAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGccc                                                                                       >  5:54526/1‑130 (MQ=255)
                                                                          agagTTCGAGCACTCGATCTCCGACCTCGAGTA‑CAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCg                                                                   >  4:284977/1‑148 (MQ=255)
                                                                             gTTCGAGCACTCGATCTCCGACCTCGAGTA‑CAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAgg                                                               <  6:153606/149‑1 (MQ=255)
                                                                                 gatcACTCGATCTCCGACCTCGAGTA‑CAACTGGTACCTGCATACCGTGTAAACGATTACGGCGCCCCTGTAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCGTCTTCGCGGGCACGCCCGCTCCCGc                                                                               <  4:96858/126‑1 (MQ=255)
                                                                                 gAGCACTCGATCTCCGACCTCGAGTA‑CAACTGGTACCTGCATACAGTCTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTCTTTGAGACCGGTGTCGCCCTCTACGCGGGCACGCCCGCTCCCGc                                                                               >  3:96858/1‑129 (MQ=255)
                                                                                   gCACTCGATCTCCGACCTCGAGTA‑CAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAg                                                         >  3:205030/1‑149 (MQ=255)
                                                                                    cACTCGATCTCCGACCTCGAGTA‑CAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGc                                                        <  2:217878/149‑1 (MQ=255)
                                                                                      cTCGATCTCCGACCTCGAGTA‑CAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTcgcg                                                                                                <  4:28498/107‑1 (MQ=255)
                                                                                      cTCGATCTCCGACCTCGAGTA‑CAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTcgcg                                                                                                >  3:28498/1‑107 (MQ=255)
                                                                                      cTCGATCTCCGACCTCGAGTA‑CAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCt                                                      >  4:61787/1‑149 (MQ=255)
                                                                                       tCGATCTCCGACCTCGAGTA‑CAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCACCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCGGAAGCCtg                                                     >  7:232610/1‑149 (MQ=255)
                                                                                         gATCTCCGACCTCGAGTA‑CAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCtgt                                                    <  4:205030/148‑1 (MQ=255)
                                                                                              ccGACCTCGAGTA‑CAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGa                                                  <  3:270999/145‑1 (MQ=255)
                                                                                              ccGACCTCGAGTA‑CAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGa                                                  >  4:270999/1‑145 (MQ=255)
                                                                                               cGACCTCGAGTA‑CAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGa                                                  <  3:8724/144‑1 (MQ=255)
                                                                                               cGACCTCGAGTA‑CAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGAc                                                 <  8:131724/145‑1 (MQ=255)
                                                                                               cGACCTCGAGTA‑CAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACtt                                               <  6:121695/147‑1 (MQ=255)
                                                                                                   cTCGAGTA‑CAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCtgt                                         <  3:274216/149‑1 (MQ=255)
                                                                                                   cTCGAGTA‑CAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCtgt                                         <  6:149639/149‑1 (MQ=255)
                                                                                                    tCGAGTA‑CAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCtgtg                                        >  7:249462/1‑149 (MQ=255)
                                                                                                    tCGAGTA‑CAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCtgtg                                        >  6:341728/1‑149 (MQ=255)
                                                                                                            cAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCCCTACCGCAGGGTAATTTGGAAGGCTGAAGCCGGTGTCTTACCGGTGGGAGCgg                                 >  7:244713/1‑149 (MQ=255)
                                                                                                             aaCTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGggc                                                                                                                                     >  7:8765/1‑48 (MQ=255)
                                                                                                             aaCTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGggc                                                                                                                                     <  8:8765/48‑1 (MQ=255)
                                                                                                             aaCTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTGGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCggg                                >  5:139030/1‑149 (MQ=255)
                                                                                                             aaCTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGAGGGAGCggg                                <  7:16623/149‑1 (MQ=255)
                                                                                                              aCTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCGGGc                               <  7:5378/149‑1 (MQ=255)
                                                                                                              aCTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCGGGc                               <  4:316255/149‑1 (MQ=255)
                                                                                                              aCTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCGGGc                               <  4:289443/149‑1 (MQ=255)
                                                                                                                  gTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCtgt                                                    <  3:28562/124‑1 (MQ=255)
                                                                                                                  gTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCtgt                                                    >  4:28562/1‑124 (MQ=255)
                                                                                                                  gTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCGGGCGt                             <  1:188387/147‑1 (MQ=255)
                                                                                                                       tGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCGGGCGTGCCcgc                       >  1:263344/1‑148 (MQ=255)
                                                                                                                                  caaaCAAAAACGCCGCCCTTGAGGGGCGGCGTTTTCGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCGGGc                               <  4:234908/128‑1 (MQ=38)
                                                                                                                                  caaaCAAAAACGCCGCCCTTGAGGGGCGGCGTTTTCGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCGGGc                               >  3:234908/2‑129 (MQ=38)
                                                                                                                                   aaaCAAAAACGCCGCCCTTGAGGGGCGGCGTTTTCGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTTGGAGCGGGCGTGCCCGCGATGACGGCGGcc          >  6:89534/1‑149 (MQ=38)
                                                                                                                                       gAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCGGTGAGGTACCTGTGGGAGCGGGCGTGCCCGCGTAGAGGGCGGCCCAgg      >  8:61722/1‑149 (MQ=255)
                                                                                                                                       gAAAACGCCGCCCCTCAAGGGCAGAGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCACCCGCAGGGGAAATTAGAAGGGGGAAGGCGGTGACATACCGGGGGGAGCGGGGGTGCCCGCGGAGAGGGGGGCCCAgg      >  7:52185/1‑149 (MQ=255)
                                                                                                                                          aaCGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGGCGGCCCAGGCt    >  5:146569/1‑148 (MQ=255)
                                                                                                                                          aaCGCCGCCCCCCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGGCGGCCCAGGCt    >  2:58476/1‑148 (MQ=255)
                                                                                                                                            cgccgcCCTTGAGGGGCGGCGTTTTCGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGGCGGCCCAGGCTTa  >  4:257882/1‑148 (MQ=255)
                                                                                                                                            cgccgcCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCAATTACCCGGCAACCCCCCAACCCCAAGGGAATTTAGAAAGCATAAACCATTGACTTTACTGTTGGAGCGGGCGTTCCCGCCAAAAAGGCCGCCCAAGCATa  >  2:64330/1‑148 (MQ=255)
                                                                                                                                             gccgccCTTGAGGGGCGGCGTTTTCGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTAc                                             >  5:150419/1‑104 (MQ=38)
                                                                                                                                             gccgccCTTGAGGGGCGGCGTTTTCGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTAc                                             <  6:150419/104‑1 (MQ=38)
                                                                                                                                                    |                                                                                                                                           
AGCGAGATCCTGAACAAGTACATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTA‑CAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGGCGGCCCAGGCTTA  >  NC_002947/5913977‑5914263

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