Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,877,475 C→G 9.4% R16G (CGC→GGC)  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,4750CG 9.4% 192.7 / 9.1 64R16G (CGC→GGC) PP_1683MarR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (27/31);  new base G (4/2);  total (31/33)
Fisher's exact test for biased strand distribution p-value = 4.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.32e-01

AGCCGATCCGCAGGCGCTTGAACCCGGACTGACGCTTCTCAATCTTTCGTGAAACTTCTTGCAACACTCCGTACGCGGGTATATACACGTATCCATGCTGACCAGTGAATGCATCTGTACCCACCTGCGTCGCGCCGCCCGCGGCGTGAGCCGGCATTACGACGAAGCCCTTGCCGGCTTCGGGGTCAATGTCGCCCAGTTTTCCCTGCTGCGGCACCTGCAGCGGCTCGACCGACCCAGCATCACAACCCTGGCCGAAGCCATGGGGCTCGAGCGCAGTACCCTGGG  >  NC_002947/1877336‑1877623
                                                                                                                                           |                                                                                                                                                      
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                      >  3:80739/1‑148 (MQ=255)
                                                                                                                        ccACCTGCGTCGCGCCGCCCGCGGCGTGAGCCGGCATTACGACGAAGCCCTTGCCGGCTTCGGGGTCAATGTCGCCCAGTTTTCCCTGCTGCGGCACCTGCAGCGGCTCGACCGACCCAGCATCACAACCCTGGCCGAAGCCATgggg                      >  8:181180/1‑148 (MQ=255)
                                                                                                                        ccACCTGCGTCGCGCCGCCCGCGGCGTGAGCCGGCATTACGACGAAGCCCTTGCCGGCTTCGGGGTCAATGTCGCCCAGTTTTCCCTGCTGCGGCACCTGCAGCGGCTCGACCGACCCAGCATCACAACCCTGGCCGAAGCCATgggg                      <  4:80739/148‑1 (MQ=255)
                                                                                                                          aCCTGCGTCGCGCCGCCCGCGGCGTGAGCCGGCATTACGACGAAGCCCTTGCCGGCTTCGGGGTCAATGTCGCCCAGTTTTCCCTGCTGCGGCACCTGCAGCGGCTCGACCGACCCAGCATCACAACCCTGGCCGa                                >  3:148245/1‑136 (MQ=255)
                                                                                                                          aCCTGCGTCGCGCCGCCCGCGGCGTGAGCCGGCATTACGACGAAGCCCTTGCCGGCTTCGGGGTCAATGTCGCCCAGTTTTCCCTGCTGCGGCACCTGCAGCGGCTCGACCGACCCAGCATCACAACCCTGGCCGa                                <  4:148245/136‑1 (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)
                                                                                                                                           cgcgGCGTGAGCCGGCATTACGACGAAGCCCTTGCCGGCTTCGGGGTCAATGTCGCCCAGTTTTCCCTGCTGCGGCACCTGCAGCGGCTCGACCGACCCAGCATCACAACCCTGGCCGAAGCCATGGGGCTCGAGCGCAGTACCCTggg  >  7:283950/1‑149 (MQ=255)
                                                                                                                                           |                                                                                                                                                      
AGCCGATCCGCAGGCGCTTGAACCCGGACTGACGCTTCTCAATCTTTCGTGAAACTTCTTGCAACACTCCGTACGCGGGTATATACACGTATCCATGCTGACCAGTGAATGCATCTGTACCCACCTGCGTCGCGCCGCCCGCGGCGTGAGCCGGCATTACGACGAAGCCCTTGCCGGCTTCGGGGTCAATGTCGCCCAGTTTTCCCTGCTGCGGCACCTGCAGCGGCTCGACCGACCCAGCATCACAACCCTGGCCGAAGCCATGGGGCTCGAGCGCAGTACCCTGGG  >  NC_002947/1877336‑1877623

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: