Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,869,859 C→A 6.8% E354* (GAG→TAG) ‡ PP_3419 ← sigma‑54 dependent transcriptional regulator/response regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,869,8590CA 6.8% 180.4 / 3.7 59E354* (GAG→TAG) ‡PP_3419sigma‑54 dependent transcriptional regulator/response regulator
Reads supporting (aligned to +/- strand):  ref base C (27/28);  new base A (2/2);  total (29/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.90e-01

AACCAAGGTCATGCACAGGTTCAACTGATGCGCTTGCACCACCTCGTTCGGCGCCAGCAACACCGTCTGCTCCAGCATGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTAGCAGGATGTCCTG‑CCCACGCGAATACAGGCGCGGCACCTTCAGGGCGATAATGCGAAGGCGGAAGAACAG  >  NC_002947/3869722‑3869997
                                                                                                                                         |                                                                                                                                           
aaCCAAGGTCATGCACAGGTTCAACTGATGCGCTTGCACCACCTCGTTCGGCGCCAGCAACACCGTCTGCTCCAGCATGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAgc                                                                                                                                  <  4:206576/149‑1 (MQ=255)
  ccAAGGTCATGCACAGGTTCAACTGATGCGCTTGCACCACCTCGTTCGGCGCCAGCAACACCGTCTGCTCCAGCATGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAgcg                                                                                                                                 >  3:180811/1‑148 (MQ=255)
  ccAAGGTCATGCACAGGTTCAACTGATGCGCTTGCACCACCTCGTTCGGCGCCAGCAACACCGTCTGCTCCAGCATGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAgcg                                                                                                                                 >  3:1715/1‑148 (MQ=255)
   cAAGGTCATGCACAGGTTCAACTGATGCGCTTGCACCACCTCGTTCGGCGCCAGCAACACCGTCTGCTCCAGCATGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAgcg                                                                                                                                 >  2:292705/1‑147 (MQ=255)
              acaGGTTCAACTGATGCGCTTGCACCACCTCGTTCGGCGCCAGCAACACCGTCTGCTCCAGCATGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTg                                                                                                                    >  7:55919/1‑149 (MQ=255)
                  gTTCAACTGATGCGCTTGCACCACCTCGTTCGGCGCCAGCAACACCGTCTGCTCCAGCATGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCa                                                                                                                 >  7:217430/1‑148 (MQ=255)
                        cTGATGCGCTTGCACCACCTCGTTCGGCGCCAGCAACACCGTCTGCTCCAGCATGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCtt                                                                                                                     >  1:282561/1‑138 (MQ=255)
                        cTGATGCGCTTGCACCACCTCGTTCGGCGCCAGCAACACCGTCTGCTCCAGCATGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCtt                                                                                                                     <  2:282561/138‑1 (MQ=255)
                         tGATGCGCTTGCACCACCTCGTTCGGCGCCAGCAACACCGTCTGCTCCAGCATGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAg                                                                                                                                   <  7:135190/123‑1 (MQ=255)
                         tGATGCGCTTGCACCACCTCGTTCGGCGCCAGCAACACCGTCTGCTCCAGCATGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAg                                                                                                                                   >  8:135190/1‑123 (MQ=255)
                             gcgcTTGCACCACCTCGTTCGGCGCCAGCAACACCGTCTGCTCCAGCATGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTa                                                                                                     >  6:95336/1‑149 (MQ=255)
                               gcTTGCACCACCTCGTTCGGCGCCAGCAACACCGTCTGCTCCAGCATGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTAtg                                                                                                   >  6:310822/1‑149 (MQ=255)
                                 ttGCACCACCTCGTTCGGCGCCAGCAACACAGTCTGCTCCAGCATGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGCTTCGGCTTGCAATAGCGCCGGCTAtggt                                                                                                 >  5:227569/1‑149 (MQ=255)
                                   gCACCACCTCGTTCGGCGCCAGCAACACCGTCTGCTCCAGCATGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTAtggt                                                                                                 <  3:17142/147‑1 (MQ=255)
                                   gCACCACCTCGTTCGGCGCCAGCAACACCGTCTGCTCCAGCATGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTAtggt                                                                                                 <  5:310822/147‑1 (MQ=255)
                                    caccacCTCGTTCGGCGACAGCAACACCGTCTGCTCCAGCATGTTGCGCAATTGACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTAtggtgg                                                                                               <  8:55919/148‑1 (MQ=255)
                                        accTGGTTCGGCGCCAGCAACACCGTCTGCGCCTGCATGTTGCGCAATTCACGCTAATTGCCAGGCCAGCTATAACCGCGCAGCATGCTTTCGGACTCGGCCGAATAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAg                                                                                          <  4:190120/149‑1 (MQ=255)
                                         ccTCGTTCGGCGCCAGCAACACCGTCTGCTCCAGCATGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTAtggtgg                                                                                               >  6:11485/1‑143 (MQ=255)
                                         ccTCGTTCGGCGCCAGCAACACCGTCTGCTCCAGCATGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTAtggtgg                                                                                               <  5:11485/143‑1 (MQ=255)
                                             gTTCGGCGCCAGCAACACCGTCTGCTCCAGCATGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAAt                                                                                     >  6:155547/1‑149 (MQ=255)
                                                  gcgcCAGCAACACCGTCTGCTCCAGCATGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCcgc                                                                                >  2:231179/1‑149 (MQ=255)
                                                  gcgcCAGCAACACCGTCTGCTCCAGCATGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCcgc                                                                                <  2:117/149‑1 (MQ=255)
                                                      cAGCAACACCGTCTGCTCCAGCATGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCgcgc                                                                              >  2:290624/1‑147 (MQ=255)
                                                        gCAACACCGTCTGCTCCAGCATGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAgcgc                                                                                                                                <  5:290347/95‑1 (MQ=255)
                                                        gCAACACCGTCTGCTCCAGCATGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAgcgc                                                                                                                                >  6:290347/1‑95 (MQ=255)
                                                         cAACACCGTCTGCTCCAGCATGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTa                                                                         <  8:217430/149‑1 (MQ=255)
                                                            cacCGTCTTCTCCAGCATGTTGCGCAATTCACGCACATTGGCAGGCCAGCTATAACCGAGCAGCAGTCTTTCGGCCTCGGACGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTAGCa                                                                      <  6:227569/149‑1 (MQ=255)
                                                               cGTCTGCTCCAGCATGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTAGCAGGa                                                                   >  1:89460/1‑149 (MQ=255)
                                                                 tCTGCTCCAGCATGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTAGCAGGATg                                                                 >  3:205646/1‑149 (MQ=255)
                                                                     cTCCAGCATGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTAGCAGGATGt                                                                <  2:284632/146‑1 (MQ=255)
                                                                     cTCCAGCATGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTAGCAGGATGt                                                                >  1:284632/1‑146 (MQ=255)
                                                                        cAGCATGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTAGCAGGATGTCCTG‑cc                                                         >  1:227347/1‑149 (MQ=255)
                                                                         aGCATGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGACCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTAGCAGGATGTCCTG‑ccc                                                        <  8:59825/149‑1 (MQ=255)
                                                                           cATGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTAGCAGGATGTCCTG‑CCCAc                                                      >  7:69076/1‑149 (MQ=255)
                                                                             tGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTAGCAGGATGTCCTG‑CCCAcgc                                                    <  1:292705/149‑1 (MQ=255)
                                                                                tGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTAGCAGGATGTCCTG‑CCCACGCGaa                                                 <  1:231179/149‑1 (MQ=255)
                                                                                  cgcAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTAGCAGGATGTCCTG‑CCCACGCGAATa                                               <  8:69076/149‑1 (MQ=255)
                                                                                                 ttGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTAGCAGGATGTCCTG‑CCCACGCGAATACAGGCGCGGCACCtt                                <  1:290624/149‑1 (MQ=255)
                                                                                                     caggccagCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTAGCAGGATGTCCTG‑CCCACGCGAATACAGGCGCGGCACCTTCAgg                            <  3:262232/149‑1 (MQ=255)
                                                                                                     caggccagCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTAGCAGGATGTCCTG‑CCCACGCGAATACAGGCGCGGCACCTTCAgg                            <  2:227347/149‑1 (MQ=255)
                                                                                                     caggccagCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTAGCAGGATGTCCTG‑CCCACGCGAATACAGGCGCGGCACCTTCAgg                            <  3:230534/149‑1 (MQ=255)
                                                                                                         ccagCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTAGCAGGATGTCCTG‑CCCACGCGAATACAGGCGCGGCACCTTCAGGGCGa                        <  6:310786/149‑1 (MQ=255)
                                                                                                          cagCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTAGCAGGATGTCCTG‑CCCACGCGAATACAGGCGCGGCACCTTCAGGGCGa                        >  3:288818/1‑148 (MQ=255)
                                                                                                            gCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTAGCAGGATGTCCTG‑CCCACGCGAATACAGGCGCGGCACCt                                 <  2:199565/137‑1 (MQ=255)
                                                                                                            gCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTAGCAGGATGTCCTG‑CCCACGCGAATACAGGCGCGGCACCt                                 >  1:199565/1‑137 (MQ=255)
                                                                                                                taACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTAGCAGGATGTCCTG‑CCCACGCGAATACAGGCGCGGCACCTTCAGGGCGATAATGc                  <  5:155547/148‑1 (MQ=255)
                                                                                                                 aaCCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTAGCAGGATGTCCTG‑CCCACGCGAATACAGGCGCGGCACCTTCAGGGCGATAATGCg                 >  6:236015/1‑148 (MQ=255)
                                                                                                                   ccGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTAGCAGGATGTCCTG‑CCCACGCGAATACAGGCGCGGCACCTTCAGGGCGATAATGCGAAg              >  2:226370/1‑149 (MQ=255)
                                                                                                                         agcagGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTAGCAGGATGTCCTG‑CCCACGCGAATACAGGCGCGGCACCTTCAGGGCGATAATGCGAAGGCGgaa        <  4:180811/149‑1 (MQ=255)
                                                                                                                          gcagGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTAGCAGGATGTCCTGCCCCACGCGAATACAGGCGCGGCACCTTCAGGGCGATAATGCGAAGGCGgaa        <  1:226370/149‑1 (MQ=255)
                                                                                                                          gcagGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTAGCAGGATGTCCTG‑CCCACGCGAATACAGGCGCGGCACCTTCAGGGCGATAATGCGAAGGCgg          >  7:159557/1‑146 (MQ=255)
                                                                                                                             gcttttCGGCCGAGGCCG‑AAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTAGCAGGATGTCCTG‑CCCACGCGAATACAGGCGCGGCACCTTCAGGGCGATAATGc                  >  1:25918/4‑134 (MQ=255)
                                                                                                                             gcttttCGGCCGAGGCCG‑AAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTAGCAGGATGTCCTG‑CCCACGCGAATACAGGCGCGGCACCTTCAGGGCGATAATGc                  >  1:30676/4‑134 (MQ=255)
                                                                                                                             gcttttCGGCCGAGGCCG‑AAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTAGCAGGATGTCCTG‑CCCACGCGAATACAGGCGCGGCACCTTCAGGGCGATAATGc                  <  2:30676/131‑1 (MQ=255)
                                                                                                                             gcttttCGGCCGAGGCCG‑AAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTAGCAGGATGTCCTG‑CCCACGCGAATACAGGCGCGGCACCTTCAGGGCGATAATGc                  <  2:25918/131‑1 (MQ=255)
                                                                                                                              gCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTAGCAGGATGTCCTG‑CCCACGCGAATACAGGCGCGGCACCTTCAGGGCGATAATGCGAAGGCGgaagaa     >  8:50224/1‑147 (MQ=255)
                                                                                                                               cTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTAGCAGGATGTCCTG‑CCCACGCGAATACAGGCGCGGCACCTTCAGGGCGATAATGCGAAGGCGGAAGAACAg  <  2:89460/149‑1 (MQ=255)
                                                                                                                                tttCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGaaa                                                                                      <  2:276030/65‑1 (MQ=255)
                                                                                                                                tttCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGaaa                                                                                      >  1:276030/1‑65 (MQ=255)
                                                                                                                                       ctcggccGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTAGCAGGATGTCCTG‑CCCACGCGAATACAGGCGCGGCACCTTCAg                             <  3:81937/114‑1 (MQ=255)
                                                                                                                                       ctcggccGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTAGCAGGATGTCCTG‑CCCACGCGAATACAGGCGCGGCACCTTCAg                             >  4:81937/1‑114 (MQ=255)
                                                                                                                                         |                                                                                                                                           
AACCAAGGTCATGCACAGGTTCAACTGATGCGCTTGCACCACCTCGTTCGGCGCCAGCAACACCGTCTGCTCCAGCATGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTAGCAGGATGTCCTG‑CCCACGCGAATACAGGCGCGGCACCTTCAGGGCGATAATGCGAAGGCGGAAGAACAG  >  NC_002947/3869722‑3869997

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