Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,833,194 C→G 6.2% intergenic (‑4/+60) pvdY ← / ← PP_4246 hydroxyproline acetylase/substrate‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,833,1940CG 6.2% 205.7 / 2.4 64intergenic (‑4/+60)pvdY/PP_4246hydroxyproline acetylase/substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (29/31);  new base G (2/2);  total (31/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

TGCAGAGAAGCGCCAGCCTCCAGTCGCAACTGGATCAGGGGGCGGCCTTCAAGGGCCAGGCTTAGCCAGCTGTCGCCCTCCTCGGCACTGAGGCTGCGCCAGCGTGGCGGTGACTGTGGCTGCGAAGCGGCATCGAACGGCATGGGCTTGACTCACGAAAATGAAGAAAAGGCACACTGCAGGGAACGTTGCCGCAACACGGGAATTTAGTCTTGCGGTGCAGCCACTTTGATTCGGTACGGCTCGAAGATGCGCGCCAGTTCACCGCTGTCACGCAGGCCACGCAGCAGG  >  NC_002947/4833048‑4833338
                                                                                                                                                  |                                                                                                                                                
tGCAGAGAAGCGCCAGCCTCCAGTCGCAACTGGATCAGGGGGCGGCCTTCAAGGGCCAGGCTTAGTCAGCTGTCGCCCTCCTCGGCACTGAGGCTGCGCCAGCGTGGCGGTGACTGTGGCTGCGAAGCGGCATCGAACGGCATGGGCtt                                                                                                                                                <  7:11316/149‑1 (MQ=255)
  cAGAGAAGCGCCAGCCTCCAGTCGCAACTGGATCAGGGGGCGGCCTTCAAGGGCCAGGCTTAGCCAGCTGTCGCCCTCCTCGGCACTGAGGCTGCGCCAGCGTGGCGGTGACTGTGGCTGCGAAGCGGCATCGAACGGCATGGGCTTGa                                                                                                                                              <  5:280538/149‑1 (MQ=255)
   agagAAGCGCCAGCCTCCAGTCGCAACTGGATCAGGGGGCGGCCTTCAAGGGCCAGGCTTAGCCAGCTGTCGCCCTCCTCGGCACTGAGGCTGCGCCAGCGTGGCGGTGACTGTGGCTGCGAAGCGGCATCGAACGGCATGGGCTTGAc                                                                                                                                             >  5:228065/1‑149 (MQ=255)
   agagAAGCGCCAGCCGCCAGGCGCAACTGGATCAGGGGGCGGCCTTCAAGGGCCAGGCGTAGCCAGCTGTCGCCCTCCTCGGCACTGAGGCTGCGCCAGCGTGGCGGTGACTGTGGCTGCGAAGCGGCATCGAACGGCATGGGCTTGAc                                                                                                                                             <  4:272230/149‑1 (MQ=255)
    gagaAGCGCCAGCCTCCAGTCGCAACTGGATCAGGGGGCGGCCTTCAAGGGCCAGGCTTAGCCAGCTGTCGCCCTCCTCGGCACTGAGGCTGCGCCAGCGTGGCGGTGACTGTGGCTGCGAAGCGGCATCGAACGGCATGGGCTTGACt                                                                                                                                            >  2:17807/1‑149 (MQ=255)
     agaAGCGCCAGCCTCCAGTCGCAACTGGATCAGGGGGCGGCCTTCAAGGGCCAGGCTTAGCCAGCTGTCGCCCTCCTCGGCACTGAGGCTGCGCCAGCGTGGCGGTGACTGTGGCTGCGAAGCGGCATCGAACGGCATGGGCTTGACt                                                                                                                                            >  8:172663/1‑148 (MQ=255)
              aGCCTCCAGTCGCAACTGGATCAGGGGGCGGCCTTCAAGGGCCAGGCTTAGCCAGCTGTCGCCCTCCTCGGCACTGAGGCTGCGCCAGCGTGGCGGTTACTGTGGCTGCGAAGCGGCATCGAACGGCATGGGCTTGACTCACGAAAATg                                                                                                                                  <  5:14854/149‑1 (MQ=255)
                gatCCAGTTGCGACTGGATCAGGGGGCGGCCTTCAAGGGCCAGGCTTAGCCAGCTGTCGCCCTCCTCGGCACTGAGGCTGCGCCAGCGTGGCGGTGACTGTGGCTGCGAAGCGGCATCGAACGGCATGGGCTTGACTCACGAAAATg                                                                                                                                  >  6:67928/3‑147 (MQ=255)
                ccTCCAGTCGCAACTGGATCAGGGGGCGGCCTTCAAGGGCCAGGCTTAGCCAGCTGTCGCCCTCCTCGGCACTGAGGCTGCGCCAGCGTGGCGGTGACTGTGGCTGCGAAGCGGCATCGAACGGCATGGGCTTGACTCACGAAAATgaa                                                                                                                                <  7:286325/149‑1 (MQ=255)
                ccTCCAGTCGCAACTGGATCAGGGGGCGGCCTTCAAGGGCCAGGCTTAGCCAGCTGTCGCCCTCCTCGGCACTGAGGCTGCGCCAGCGTGGCGGTGACTGTGGCTGCGAAGCGGCATCGAACGGCATGGGCTTGACTCACGAAAATgaa                                                                                                                                <  4:237837/149‑1 (MQ=255)
                ccTCCAGTCGCAACTGGATCAGGGGGCGGCCTTCAAGGGCCAGGCTTAGCCAGCTGTCGCCCTCCTCGGCACTGAGGCTGCGCCAGCGTGGCGGTGACTGTGGCTGCGAAGCGGCATCGAACGGCATGGGCTTGACTCACGAAAATgaa                                                                                                                                <  6:228065/149‑1 (MQ=255)
                ccTCCAGTCGCAACTGGATCAGGGGGCGGCCTTCAAGGGCCAGGCTTAGCCAGCTGTCGCCCTCCTCGGCACTGAGGCTGCGCCAGCGTGGCGGTGACTGTGGCTGCGAAGCGGCATCGAACGGCATGGGCTTGACTCACGAAAATg                                                                                                                                  >  2:303505/1‑147 (MQ=255)
                 cTCCAGTCGCAACTGGATCAGGGGGCGGCCTTCAAGGGCCAGGCTTAGCCAGCTGTCGCCCTCCTCGGCACTGAGGCTGCGCCAGCGTGGCGGTGACTGTGGCTGCGAAGCGGCATCGAACGGCATGGGCTTGACTCACGAAAATgaa                                                                                                                                <  1:134303/148‑1 (MQ=255)
                     aGTCGCAACTGGATCAGGGGGCGGCCTTCAAGGGCCAGGCTTAGCCAGCTGTCGCCCTCCTCGGCACTGAGGCTGCGCCAGCGTGGCGGTGACTGTGGCTGCGAAGCGGCATCGAACGGCATGGGCTTGACTCACGAAAATgaagaa                                                                                                                             >  8:172042/1‑147 (MQ=255)
                         gCAACTGGATCAGGGGGCGGCCTTCAAGGGCCAGGCTTAGTGCACTGTCGCCCTCCTCGGCACTGAGGCTGCGCCAGCGTGGCGGTGACTGTGGCTGCGAAGCGGCATCGAACGGCATGGGCTTGACTCACGAAAATGAAGAAAAGGc                                                                                                                        <  1:316250/148‑1 (MQ=255)
                                gATCAGGGGGCGGCCTTCAAGGGCCAGGCTTAGCCAGCTGTCGCCCTCCTCGGCACTGAGGCTGCGCCAGCGTGGCGGTGACTGTGGCTGCGAAGCGGCATCGAACGGCATGGGCTTGACTCACGAAAATGAAGAAAAGGCACACTGc                                                                                                                 <  5:67928/148‑1 (MQ=255)
                                gATCAGGGGGCGGCCTTCAAGGGCCAGGCTTAGCCAGCTGTCGCCCTCCTCGGCACTGAGGCTGCGCCAGCGTGGCGGTGACTGTGGCTGCGAAGCGGCATCGAACGGCATGGGCTTGACTCACGAAAATGAAGAAAAGGCACACTGc                                                                                                                 <  5:324523/148‑1 (MQ=255)
                                    aGGGGGCGGCCTTCAAGGGCCAGGCTTAGCCAGCTGTCGCCCTCCTCGGCACTGAGGCTGCGCCAGCGTGGCGGTGACTGTGGCTGCGAAGCGGCATCGAACGGCATGGGCTTGACTCACGAAAATGAAGAAAAGGCACACTGCAggg                                                                                                             >  8:182853/1‑148 (MQ=255)
                                     gggggCGGCCTTCAAGGGCCAGGCTTAGGCAGCTGTCGCCCTCCTCGGCACTGAGGCTGCGCCAGCGTGGCGGTGACTGTGGCTGCGAAGCGGCATCGAACGGCATGGGCTTGACTCACGAAAATGAAGAAAAGGCACACTGCAGGGaa                                                                                                           >  1:214362/1‑149 (MQ=255)
                                       gggCGGCCTTCAAGGGCCAGGCTTAGCCAGCTGTCGCCCTCCTCGGCACTGAGGCTGCGCCAGCGTGGCGGTGACTGTGGCTGCGAAGCGGCATCGAACGGCATGGGCTTGACTCACGAAAATGAAGAAAAGGCACACTGCAGGGAACg                                                                                                         >  4:366496/1‑149 (MQ=255)
                                           ggcCTTCAAGGGCCAGGCTTAGCCAGCTGTCGCCCTCCTCGGCACTGAGGCTGCGCCAGCGTGGCGGTGACTGTGGCTGCGAAGCGGCATCGAACGGCATGGGCTTGACTCACGAAAATGAAGAAAAGGCACACTGCAGGGAACGTTGc                                                                                                     <  8:29090/149‑1 (MQ=255)
                                            gcCTTCAAGGGCCAGGCTTAGCCAGCTGTCGCCCTCCTCGGCACTGAGGCTGCGCCAGCGTGGCGGTGACTGTGGCTGCGAAGCGGCATCGAACGGCATGGGCTGGACTCACGAAAATGAAGAAAAGGCACACTGCAGGGAACGTTGcc                                                                                                    <  1:268097/149‑1 (MQ=255)
                                             ccTTCAAGGGCCAGGCTTAGCCAGCTGTCGCCCTCCTCGGCACTGAGGCTGCGCCAGCGTGGCGGTGACTGTGGCTGCGAAGCGGCATCGAACGGCATGGGCTTGACTCACGAAAATGAAGAAAAGGCACACTGCAGGGAACGTTGc                                                                                                     <  7:172663/147‑1 (MQ=255)
                                                 cAAGGGCCAGGCTTAGCCAGCTGTCGCCCTCCTCGGCACTGAGGCTGCGCCAGCGTGGCGGTGACTGTGGCTGCGAAGCGGCATCGAACGGCATGGGc                                                                                                                                                  <  8:175538/98‑1 (MQ=255)
                                                 cAAGGGCCAGGCTTAGCCAGCTGTCGCCCTCCTCGGCACTGAGGCTGCGCCAGCGTGGCGGTGACTGTGGCTGCGAAGCGGCATCGAACGGCATGGGc                                                                                                                                                  >  7:175538/1‑98 (MQ=255)
                                                 cAAGGGCCAGGCTTAGCCAGCTGTCGCCCTCCTCGGCACTGAGGCTGCGCCAGCGTGGCGGTGACTGTGGCTGCGAAGCGGCATCGAACGGCATGGGCTTGACTCACGAAAATGAAGAAAAGGCACACTGCAGGGAACGTTGCCGCaa                                                                                                >  7:109668/1‑148 (MQ=255)
                                                     ggCCAGGCTTAGCCAGCTGTCGCCCTCCTCGGCACTGAGGCTGCGCCAGCGTGGCGGTGACTGTGGCTGCGAAGCGGCATCGAACGGCATGGGCTTGACTCACGAAAATGAAGAAAAGGCCCCCTGCAGGGAACGTTGCCGCAACACgg                                                                                           >  3:125579/1‑149 (MQ=255)
                                                       ccAGGCTTAGCCAGCTGTCGCCCTCCTCGGCACTGAGGCTGCGCCAGCGTGGCGGTGACTGTGGCTGCGAAGCGGCATCGAACGGCATGGGCTTGACTCACGAAAATGAAGAAAAGGCACACTGCAGGGAAc                                                                                                          >  8:240839/1‑132 (MQ=255)
                                                       ccAGGCTTAGCCAGCTGTCGCCCTCCTCGGCACTGAGGCTGCGCCAGCGTGGCGGTGACTGTGGCTGCGAAGCGGCATCGAACGGCATGGGCTTGACTCACGAAAATGAAGAAAAGGCACACTGCAGGGAAc                                                                                                          <  7:240839/132‑1 (MQ=255)
                                                       ccAGGCTTAGCCAGCTGTCGCCCTCCTCGGCACTGAGGCTGCGCCAGCGTGGCGGTGACTGTGGCTGCGAAGCGGCATCGAACGGCATGGGCTTGACTCACGAAAATGAAGAAAAGGCACACTGCAGGGAACGTTGCCGCAACACGGGa                                                                                         <  1:17807/149‑1 (MQ=255)
                                                           gCTTAGCCAGCTGTCGCCCTCCTCGGCACTGAGGCTGCGCCAGCGTGGCGGTGACTGTGGCTGCGAAGCGGCATCGAACGGCATGGGCTTGACTCACGAAAATGAAGAAAAGGCACACTGCAGGGAACGTTGCCGCAACACGGGAAtt                                                                                      >  3:105377/1‑148 (MQ=255)
                                                            cTTAGCCAGCTGTCGCCCTCCTCGGCACTGAGGCTGCGCCAGCGTGGCGGTGACTGTGGCTGCGAAGCGGCATCGAACGGCATGGGCTTGACTCACGAAAATGAAGAAAAGGCACACt                                                                                                                   >  5:341842/1‑118 (MQ=255)
                                                            cTTAGCCAGCTGTCGCCCTCCTCGGCACTGAGGCTGCGCCAGCGTGGCGGTGACTGTGGCTGCGAAGCGGCATCGAACGGCATGGGCTTGACTCACGAAAATGAAGAAAAGGCACACt                                                                                                                   <  6:341842/118‑1 (MQ=255)
                                                            cTTAGCCAGCTGTCGCCCTCCTCGGCACTGAGGCTGCGCCAGCGTGGCGGTGACTGTGGCTGCGAAGCGGCATCGAACGGCATGGGCTTGACTCACGAAAATGAAGAAAAGGCACACTGCAGGGAACGTTGCCGCAACACGTGAATTTa                                                                                    <  4:105377/149‑1 (MQ=255)
                                                               aGCCAGCTGTCGCCCTCCTCGGCACTGAGGCTGCGCCAGCGTGGCGGTGACTGTGGCTGCGAAGCGGCATCGAACGGCATGGGCTTGACTCACGAAAATGAAGAAAAGGCACACTGCAGTGAACGTTGCCGCAACACGGGAATTTAGt                                                                                  <  1:22662/148‑1 (MQ=255)
                                                                 ccAGCTGTCGCCCTCCTCGGCACTGAGGCTGCGCCAGCGTGGCGGTGACTGTGGCTGCGAAGCGGCATCGAACGGCATGGGCTTGACTCACGAAAATGAAGAAAAGGCACACTGCAGGGAACGTTGCCGCAACACGGGAATTTAGTCtt                                                                               <  4:206906/149‑1 (MQ=255)
                                                                 ccAGCTGTCGCCCTCCTCGGCACTGAGGCTGCGCCAGCGTGGCGGTGACTGTGGCTGCGAAGCGGCATCGAACGGCATGGGCTTGACTCACGAAAATGAAGAAAAGGCACACTGCAGGGAACGTTGCCGCAACACGGGAATTTAGTCtt                                                                               >  4:222318/1‑149 (MQ=255)
                                                                        tCGCCCTCCTCGGCACTGAGGCTGCGCCAGCGTGGCGGTGACTGTGGCTGCGAAGCGGCATCGAACGGCATGGGCTTGACTCACGAAAATGAAGAAAAGGCACACTGCAGGGAACGTTGCCGCAACACGGGAATTTAGTCTTGCGGTg                                                                         >  1:43671/1‑148 (MQ=255)
                                                                           cccTCCTCGGCACTGAGGCTGCGCCAGCGTGGCGGTGACTGTGGCTGCGAAGCGGCATCGAACGGCATGGGCTTGACTCACGAAAATGAAGAAAAGGCACACTGCAGGGAACGTTGCCGCAACACGGGAATTTAGTCTTGCGGTGCAGc                                                                     >  2:58000/1‑149 (MQ=255)
                                                                           cccTCCTCGGCACTGAGGCTGCGCCAGCGTGGCGGTGACTGTGGCTGCGAAGCGGCATCGAACGGCATGGGCTTGACTCACGAAAATGAAGAAAAGGCACACTGCAGGGAACGTTGCCGCAACACGGGAATTTAGTCTTGCGGTGCAGc                                                                     <  7:352359/149‑1 (MQ=255)
                                                                              tcctcGGCACTGAGGCTGCGCCAGCGTGGCGGTGACTGTGGCTGCGAAGCGGCATCGAACGGCATGGGGTTGACTCACGAAAATGAAGAAAAGGCACACTGCAGGGAACGTTGCCGCAACACGGGAATTTAGTCTTGCGGTGCAGCCAc                                                                  >  8:38931/1‑149 (MQ=255)
                                                                              tcctcGGCACTGAGGCTGCGCCAGCGTGGCGGTGACTGTGGCTGCGAAGCGGCATCGAACGGCATGGGGTTGACTCACGAAAATGAAGAAAAGGCACACTGCAGGGAACGTTGCCGCAACACGGGAATTTAGTCTTGCGGTGCAGCCAc                                                                  >  8:44791/1‑149 (MQ=255)
                                                                                ctcGGCACTGAGGCTGCGCCAGCGTGGCGGTGACTGTGGCTGCGAAGCGGCATCGAACGGCATGGGCTTGACTCACGAAAATGAAGAAAAGGCACACTGCAGGGAACGTTGCCGCAACACGGGAATTTAGTCTTGCGGTGCAGCCAc                                                                  <  1:38865/147‑1 (MQ=255)
                                                                                ctcGGCACTGAGGCTGCGCCAGCGTGGCGGGGACTGTGGCTGCGAAGCGGCATCGAACGGCATGGGCTTGACTCACGAAAATGAAGAAAAGGCACACTGCAGGGAACGTTGCCGCAACACGGGAATTTAGTCTTGCGGTGCAGCCAc                                                                  <  6:213835/147‑1 (MQ=255)
                                                                                 tcGGCACTGAGGCTGCGCCAGCGTGGCGGTGACTGTGGCTGCGAAGCGGCATCGAACGGCATGGGCTTGACTCACGAAAATGAAGAAAAGGCACACTGCAGGGAACGTTGCCGCAACACGGGAATTTAGTCTTGCGGTGCAGCCACttt                                                               >  7:32424/1‑149 (MQ=255)
                                                                                       cTGAGGCTGCGCCAGCGTGGCGGTGACTGTGGCTGCGAAGCGGCATCGAACGGCATGGGc                                                                                                                                                  <  6:25778/60‑1 (MQ=255)
                                                                                       cTGAGGCTGCGCCAGCGTGGCGGTGACTGTGGCTGCGAAGCGGCATCGAACGGCATGGGc                                                                                                                                                  >  5:25778/1‑60 (MQ=255)
                                                                                          aGGCTGCGCCAGCGTGGCGGTGACTGTGGCTGCGAAGCGGCATCGAACGGCATGGGCTTGACTCACGAAAATGAAGAAAAGGCACACTGCAGGGAACGTTGCCGCAACACGGGAATTTAGTCTTGCGGTGCAGCCACTTTGATTCGGTa                                                      <  3:222318/149‑1 (MQ=255)
                                                                                             cTGCGCCAGCGTGGCGGTGACTGTGGCTGCGAAGCGGCATCGAACGGCATGGGCTTGACTCACGAAAATGAAGAAAAGGCACACTGCAGGGAACGTTGCCGCAACACGGGAATTTAGTCTTGCGGTGCAGCCACTTTGATTCGGTACgg                                                   <  8:109668/149‑1 (MQ=255)
                                                                                              tGCGCCAGCGTGGCGGTGACTGTGGCTGCGAAGCGGCATCGAACGGCATGGGCTTGACTCACGAAAATGAAGAAAAGGCACACTGCAGGGAACGTTGCCGCAACACGGGAATTTAGTCTTGCGGTGCAGCCACTTTGATTCGGTACGGc                                                  >  1:71278/1‑149 (MQ=255)
                                                                                                  ccAGCGTGGCGGTGACTGTGGCTGCGAAGCGGCATCGAACGGCATGGGCTTGACTCACGAAAATGAAGAAAAGGCACACTGCAGGGAACGTTGCCGCAACACGGGAATTTAGTCTTGCGGTGCAGCCACTTTGATTCGGTACGGCTCg                                               >  7:246489/1‑148 (MQ=255)
                                                                                                   cAGCGTGGCGGTGACTGTGGCTGCGAAGCGGCATCGAACGGCATGGGGTTGACTCACGAAAATGAAGAAAAGGCACACTGCAGGGAACGTTGCCGCAACACGGGAATTTAGTCTTGCGGTGCAGCCACTTTGATTCGGTACGGCTCGaa                                             <  7:38931/149‑1 (MQ=255)
                                                                                                   cAGCGTGGCGGTGACTGTGGCTGCGAAGCGGCATCGAACGGCATGGGGTTGACTCACGAAAATGAAGAAAAGGCACACTGCAGGGAACGTTGCCGCAACACGGGAATTTAGTCTTGCGGTGCAGCCACTTTGATTCGGTACGGCTCGaa                                             <  7:44791/149‑1 (MQ=255)
                                                                                                          gCGGTGACTGTGGCTGCGAAGCGGCATCGAACGGCATGGGCTTGACTCACGAAAATGAAGAAAAGGCACACTGCAGGGAACGTTGCCGCAACACGGGAATTTAGTCTTGCGGTGCAGCCACTTTGATTCGGTACGGCTCGAAGATgcgc                                      >  7:138419/1‑149 (MQ=255)
                                                                                                                aCTGTGGCTGCGAAGCGGCATCGAACGGCATGGGCTTGACTCACGAAAATGAAGAAAAGGCACACTGCAGGGAACGTTGCCGCAACACGGGAATTTAGTCTTGCGGTGCAGCCACTTTGATTCGGTACGGCTCGAAGATGCGCGCCAg                                 <  1:58000/148‑1 (MQ=255)
                                                                                                                      gCTGCGAAGCGGCATCGAACGGCATGGGCTTGACTCACGAAAATGAAGAAAAGGCACACTGCAGGGAACGTTGCCGCAACACGGGAATTTAGTCTTGCGGTGCAGCCACTTTGATTCGGTACGGCTCGAAGATGCGCGCCAGTTCACCg                          <  1:303505/149‑1 (MQ=255)
                                                                                                                          cGAAGCGGCATCGAACGGCATGGGCTTGACTCACGAAAATGAAGAAAAGGCACACTGCAGGGAACGTTGCCGCAACACGGGAATTTAGTCTTGCGGTGCAGCCACTTTGATTCGGTACGGCTCGAAGATGCGCGCCAGTTCACCGCt                        <  7:182853/147‑1 (MQ=255)
                                                                                                                          cGAAGCGGCATCGAACGGCATGGGCTTGACTCACGAAAATGAAGAAAAGGCACACTGCAGGGAACGTTGCCGCAACACGGGAATTTAGTCTTGCGGTGCAGCCACTTTGATTCGGTACGGCTCGAAGATGCGCGCCAGTTCACCGCTGt                      >  5:86327/1‑149 (MQ=255)
                                                                                                                          cGAAGCGGCATCGAACGGCATGGGCTTGACTCACGAAAATGAAGAAAAGGCACACTGCAGGGAACGTTGCCGCAACACGGGAATTTAGTCTTGCGGTGCAGCCACTTTGATTCGGTACGGCTCGAAGATGCGCGCCAGTTCACCGCTGt                      >  2:119673/1‑149 (MQ=255)
                                                                                                                           gAAGCGGCATCGAACGGCATGGGCTTGACTCACGAAAATGAAGAAAAGGCACACTGCAGGGAACGTTGCCGCAACACGGGAATTTAGTCTTGCGGTGCAGCCACTTTGATTCGGTACGGCTCGAAGATGCGCGCCAGTTCACCGCTGt                      >  7:304258/1‑148 (MQ=255)
                                                                                                                                   aTCGAACGGCATGGGCTTGACTCACGAAAATGAAGAAAAGGCACACTGCAGGGAACGTTGCCGCAACACGGGAATTTAGTCTTGCGGTGCAGCCACTTTGATTCGGTACGGCTCGAAGATGCGCGCCAGTTCACCGCTGTCACGCAGGc             >  8:24903/1‑149 (MQ=255)
                                                                                                                                      gAACGGCATGGGCTTGACTCACGAAAATGAAGAAAAGGCACACTGCAGGGAACGTTGCCGCAACACGGGAATTTAGTCTTGCGGTGCAGCCACTTTGATTCGGTACGGCTCGAAGATGCGCGCCAGTTCACCGCTGTCACGCAGGCCa           <  7:172042/148‑1 (MQ=255)
                                                                                                                                       aaCGGCATGGGCTTGACTCACGAAAATGAAGAAAAGGCACACTGCAGGGAACGTTGCCGCAACACGGGAATTTAGTCTTGCGGTGCAGCCACTTTGATTCGGTACGGCTCGAAGATGCGCGCCAGTTCACCGCTGTCACGCAGGCCACg         >  5:78488/1‑149 (MQ=255)
                                                                                                                                         cGGCATGGGCTTGACTCACGAAAATGAAGAAAAGGCACACTGCAGGGAACGTTGCCGCAACACGGGAATTTAGTCTTGCGGTGCAGCCACTTTGATTCGGTACGGCTCGAAGATGCGCGCCAGTTCACCGCTGTCACGCAGGCCACg         >  6:5068/1‑147 (MQ=255)
                                                                                                                                           gCATGGGCTTGACTCACGAAAATGAAGAAAAGGCACACTGCAGGGAACGTTGCCGCAACACGGGAATTTAGTCTTGCGGTGCAGCCACTTTGATTCGGTACGGCTCGAAGATGCGCGCCAGTTCACCGCTGTCACGCAg               <  7:77077/139‑1 (MQ=255)
                                                                                                                                           gCATGGGCTTGACTCACGAAAATGAAGAAAAGGCACACTGCAGGGAACGTTGCCGCAACACGGGAATTTAGTCTTGCGGTGCAGCCACTTTGATTCGGTACGGCTCGAAGATGCGCGCCAGTTCACCGCTGTCACGCAGGCCACgcag      >  5:319298/1‑148 (MQ=255)
                                                                                                                                           gCATGGGCTTGACTCACGAAAATGAAGAAAAGGCACACTGCAGGGAACGTAGCCGCAACACGGGAATTTAGTCTTGCGGTGCAGCCACTTTGATTCGGTACGGCTCGAAGATGCGCGCCAGTTCACCGCTGTCACGCAg               >  8:77077/1‑139 (MQ=255)
                                                                                                                                              tGGGCTTGACTCACGAAAATGAAGAAAAGGCACACTGCAGGGAACGTTGCCGCAACACGGGAATTTAGTCTTGCGGTGCAGCCACTTTGATTCGGTACGGCTCGAAGATGCGCGCCAGTTCACCGCTGTCACGCAGGCCACGCAGCAgg  <  2:43671/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
TGCAGAGAAGCGCCAGCCTCCAGTCGCAACTGGATCAGGGGGCGGCCTTCAAGGGCCAGGCTTAGCCAGCTGTCGCCCTCCTCGGCACTGAGGCTGCGCCAGCGTGGCGGTGACTGTGGCTGCGAAGCGGCATCGAACGGCATGGGCTTGACTCACGAAAATGAAGAAAAGGCACACTGCAGGGAACGTTGCCGCAACACGGGAATTTAGTCTTGCGGTGCAGCCACTTTGATTCGGTACGGCTCGAAGATGCGCGCCAGTTCACCGCTGTCACGCAGGCCACGCAGCAGG  >  NC_002947/4833048‑4833338

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: ‑