Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,070,246:1 +G 88.1% intergenic (+46/+96) aroP‑I → / ← PP_0928 aromatic amino acid transport protein/K+‑dependent Na+/Ca+ exchanger‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,070,2461.G88.1% 171.6 / 18.9 59intergenic (+46/+96)aroP‑I/PP_0928aromatic amino acid transport protein/K+‑dependent Na+/Ca+ exchanger‑like protein
Reads supporting (aligned to +/- strand):  ref base . (2/5);  new base G (26/24);  total (29/30)
Fisher's exact test for biased strand distribution p-value = 4.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.30e-01

GCCGCCGATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGT‑‑GATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTA  >  NC_002947/1070104‑1070377
                                                                                                                                                 |                                                                                                                                    
gccgccGATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGT‑‑Gagatc                                                                                                                               <  2:373308/149‑5 (MQ=255)
  cgccgATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGT‑‑Gagatc                                                                                                                               >  8:220714/1‑143 (MQ=255)
   gccgATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGT‑‑Gagatcgag                                                                                                                            <  6:72512/149‑8 (MQ=255)
   gccgATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGT‑‑Gagatcgag                                                                                                                            <  8:164251/149‑8 (MQ=255)
         tCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGGTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGATCGAGcgccg                                                                                                                         <  5:207272/148‑1 (MQ=255)
              gatcGGGAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGGGGCTTTCAAGCTGAAGAAGGCGCGCCAGGCCAAGGGCTGATCGGGTGTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGATCGAGCGCCGCgcggg                                                                                                                   <  8:143157/147‑1 (MQ=255)
                  ggTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCCGTGA‑‑GATCGAGCGCCGCGCGggcggc                                                                                                               >  2:303939/1‑149 (MQ=255)
                         cTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGAt                                                                                                        >  6:337199/1‑149 (MQ=255)
                            aTCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTc                                                                                                           <  4:276101/143‑1 (MQ=255)
                            aTCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTc                                                                                                           >  3:276101/1‑143 (MQ=255)
                            aTCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGATCGAGCGCCGCCCGCGCGgcgc                                                                                                             <  8:121973/141‑1 (MQ=38)
                            aTCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGATCGAGCGCCGCCCGCGCGgcgc                                                                                                             >  7:121973/1‑141 (MQ=38)
                                cGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGGTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAg                                                                                                 <  3:61208/149‑1 (MQ=255)
                                     tGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTcac                                                                                                   >  4:328208/1‑142 (MQ=255)
                                     tGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTcac                                                                                                   <  3:328208/142‑1 (MQ=255)
                                           gctgtgctgtGGGTGGCTTTCAAGCTGAAGAAGGCGCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGc                                                                                      <  8:342015/149‑1 (MQ=255)
                                               tgctgtgGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTa                                                                                   <  8:46910/148‑1 (MQ=255)
                                                 ctgtgGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAg                                                                                <  7:220714/149‑1 (MQ=255)
                                                 ctgtgGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAg                                                                                >  1:63157/1‑149 (MQ=255)
                                                            tttCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTcag                                                                     >  7:336136/1‑147 (MQ=255)
                                                            tttCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTg                                                                     >  7:244707/1‑149 (MQ=255)
                                                            tttCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTg                                                                     >  8:273762/1‑149 (MQ=255)
                                                            tttCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTg                                                                     >  4:185696/1‑149 (MQ=255)
                                                                  gCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGcc                                                               >  1:38312/1‑149 (MQ=255)
                                                                     gaagaaGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCt                                                             >  4:361270/1‑148 (MQ=255)
                                                                      aagaagGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGa                                                                                                         >  4:220352/1‑103 (MQ=255)
                                                                      aagaagGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGa                                                                                                         <  3:220352/103‑1 (MQ=255)
                                                                       agaagGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGa                                                                                                         >  6:423294/1‑102 (MQ=255)
                                                                       agaagGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGa                                                                                                         <  5:423294/102‑1 (MQ=255)
                                                                       agaagGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTACGCCGTTAGGCTTTCATCGCC‑‑tgt                                                                                                                                       >  4:178829/1‑72 (MQ=38)
                                                                       agaagGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTACGCCGTTAGGCTTTCATCGCC‑‑tgt                                                                                                                                       <  3:178829/72‑1 (MQ=38)
                                                                         aagGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGAcgc                                                        <  2:63157/149‑1 (MQ=255)
                                                                         aagGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGATCGAGCGCCGCGCCGGCGGCGCGCCATCTCACAGGCGATTAAAGCCTTAAGGCAAGCGTCTGGGGGCCCTggagc                                                        >  1:186168/1‑145 (MQ=37)
                                                                            gCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGtt                                                     >  4:326462/1‑149 (MQ=255)
                                                                               cGCCAGGCCAAGGGCTGATCGGGTTGGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCGGTGA‑‑GATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTccc                                                  <  5:363818/149‑1 (MQ=38)
                                                                                gCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCCGTGA‑‑GATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACTGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTccc                                                  <  1:303939/148‑1 (MQ=38)
                                                                                gCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGATCGAGCGCCGCCCGCGCGGCGCTc                                                                                                           <  2:158840/91‑1 (MQ=25)
                                                                                gCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGATCGAGCGCCGCCCGCGCGGCGCTc                                                                                                           >  1:158840/1‑91 (MQ=25)
                                                                                  cAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCtt                                               <  5:337199/149‑1 (MQ=255)
                                                                                  cAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCtt                                               <  5:35426/149‑1 (MQ=255)
                                                                                   aGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTc                                              >  4:115491/1‑149 (MQ=255)
                                                                                   aGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTc                                              <  2:38312/149‑1 (MQ=255)
                                                                                      ccAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCAt                                           >  7:313825/1‑149 (MQ=255)
                                                                                        aaGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGc                                         <  3:77835/149‑1 (MQ=255)
                                                                                          gggCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGccc                                       <  8:244707/149‑1 (MQ=255)
                                                                                           ggCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGcccc                                      >  5:93396/1‑149 (MQ=255)
                                                                                                 tCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGGTGGATCTGACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGGCCTGACGCGGTGCCCATCCAATCCCCTCaaag                                >  5:284528/1‑147 (MQ=37)
                                                                                                   gggTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAAtgct                              <  3:361270/149‑1 (MQ=255)
                                                                                                   gggTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAAtgct                              >  6:113459/1‑149 (MQ=255)
                                                                                                   gggTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAAtgct                              <  5:113459/149‑1 (MQ=255)
                                                                                                          ttCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTcc                        <  3:115491/148‑1 (MQ=255)
                                                                                                              tAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCaa                   <  1:8564/149‑1 (MQ=255)
                                                                                                                  gggTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGcccc                                      <  7:392144/126‑1 (MQ=255)
                                                                                                                  gggTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGcccc                                      >  8:392144/1‑126 (MQ=255)
                                                                                                                  gggTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAgcac               <  7:273762/149‑1 (MQ=255)
                                                                                                                   ggTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAgcacg              >  4:240392/1‑149 (MQ=255)
                                                                                                                        gCCGTTAGGTTTTCGGCGCC‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCGCAAtgctgt                            >  4:330611/1‑129 (MQ=38)
                                                                                                                          cGTTAGGTTTTCGGCGCC‑‑TGTCAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAgcac               >  7:4496/1‑141 (MQ=255)
                                                                                                                          cGTTAGGTTTTCGGCGCC‑‑TGTCAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAgcac               <  8:4496/141‑1 (MQ=255)
                                                                                                                                 ttttCGGCGCC‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTa  >  2:314610/1‑147 (MQ=255)
                                                                                                                                 ttttCGGCGCC‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTa  >  3:357394/1‑147 (MQ=255)
                                                                                                                                                 |                                                                                                                                    
GCCGCCGATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGT‑‑GATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTA  >  NC_002947/1070104‑1070377

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: