Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,189,035 A→G 8.8% intergenic (+25/+17) PP_1038 → / ← tatC‑I hypothetical protein/Sec‑independent protein translocase protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,189,0350AG 8.8% 156.0 / 5.5 58intergenic (+25/+17)PP_1038/tatC‑Ihypothetical protein/Sec‑independent protein translocase protein
Reads supporting (aligned to +/- strand):  ref base A (26/26);  new base G (2/3);  total (28/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.68e-01

GGAAGAGTGGAAGGTGGAGCTGGTGGTGGCTGACGAGCTGATTGTGCAGGTCGTCGCTGCGCTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAACCCCATCTCGAACAGCACCCACATCGGCACTGCCAGCATCGTCTGCGAAAACACATCCGGCGGTGTCAG  >  NC_002947/1188887‑1189183
                                                                                                                                                    |                                                                                                                                                    
ggAAGAGTGGAAGGTGGAGCTGGTGGTGGCTGACGAGCTGATTGTGCAGGTCGTCGCTGCGCTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCa                                                                                                                                                      >  1:234812/1‑149 (MQ=255)
       tGGAAGGTGGAGCTGGTGGTGGCTGACGAGCTGATTGTGCAGGTCGTCGCTGCGCTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCcgc                                                                                                                                               <  8:180369/149‑1 (MQ=255)
         gAAGGTGGAGCTGGTGGTGGCTGACGAGCTGATTGTGCAGGTCGTCGCTGCGCTCAAGCAGAGACACCAGTAAGAGACGCAAGCCTATGAAGTCTGGCTGATCGCCGAGTTCTAGAACGAATCACGGCTATTCGCGGGCACGCCcgcg                                                                                                                                              >  7:156869/1‑148 (MQ=255)
           aGGTGGAGCTGGTGGTGGCTGACGAGCTGATTGTGCAGGTCGTCGCTGCGCTCAAGCAGAGGCAGCCGTACGAGACGCCAGCCTATGAATTCTGTCGGCTCGCCGAGTTCTAGAACTCATCACGGCTATTCGCGGGCGTGCCCACGaa                                                                                                                                            <  4:143210/148‑1 (MQ=255)
            ggTGGAGCTGGTGGTGGCTGACGAGCTGAATGTGCAGGTCGTCGCTGCGCTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCGTGCCCGCGaata                                                                                                                                          <  5:388331/149‑3 (MQ=255)
                gAGCTGGTGGTGGCTGACGAGCTGATTGTGCAGGTCGTCGCTGCGCTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGcc                                                                                                                                      <  1:144103/149‑1 (MQ=255)
                 aGCTGGTGGTGGCTGACGAGCTGATTGTGCAGGTCGTCGCTGCGCTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGccc                                                                                                                                     >  2:86618/1‑149 (MQ=255)
                   cTGGTGGTGGCTGACGAGCTGATTGTGCAGGTCGTCGCTGCGCTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTa                                                                                                                                   <  3:38463/149‑1 (MQ=255)
                         gtggCTGACGAGCTGATTGTGCAGGTCGTCGCTGCGCTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATc                                                                                                                              <  5:41488/148‑1 (MQ=255)
                         gtggCTGACGAGCTGATTGTGCAGGTCGTCGCTGCGCTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATc                                                                                                                              >  4:275537/1‑148 (MQ=255)
                          tggCTGACGAGCTGATTGTGCAGGTCGTCGCTGCGCTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCa                                                                                                                            <  3:306857/149‑1 (MQ=255)
                            ctatctgAGCTGATTGTGCAGGTCGTCGCTGCGCTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATTAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAg                                                                                                                           <  5:287760/142‑1 (MQ=255)
                              tGACGAGCTGATTGTGCAGGTCGTCGCTGCGCTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCca                                                                                                                        >  5:222813/1‑149 (MQ=255)
                                  gAGCTGATTGTGCAGGTCGTCGCTGCGCTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAAt                                                                                                                    >  1:175417/1‑149 (MQ=255)
                                  gAGCTGATTGTGCAGGTCGTCGCTGCGCTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCCAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAAt                                                                                                                    >  4:131347/1‑149 (MQ=255)
                                   aGCTGATTGTGCAGGTCGTCGCTGCGCTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAg                                                                                                                           >  6:287760/1‑141 (MQ=255)
                                        aTTGGGCCGGTCGTCGCGGCGCGCAAGCAGCGCCACCCGGACGAGACGCCAGCCTATGAAGGCGGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTcc                                                                                                                <  8:158494/147‑1 (MQ=255)
                                                    gtcgCTGCGCTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCt                                                                                                          >  2:247787/1‑141 (MQ=255)
                                                    gtcgCTGCGCTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCt                                                                                                          <  1:247787/141‑1 (MQ=255)
                                                    gtcgCTGCGCTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCg                                                                                                   <  7:113542/148‑1 (MQ=255)
                                                           cgcTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCt                                                                                             <  2:175417/147‑1 (MQ=255)
                                                            gcTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCGTGCCCGCGaat                                                                                                                                           <  3:174717/100‑2 (MQ=255)
                                                            gcTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCGTGCCCGCGaat                                                                                                                                           >  4:174717/1‑99 (MQ=255)
                                                             cTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCaa                                                                                         >  5:251184/1‑149 (MQ=255)
                                                                   cAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTcc                                                                                   <  2:234812/149‑1 (MQ=255)
                                                                       gCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACaa                                                                               <  6:222813/149‑1 (MQ=255)
                                                                        ccACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAg                                                                              >  2:108522/1‑149 (MQ=255)
                                                                            ccGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCca                                                                                  >  2:247378/1‑141 (MQ=255)
                                                                            ccGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCca                                                                                  <  1:247378/141‑1 (MQ=255)
                                                                            ccGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAg                                                                          <  3:131347/149‑1 (MQ=255)
                                                                             cGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACATTCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAg                                                                          >  4:180964/1‑148 (MQ=255)
                                                                              gTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGc                                                                         >  8:104358/1‑148 (MQ=255)
                                                                                  gagaCGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAAccc                                                                    >  7:226389/1‑149 (MQ=255)
                                                                                   agaCGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTCGCTCCCGCCGCTTCAACATTCCACAAGCCAGCAAccc                                                                    >  5:385714/1‑148 (MQ=255)
                                                                                      cGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAACCCCAtc                                                                <  7:104358/149‑1 (MQ=255)
                                                                                         cAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAACCCCAtctc                                                              >  3:116806/1‑148 (MQ=255)
                                                                                            ccTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAACCCCATCTCGaa                                                           >  4:69421/1‑148 (MQ=255)
                                                                                            ccTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATACAGCCACAATCTCCTTGCCTGGCTCCGGCCCCTTCAACATTCCACAAGCCAGCAACCCCATCTCGaa                                                           >  4:400819/1‑148 (MQ=255)
                                                                                               aTGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAACCCCATCTCGAACAg                                                        >  7:302446/1‑148 (MQ=255)
                                                                                                 gAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAACCCCATCTCGAACAGCAc                                                     >  3:402707/1‑149 (MQ=255)
                                                                                                  aaGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAACCCCATCTCGAACAGCAc                                                     <  8:347850/148‑1 (MQ=255)
                                                                                                         gcggcTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAACCCCATCTCGAACAGCACCCACATc                                              <  5:3214/148‑1 (MQ=255)
                                                                                                                gCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAACCCCATCTCGAACAGCACCCACATCGGCACt                                        <  3:69421/147‑1 (MQ=255)
                                                                                                                gCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAACCCCATCTCGAACAGCACCCACATCGGCACTGc                                      >  7:124806/1‑149 (MQ=255)
                                                                                                                 ccGATTTCTAGAACTCATCACTGCTAATCGCGGTCACGCCCGCGAAAGGTCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGTCCGCTTCAACATTCCACAAGCCAGCAACCCCATCTCGAACAGCACCCACATCGGCACTGc                                      <  8:156869/148‑1 (MQ=255)
                                                                                                                 ccGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAACCCCATCTCGAACAGCACCCACATCgg                                            >  1:178766/1‑142 (MQ=255)
                                                                                                                 ccGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAACCCCATCTCGAACAGCACCCACATCgg                                            <  2:178766/142‑1 (MQ=255)
                                                                                                                 ccGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAACCCCATCTCGAACAGCACCCACATCGGCACTGcc                                     <  1:86618/149‑1 (MQ=255)
                                                                                                                  cGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAACCCCATCTCGAACAGCACCCACATCGGCACTGCCa                                    <  8:124806/149‑1 (MQ=255)
                                                                                                                              cGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCACTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAACCCCATCTCGAACAGCACCCACATCGGCACTGCCAGCATCGTCTGc                         >  1:72276/1‑148 (MQ=255)
                                                                                                                                  tCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAACCCCATCTCGAACAGCACCCACATCGGCACTGCCAGCATCGTCTGCGaaaa                    <  3:275537/149‑1 (MQ=255)
                                                                                                                                    aCGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAACCCCATCTCGAACAGCACCCACATCGGCACTGCCAGCATCGTCTGCGAAAac                   <  8:302446/148‑1 (MQ=255)
                                                                                                                                     cGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAACCCCATCTCGAACAGCACCCACATCGGCACTGCCAGCATCGTCTGCGaaaa                    >  7:199544/1‑146 (MQ=255)
                                                                                                                                     cGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAACCCCATCTCGAACAGCACCCACATCGGCACTGCCAGCATCGTCTGCGaaaa                    <  8:199544/146‑1 (MQ=255)
                                                                                                                                         tATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAACCCCATCTCGAACAGCACCCACATCGGCACTGCCAGCATCGTCTGCGAAAACACATcc             >  1:628/1‑149 (MQ=255)
                                                                                                                                         ctttCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAACCCCATCTCGAACAGCACCCACATCGGCACTg                                       >  6:29486/3‑123 (MQ=255)
                                                                                                                                         ctttCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAACCCCATCTCGAACAGCACCCACATCGGCACTg                                       <  5:29486/121‑1 (MQ=255)
                                                                                                                                           ttCGCGGGCGTGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAACCCCATCTCGAACAGCACCCACATCGGCACTGCCAGCATCGTCTGCGAAAACACATCcgg           >  1:190151/1‑149 (MQ=255)
                                                                                                                                               cgGGCTCGCCCGCGAAAGTGCCCTACAATCCAGCCAGAATCTCCTTGCCTGGCTCCGGCCGTTTCAAGATTCCACAAGCCAGCAACCCCATCTCGAACATCACCCACATCGGCACTGCCAGCATCGTCTGCGAAAACACATCCGGCGgt       <  7:179393/149‑1 (MQ=255)
                                                                                                                                                    aCGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAACCCCATCTCGCACAGCACCCACATCGGCACTGCCAGCATCGTCTGCGAAAACACATCCGGCGGTGTCAg  >  3:326691/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GGAAGAGTGGAAGGTGGAGCTGGTGGTGGCTGACGAGCTGATTGTGCAGGTCGTCGCTGCGCTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAACCCCATCTCGAACAGCACCCACATCGGCACTGCCAGCATCGTCTGCGAAAACACATCCGGCGGTGTCAG  >  NC_002947/1188887‑1189183

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: