Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,869,858 T→G 6.8% E354A (GAG→GCG) ‡ 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,8580TG 6.8% 171.8 / 4.1 59E354A (GAG→GCG) ‡PP_3419sigma‑54 dependent transcriptional regulator/response regulator
Reads supporting (aligned to +/- strand):  ref base T (27/28);  new base G (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 = 3.03e-01

CAACGGCTCGTCAACCAAGGTCATGCACAGGTTCAACTGATGCGCTTGCACCACCTCGTTCGGCGCCAGCAACACCGTCTGCTCCAGCATGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTAGCAGGATGTCCTG‑CCCACGCGAATACAGGCGCGGCACCTTCAGGGCGATAATGCGAAGGCGGAAGAACAG  >  NC_002947/3869710‑3869997
                                                                                                                                                    |                                                                                                                                            
cAACGGCTCGTCAACCAAGGTCATGCACAGGTTCAACTGATGCGCTTGCACCACCTCGTTCGGCGCCAGCAACACCGTCTGCTCCAGCATGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTtcggcct                                                                                                                                              >  3:206576/1‑149 (MQ=255)
            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                                                                                               <  5:11485/143‑1 (MQ=255)
                                                     ccTCGTTCGGCGCCAGCAACACCGTCTGCTCCAGCATGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTAtggtgg                                                                                               >  6:11485/1‑143 (MQ=255)
                                                         gTTCGGCGCCAGCAACACCGTCTGCTCCAGCATGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAAt                                                                                     >  6:155547/1‑149 (MQ=255)
                                                              gcgcCAGCAACACCGTCTGCTCCAGCATGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCcgc                                                                                <  2:117/149‑1 (MQ=255)
                                                              gcgcCAGCAACACCGTCTGCTCCAGCATGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCcgc                                                                                >  2:231179/1‑149 (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)
                                                                                                                                                    |                                                                                                                                            
CAACGGCTCGTCAACCAAGGTCATGCACAGGTTCAACTGATGCGCTTGCACCACCTCGTTCGGCGCCAGCAACACCGTCTGCTCCAGCATGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTAGCAGGATGTCCTG‑CCCACGCGAATACAGGCGCGGCACCTTCAGGGCGATAATGCGAAGGCGGAAGAACAG  >  NC_002947/3869710‑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: ‑