Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,526,736 C→A 6.7% A308E (GCA→GAA)  ddlB → D‑alanine‑‑D‑alanine ligase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,526,7360CA 6.7% 190.6 / 3.7 60A308E (GCA→GAA) ddlBD‑alanine‑‑D‑alanine ligase B
Reads supporting (aligned to +/- strand):  ref base C (29/27);  new base A (2/2);  total (31/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.85e-01

GGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCGATGATACGTCAGCAGCAACCCGTTACCGGCCGTAGCAAGCCGGTGCCGCGTGGTGCCAGCCG‑GCTGGTGGCCGACGAGCCCGTATCGGCGCGCCTG  >  NC_002947/1526591‑1526883
                                                                                                                                                 |                                                                                                                                                    
ggCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAAt                                                                                                                                                   >  7:42959/1‑149 (MQ=255)
      gggCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGc                                                                                                                                             <  7:148532/149‑1 (MQ=255)
       ggCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGcc                                                                                                                                            <  2:172406/149‑1 (MQ=255)
                aCGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTg                                                                                                                                   >  2:9377/1‑149 (MQ=255)
                    gATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGGTCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGa                                                                                                                               >  2:102190/1‑149 (MQ=255)
                     aTGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGa                                                                                                                               <  3:317557/148‑1 (MQ=255)
                       gCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGAcgc                                                                                                                            >  4:20350/1‑149 (MQ=255)
                       gCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGAcgc                                                                                                                            >  3:204/1‑149 (MQ=255)
                        cAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCCCCCGGCCTGACAGATCCTAGCCTGGGGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGAcgcg                                                                                                                           <  5:103821/149‑1 (MQ=255)
                          ggACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAg                                                                                                                         <  7:121082/149‑1 (MQ=255)
                          ggACGAGCAGGGCCCGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATCGCGCCCCGTGCGGATGGCCTGGACTTCCAGCAACTGCTGCTGGCAATCCTGGCCGAAAGCGTGGCGAAGCGAg                                                                                                                         >  6:332352/1‑149 (MQ=255)
                           gACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAgg                                                                                                                        >  6:289538/1‑149 (MQ=255)
                                gCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTtaat                                                                                                                   <  6:10065/149‑1 (MQ=255)
                                 cAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTtaat                                                                                                                   >  6:182799/1‑148 (MQ=255)
                                     gCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCa                                                                                                              <  4:165277/149‑1 (MQ=255)
                                                  gctCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCgat                                                                                                 >  1:227904/1‑149 (MQ=255)
                                                    tCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCgatga                                                                                               >  3:54845/1‑149 (MQ=255)
                                                     cGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCgatga                                                                                               >  1:136906/1‑148 (MQ=255)
                                                      gAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCgatgat                                                                                              <  1:9377/148‑1 (MQ=255)
                                                       aaGTCCACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGCCTTCCAGCAACTGGGGCTGGCAATCCGGGCCGCAAGCGTg                                                                                                                                   >  5:188704/1‑110 (MQ=255)
                                                          tCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCGATGATACGt                                                                                          >  5:132687/1‑148 (MQ=255)
                                                            aaCACCGCACCAGGCATGACAGGTCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCGATGATACGTca                                                                                        <  1:102190/148‑1 (MQ=255)
                                                            aaCACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCGATGATACGTcag                                                                                       >  5:49170/1‑149 (MQ=255)
                                                               acCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCGATGATACGTcagcag                                                                                    >  7:236973/1‑149 (MQ=255)
                                                                     ccAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCGATGATACGTCAGCAGCAAccc                                                                              <  4:204/149‑1 (MQ=255)
                                                                             gACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCGATGATACGTCAGCAGCAACCCGTTACCgg                                                                      >  2:48285/1‑149 (MQ=255)
                                                                             gACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCGATGATACGTCAGCAGCAACCCGTTACCgg                                                                      >  2:43030/1‑149 (MQ=255)
                                                                               cAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCGATGATACGTCAGCAGCAACCCGTTACCGGcc                                                                    >  1:109269/1‑149 (MQ=255)
                                                                                aGATCATAGCCTTGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCTGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTTATAACCATGCAAGGCGCGATGATACGTCAGCAGGAACCCGTTACCGGCCg                                                                   <  5:332352/149‑1 (MQ=255)
                                                                                aGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCGATGATACGTCAGCAGCAAccc                                                                              >  3:53615/1‑138 (MQ=255)
                                                                                aGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCGATGATACGTCAGCAGCAAccc                                                                              <  4:53615/138‑1 (MQ=255)
                                                                                aGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCGATGATACGTCAGCAGCAACCCGTTACCGGCCg                                                                   <  3:20350/149‑1 (MQ=255)
                                                                                  aTCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCGATGATACGTCAGCAGCAACCCGTTACCGGCCg                                                                   >  8:131477/1‑147 (MQ=255)
                                                                                   tCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTTATAACCATGCAAGGCGCGATGATACGTCAGCAGCAACCCGTTACCGGCCGTAg                                                                >  1:293899/1‑149 (MQ=255)
                                                                                          cTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCGATGATACGTCAGCAGCAACCCGTTACCGGCCGTAGCAAGCCg                                                         <  8:42959/149‑1 (MQ=255)
                                                                                             gTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCGATGATACGTCAGCAGCAACCCGTTACCGGCCGTAGCAAGCCGGt                                                       <  2:136906/148‑1 (MQ=255)
                                                                                              tGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCGATGATACGTCAGCAGCAACCCGTTACCGGCCGTAGCAAGCCGGTGc                                                     <  8:236973/149‑1 (MQ=255)
                                                                                                 ccATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGAGACGCGAGGTTAATAACCATGCAAGGCGCGATGATACGTCAGCAGCAACCCGTTACCGGCCGTAGCAAGCCGGTGCcgc                                                  >  1:40189/1‑149 (MQ=255)
                                                                                                   aTGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCGATGATACGTCAGCAGCAACCCGTTACCGGCCGTAGCAAGCCGGTGCCGCgt                                                <  2:227904/149‑1 (MQ=255)
                                                                                                   aTGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCAAGCAAGGCGCGATGATACGTCAGCAGCAACCCGTTACCGGCCGTAGCAAGCCGGTGCCGCgt                                                >  8:169328/1‑149 (MQ=255)
                                                                                                    tGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCGATGATACGTCAGCAGCAACCCGTTACCGGCCGTAGCAAGCCGGTGCCGCgtg                                               <  5:289538/149‑1 (MQ=255)
                                                                                                      gcggcCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCGATGATACGTCAGCAGCAACCCGTTACCGGCCGTAGCAAGCCGGTGCCGCgtgg                                              <  6:49170/148‑1 (MQ=255)
                                                                                                        ggcCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCGATGATACGTCAGCAGCAACCCGTTACCGGCCGTAGCAAGCCGGTGCCGCGTGGTGc                                           <  4:59957/149‑1 (MQ=255)
                                                                                                                cGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCGATGATACGTCAGCAGCAACCCGTTACCGGCCGTAGCAAGCCGGTGCCGCGTGGTGCCAGCCG‑Gc                                  <  8:233501/149‑1 (MQ=255)
                                                                                                                   cTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCGATGATACGTCAGCAGCAACCCGTTACCGGCCGTAGCAAGCCGGTGCCGCGTGGTGCCAGCCG‑GCtgg                               <  3:318793/149‑1 (MQ=255)
                                                                                                                     ggCCAGGATTGCCAGCACCAGTTGCTGGAAGTCCAGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCGATGATACGTCAGCAGCAACCCGTTACCGGCCGTAg                                                                >  7:37293/1‑115 (MQ=25)
                                                                                                                     ggCCAGGATTGCCAGCACCAGTTGCTGGAAGTCCAGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCGATGATACGTCAGCAGCAACCCGTTACCGGCCGTAg                                                                <  8:37293/115‑1 (MQ=25)
                                                                                                                     ggCCAGGATTGCCAGCACCAGTTGCTGGAAGTCCAGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCGATGATACGTCAGCAGCAACCCGTTACCGGCCGTAGCAAGCCGGTGCCGCGTGGTGCCAGCCG‑GCtggt                              >  1:205562/1‑148 (MQ=38)
                                                                                                                     ggCCAGGATTGCCAGCACCAGTTGCTGGAAGTCCAGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCGATGATACGTCAGCAGCAACCCGTTACCGGCCGTAGCAAGCCGGTGCCGCGTGGTGCCAGCCG‑GCtggt                              <  2:205562/148‑1 (MQ=38)
                                                                                                                            attgccAGCACCAGTTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCGATGATACGTCAGCAGTAACCCGTTACCGGCCGTAGCAAGCCGGTGCCGCGTGGTGCCAGCCG‑GCTGGTGGCcgac                       >  6:242015/5‑148 (MQ=255)
                                                                                                                                ccAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCGATGATACGTCAGCAGCAACCCGTTACCGGCCGTAGCAAGCCGGTGCCGCgtggt                                             >  6:299847/1‑123 (MQ=255)
                                                                                                                                ccAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCGATGATACGTCAGCAGCAACCCGTTACCGGCCGTAGCAAGCCGGTGCCGCgtggt                                             <  5:299847/123‑1 (MQ=255)
                                                                                                                                 cAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCGATGATACGTCAGCAGCAACCCGTTACCGGCCGTAGCAAGCCGGTGCCGCGTGGTGCCAGc                                       <  7:209254/128‑1 (MQ=255)
                                                                                                                                 cAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCGATGATACGTCAGCAGCAACCCGTTACCGGCCGTAGCAAGCCGGTGCCGCGTGGTGCCAGc                                       >  8:209254/1‑128 (MQ=255)
                                                                                                                                           tGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCAAGCAAGGCGCGATGATACGTCAGCAGCAACCCGTTACCGGCCGTAGCAAGCCGGTGCCGCGTGGTGCCAGCCG‑GCTGGTGGCCGACGAGCCCGTATCGgcgc       <  7:169328/149‑1 (MQ=255)
                                                                                                                                            gCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCGATGATACGTCAGCAGCAACCCGTTACCGGCCGTAGCAAGCCGGTGCCGCGTGGTGCCAGCCG‑GCT‑GTGGCCGACGAGCCCGTATCGgcgcgc     >  6:325909/1‑149 (MQ=255)
                                                                                                                                                gCAATCCTGGCCGAAAGCGTGGCGATGCGAGGTTAATAACCATGCAAGGCGCGATGATACGTCAGCAGCAACCCGTTACCGGCCGTAGCAAGCCGGTGCCGCGTGGTGCCAGCCG‑GCTGGTGGCCGAGGAGCCCGTATCGGCGCGCCTg  >  8:121224/1‑149 (MQ=255)
                                                                                                                                                gCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCGATGATACGTCAGCAGCAACCCGTTACCGGCCGTAGCAAGCCGGTGCCGCGTGGTGCCAGCCGCGCTGGTGGCCGACGAGCCCGTATCGGCGCGCCt   <  1:43030/149‑1 (MQ=255)
                                                                                                                                                gCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCGATGATACGTCAGCAGCAACCCGTTACCGGCCGTAGCAAGCCGGTGCCGCGTGGTGCCAGCCGCGCTGGTGGCCGACGAGCCCGTATCGGCGCGCCt   <  1:48285/149‑1 (MQ=255)
                                                                                                                                                gCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCGATGATACGTCAGCAGCAACCCGTTACCGGCCGTAGCAAGCCGGTGCCGCGTGGTGCCAGCCG‑GCTGGTGGCCGACGAGCCCGTATCGGCGCGCCTg  >  1:185677/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
GGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCGATGATACGTCAGCAGCAACCCGTTACCGGCCGTAGCAAGCCGGTGCCGCGTGGTGCCAGCCG‑GCTGGTGGCCGACGAGCCCGTATCGGCGCGCCTG  >  NC_002947/1526591‑1526883

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