Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,526,738 A→G 6.6% I309V (ATC→GTC)  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,7380AG 6.6% 172.3 / 2.3 60I309V (ATC→GTC) ddlBD‑alanine‑‑D‑alanine ligase B
Reads supporting (aligned to +/- strand):  ref base A (29/27);  new base G (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 = 6.27e-01

GGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCGATGATACGTCAGCAGCAACCCGTTACCGGCCGTAGCAAGCCGGTGCCGCGTGGTGCCAGCCG‑GCTGGTGGCCGACGAGCCCGTATCGGCGCGCCTGCCC  >  NC_002947/1526591‑1526886
                                                                                                                                                   |                                                                                                                                                     
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                                                                   <  8:37293/115‑1 (MQ=25)
                                                                                                                     ggCCAGGATTGCCAGCACCAGTTGCTGGAAGTCCAGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCGATGATACGTCAGCAGCAACCCGTTACCGGCCGTAg                                                                   >  7:37293/1‑115 (MQ=25)
                                                                                                                     ggCCAGGATTGCCAGCACCAGTTGCTGGAAGTCCAGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCGATGATACGTCAGCAGCAACCCGTTACCGGCCGTAGCAAGCCGGTGCCGCGTGGTGCCAGCCG‑GCtggt                                 <  2:205562/148‑1 (MQ=38)
                                                                                                                     ggCCAGGATTGCCAGCACCAGTTGCTGGAAGTCCAGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCGATGATACGTCAGCAGCAACCCGTTACCGGCCGTAGCAAGCCGGTGCCGCGTGGTGCCAGCCG‑GCtggt                                 >  1:205562/1‑148 (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)
                                                                                                                                                   aTCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCGATGATACGTCAGCAGCAACCCGTTACCGGCCGTAGCAAGCCGGTGCCGCGTGGTGCCAGCCG‑GCTGGTGGCCGACGAGCCCGTATCGGCGCGCCTGccc  <  2:293899/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCGATGATACGTCAGCAGCAACCCGTTACCGGCCGTAGCAAGCCGGTGCCGCGTGGTGCCAGCCG‑GCTGGTGGCCGACGAGCCCGTATCGGCGCGCCTGCCC  >  NC_002947/1526591‑1526886

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