Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,956,707 G→C 7.8% R37R (CGG→CGC)  cyaA → adenylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,956,7070GC 7.8% 183.4 / 4.6 61R37R (CGG→CGC) cyaAadenylate cyclase
Reads supporting (aligned to +/- strand):  ref base G (28/28);  new base C (3/2);  total (31/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 = 9.91e-01

CTCGCCACAGCTTCAATGGCCTGCCCGAGCGCGCAATGAACCACCCTCACGAAATCCGCCCTGACCTGGACGAAGGCATCGACCGCAAGGTCCTGGCGACATTGCGTGCGCGCTTCCTGCATCTGAACCAGGGCCGACTGCAGCGGGCTAAGGAAGGCCTGTCAACGCGCCAGCAACAGGTACTGTCATTGCTGCCGCTGCTGTTTCACGTGAACCACCCGCTATTGCCGGGTTACGTCTCCGGCAGCACCCCCGCGGGCGTATCGGCTTACGAGCCAAGCGCCGAGCT  >  NC_002947/5956562‑5956850
                                                                                                                                                 |                                                                                                                                               
cTCGCCACAGCTTCAATGGCCTGCCCGAGCGCGCAATGAACCACCCTCACGAAATCCGCCCTGACCTGGACGAAGGCATCGACCGCAAGGTCCTGGCGACATTGCGTGCGCGCTTCCTGCATCTGAACCAGGGCCGACTGCAGCGGGc                                                                                                                                               >  8:214481/1‑148 (MQ=255)
   gCCACAGCTTCAATGGCCTGCCCGAGCGCGCAATGAACCACCCTCACGAAATCCGCCCTGACCTGGACGAAGGCATCGACCGCAAGGTCCTGGCGACATTGCGTGCGCGCTTCCTGCATCTGAACCAGGGCCGACTGCAGCGGGCTaag                                                                                                                                           <  2:20164/149‑1 (MQ=255)
      acaGATTCAATGGCCTGCCCGAGCGCGCAATGAACCACCCTCACGAAATCCGCCCTGACCTGGACGGAGGCATCGACCGCAAGGTCCTGGCGACATTGCGTGCGCGCTTCCTGCATCTGAACCAGGGCCGACTGCAGCGGGCTaaggaa                                                                                                                                        <  5:222312/149‑1 (MQ=255)
         gCTTCAATGGCCTGCCCGAGCGCGCAATGAACCACCCTCACGAAATCCGCCCTGACCTGGACGAAGGCATCGACCGCAAGGTCCTGGCGACATTGCGTGCGCGCTTCCTGCATCTGAACCAGGGCCGACTGCAGCGGGCTAAGGAAGGc                                                                                                                                     <  6:206368/149‑1 (MQ=255)
          cTTCAATGGCCTGCCCGAGCGCGCAATGAACCACCCTCACGAAATCCGCCCTGACCTGGACGAAGGCATCGACCGCAAGGTCCTGGCGACATTGCGTGCGCGCTTCCTGCATCTGAACCAGGGCCGACTGCAGCGGGCTAAGGAAGGcc                                                                                                                                    >  3:226508/1‑149 (MQ=255)
             cAATGGCCTGCCCGAGCGCGCAATGAACCACCCTCACGAAATCCGCCCTGACCTGGACGAAGGCATCGACCGCAAGGTCCTGGCGACATTGCGTGCGCGCTTCCTGCATCTGAACCAGGGCCGACTGCAGCGGGCTAAGGAAGGCCTGt                                                                                                                                 >  3:245151/1‑149 (MQ=255)
              aaTGGCCTGCCCGAGCGCGCAATGAACCACCCTCACGAAATCCGCCCTGACCTGGACGAAGGCATCGACCGCAAGGTCCTGGCGACATTGCGTGCGCGCTTCCTGCATCTGAACCAGGGCCGACTGCAGCGGGCTAAGGAAGGCCTGTc                                                                                                                                <  6:35475/149‑1 (MQ=255)
               aTGGCCTGCCCGAGCGCGCAATGAACCACCCTCACGAAATCCGCCCTGACCTGGACGAAGGCATCGACCGCAAGGTCCTGGCGACATTGCGTGCGCGCTTCCTGCATCTGAACCAGGGCCGACTGCAGCGCGCTAAGGAAGGCCTGTc                                                                                                                                >  4:81945/1‑148 (MQ=255)
               aTGGCCTGCCCGAGCGCGCAATGAACCACCCTCACGAAATCCGCCCTGACCTGGAAGAAGGCATCGACCGCAAGGTCCTGGCGACATTGCGTGCGCGCTTCCTGCATCTGAACCAGGGCCGACTGCAGCGCGCTAAGGAAGGCCTGTc                                                                                                                                >  4:77440/1‑148 (MQ=255)
                      gCCCGAGCGCGCAATGAACCACCCTCACGAAATCCGCCCTGACCTGGACGAAGGCATCGACCGCAAGGTCCTGGCGACATTGCGTGCGCGCTTCCTGCATCTGAACCAGGGCCGACTGCAGCGGGCTAAGGAAGGc                                                                                                                                     >  2:149763/1‑136 (MQ=255)
                      gCCCGAGCGCGCAATGAACCACCCTCACGAAATCCGCCCTGACCTGGACGAAGGCATCGACCGCAAGGTCCTGGCGACATTGCGTGCGCGCTTCCTGCATCTGAACCAGGGCCGACTGCAGCGGGCTAAGGAAGGc                                                                                                                                     <  1:149763/136‑1 (MQ=255)
                       cccGAGCGCGCAATGAACCACCCTCACGAAATCCGCCCTGACCTGGACGAAGGCATCGACCGCAAGGTCCTGGCGACATTGCGTGCGCGCTTCCTGCATCTGAACCAGGGCCGACTGCAGCGGGCTAAGGAAGGCCTGTCAACgcgc                                                                                                                         <  7:214481/147‑1 (MQ=255)
                       cccGAGCGCGCAATGAACCACCCTCACGAAATCCGCCCTGACCTGGACGAAGGCATCGACCGCAAGGTCCTGGCGACATTGCGTGCGCGCTTCCTGCATCTGAACCAGGGCCGACTGCAGCGGGCTAAGGAAGGCCTGTCAACgcgc                                                                                                                         >  6:155185/1‑147 (MQ=255)
                       cccGAGCGCGCAATGAACCACCCTCACGAAATCCGCCCTGACCTGGACGAAGGCATCGACCGCAAGGTCCTGGCGACATTGCGTGCGCGCTTCCTGCATCTGAACCAGGGCCGACTGCAGCGGGCTAAGGAAGGCCTGTCAACgcgc                                                                                                                         >  3:38257/1‑147 (MQ=255)
                         cGAGCGCGCAATGAACCACCCTCACGAAATCCGCCCTGACCTGGACGAAGGCATCGACCGCAAGGTCCTGGCGACATTGCGTGCGCGCTTCCTGCATCTGAACCAGGGCCGACTGCAGCGGGCTAAGGAAGGCCTGTCAACGCGCCAGc                                                                                                                     <  2:315968/149‑1 (MQ=255)
                            gcgcgcAATGAACCACCCTCACGAAATCCGCCCTGACCTGGACGAAGGCATCGACCGCAAGGTCCTGGCGACATTGCGTGCGCGCTTCCTGCATCTGAACCAGGGCCGACTGCAGCGGGCTAAGGAAGGCCTGTCAACGCGCCAGCAAc                                                                                                                  <  8:144353/149‑1 (MQ=255)
                                gcAATGAACCACCCTCACGAAATCCGCCCTGACCTGGACGAAGGCATCGACCGCAAGGTCCTGGCGACATTGCGTGCGCGCTTCCTGCATCTGAACCAGGGCCGACTGCAGCGGGCTAAGGAAGGCCTGTCAACGCGCCAGCAACAgg                                                                                                               >  6:59442/1‑148 (MQ=255)
                                gcAATGAACCACCCTCACGAAATCCGCCCTGACCTGGACGAAGGCATCGACCGCAAGGTCCTGGCGACATTGCGTGCGCGCTTCCTGCATCTGAACCAGGGCCGACTGCAGCGGGCTAAGGAAGGCCTGTCAACGCGCCAGCAACAGGt                                                                                                              <  3:173966/149‑1 (MQ=255)
                                    tGAACCACCCTCACGAAATCCGCCCTGACCTGGACGAAGGCATCGACCGCAAGGTCCTGGCGACATTGCGTGCGCGCTTCCTGCATCTGAACCAGGGCCGACTGCAGCGGGCTAAGGAAGGCCTGTCAACGCGCCAGCAACAGGTACTg                                                                                                          <  4:226508/149‑1 (MQ=255)
                                     gccccaccCTCCCGCCATCCGCGCTGCCCCGGACGCAGGCATCGACCGCAAGGTCCTGGCGACATTGCGTGCGCGCTTCCTGCATCTGACCCAGCGCCGACTGCAGCGGGCTAAGGAAGGCCTGTCAACGCGCCAGCAACAGGTACTGt                                                                                                         <  4:245151/146‑1 (MQ=255)
                                      aaCCACCCTCACGAAATCCGCCCTGACCTGGACGAAGGCATCGACCGCAAGGTCCTGGCGACATTGCGTGCGCGCTTCCTGCATCTGAACCAGGGCCGACTGCAGCGGGCTAAGGAAGGCCTGTCAACGCGCCAGCAACAGGTACTGt                                                                                                         >  2:161038/1‑148 (MQ=255)
                                       accaccCTCACGAAATCCGCCCTGACCTGGACGAAGGCATCGACCGCAAGGTCCTGGCGACATTGCGTGCGCGCTTCCTGCATCTGAACCAGGGCCGACTGCAGCGGGCTAAGGAAGGCCTGTCAACGCGCCAGCAACAGGTACt                                                                                                           <  1:36271/145‑1 (MQ=255)
                                       accaccCTCACGAAATCCGCCCTGACCTGGACGAAGGCATCGACCGCAAGGTCCTGGCGACATTGCGTGCGCGCTTCCTGCATCTGAACCAGGGCCGACTGCAGCGGGCTAAGGAAGGCCTGTCAACGCGCCAGCAACAGGTACt                                                                                                           >  2:36271/1‑145 (MQ=255)
                                        ccaccCTCACGAACTCCGCCCTGACCTGGACGAAGGCATCGACCGCAAGGTCCTGGCGACATTGCGTGCGCGCTTCCTGCATCTGAACCAGGGCCGACTGCAGCGGGCTAAGGAAGGCCTGTCAACGCGCCAGCAACAGGTACTGTCAt                                                                                                      <  8:139257/149‑1 (MQ=255)
                                             cTCACGAAATCCGCCCTGACCTGGACGAAGGCAACGACCGCAAGGTCCTGGCGACATTGCGTGCGCGCTTCCTGCATCTGAACCAGGGCCGACTGCAGCGGGCTAAGGAAGGCCTGTCAACGCGCCAGCAACAGGTACTGTCATtgctg                                                                                                 >  1:208795/1‑149 (MQ=255)
                                              tCACGAAATCCGCCCTGACCTGGACGAAGGCATCGACCGCAAGGTCCTGGCGACATTGCGTGCGCGCTTCCTGCATCTGAACCAGGGCCGACTGCAGCGCGCTAAGGAAGGCCTGTCAACGCGCCAGCAACAGGTACTGTCATtgctgc                                                                                                <  3:77440/149‑1 (MQ=255)
                                                aCGAAATCCGCCCTGACCTGGACGAAGGCATCGACCGCAAGGTCCTGGCGACATTGCGTGCGCGCTTCCTGCATCTGAACCAGGGCCGACTGCAGCGGGGTAAGGAGGGGCTGTCAACGCGGCAGCAAGAGGTACGGTCATTGCTGCCg                                                                                              >  6:186752/1‑149 (MQ=255)
                                                aCGAAATCCGCCCTGACCTGGACGAAGGCATCGACCGCAAGGTCCTGGCGACATTGCGTGCGCGCTTCCTGCATCTGAACCAGGGCCGACTGCAGCGGGCTAAGGAAGGCCTGTCAACGCGCCAGCAACAGGTACTGTCATtgctgc                                                                                                <  2:6663/147‑1 (MQ=255)
                                                aCGAAATCCGCCCTGACCTGGACGAAGGCATCGACCGCAAGGTCCTGGCGACATTGCGTGCGCGCTTCCTGCATCTGAACCAGGGCCGACTGCAGCGCGCTAAGGAAGGCCTGTCAACGCGCCAGCAACAGGTACTGTCATtgctgc                                                                                                <  3:81945/147‑1 (MQ=255)
                                                  gAAATCCGCCCTGACCTGGACGAAGGCATCGACCGCAAGGTCCTGGCGACATTGCGTGCGCGCTTCCTGCATCTGAACCAGGGCCGACTGCAGCGGGCTAAGGAAGGGCTGTCAACGCGCCAGCAACAGGTACTGTCATTGCTGCCgct                                                                                            >  5:276820/1‑149 (MQ=255)
                                                          cccTGACCTGGACGAAGGCATCGACCGCAAGGTCCTGGCGACATTGCGTGCGCGCTTCCTGCATCTGAACCAGGGCCGACTGCAGCGGGCTAAGGAAGGCCTGTCAACGCGCCAGCAACAGGTACTGTCATTGCTGCCGCTGCTGTTTc                                                                                    <  6:302268/149‑1 (MQ=255)
                                                            cTGACCTGGACGAAGGCATCGACCGCAAGGTCCTGGCGACATTGCGTGCGCGCTTCCTGCATCTGAACCAGGGCCGACTGCAGCGGGCTAAGGAAGGCCTGTCAACGCGCCAGCAACAGGTACTGTCATTGCTGCCGCTGCTGTTTCAc                                                                                  <  5:155185/149‑1 (MQ=255)
                                                               acatGGACGAAGGCATCGACCGCAAGGTCCTGGCGACATTGCGTGCGCACTTCCTGCATCTGAACCAGGGCCGACTGCAGCGGGCTAAGGAAGGCCTGTCAACGAGCCAGCAACAGGTACTGTCATTGCTGCCGCTGCTGTTTCACGt                                                                                >  6:29705/4‑148 (MQ=255)
                                                                ccTGGACGAAGGCATCGACCGCAAGGTCCTGGCGACATTGCGTGCGCGCTTCCTGCATCTGAACCAGGGCCGACTGCAGCGGGCTAAGGAAGGCCTGTCAACGCGCCAGCAACAGGTACTGTCATTGCTGCCGCTGCTGTTTCACGTg                                                                               <  4:38257/148‑1 (MQ=255)
                                                                    gACGAAGGCATCGACCGCAAGGTCCTGGCGACATTGCGTGCGCGCTTCCTGCATCTGAACCAGGGCCGACTGCAGCGGGCTAAGGAAGGCCTGTCAACGCGCCAGCAACAGGTACTGTCATTGCTGCCGCTGCTGTTTCACGTGAacca                                                                          >  5:291322/1‑149 (MQ=255)
                                                                    gACGAAGGCATCGACCGCAAGGTCCTGGCGACATTGCGTGCGAGCTTCATGCATCTGAACCAGCGCCGACTGCAGCGGGCTAAGGAAGCCCTGTCAACTCGCCAGCAACAGGTAGAGTCAATGCTGACGCTGCTGATTCACGTGAacca                                                                          >  2:229348/1‑149 (MQ=255)
                                                                                       aGGTCCTGGCGACATTGCGTGCGCGCTTCCTGCATCTGAACCAGGGCCGACTGCAGCGGGCTAAGGAAGGCCTGTCAACGCGCCAGCAACAGGTACTGTCATTGCTGCCGCTGCTGTTTCACGTGAACCACCCGCTATTGCCGGGTTAc                                                       <  6:291322/149‑1 (MQ=255)
                                                                                             ccaggACCTTGCGTGCGCGCTTCCTGCATCTGAACCAGGGCCGACTGCAGCGGGCTAAGGAAGGCCTGTCAACGCGCCAGCAACAGGTACTGTCATTGCTGCCGCTGCTGTTTCACGTGAACCACCCGCTATTGCCGGGTTACGTCTcc                                                 <  1:145843/145‑1 (MQ=255)
                                                                                               gCGACATTGCGTGCGCGCTTCCTGCATCTGAACCAGGGCCGACTGCAGCGGGCTAAGGAAGGCCTGTCAACGCGCCAGCAACAGGTACTGTCATTGCTGCCGCTGCTGTTTCACGTGAACCACCCGCTATTGCCGGGTTACGTCTCCgg                                               <  2:208795/149‑1 (MQ=255)
                                                                                                 gACATTGCGTGCGCGCTTCCTGCATCTGAACCAGGGCCGACTGCAGCGGGCTAAGGAAGGCCTGTCAACGCGCCAGCAACAGGTACTGTCATTGCTGCCGCTGCTGTTTCACGTGAACCACCCGCTATTGCCGGGTTACGTCTCCGgca                                             <  1:161038/149‑1 (MQ=255)
                                                                                                            cgcgcTTCCTGCATCTGAACCAGGGCCGACTGCAGCGGGCTAAGGAAGGCCTGTCAACGCGCCAGCAACAGGTACTGTCATTGCTGCCGCTGCTGTTTCACGTGAACCACCCGCTATTGCCGGGTTACGTCTCCGGCAGCACCCCcgcg                                  >  4:81511/1‑149 (MQ=255)
                                                                                                               gcTTCCTGCAGCTGCACCCGGGCCGACTGCAGCGGGCGAAGGAAGGCCTGTCAACGCGCCAGCAACAGGTACTGTCATTGCTGCCGCTGCTGTTTCACGTGAACCACCCGCTATTGCCGGGTTACGTCTCCGGCAGCACCCCCGCGGGc                               <  7:303390/149‑1 (MQ=255)
                                                                                                                 ttCCTGCATCTGAACCAGGGCCGACTGCAGCGGGCTAAGGAAGGCCTGTCAACGCGCCAGCAACAGGTACTGTCATTGCTGCCGCTGCTGTTTCACGTGAACAACCCGCTATTGCCGGGTTACGTCTCCGGCAGCACCCCCCCGGGCg                              >  1:123062/1‑148 (MQ=255)
                                                                                                                    cTGCATCTGAACCAGGGCCGACTGCAGCGGGCTAAGGAAGGCCTGTCAACGCGCCAGCAACAGGTACTGTCATTGCTGCCGCTGCTGTTTCACGTGAACCACCCGCTATTGCCGGGTTACGTCTCCGGCAGCACCCCCGCGGGCGTATc                          >  3:16713/1‑149 (MQ=255)
                                                                                                                     tGCATCTGAACCAGGGCCGACTGCAGCGGGCTAAGGAAGGCCTGTCAACGCGCCAGCAACAGGTACTGTCATTGCTGCCGCTGCTGTTTCACGTGACCCACCCGCTATTGCCGGGTTACGTCTCCGGCAGCACCCACGCGGGCGTATCg                         >  8:112481/1‑149 (MQ=255)
                                                                                                                     tGCATCTGAACCAGGGCCGACTGCAGCGGGCTAAGGAAGGCCTGTCAACGCGCCAGCAACAGGTACTGTCATTGCTGCCGCTGCTGTTTCACGTGAACCACCCGCTATTGCCGGGTTACGTCTCCGGCAGCACCCCCGCGGGCGTATCg                         >  5:76685/1‑149 (MQ=255)
                                                                                                                      gCATCTGAACCAGGGCCGACTGCAGCGGGCTAAGGAAGGCCTGTCAATGCGCCAGCAACAGGTACTGTCATTGCTGCCGCTGCTGTTTCACGTGAACCACCCGCTATTGCCGGGTTACGTCTCCGGCAGCACCCCCGCGGGCGTATCgg                        >  6:167488/1‑149 (MQ=255)
                                                                                                                       cATCTGAACCAGGGCCGACTGCAGCGGGCTAAGGAAGGCCTGTCAACGCGCCAGCAACAGGTACTGTCATTGCTGCCGCTGCTGTTTCACGTGAACCACCCGCTATTGCCGGGTTACGTCTCCGGCAGCACCCCCGCGGGCGTATCgg                        <  3:81511/148‑1 (MQ=255)
                                                                                                                        aTCTGAACCAGGGCCGACTGCAGCGGGCTAAGGAAGGCCTGTCAACGCGCCAGCAACAGGTACTGTCATTGCTGCCGCTGCTGTTTCACGTGAACCACCCGCTATTGCCGGGTTACGTCTcc                                                 <  1:291198/122‑1 (MQ=255)
                                                                                                                        aTCTGAACCAGGGCCGACTGCAGCGGGCTAAGGAAGGCCTGTCAACGCGCCAGCAACAGGTACTGTCATTGCTGCCGCTGCTGTTTCACGTGAACCACCCGCTATTGCCGGGTTACGTCTcc                                                 >  2:291198/1‑122 (MQ=255)
                                                                                                                         tCTGAACCAGGGCCGACTGCAGCGGGCTAAGGAAGGCCTGTCAACGCGCCAGCAACAGGTACTGTCATTGCTGCCGCTGCTGTTTCACGTGAACCACCCGCTATTGCCGGGTTACGTCTCCGGCAGCACCCCCGCgg                                 <  1:310548/137‑1 (MQ=255)
                                                                                                                         tCTGAACCAGGGCCGACTGCAGCGGGCTAAGGAAGGCCTGTCAACGCGCCAGCAACAGGTACTGTCATTGCTGCCGCTGCTGTTTCACGTGAACCACCCGCTATTGCCGGGTTACGTCTCCGGCAGCACCCCCGCgg                                 >  2:310548/1‑137 (MQ=255)
                                                                                                                             aaacAGGGCCGACTGCAGCGCGCTAAGGAAGGCCTCTCAACGCGCCAGCAACAGGTACTGTCATTGCTGCCGCTGCTGTTTCACGTGCACCACCCGCTATTGCCGGGTTACGTCTCCGGCAGCACCCCCGCGGGCGTATCGGCTTAc                   >  2:133550/4‑147 (MQ=255)
                                                                                                                              aCCAGGGCCGACTGCAGCGGGCTAAGGAAGGCCTGTCAACGCGCCAGCAACAGGTACTGTCATTGCTGCCGCTGCTGTTTCACGTGAACCACCCGCTATTGCCGGGTTACGTCTCCGGCAGCACCCCCGCGGGCGTATCGGCTTACGa                 <  5:59442/148‑1 (MQ=255)
                                                                                                                               ccAGGGCCGACTGCAGCGGGCTAAGGAAGGCCTGTCAACGCGCCAGCAACAGGTACTGTCATTGCTGCCGCTGCTGTTTCACGTGAACCACCCGCTATTGCCGGGTTACGTCTCCGGCAGCACCCCCGCGGGCGTATCGGCTTACGAGc               >  3:79476/1‑149 (MQ=255)
                                                                                                                                cAGGGCCGACTGCAGCGGGCTAAGGAAGGCCTGTCAACGCGCCAGCAACAGGTACTGTCATTGCTGCCGCTGCTGTTTCACGTGAACCACCCGCTATTGCCGGGTTACGTCTCCGGCAGCACCCCCGCGGGCGTATCGGCTTACGAGc               >  1:179365/1‑148 (MQ=255)
                                                                                                                                    gCCGAGTGCAGCGGGCTAAGGAAGGCCTGTCAACGCGCCTGCAAGAGGTAGTGTCATTGCTGCCGCTGCTGTTTCTCGTGATCCTCGCTCTATTGCCGGGTTACGTCTCCGGCAGCACCCGGTCTGGCGTATCGGGTTACGTGCCTAgc          <  1:229348/149‑1 (MQ=255)
                                                                                                                                         cTGCAGCGGGCTAAGGAAGGCCTGTCAACGCGCCAGCAACAGGTACTGTCATTGCTGCCGCTGCTGTTTCACGTGAACCACCCGCTATTGCCGGGTTACGTCTCCGGCAGCACCCCCGCGGGCGTATCGGCTTACGAGCCAAGCGCCGa     >  8:200293/1‑149 (MQ=255)
                                                                                                                                          tGCAGCGGGCTAAGGAAGGCCTGTCAACGCGCCAGCAACAGGTACTGTCATTGCTGCCGCTGCTGTTTCACGTGAACCACCCGCTATTGCCGGGTTACGTCTCCGGCAGCACCCCCGCGGGCGTATCGGCTTACGAGCCAAGCGCCGAg    >  8:115229/1‑149 (MQ=255)
                                                                                                                                            cAGCGGGCTAAGGAAGGCCTGTCAACGCGCCAGCAACAGGTACTGTCATTGCTGCCGCTGCTGTTTCACGTGAACCACCCGCTATTGCCGGGTTACGTCTCCGGCAGCACCCCCGCGGGCGTATCGGCTTACGAGCCAAGCGCCGAGc   <  6:276820/148‑1 (MQ=255)
                                                                                                                                            cAGCGGGCTAAGGAAGGCCTGTCAACGCGCCAGCAACAGGTACTGTCATTGCTGCCGCTGCTGTTTCACGTGAACCACCCGCTATTGCCGGGTTACGTCTCCGGCAGCACCCCCGCGGGCGTATCGGCTTACGAGCCAAGCGCCGAGCt  >  5:62688/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
CTCGCCACAGCTTCAATGGCCTGCCCGAGCGCGCAATGAACCACCCTCACGAAATCCGCCCTGACCTGGACGAAGGCATCGACCGCAAGGTCCTGGCGACATTGCGTGCGCGCTTCCTGCATCTGAACCAGGGCCGACTGCAGCGGGCTAAGGAAGGCCTGTCAACGCGCCAGCAACAGGTACTGTCATTGCTGCCGCTGCTGTTTCACGTGAACCACCCGCTATTGCCGGGTTACGTCTCCGGCAGCACCCCCGCGGGCGTATCGGCTTACGAGCCAAGCGCCGAGCT  >  NC_002947/5956562‑5956850

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