Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,877,474 C→G 7.9% A15A (GCC→GCG PP_1683 → MarR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,877,4740CG 7.9% 198.0 / 5.3 63A15A (GCC→GCGPP_1683MarR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (27/31);  new base G (3/2);  total (30/33)
Fisher's exact test for biased strand distribution p-value = 6.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.01e-01

AGCCGATCCGCAGGCGCTTGAACCCGGACTGACGCTTCTCAATCTTTCGTGAAACTTCTTGCAACACTCCGTACGCGGGTATATACACGTATCCATGCTGACCAGTGAATGCATCTGTACCCACCTGCGTCGCGCCGCCCGCGGCGTGAGCCGGCATTACGACGAAGCCCTTGCCGGCTTCGGGGTCAATGTCGCCCAGTTTTCCCTGCTGCGGCACCTGCAGCGGCTCGACCGACCCAGCATCACAACCCTGGCCGAAGCCATGGGGCTCGAGCGCAGTACCCTG  >  NC_002947/1877336‑1877621
                                                                                                                                          |                                                                                                                                                     
aGCCGATCCGCAGGCGCTTGAACCCGGACTGACGCTTCTCAATCTTTCGTGAAACTTCTTGCAACACTCCGTACGCGGGTATATACACGTATCCATGCTGACCAGTGAATGCATCTGTACCCACCTGCGTCGCCCCCCCGgccgcgcgc                                                                                                                                             >  4:106873/1‑146 (MQ=255)
  ccGATCCGCAGGCGCTTGAACCCGGACTGACGCTTCTCAATCTTTCGTGAAACTTCTTGCAACACTCCGTACGCGGGTATATACACGTATCCATGCTGACCAGTGAATGCATCTGTACCCACCTGCGTCGCGCCGCCCGCGGCGTGAGc                                                                                                                                           <  3:117792/149‑1 (MQ=255)
     aTCCGCAGGCGCTTGAACCCGGACTGACGCTTCTCAATCTTTCGTGAAACTTCTTGCAACACTCCGTACGCGGGTATATACACGTATCCATGCTGACCAGTGAATGCATCTGTACCCACCTGCGTCGCGCCGCGGGCGGCGCGacgcag                                                                                                                                        <  3:106873/149‑6 (MQ=255)
     aTCCGCAGGCGCTTGAACCCGGACTGACGCTTCTCAATCTTTCGTGAAACTTCTTGCAACACTCCGTACGCGGGTATATACACGTATCCATGCTGACCAGTGAATGCATCTGTACCCACCTGCGTCGCGCCGCCCGCGGCGTGAGCCgg                                                                                                                                        <  1:187541/149‑1 (MQ=255)
       ccGCAGGCGCTTGAACCCGGACTGACGCTTCTCAATCTTTCGTGAAACTTCTTGCAACACTCCGTACGCGGGTATATACACGTATCCATGCTGACCAGTGAATGCATTTGTACCCACCTGCGTCGCGCCGCCCGCGGCGTGAGCCGGCa                                                                                                                                      >  1:274838/1‑149 (MQ=255)
         gCAGGCGCTTGAACCCGGACTGACGCTTCTCAATCTTTCGTGAAACTTCTTGCAACACTCCGTACGCGGGTATATACACGTATCCATGCTGACCAGTGAATGCATCTGTACCCACCTGCGTCGCGCCGCCCGCGGCGTGAGCCGGCAtt                                                                                                                                    <  8:7363/149‑1 (MQ=255)
         gCAGGCGCTTGAACCCGGACTGACGCTTCTCAATCTTTCGTGAAACTTCTTGCAACACTCCGTACGCGGGTATATACACGTATCCATGCTGACCAGTGAATGCATCTGTACCCACCTGCGTCGCGCCGCCCGCGGCGTGAGCCGGCAt                                                                                                                                     >  4:69382/1‑148 (MQ=255)
          cAGGCGCTTGAACCCGGACTGACGCTTCTCAATCTTTCGTGAAACTTCTTGCAACACTCCGTACGCGGGTATATACACGTATCCATGCTGACCAGTGAATGCATCTGTACCCACCTGCGTCGCGCCGCCCGCGGCGTGAGCCGGCATTa                                                                                                                                   <  7:188506/149‑1 (MQ=255)
               gcTTGAACCCGGACTGACGCTTCTCAATCTTTCGTGAAACTTCTTGCAACACTCCGTACGCGGGTATATACACGTATCCATGCTGACCAGTGAATGCATCTGTACCCACCTGCGTCGCGCCGCCCGCGGCGTGAGCCGGCATTacgac                                                                                                                               >  5:60122/1‑148 (MQ=255)
                         ggACTGACGCTTCTCAATCTTTCGTGAAACTTCTTGCAACACTCCGTACGCGGGTATATACACGTATCCATGCTGACCAGTGAATGCATCTGTACCCACCTGCGTCGCGCCGCCCGCGGCGTGAGCCGGCATTACGACGAAGCCCTTGc                                                                                                                    >  4:49694/1‑149 (MQ=255)
                               aCGCTTCTCAATCTTTCGTGAAACTTCTTGCAACACTCCGTACGCGGGTATATACACGTATCCATGCTGACCAGTGAATGCATCTGTACCCACCTGCGTCGCGCCGCCCGCGGCGTGAGCCGGCATTACGACGAAGCCCTTGCCGGCtt                                                                                                            >  2:210125/1‑149 (MQ=255)
                                 gCTTCTCAATCTTTCGTGAAACTTCTTGCAACACTCCGTACGCGGGTATATACACGTATCCATGCTGACCAGTGAATGCATCTGTACCCACCTGCGTCGCGCCGCCCGCGGCGTGAGCCGGCATTACGACGAAGCCCTTGCCGGCTTCg                                                                                                          <  6:295946/149‑1 (MQ=255)
                                   ttCTCAATCTTTCGTGAAACTTCTTGCAACACTCCGTACGCGGGTATATACACGTATCCATGCTGACCAGTGAATGCATCTGTACCCACCTGCGTCGCGCCGCCCGCGGCGTGAGCCGGCATTACGACGAAGCCCTTGCCGGCTTCggg                                                                                                        >  6:75779/1‑149 (MQ=255)
                                     ctcAATCTTTCGTGAAACTTCTTGCAACACTCCGTACGCGGGTATATACACGTATCCATGCTGACCAGTGAATGCATCTGTACCCACCTGCGTCGCGCCGCCCGCGGCGTGAGCCGGCATTACGACGAAGCCCTTGCCGGCTTCggg                                                                                                        <  5:75779/147‑1 (MQ=255)
                                     ctcAATCTTTCGTGAAACTTCTTGCAACACTCCGTACGCGGGTATATACACGTATCCATGCTGACCAGTGAATGCATCTGTACCCACCTGCGTCGCGCCGCCCGCGGCGTGAGCCGGCATTACGACGAAGCCCTTGCCGGCTTCggg                                                                                                        <  6:218534/147‑1 (MQ=255)
                                     ctcAATCTTTCGTGAAACTTCTTGCAACACTCCGTACGCGGGTATATACACGTATCCATGCTGACCAGTGAATGCATCTGTACCCACCTGCGTCGCGCCGCCCGCGGCGTGAGCCGGCATTACGACGAAGCCCTTGCCGGCTTCggg                                                                                                        <  6:60122/147‑1 (MQ=255)
                                         aTCTTTCGTGAAACTTCTTGCAACACTCCGTACGCGGGTATATACACGTATCCATGCTGACCAGTGAATGCATCTGTACCCACCTGCGTCGCGCCGCCCGCGGCGTGAGCCGGCATTACGACGAAGCCCTTGCCGGCTTCGGGGTCAAt                                                                                                  <  8:271435/149‑1 (MQ=255)
                                          tCTTTCGTGAAACTTCTTGCAACACTCCGTACGCGGGTATATACACGTATCCATGCTGACCAGTGAATGCATCTGTACCCACCTGCGTCGCGCCGCCCGCGGCGTGAGCCGGCATTACGACGAAGCCCTTGCCGGCTTCGGGGTCAAt                                                                                                  >  3:134148/1‑148 (MQ=255)
                                          tCTTTCGTGAAACTTCTTGCAACACTCCGTACGCGGGTATATACACGTATCCATGCTGACCAGTGAATGCATCTGTACCCACCTGCGTCGCGCCGCCCGCGGCGTGAGCCGGCATTACGACGAAGCCCTTGCCGGCTTCGGGGTCAAt                                                                                                  <  4:134148/148‑1 (MQ=255)
                                            tttCGTGAAACTTCTTGCAACACTCCGTACGCGGGTATATACACGTATCCATGCTGACCAGTGAATGCATCTGTACCCACCTGCGTCGCGCCGCCCGCGGCGTGAGCCGGCATTACGACGAAGCCCTTGCCGGCTTCGGGGTCAATGTc                                                                                               >  7:155646/1‑149 (MQ=255)
                                                gTGAAACTTCTTGCAACACTCCGTACGCGGGTATATACACGTATCCATGCTGACCAGTGAATGCATCTGTACCCACCTGCGTCGCGCCGCCCGCGGCGTGAGCCGGCATTACGACGAAGCCCTTGCCGGCTTCGGGGTCAATGTCGccc                                                                                           <  8:155646/149‑1 (MQ=255)
                                                gTGAAACTTCTTGCAACACTCCGTACGCGGGTATATACACGTATCCATGCTGACCAGTGAATGCATCTGTACCCACCTGCGTCGCGCCGCCCGCGGCGTGAGCCGGCATTACGACGAAGCCCTTGCCGGCTTCGGGGTCAATGTCGcc                                                                                            >  7:84962/1‑148 (MQ=255)
                                                 tGAAACTTCTTGCAACACTCCGTACGCGGGTATATACACGTATCCATGCTGACCAGTGAATGCATCTGTACCCACCTGCGTCGCGCCGCCCGCGGCGTGAGCCGGCATTACGACGAAGCCCTTGCCGGCTTCGGGGTCAATGTCGcc                                                                                            >  5:168183/1‑147 (MQ=255)
                                                 tGAAACTTCTTGCAACACTCCGTACGCGGGTATATACACGTATCCATGCTGACCAGTGAATGCATCTCTACCCACCTGCGTCGCGCCGCCCGCGGCGTGAGCCGGCATTACGACGAAGc                                                                                                                          <  7:334980/119‑1 (MQ=255)
                                                 tGAAACTTCTTGCAACACTCCGTACGCGGGTATATACACGTATCCATGCTGACCAGTGAATGCATCTATACCCACCTGCGTCGCGCCGCCCGCGGCGTGAGCCGGCATTACGACGAAGc                                                                                                                          >  8:334980/1‑119 (MQ=255)
                                                  gAAACTTCTTGCAACACTCCGTACGCGGGTATATACACGTATCCATGCTGACCAGTGAATGCATCTGTACCCACCTGCGTCGCGCCGCCCGCGGCGTGAGCCGGCATTACGACGAAg                                                                                                                           <  3:262997/117‑1 (MQ=255)
                                                  gAAACTTCTTGCAACACTCCGTACGCGGGTATATACACGTATCCATGCTGACCAGTGAATGCATCTGTACCCACCTGCGTCGCGCCGCCCGCGGCGTGAGCCGGCATTACGACGAAg                                                                                                                           >  4:262997/1‑117 (MQ=255)
                                                         cttGCAACACTCCGTACGCGGGTATATACACGTATCCATGCTGACCAGTGAATGCATTTGTACCCACCTGCGTCGCGCCGCCCGCGGCGTGAGCCGGCATTACGACGAAGCCCTTGCCGGCTTCGGGGTCAATGTCGCCCAGTTTTcc                                                                                   <  2:274838/148‑1 (MQ=255)
                                                         cttGCAACACTCCGTACGCGGGTATATACACGTATCCATGCTGACCAGTGAATGCATCTGTACCCACCTGCGTCGCGCCGCCCGCGGCGTGAGCCGGCATTACGACGAAGCCCTTGCCGGCTTCGGGGTCAATGTCGCCCAGTTTTc                                                                                    <  3:49694/147‑1 (MQ=255)
                                                            gCAACACTCCGTACGCGGGTATATACACGTATCCATGCTGACCAGTGAATGCATCTGTACCCACCTGCGTCGCGCCGCCCGCGGCGTGAGCCGGCATTACGACGAAGCCCTTGCCGGCTTCGGGGTCAAt                                                                                                  <  7:32503/130‑1 (MQ=255)
                                                            gCAACACTCCGTACGCGGGTATATACACGTATCCATGCTGACCAGTGAATGCATCTGTACCCACCTGCGTCGCGCCGCCCGCGGCGTGAGCCGGCATTACGACGAAGCCCTTGCCGGCTTCGGGGTCAAt                                                                                                  >  8:32503/1‑130 (MQ=255)
                                                             cAACACTCCGTACGCGGGTATATACACGTATCCATGCTGACCAGTGAATGCATCTGTACCCACCTGCGTCGCGCCGCCCGCGGCGTGAGCCGGCATTACGACGAAGCCCTTGCCGGCTTCGGGGTCAATGTCGCCCAGTTTTCCctgc                                                                               <  6:168183/148‑1 (MQ=255)
                                                              aaCACTCCGTACGCGGGTATATACACGTATCCATGCTGACCAGTGAATGCATCTGTACCCACCTGCGTCGCGCCGCCCGCGGCGTGAGCCGGCATTACGACGAAGCCCTTGCCGGCTTCGGGGTCAATGTCGcc                                                                                            <  6:4085/134‑1 (MQ=255)
                                                              aaCACTCCGTACGCGGGTATATACACGTATCCATGCTGACCAGTGAATGCATCTGTACCCACCTGCGTCGCGCCGCCCGCGGCGTGAGCCGGCATTACGACGAAGCCCTTGCCGGCTTCGGGGTCAATGTCGcc                                                                                            >  5:4085/1‑134 (MQ=255)
                                                               acacTCCGTACGCGGGTATATACACGTATCCATGCTGACCAGTGAATGCATCTGTACCCACCTGCGTCGCGCCGCCCGCGGCGTGAGCCGGCATTACGACGAAGCCCTTGCCGGCTTCGGGGTCAATGTCGCCCAGTTTTCCCtgctgc                                                                            <  5:273697/149‑1 (MQ=255)
                                                                cacTCCGTACGCGGGTATATACACGTATCCATGCTGACCAGTGAATGCATCTGTACCCACCTGCGTCGCGCCGCCCGCGGCGTGAGCCGGCATTACGACGAAGCCCTTGCCGGCTTCGGGGTCAATGTCGCCCAGTTTTCCCTGCTGCg                                                                           <  8:84962/149‑1 (MQ=255)
                                                                   tCCGTACGCGGGTATATACACGTATCCATGCTGACCAGTGAATGCATCTGTACCCACCTGCGTCGCGCCGCCCGCGGCGTGAGCCGGCATTACGACGAAGCCCTTGCCGGCTTCGGGGTCAATGTCGCCCAGTTTTCCCTGCTGCGGCa                                                                        >  2:100126/1‑149 (MQ=255)
                                                                         cgcgGGTATATACACGTATCCATGCTGACCAGTGAATGCATCTGTACCCACCTGCGTCGCGCCGCCCGCGGCGTGAGCCGGCATTACGACGAAGCCCTTGCCGGCTTCGGGGTCAATGTCGCCCAGTTTTCCCTGCTGCGGCACCTGCa                                                                  <  3:69382/149‑1 (MQ=255)
                                                                                 tataCACGTATCCATGCTGACCAGTGAATGCATCTGTACCCACCTGCGTCGCGCCGCCCGCGGCGTGAGCCGGCATTACGACGAAGCCCTTGCCGGCTTCGGGGTCAATGTCGCCCAGTTTTCCCTGCTGCGGCACCTGCAGCGGCTc                                                           <  1:210125/148‑1 (MQ=255)
                                                                                      acGTATCCATGCTGACCAGTGAATGCATCTGTACCCACCTGCGTCGCGCCGCCCGCGGCGTGAGCCGGCATTACGACGAAGCCCTTGCCGGCTTCGGGGTCAATGTCGc                                                                                             >  1:218409/1‑109 (MQ=255)
                                                                                      acGTATCCATGCTGACCAGTGAATGCATCTGTACCCACCTGCGTCGCGCCGCCCGCGGCGTGAGCCGGCATTACGACGAAGCCCTTGCCGGCTTCGGGGTCAATGTCGc                                                                                             <  2:218409/109‑1 (MQ=255)
                                                                                         tATCCATGCTGACCAGTGAATGCATCTGTACCCACCTGCGTCGCGCCGCCCGCGGCGTGAGCCGGCATTACGACGAAGCCCTTGCCGGCTCTCGGGGTCAATGTCGCCCAGTTTTCCCTGCTGCGGCACCTGCAGCGGCTCgaccgacc                                                   <  1:100126/149‑1 (MQ=255)
                                                                                                      cAGTGAATGCATCTGTACCCACCTGCGTCGCGCCGCCCGCGGCGTGAGCCGGCATTACGACGAAGCCCTTGCCGGCTTCGGGGTCAATGTCGCCCAGTTTTCCCTGCTGCGGCACCTGCAGCGGCTCGACCGACCCAGCATCACAAccc                                     >  5:239712/1‑149 (MQ=255)
                                                                                                      cAGTGAATGCATCTGTACCCACCTGCGTCGCGCCGCCCGCGGCGTGAGCCGGCATTACGACGAAGCCCTTGCCGGCTTCGGGGTCAATGTCGCCCAGTTTTCCCTGCTGCGGCACCTGCAGCGGCTCGACCGACCCAGCATCACAAccc                                     >  6:145835/1‑149 (MQ=255)
                                                                                                       aGTGAATGCATCTGTACCCACCTGCGTCGCGCCGCCCGCGGCGTGAGCCGGCATTACGACGAAGCCCTTGCCGGCTTCGGGGTCAATGTCGCCCAGTTTTCCCTGCTGCGGCACCTGCAGCGGCTCGACCGACCCAGCATCACAACCCt                                    >  4:44814/1‑149 (MQ=255)
                                                                                                        gTGAATGCATCTGTACCCACCTGCGTCGCGCCGCCCGCGGCGTGAGCCGGCATTACGACGAAGCCCTTGCCGGCTTCGGGGTCAATGTCGCCCAGTTTTCCCTGCTGCGGCACCTGCAGCGGCTCGACCGACCCAGCATCACAACCCTg                                   <  5:145835/149‑1 (MQ=255)
                                                                                                         tGAATGCATCTGTACCCACCTGCGTCGCGCCGCCCGCGGCGTGAGCCGGCATTACGACGAAGCCCTTGCCGGCTTCGGGGTCAATGTCGCCCAGTTTTCCCTGCTGCGGCACCTGCAGCGGCTCGACCGACCCAGCATCACAACCCTgg                                  <  2:194248/149‑1 (MQ=255)
                                                                                                           aaTGCATCTGTACCCACCTGCGTCGCGCCGCCCGCGGCGTGAGCCGGCATTACGACGAAGCCCTTGCCGGCTTCGGGGTCAATGTCGCCCAGTTTTCCCTGCTGCGGCACCTGCAGCGGCTCGACCGACCCAGCATCACAACCCTGGcc                                >  6:165142/1‑149 (MQ=255)
                                                                                                              gCATCTGTACCCACCTGCGTCGCGCCGCCCGCGGCGTGAGCCGGCATTACGACGAAGCCCTTGCCGGCTTCGGGGTCAATGTCGCCCAGTTTTCCCTGCTGCGGCACCTGCAGCGGCTCGACCGAccc                                                  >  1:303409/1‑128 (MQ=255)
                                                                                                              gCATCTGTACCCACCTGCGTCGCGCCGCCCGCGGCGTGAGCCGGCATTACGACGAAGCCCTTGCCGGCTTCGGGGTCAATGTCGCCCAGTTTTCCCTGCTGCGGCACCTGCAGCGGCTCGACCGAccc                                                  <  2:303409/128‑1 (MQ=255)
                                                                                                              gCATCTGTACCCACCTGCGTCGCGCCGCCCGCGGCGTGAGCCGGCATTACGACGAAGCCCTTGCCGGCTTCGGGGTCAATGTCGCCCAGTTTTCCCTGCTGCGGCACCTGCAGCGGCTCGACCGACCCAGCATCACAACCCTGGCCg                               >  5:266517/1‑147 (MQ=255)
                                                                                                               cATCTGTACCCACCTGCGTCGCGCCGCCCGCGGCGTGAGCCGGCATTACGACGAAGCCCTTGCCGGCTTCGGGGTCAATGTCGCCCAGTTTTCCCTGCTGCGGCACCTGCAGCGGCTCGACCGACCCAGCATCACAACCCTGGCCGaa                             <  6:239712/148‑1 (MQ=255)
                                                                                                                     tACCCACCTGCGTCGCGCCGCCCGCGGCGTGAGCCGGCATTACGACGAAGCCCTTGCCGGCTTCGGGGTCAATGTCGCCCAGTTTTCCCTGCTGCGGCACCTGCAGCGGCTCGACCGACCCAGCATCACAACCCTGGCCGAAGCCATgg                      <  5:165142/149‑1 (MQ=255)
                                                                                                                       cccACCTGCGTCGCGCCGCCCGCGGCGTGAGCCGGCATTACGACGAAGCCCTTGCCGGCTTCGGGGTCAATGTCGCCCAGTTTTCCCTGCTGCGGCACCTGCAGCGGCTCGACCGACCCAGCATCACAACCCTGGCCGAAGCCATgggg                    >  7:72670/1‑149 (MQ=255)
                                                                                                                        ccACCTGCGTCGCGCCGCCCGCGGCGTGAGCCGGCATTACGACGAAGCCCTTGCCGGCTTCGGGGTCAATGTCGCCCAGTTTTCCCTGCTGCGGCACCTGCAGCGGCTCGACCGACCCAGCATCACAACCCTGGCCGAAGCCATgggg                    <  4:80739/148‑1 (MQ=255)
                                                                                                                        ccACCTGCGTCGCGCCGCCCGCGGCGTGAGCCGGCATTACGACGAAGCCCTTGCCGGCTTCGGGGTCAATGTCGCCCAGTTTTCCCTGCTGCGGCACCTGCAGCGGCTCGACCGACCCAGCATCACAACCCTGGCCGAAGCCATgggg                    >  8:181180/1‑148 (MQ=255)
                                                                                                                        ccACCTGCGTCGCGCCGCCCGCGGCGTGAGCCGGCATTACGACGAAGCCCTTGCCGGCTTCGGGGTCAATGTCGCCCAGTTTTCCCTGCTGCGGCACCTGCAGCGGCTCGACCGACCCAGCATCACAACCCTGGCCGAAGCCATgggg                    >  3:80739/1‑148 (MQ=255)
                                                                                                                          aCCTGCGTCGCGCCGCCCGCGGCGTGAGCCGGCATTACGACGAAGCCCTTGCCGGCTTCGGGGTCAATGTCGCCCAGTTTTCCCTGCTGCGGCACCTGCAGCGGCTCGACCGACCCAGCATCACAACCCTGGCCGa                              <  4:148245/136‑1 (MQ=255)
                                                                                                                          aCCTGCGTCGCGCCGCCCGCGGCGTGAGCCGGCATTACGACGAAGCCCTTGCCGGCTTCGGGGTCAATGTCGCCCAGTTTTCCCTGCTGCGGCACCTGCAGCGGCTCGACCGACCCAGCATCACAACCCTGGCCGa                              >  3:148245/1‑136 (MQ=255)
                                                                                                                             cggctcacgccgcGGGCGGCGCGAGCCGGCATTACGACGAAGCCCTTGCCGGCTTCGGGGTCAATGTCGCCCAGTTTTCCCTGCTGCGGCACCTGCAGCGGCTCGACCGACCCAGCATCACAACCCTGGCCGAAGCCATGGGGCTCGAg              >  2:192537/8‑149 (MQ=255)
                                                                                                                              ggctcacgccgcGGGCGGCGCGAGCCGGCATTACGACGAAGCCCTTGCCGGCTTCGGGGTCAATGTCGCCCAGTTTTCCCTGCTGCGGCACCTGCAGCGGCTCGACCGACCCAGCATCACAACCCTGGCCGAAGCCATGGGGCTCGAgc             <  1:192537/143‑1 (MQ=255)
                                                                                                                               gctcacgccgcGGGCGGCGCGAGCCGGCATTACGACGAAGCCCTTGCCGGCTTCGGGGTCAATGTCGCCCAGTTTTCCCTGCTGCGGCACCTGCAGCGGCTCGACCGACCCAGCATCACAACCCTGGCCGAAGCCATGGGGCTCGAgcg            >  7:219250/6‑149 (MQ=255)
                                                                                                                                ctcacgccgcGGGCGGCGTGAGCCGGCATTACGACGAAGCCCTTGCCGGCTTCAGGGTCAATGTCGCCCAGTTTTCCCTGCTGCGGCACCTGCAGCGGCTCGACCGACCCAGCTTCACAACCCTGGCCGAAGCCATGGGGCTCGAgcg            >  5:101701/5‑148 (MQ=255)
                                                                                                                                     gccgccCGCGGCGTGAGCCGGCATTACGACGAAGCCCTTGCCGGGCTTCGGGGTCAATGTCGCCCAGTTTTCCCTGCTGCGGCCCCTGCAGCGGCTCGACCGCCCCAGCAGCACAACCCTGGCCGAAGCCATGGGGCTCGAGCGCAGtc       >  1:242392/1‑148 (MQ=255)
                                                                                                                                         cccgcGGCGTGACCCGGCATTACGACGAAGCCCTTGCCGGCTTCGGGGTCAATGTCGCCCAGTTTTCCCTGCTGCGGCACCTGCAGCGGCTCGACCGACCCAGCATCACAACCCTGGCCGAAGCCATGGGGCTCGAGCGCAGTACCCTg  >  6:76803/1‑149 (MQ=255)
                                                                                                                                          |                                                                                                                                                     
AGCCGATCCGCAGGCGCTTGAACCCGGACTGACGCTTCTCAATCTTTCGTGAAACTTCTTGCAACACTCCGTACGCGGGTATATACACGTATCCATGCTGACCAGTGAATGCATCTGTACCCACCTGCGTCGCGCCGCCCGCGGCGTGAGCCGGCATTACGACGAAGCCCTTGCCGGCTTCGGGGTCAATGTCGCCCAGTTTTCCCTGCTGCGGCACCTGCAGCGGCTCGACCGACCCAGCATCACAACCCTGGCCGAAGCCATGGGGCTCGAGCGCAGTACCCTG  >  NC_002947/1877336‑1877621

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: