Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 438,176 C→T 6.2% A62T (GCC→ACC)  rdoA ← serine/threonine protein kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947438,1760CT 6.2% 189.7 / 2.0 65A62T (GCC→ACC) rdoAserine/threonine protein kinase
Reads supporting (aligned to +/- strand):  ref base C (29/31);  new base T (2/2);  total (31/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.35e-01

GAAGCCCTTGTGTTCGAACAAGGTGCGGTCATCGTGCTGCATGGGCGCAACCACCGGGACTTCGCAGTCGGCGAGTTCGGCCGTGAAGGCGTGCTCTTCGAGGATGGCTGCATCGCTCCAGCGGCCTGGGCGGTAGAACTTGGCAATCAACGGCTGCGCGTCTTCGATGCCCACCTGGTACACGCGGTTCTCGTAGCTGTTCAGCGCCAGCACCCGAGCGTCGCTGAGAAAGCCCAGGCTTTCCACTGCGTCCAGGACCAGGTCGGGGGTGAGGGTGTCGA  >  NC_002947/438033‑438313
                                                                                                                                               |                                                                                                                                         
gAAGCCCTTGTGTTCGAACAAGGTGCGGTCATCGTGCTGCATGGGCGCAACCACCGGGACTTCGCAGTCGGCGAGTTCGGCCGTGAAGGCGTGCTCTTCGAGGATGGCTGCATCGCTCCAGCGGCCTGGGCGGTAGAACTTGGCAATc                                                                                                                                       >  4:17035/1‑148 (MQ=255)
gAAGCCCTTGTGTTCGAACAAGGTGCGGTCATCGTGCTGCATGGGCGCAACCACCGGGACTTCGCAGTCGGCGAGTTCGGCCGTGAAGGCGTGCTCTTCGAGGATGGCTGCATCGCTCCAGCGGCCTGGGCGGTAGAACTTGGCAATc                                                                                                                                       >  2:4706/1‑148 (MQ=255)
     ccTTGTGTTCGAACAAGGTGCGGTCATCGTGCTGCATGGGCGCAACCACCGGGACTTCGCAGTCGGCGAGTTCGGCCGTGAAGGCGTGCTCTTCGAGGATGGCTGCATCGCTCCAGCGGCCTGGGCGGTAGAACTTGGCAATCAACGGc                                                                                                                                 >  1:42446/1‑149 (MQ=255)
        tgtgTTCGAACAAGGTGCGGTCATCGTGCTGCATGGGCGCAACCACCGGGACTTCGCAGTCGGCTAGTTCGGCCGTGAAGGCGTGCTCTTCGAGGATGGCTGCATCGCTCCAGCGGCCTGGGCGGTAGAACTTGGCAATCAACGGCTgc                                                                                                                              >  5:304002/1‑149 (MQ=255)
            ttCGAACAAGGTGCGGTCATCGTGCTGCATGGGCGCAACCACCGGGACTTCGCAGTCGGCGAGTTCGGCCGTGAAGGCGTGCTCTTCGAGGATGTCTGCATCGCTCCAGCGGCCTGGGCGGTAGAACTTGGCAATCAACGGCTgcgc                                                                                                                            <  8:221396/147‑1 (MQ=255)
            ttCGAACAAGGTGCGGTCATCGTGCTGCATGGGCGCAACCACCGGGACTTCGCAGTCGGCGAGTTCGGCCGTGAAGGCGTGCTCTTCGAGGATGGCTGCATCGCTCCAGCGGCCTGGGCGGTAGAACTTGGCAATCAACGGCTgcgc                                                                                                                            >  7:221396/1‑147 (MQ=255)
            ttCGAACAAGGTGCGGTCATCGTGCTGCATGGGCGCAACCACCGGGACTTCGCAGTCGGCGAGTTCGGCCGTGAAGGCGTGCTCTTCGAGGATGGCTGCATCGCTCCAGCGGCCTGGGCGGTAGAACTTGGCAATCAACGGCTGCGCGt                                                                                                                          >  3:348917/1‑149 (MQ=255)
            ttCGAACAAGGTGCGGTCATCGTGCTGCATGGGCGCAACCACCGGGACTTCGCAGTCGGCGAGTTCGGCCGTGAAGGCGTGCTCTTCGAGGATGGCTGCATCGCTCCAGCGGCCTGGGCGGTAGAACTTGGCAATCAACGGCTGCGCGt                                                                                                                          >  3:293050/1‑149 (MQ=255)
              cGAACAAGGTGCGGTCATCGTGCTGCATGGGCGCAACCACCGGGACTTCGCAGTCGGCGAGTTCGGCCGTGAAGGCGTGCTCTTCGAGGATGGCTGCATCGCTCCAGCGGCCTGGGCGGTAGAACTTGGCAATCAACGGCTGCGCGTc                                                                                                                         <  3:43490/148‑1 (MQ=255)
              cGAACAAGGTGCGGTCATCGTGCTGCATGGGCGCAACCACCGGGACTTCGCAGTCGGCGAGTTCGGCCGTGAAGGCGTGCTCTTCGAGGATGGCTGCATCGCTCCAGCGGCCTGGGCGGTAGAACTTGGCAATCAACGGCTGCGCGTc                                                                                                                         <  3:38593/148‑1 (MQ=255)
              cGAACAAGGTGCGGTCATCGTGCTGCATGGGCGCAACCACCGGGACTTCGCAGTCGGCGAGTTCGGCCGTGAAGGCGTGCTCTTCGAGGATGGCTGCATCGCTCCAGCGGCCTGGGCGGTAGAACTTGGCAATCAACGGCTGCGCGTCt                                                                                                                        >  7:224534/1‑149 (MQ=255)
               gAACAAGGTGCGGTCATCGTGCTGCATGGGCGCAACCACCGGGACTTCGCAGTCGGCGAGTTCGGCCGTGAAGGCGTGCTCTTCGAGGATGGCTGCATCGCTCCAGCGGCCTGGGCGGTAGAACTTGGCAATCAACGGCTGCGCGTCtt                                                                                                                       <  2:183741/149‑1 (MQ=255)
                      gTGCGGTCATCGTGCTGCATGGGCGCAACCACCGGGACTTCGCAGTCGGCGAGTTCGGCCGTGAAGGCGTGCTCTTCGAGGATGGCTGCATCGCTCCAGCGGCCTGGGCGGTAGAACTTGGCAATCAACGGCTGCGCGTCTTCGATGcc                                                                                                                >  6:246012/1‑149 (MQ=255)
                         gggtgttCGTGCTGGATGGGCGCAACGACCGGGAGTTCGCAGTCGGCGAGTTGGGCCGTGAAGGCGTGCTCTTCGAGGATTGCTGCATCGCTCCAGCGGCATGGGCTGTATAACTTGGGAATCAACGTCTGCGCTTATTCGATGCCCAc                                                                                                             <  8:201310/143‑1 (MQ=255)
                          ggTCATCGTGCTGCATGGGCGCAACCACCGGGACTTCGCAGTCGGCGAGTTCGGCCGTGAAGGCGTGCTCTTCGAGGATGGCTGCATCGCTCCAGCGGCCTGGGCGGTAGAACTTGGCAATCAACGGCTGCGCGTCATCGATGCCCAcc                                                                                                            >  6:231034/1‑149 (MQ=255)
                              aTCGTGCTGCATGGGCGCAACCACCGGGACTTCGCAGTCGGCGAGTTCGGCCGTGAAGGCGTGCTCTTCGAGGATGGCTGCATCGCTCCAGCGGCCTGGGCGGTAGAACTTGGCAATCAACGGCTGCGCGTCTTCGATGCCCACCTgg                                                                                                         >  5:295918/1‑148 (MQ=255)
                                  tgctgcATGGGCGCAACCACCGGGACTTCGCAGTCGGCGAGTTCGGCCGTGAAGGCGTGCTCTTCGAGGATGGCTGCATCGCTCCAGCGGCCTGGGCGGTAGAACTTGGCAATCAACGGCTGCGCGTCTTCGATGCCCACCTGGTacac                                                                                                    >  5:387913/1‑149 (MQ=255)
                                  tgctgcATGGGCGCAACCACCGGGACTTCGCAGTCGGCGAGTTCGGCCGTGAAGGCGTGCTCTTCGAGGATGGCTGCATCGCTCCAGCGGCCTGGGCGGTAGAACTTGGCAATCAACGGCTGCGCGTCTTCGATGCCCACCTGGTacac                                                                                                    >  7:266476/1‑149 (MQ=255)
                                       cATGGGCGCAACCACCGGGACTTCGCAGTCGGCGAGTTCGGCCGTGAAGGCGTGCTCTTCGAGGATGGCTGCATCGCTCCAGCGGCCTGGGCGGTAGAACTTGGCAATCAACGGCTGCGCGTCTTCGATGCCCACCTGGTACACGCgg                                                                                                <  4:348917/148‑1 (MQ=255)
                                       cATGGGCGCAACCACCGGGACTTCGCAGTCGGCGAGTTCGGCCGTGAAGGCGTGCTCTTCGAGGATGGCTGCATCGCTCCAGCGGCCTGGGCGGTAGAACTTGGCAATCAACGGCTGCGCGTCTTCGATGCCCACCTGGTACACGCGGt                                                                                               <  4:293671/149‑1 (MQ=255)
                                          gggCGCAACCACCGGGACTTCGCAGTCGGCGAGTTCGGCCGTGAAGGCGTGCTCTTCGAGGATGGCTGCATCGCTCCAGCGGCCTGGGCGGTAGAACTTGGCAATCAACGGCTGCGCGTCTTCGATGCCCACCTGGTACACGCGGTtct                                                                                            >  3:351326/1‑149 (MQ=255)
                                            gcgcAACCACCGGGACTTCGCAGTCGGCGAGTTCGGCCGTGAAGGCGTGCTCTTCGAGGATGGCTGCATCGCTCCAGCGGCCTGGGCGGTAGAACTTGGCAATCAACGGCTGCGCGTCTTCGATGCCCACCTGGTACACGCGGTTCTCg                                                                                          >  2:337391/1‑149 (MQ=255)
                                              gcAACCACCGGGACTTCGCAGTCGGCGAGTTCGGCCGTGAAGGCGTGCTCTTCGAGGATGGCTGCATCGCTCCAGCGGCCTGGGCGGTAGAACTTGGCAATCAACGGCTGCGCGTCTTCGATGCCCACCTGGTACACGCGGTTCTCGTa                                                                                        <  7:147204/149‑1 (MQ=255)
                                              gcAACCACCGGGACTTCGCAGTCGGCGAGTTCGGCCGTGAAGGCGTGCTCTTCGAGGATGGCTGCATCGCTCCAGCGGCCTGGGCGGTAGAACTAGGCAATAAACGGCTGCGCGTCTTCGATGCCCACCTGGTACACGCGGTTCTCg                                                                                          >  3:30481/1‑147 (MQ=255)
                                               cAACCACCGGGACTTCGCAGTCGGCGAGTTCGGCCGTGAAGGCGTGCTCTTCGAGGATGGCTGCATCGCTCCAGCGGCCTGGGCGGTAGAACTTGGCAATCAACGGCTGCGCGTCTTCGATGCCCACCTGGTACACGCGGTTCTCGTAg                                                                                       <  3:17035/149‑1 (MQ=255)
                                               cAACCACCGGGACTTCGCAGTCGGCGAGTTCGGCCGTGAAGGCGTGCTCTTCGAGGATGGCTGCATCGCTCCAGCGGCCTGGGCGGTAGAACTTGGCAATCAACGGCTGCGCGTCTTCGATGCCCACCTGGTACACGCGGTTCTCGTAg                                                                                       <  8:224534/149‑1 (MQ=255)
                                                aaCCACCGGGACTTCGCAGTCGGCGAGTTCGGCCGTGAAGGCGTGCTCTTCGAGGATGGCTGCATCGCTCCAGCGGCCTGGGCGGTAGAACTTGGTAATCAACGGCTGCGCGTCTTCGATGCCCACCTGGTACACGCGGTTCTCGTAGc                                                                                      >  4:148421/1‑149 (MQ=255)
                                                 accaccGGGACTTCGCAGTCGGCGAGTTCGGCCGTGAAGGCGTGCTCTTCGAGGATGGCTGCATCGCTCCAGCGGCCTGGGCGGTAGAACTTGGCAATCAACGGCTGCGCGTCTTCGATGCCCACCTGGTACACGCGGTTCTCGTAGCt                                                                                     <  6:295918/149‑1 (MQ=255)
                                                       gggACTTCGCAGTCGGCGAGTTCGGCCGTGAAGGCGTGCTCTTCGAGGATGGCTGCATCGCTCCAGCGGCCTGGGCGGTAGAACTTGGCAATCAACGGCTGCGCGTCTTCGATGCCCACCTg                                                                                                          >  4:295389/1‑122 (MQ=255)
                                                       gggACTTCGCAGTCGGCGAGTTCGGCCGTGAAGGCGTGCTCTTCGAGGATGGCTGCATCGCTCCAGCGGCCTGGGCGGTAGAACTTGGCAATCAACGGCTGCGCGTCTTCGATGCCCACCTg                                                                                                          <  3:295389/122‑1 (MQ=255)
                                                            ttCGCAGTCGGCGAGTTCGGCCGTGAAGGCGTGCTCTTCGAGGATGGCTGCATCGCTCCAGCGGCCTGGGCGGTAGAACTTGGCAATCAACGGCTGCGCGTCTTc                                                                                                                      >  5:43879/1‑105 (MQ=255)
                                                            ttCGCAGTCGGCGAGTTCGGCCGTGAAGGCGTGCTCTTCGAGGATGGCTGCATCGCTCCAGCGGCCTGGGCGGTAGAACTTGGCAATCAACGGCTGCGCGTCTTc                                                                                                                      <  6:43879/105‑1 (MQ=255)
                                                              cGCAGTCGGCGAGTTCGGCCGTGAAGGCGTGCTCTTCGAGGATGGCTGCATCGCTCCAGCGGCCTGGGCGGTAGAACTTGGCAATCAACGGCTGCGCGTCTTCGATGCCCACCTGGTACACGCGGTTCTCGTAGCTGTTCAGCGCCAg                                                                         >  3:280367/1‑148 (MQ=255)
                                                               gCAGTCGGCGAGTTCGGCCGTGAAGGCGTGCTCTTCGAGGATGGCTGCATCGCTCCAGCGGCCTGGGCGGTAGAACTTGGCAATCAACGGCTGCGCGTCTTCGATGCCCACCTGGTACACGCGGTTCTCGTAGCTGTTCAGCGCCAGc                                                                        <  2:42446/148‑1 (MQ=255)
                                                                  gTCGGCGAGTTCGGCCGTGAAGGCGTGCTCTTCGAGGATGGCTGCATCGCTCCAGCGGCCTGGGCGGTAGAACTTGGCAATCAACGGCTGCGCGTCTTCGATGCCCACCTGGTACACGCGGTtctc                                                                                           >  2:132656/1‑126 (MQ=255)
                                                                  gTCGGCGAGTTCGGCCGTGAAGGCGTGCTCTTCGAGGATGGCTGCATCGCTCCAGCGGCCTGGGCGGTAGAACTTGGCAATCAACGGCTGCGCGTCTTCGATGCCCACCTGGTACACGCGGTtctc                                                                                           <  1:132656/126‑1 (MQ=255)
                                                                   tCGGCGAGGTCGGCGGTGTAGGCGTGCTCTTCGAGGATGGCTGCATCGCTCCAGCGGCCTTGGCGGTAGAACTTGGCAATCAACGGCTGCGCGTCTTCGATGCCCACCTGGTACACGCGGTTCTCGTAGCTGTTCAGCGCCAGCAccc                                                                    <  4:30481/148‑1 (MQ=255)
                                                                      gCGAGTTCGGCCGTGAAGGCGTGCTCTTCGAGGATGGCTGCATCGCTCCAGCGGCCTGGGCGGTAGAACTTGGCAATCAACGGCTGCGCGTCTTCGATGCCCACCTGGTACACGCGGTTCTCGTAGCTGTTCAGCGCCAGCACCCGAGc                                                                <  6:387913/149‑1 (MQ=255)
                                                                      gCGAGTTCGGCCGTGAAGGCGTGCTCTTCGAGGATGGCTGCATCGCTCCAGCGGCCTGGGCGGTAGAACTTGGCAATCAACGGCTGCGCGTCTTCGATGCCCACCTGGTACACGCGGTTCTCGTAGCTGTTCAGCGCCAGCACCCGAGc                                                                <  8:266476/149‑1 (MQ=255)
                                                                       cGAGTTCGGCCGTGAAGGCGTGCTCTTCGAGGATGGCTGCATCGCTCCAGCGGCCTGGGCGGTAGAACTTGGCAATCAACGGCTGCGCGTCTTCGATGCCCACCTGGTACACGCGGTTCTCGTAGCTGTTCAGCGCCAGCACCCGAGCg                                                               <  5:311812/149‑1 (MQ=255)
                                                                           ttCGGCCGTGAAGGCGTGCTCTTCGAGGATGGCTGCATCGCTCCAGCGGCCTGGGCGGTAGAACTTGGCAATCAACGGCTGCGCGTCTTCGATGCCCACCTGGTACACGCGGTTCTCGTAGCTGTTCAGCGCCAGCACCCGAGCGTCg                                                            <  6:304002/148‑1 (MQ=255)
                                                                           ttCGGCCGTGAAGGCGTGCTCTTCGAGGATGGCTGCATCGCTCCAGCGGCCTGGGCGGTAGAACTTGGCAATCAACGGCTGCGCGTCTTCGATGCCCACCTGGTACACGCGGTTCTCGTAGCTGTTCAGCGCCAGCACCCGAGCGTCGc                                                           >  1:179798/1‑149 (MQ=255)
                                                                              ggCCGTGAAGGCGTGCTCTTCGAGGATGGCTGCATCGCTCCAGCGGCCTGGGCGGTAGAACTTGGCAATCAACGGCTGCGCGTCTTCGATGCCCACCTGGTACACGCGGTTCTCGTAGCTGTTCAGCGCCAGCACCCGAGCGTCGCt                                                          >  6:146916/1‑147 (MQ=255)
                                                                              ggCCGTGAAGGCGTGCTCTTCGAGGATGGCTGCATCGCTCCAGCGGCCTGGGCGGTAGAACTTGGCAATCAACGGCTGCGCGTCTTCGATGCCCACCTGGTACACGCGGTTCTCGTAGCTGTTCAGCGCCAGCACCCGAGCGTCGCTga                                                        <  1:337391/149‑1 (MQ=255)
                                                                                        gCGTGCTCTTCGAGGATGGCTGCATCGCTCCAGCGGCCTGGGCGGTAGAACTTGGCAATCAACGGCTGCGCGTCTTCGATGCCCACCTGGTACACGCGGTTCTCGTAGCTGTTCAGCGCCAGCACCCGAGCGTCGCTGAGAAAGCCCAg                                              <  5:246012/149‑1 (MQ=255)
                                                                                         cGTGCTCTTCGAGGATGGCTGCATCGCTCCAGCGGCCTGGGCGGTAGAACTTGGCAATCAACGGCTGCGCGTCTTCGATGCCCACCTGGTACACGCGGTTCTCGTAGCTGTTCAGCGCCAGCACCCGAGCGTCGCTGAGAAAGCCCAgg                                             <  6:18610/149‑1 (MQ=255)
                                                                                                ttCGAGGATGGCTGCATCGCTCCAGCGGCCTGGGCGGTAGAACTTGGCAATCAACGGCTGCGCGTCTTCGATGCCCACCTGGTACACGCGGTTCTCGTAGCTGTTCAGCGCCAGCACCCGAGCGTCGCTGAGAAAGCCCAGGCTTTCCa                                      >  3:260555/1‑149 (MQ=255)
                                                                                                  cGAGGATGGCTGCATCGCTCCAGCGGCCTGGGCGGTAGAACTTGGCAATCAACGGCTGCGCGTCTTCGATGCCCACCTGGTACACGCGGTTCTCGTAGCTGTTCAGCGCCAGCACCCGAGCGTCGCTGAGAAAGCCCAGGCTTTCCACt                                    >  2:166334/1‑149 (MQ=255)
                                                                                                  cGAGGATGGCTGCATCGCTCCAGCGGCCTGGGCGGTAGAACTTGGCAATCAACGGCTGCGCGTCTTCGAAGCCCACCTGGTACACGCGGTTCTCGTAGCTGTTCAGCGCCAGCACCCGAGCGTCGCTGAGAAAGCCCAgg                                             <  1:362502/140‑1 (MQ=255)
                                                                                                  cGAGGATGGCTGCATCGCTCCAGCGGCCTGGGCGGTAGAACTTGGCAATCAACGGCTGCGCGTCTTCGAAGCCCACCTGGTACACGCGGTTCTCGTAGCTGTTCAGCGCCAGCACCCGAGCGTCGCTGAGAAAGCCCAgg                                             >  2:362502/1‑140 (MQ=255)
                                                                                                   gAGGATGGCTGCATCGCTCCAGCGGCCTGGGCGGTAGAACTTGGCAATCAACTGCTGCGCGTCTTCGATGCCCACCTGGTACACGCGGTTCTCGTAGCTGTTCAGCGCCAGCGCCCGAGCGTCGCTGAGAAAGCCAAGGCTTTCCACTg                                   >  3:292412/1‑149 (MQ=255)
                                                                                                   gAGGATGGCTGCATCGCTCCAGCGGCCTGGGCGGTAGAACTTGGCAATCAACGGCTGCGCGTCTTCGATGCCCACCTGGTACACGCGGTTCTCGTAGCTGTTCAGCGCCAGCACCCGAGCGTCGCTGAGAAAGCCTGGGCTTTctcag                                    <  2:179798/148‑5 (MQ=255)
                                                                                                    aGGATGGCTGCATCGCTCCAGCGGCCTGGGCGGTAGAACTTGGCAATCAACGGCTGCGCGTCTTCGATGCCCACCTGGTACACGCGGTTCTCGTAGCTGTTCAGCGCCAGCACCCGAGCGTCGCTGAGAAAGCCCa                                               >  7:128470/1‑136 (MQ=255)
                                                                                                              cATCGCTCCAGCGGCCTGGGCGGTAGAACTTGGCAATCAACGGCTGCGCGTCTTCGATGCCCACCTGGTACACGCGGTTCTCGTAGCTGTTCAGCGCCAGCACCCGAGCGTCGCTGAGAAAGCCCAGGCTTTCCACTGCGTCCAGGa                          <  8:215949/147‑1 (MQ=255)
                                                                                                                tCGCTCCAGCGGCCTGGGCGGTAGAACTTGGCAATCAACGGCTGCGCGTCTTCGATGCCCACCTGGTACACGCGGTTCTCGTAGCTGTTCAGCGCCAGCACCCGAGCGTCGCTGAGAAAGCCCAGGCTTTCCACTGCGTCCAGGACCAg                      >  8:1254/1‑149 (MQ=255)
                                                                                                                  gCTCCAGCGGCCTGGGCGGTAGAACTTGGCAATCAACGGCTGCGCGTCTTCGATGCCCACCTGGTACACGCGGTTCTCGTAGCTGTTCAGCGCCAGCACCCGAGCGTCGCTGAGAAAGCCCAGGCTTTCCACTGCGTCCAGGACCAGGt                    <  4:280367/149‑1 (MQ=255)
                                                                                                                      cAGCGGCCTGGGCGGTAGAACTTGGTAATCAACGGCTGCGCGTCTTCGATGCCCACCTGGTACACGCGGTTCTCGTAGCTGTTCAGCGCCAGCACCCGAGCGTCGCTGAGAAAGCCCAGGCTTTCCACTGCGTCCAGGACCAGGTCgg                 <  3:148421/148‑1 (MQ=255)
                                                                                                                       aGCGGCCTGGGCGGTAGAACTTGGCAATCAACGGCTGCGCGTCTTCGATGCCCACCTGGTACACGCGGTTCTCGTAGCTGTTCAGCGCCAGCACCCGAGCGTCGCTGAGAAAGCCCAGGCTTTCCACTGCGTCCAGGACCAGGTCgggg               <  4:293050/149‑1 (MQ=255)
                                                                                                                         cGGCCTGGGCGGTAGAACTTGGCAATCAACGGCTGCGCGTCTTCGATGCCCACCTGGTACACGCGGTTCTCGTAGCTGTTCAGCGCCAGCACCCGAGCGTCGCTGAGAAAGCCCAGGCTTTCCACTGCGTCCAGGACCAGGTCGGGGGt             >  7:132683/1‑149 (MQ=255)
                                                                                                                            ccTGGGCGGTAGAACTTGGCAATCAACGGCTGCGCGTCTTCGATGCCCACCTGGTACACGCGGTTCTCGTAGCTGTTCAGCGCCAGCACCCGAGCGTCGCTGAGAAAGCCCAGGCTTTCCACTGCGTCCAGGACCAGGTCGGGGGTGAg          <  8:132683/149‑1 (MQ=255)
                                                                                                                             cTGGGCGGTAGAACTTGGTAATCAACGGCTGCGCGTCTTCGATGCCCACCTGGTACACGCGGTTCTCGTAGCTGTTCAGCGCCAGCACCCGAGCGTCGCTGAGAAAGCCCAGGCTTTCCACt                                    <  6:131019/122‑1 (MQ=255)
                                                                                                                             cTGGGCGGTAGAACTTGGTAATCAACGGCTGCGCGTCTTCGATGCCCACCTGGTACACGCGGTTCTCGTAGCTGTTCAGCGCCAGCACCCGAGCGTCGCTGAGAAAGCCCAGGCTTTCCACt                                    >  5:131019/1‑122 (MQ=255)
                                                                                                                             cTGGGCGGTAGAACTTGGCAATCAACGGCTGCGCGGCTTCGAGGCCCACCTGGTACACGCGGTTCTCGTAGCTGTTCAGCGCCAGCACCCGAGCGTCGCTGAGAAAGCCCAGGCTTTCCACTGCGTCCAGGACCAGGTCGGGGGTGAgg         <  3:49096/149‑1 (MQ=255)
                                                                                                                                ggCGGTAGAACTTGGCAATCAACGGCTGCGCGTCTTCGATGCCCACCTGGTACACGCGGTTCTCGTAGCTGTTCAGCGCCAGCACCCGAGCGTCGCTGAGAAAGCCCAGGCTTTCCACTGCGTCCAGGACCAGGTCGGGGGTGAGGgt       <  1:166334/148‑1 (MQ=255)
                                                                                                                                    gTAGAACTTGGCAATCAACGGCTGCGCGTCTTCGATGCCCACCTGGTACACGCGGTTCTCGTAGCTGTTCAGCGCCAGCACCCGAGCGTCGCTGAGAAAGCCCAGGCTTTCCACTGCGTCCAGGACCAGGTCGGGGGTGAGGGTGTCGa  <  5:231034/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                         
GAAGCCCTTGTGTTCGAACAAGGTGCGGTCATCGTGCTGCATGGGCGCAACCACCGGGACTTCGCAGTCGGCGAGTTCGGCCGTGAAGGCGTGCTCTTCGAGGATGGCTGCATCGCTCCAGCGGCCTGGGCGGTAGAACTTGGCAATCAACGGCTGCGCGTCTTCGATGCCCACCTGGTACACGCGGTTCTCGTAGCTGTTCAGCGCCAGCACCCGAGCGTCGCTGAGAAAGCCCAGGCTTTCCACTGCGTCCAGGACCAGGTCGGGGGTGAGGGTGTCGA  >  NC_002947/438033‑438313

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: ‑