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

GCGAGATCCTGAACAAGTACATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCC‑‑TCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTC‑CCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCGGGCGTGCCCG  >  NC_002947/5913978‑5914241
                                                                                                                                                      |                                                                                                                    
gCGAGATCCTGAACAAGTACATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCC‑‑Tc                                                                                                                      <  7:111096/149‑1 (MQ=255)
  gagaTCCTGAACAAGTACATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCC‑‑TCaa                                                                                                                    <  6:302990/149‑1 (MQ=255)
   agaTCCTGAACAAGTACATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCC‑‑TCAAg                                                                                                                   >  8:106391/1‑149 (MQ=255)
   agaTCCTGAACAAGTACATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCC‑‑TCAAg                                                                                                                   <  2:238114/149‑1 (MQ=255)
         tGAACAAGTACATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCC‑‑TCAAGggcggc                                                                                                             >  1:295346/1‑149 (MQ=255)
         tGAACAAGTACATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCC‑‑TCAAGggcggc                                                                                                             >  8:263948/1‑149 (MQ=255)
          gAACAAGTACATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCC‑‑TCAAGGgcggcg                                                                                                            >  1:69044/1‑149 (MQ=255)
           aaCAAGTACATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCC‑‑TCAAGGGCGGCGt                                                                                                           >  4:35932/1‑149 (MQ=255)
              aaGTACATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCC‑‑TCAAGGGCGGCGtttt                                                                                                        >  7:308868/1‑149 (MQ=255)
                   cATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCC‑‑TCAAGGGCGGCGTTtttgttt                                                                                                   <  2:69044/149‑1 (MQ=255)
                     tCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCC‑‑TCAAGGGCGGCGTTttcgtttac                                                                                                 <  7:94949/149‑3 (MQ=255)
                     tCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCC‑‑TCAAGGGCGGCGTTTTTGTTTga                                                                                                 <  2:291953/149‑1 (MQ=255)
                        aTCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCC‑‑TCAAGGGCGGCGTTtttgttt                                                                                                   <  6:316062/144‑1 (MQ=255)
                        aTCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCC‑‑TCAAGGGCGGCGTTtttgttt                                                                                                   >  5:316062/1‑144 (MQ=255)
                        aTCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCC‑‑TCAAGGGCGGCGTTTTTGTTTGAGAc                                                                                              <  3:58911/149‑1 (MQ=255)
                        aTCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCC‑‑TCAAGGGCGGCGTTTTTGTTTGAGAc                                                                                              <  8:1485/149‑1 (MQ=255)
                         tCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCC‑‑TCAAGGGCGGCGTTTTTGTTTGAGAcc                                                                                             >  4:286312/1‑149 (MQ=255)
                             aaGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACAAAAACGCCGCCCT‑‑TGAGGGGCGGCGttttc                                                                                                       <  2:119246/135‑2 (MQ=255)
                             aaGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACAAAAACGCCGCCCT‑‑TGAGGGGCGGCGttttc                                                                                                       >  1:119246/1‑134 (MQ=255)
                               gTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCC‑‑TCAAGGGCGGCGttt                                                                                                         <  2:325843/131‑1 (MQ=255)
                               gTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCC‑‑TCAAGGGCGGCGttt                                                                                                         >  1:325843/1‑131 (MQ=255)
                                 aCATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCC‑‑TCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGc                                                                                     >  5:155938/1‑149 (MQ=255)
                                      gACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCC‑‑TCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCctct                                                                                <  3:35932/149‑1 (MQ=255)
                                         aTCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGACTTCGAGCACTCGATCTCAGACATCGAGTACAACTGGTAACTGCATACCGTGTAAACGAAAACGCCGACCA‑‑TCAAGGGCGGCGTTTTTGTTAGAGACAGGTGTCGCACTCTTcg                                                                             >  8:78028/1‑149 (MQ=255)
                                              cgtcgCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCC‑‑TCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCgg                                                                          <  2:73923/147‑1 (MQ=255)
                                               gtcgCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCC‑‑TCAAGGGCGGCGTTtttgttt                                                                                                   >  8:206932/1‑121 (MQ=255)
                                               gtcgCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCC‑‑TCAAGGGCGGCGTTtttgttt                                                                                                   <  7:206932/121‑1 (MQ=255)
                                                   cgTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCC‑‑TCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGcc                                                                   >  4:285583/1‑149 (MQ=255)
                                                     tGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTAGAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCC‑‑TCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCg                                                                 <  1:277021/149‑1 (MQ=255)
                                                        aaGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCC‑‑TCCAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTc                                                              <  3:50200/149‑1 (MQ=255)
                                                        aaGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCC‑‑TCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTc                                                              >  7:36625/1‑149 (MQ=255)
                                                         aGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCC‑‑TCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTc‑c                                                            <  7:106391/149‑1 (MQ=255)
                                                         aGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCC‑‑TCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTc‑c                                                            >  6:15900/1‑149 (MQ=255)
                                                          ggAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCC‑‑TCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTc‑cc                                                           >  5:167340/1‑149 (MQ=255)
                                                             aaGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCC‑‑TCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTC‑CCGCa                                                        <  5:15900/149‑1 (MQ=255)
                                                              aGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCC‑‑TCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTC‑CCACAg                                                       <  2:246294/149‑1 (MQ=255)
                                                                cGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCC‑‑TCAAGGGCGGCGTTTTTGTTTGAGACCGGTGCCGCCCTCTTCGCGGGCACGCCCGCTC‑CCACAgg                                                      >  5:231785/1‑148 (MQ=255)
                                                                   gCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCC‑‑TCAAGGGCGGCGTTTTTGTTTGAGACCGgtg                                                                                         <  4:107819/111‑1 (MQ=255)
                                                                   gCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCC‑‑TCAAGGGCGGCGTTTTTGTTTGAGACCGgtg                                                                                         >  3:107819/1‑111 (MQ=255)
                                                                            gTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCC‑‑TCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTC‑CCGCAGGGTAATTTCGAAgg                                         <  6:167340/149‑1 (MQ=255)
                                                                                        gATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCGATCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTC‑CCGCAGGGTAATTTCGAAGGCTGAAGCCtg                               >  1:97988/1‑149 (MQ=255)
                                                                                        gATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCC‑‑TCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCTCCGCAGGGTAATTTCGAAGGCTGAAGCCt                                <  1:26422/147‑1 (MQ=255)
                                                                                        gATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCC‑‑TCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTC‑CCGCAGGGTAATTTCGAAGGCTGAAGCCtg                               <  3:286312/147‑1 (MQ=255)
                                                                                             ccGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCC‑‑TCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTC‑CCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTa                        >  3:309970/1‑149 (MQ=255)
                                                                                              cGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCC‑‑TCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTC‑CCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGa                            <  2:295346/144‑1 (MQ=255)
                                                                                              cGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCC‑‑TCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTC‑CCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGa                            <  1:316320/144‑1 (MQ=255)
                                                                                              cGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCC‑‑TCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTC‑CCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGa                            <  6:155938/144‑1 (MQ=255)
                                                                                                     gAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCC‑‑TCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTC‑CCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTgg                 <  7:263948/148‑1 (MQ=255)
                                                                                                     gAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCC‑‑TCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTC‑CCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTgg                 <  8:308868/148‑1 (MQ=255)
                                                                                                      aGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCACC‑‑TCAAGGGCTGCGTTTTTGTTTGAGAGCTGTGTCGCCCTCTTCGCTGGCACGCCCTCTC‑CCGAAGTTTAATTTCGAAGGCTGAAGCCTGTGACTTACCTTTGGGa               <  7:78028/149‑1 (MQ=255)
                                                                                                          cAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCC‑‑TCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTG‑CCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCgg           <  2:97988/149‑1 (MQ=255)
                                                                                                            aCTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCC‑‑TCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTC‑CCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCGGGc         <  3:285583/149‑1 (MQ=255)
                                                                                                                gTACCTGCATACCGTGTAAACGAAAACGCCGCCCC‑‑TCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTC‑CCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCGGGCGTGc     <  4:309970/149‑1 (MQ=255)
                                                                                                                   ccTGCATACCGTGTAAACGAAAACGCCGCCCC‑‑TCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTC‑CCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCGGGCGTGCCcg  <  8:36625/149‑1 (MQ=255)
                                                                                                                    cTGCATACCGTGTAAACGAAAACGCCGCCCC‑‑TCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTC‑CCGCAgg                                                      >  4:140788/1‑96 (MQ=255)
                                                                                                                    cTGCATACCGTGTAAACGAAAACGCCGCCCC‑‑TCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTC‑CCGCAgg                                                      <  3:140788/96‑1 (MQ=255)
                                                                                                                      gCATACCGTGTAAACGAAAACGCCGCCCC‑‑TCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTC‑CCGCAGGGTAATTTCGAAGGCt                                       <  2:228055/109‑1 (MQ=255)
                                                                                                                      gCATACCGTGTAAACGAAAACGCCGCCCC‑‑TCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTC‑CCGCAGGGTAATTTCGAAGGCt                                       >  1:228055/1‑109 (MQ=255)
                                                                                                                                     gAAAACGCCGCCCC‑‑TCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTC‑CCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCtgt                   <  4:136493/114‑1 (MQ=255)
                                                                                                                                     gAAAACGCCGCCCC‑‑TCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTC‑CCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCtgt                   >  3:136493/1‑114 (MQ=255)
                                                                                                                                           gccgccCT‑‑TGAGGGGCGGCGTTTTCGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTC‑CCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCGGGc         >  4:88908/1‑118 (MQ=39)
                                                                                                                                           gccgccCT‑‑TGAGGGGCGGCGTTTTCGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTC‑CCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCGGGc         <  3:88908/118‑1 (MQ=39)
                                                                                                                                                      |                                                                                                                    
GCGAGATCCTGAACAAGTACATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCC‑‑TCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTC‑CCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCGGGCGTGCCCG  >  NC_002947/5913978‑5914241

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