Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 399,993 C→A 11.0% V311L (GTG→TTG)  ddlA ← D‑alanine‑D‑alanine ligase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913399,9930CA11.0% 183.5 / 2.4 61V311L (GTG→TTG) ddlAD‑alanine‑D‑alanine ligase A
Reads supporting (aligned to +/- strand):  ref base C (21/32);  new base A (6/2);  total (27/34)
Fisher's exact test for biased strand distribution p-value = 1.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGGTGTAACCCAGACCGCTGGCTTGCCACAGCTTCGGATACATACTGATATTAGTGAAGCCAGGCAGGGTGTTGATCTCGTTGATCACCACTTCGTTCTCTGGGGTTAAAAACACGTCTACACGCGCCATGCCTGCGCATCCCAACGTTTGATAAGCCTGAACGGCAATCGCCCGGATCTTATCGTTGATTTCTGGCGCAATGGCTGCCGGAACTACCACTTTCGCGCCATCTTCG  >  NC_000913/399846‑400138
                                                                                                                                                   |                                                                                                                                                 
aTGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGGTGTAACCCAGACCGCTGGCTTGCCACAGCTTCGGATACATACTGATATTAGTGAAGCCAGGCAGGGTGTTGATCTCGTTGATCACCACt                                                                                                                                                  >  7:13450/1‑149 (MQ=255)
   cggtataCGCAGCGTGACGCTCCAGCGGCAGTTCAATCAAACGTGTGATGAGATCGGTGTAACCCATACCGCTGGCTTGCAACAGCTTCGGATACATACTGATATTAGTGAAGCCAGGCAGGGTGTTGATCTCGTTGATAACCACTTCg                                                                                                                                               <  6:323793/146‑1 (MQ=255)
   cGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGGTGTAACCCAGACCGCTGGCTTGCCACAGCTTCGGATACATACTGATATTAGTGAAGCCAGGCAGGGTGTTGATCTCGTTGATCACCACTTc                                                                                                                                                <  5:180375/148‑1 (MQ=255)
     ttATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGGTGTAACCCAGACCGCTGGCTTGCCACAGCTTCGGATACATACTGATATTAGTGAAGCCAGGCAGGGTGTTGATCTCGTTGATCACCACTTCGtt                                                                                                                                             <  2:150044/149‑1 (MQ=255)
          cGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGGTGTAACCCAGACCGCTCGCTTGCCACAGCTTCGGATACATACTGATATTAGTGAAGCCAGGCAGGGTGTTGATCTCGTTGATCACCACTTCGTTctct                                                                                                                                         <  5:139240/148‑1 (MQ=255)
           gCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGGTGTAACCCAGACCGCTGGCTTGCCACAGCTTCGGATACATACTGATATTAGTGAAGCCAGGCAGGGTGTTGATCTCGTTGATCACCACTTCGTTCTCTgg                                                                                                                                       <  5:352998/149‑1 (MQ=255)
           gCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGGTGTAACCCAGACCGCTGGCTTGCCACAGCTTCGGATACATACTGATATTAGTGAAGCCAGGCAGGGTGTTGATCTCGTTGATCACCACTTCGTTCTCTgg                                                                                                                                       <  5:74626/149‑1 (MQ=255)
                       tCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGGTGTAACCCAGACCGCTGGCTTGCCACAGCTTCGGATACATACTGATATTAGTGAAGCCAGGCAGGGTGTTGATCTCGTTGATCACCACTTCGTTCTCTGGGGTTAAAAacac                                                                                                                           <  4:344662/149‑1 (MQ=255)
                         cAGCGCCAGTTCAATCAAACGTGTGATCAGATCGGTGTAACCCAGACCGCTGGCTTGCCACAGCTTCGGATACATACTGATATTAGTGAAGCCAGGCAGGGTGTTGATCTCGTTGATCACCACTTCGTTCTCTGGGGTTAAAAacac                                                                                                                           <  8:157678/147‑1 (MQ=255)
                               cAGTTCAATCAAACGTGTGATCAGATCGGTGTAACCCAGACCGCTGGCTTGCCACAGCTTCGGATACATACTGATATTAGTGAAGCCAGGCAGGGTGTTGATCTCGTTGATCACCACTTCGTTCGCTGGGGTTAAAAACCCGTCTacac                                                                                                                   >  7:167061/1‑149 (MQ=255)
                                aGTTCAATCAAACGTGTGATCAGATCGGTGTAACCCAGACCGCTGGCTTGCCACAGCTTCGGATACATACTGATATTAGTGAAGCCAGGCAGGGTGTTGATCTCGTTGATCACCACTTCGTTATATCTGGTTAAAAACACGTCTAAAcg                                                                                                                  >  2:64013/1‑149 (MQ=255)
                                 gTTCAATCAAACGTGTGATCAGATCGGTGTAACCCAGACCGCTGGCTTGCCACAGCTTCGGATACATACTGATATTAGTGAAGCCAGGCAGGGTGTTGATCTCGTTGATCACCACTTCGTTCTCTGGGGTTAAAAACACGTCTACAcgc                                                                                                                 >  3:62913/1‑149 (MQ=255)
                                 gTTCAATCAAACGTGTGATCAGATCGGTGTAACCCAGACCGCTGGCTTGCCACAGCTTCGGATACATACTGATATTAGTGAAGCCAGGCAGGGTGTTGATCTCGTTGATCACCACTTCGTTCTCTGGGGTTAAAAACACGTCTACAcgc                                                                                                                 >  5:50400/1‑149 (MQ=255)
                                 gTTCAATCAAACGTGTGATCAGATCGGTGTAACCCAGACCGCTGGCTTGCCACAGCTTCGGATACATACTGATATTAGTGAAGCCAGGCAGGGTGTTGATCTCGTTGATCACCACATCCTTATCTGATGTTAAAAACACGTCTACAcgc                                                                                                                 >  2:372158/1‑149 (MQ=255)
                                       tcaaACGTGTGATCAGATCGGTGTAACCCAGACCGCTGGCTTGCCACAGCTTCGGATACATACTGATATTAGTGAAGCCAGGCAGGGTGTTGATCTCGTTGATCACCACTTCGTTCTCTGGGGTTAAAAACACGTCTACACGCGCCATg                                                                                                           >  7:94558/1‑149 (MQ=255)
                                         aaaCGTGTGATCAGATCGGTGTAACCCAGACCGCTGGCTTGCCACAGCTTCGGATACATACTGATATTAGTGAAGCCAGGCAGGGTGTTGATCTCGTTGATCACCACTTCGTTCTCTGGGGTTAAAAACACGTCTACACGCGCCATGcc                                                                                                         >  7:338881/1‑149 (MQ=255)
                                            cGTGTGATCAGATCGGTGTAACCCAGACCGCTGGCTTGCCACAGCTTCGGATACATACTGATATTAGTGAAGCCAGGCAGGGTGTTGATCTCGTTGATCACCACTTCGTTCTCTGGGGTTAAAAACACGTCTACACGCGCCATGCCTgc                                                                                                      <  8:161086/149‑1 (MQ=255)
                                             gtgtGATCAGATCGGTGTAACCCAGACCGCTGGCTTGCCACAGCTTCGGATACATACTGATATTAGTGAAGCCAGGCAGGGTGTTGATCTCGTTGATCACCACTTCGTTCTCTGGGGTTAAAAACACGTCTACACGCGCCATGCCTgcg                                                                                                     >  7:272912/1‑149 (MQ=255)
                                                       aTCGGTGTAACCCAGACCGCTGGCTTGCCACAGCTTCGGATACATACTGATATTAGTGAAGCCAGGCAGGGTGTTGATCTCGTTGATCACCACTTCGTTCTCTGGGGTTAAAAACACGTCTACACGCGCCATGCCTGCGCATCCCAACg                                                                                           >  3:241399/1‑149 (MQ=255)
                                                            tgtAACCCAGACCGCTGGCTTGCCACAGCTTCGGATACATACTGATATTAGTGAAGCCAGGCAGGGTGTTGATCTCGTTGATCACCACTTCGTTCTCTGGGGTTAAAAACACGTCTACACGCGCCATGCCTGCGCATCCCAACGTTTGa                                                                                      <  8:13450/149‑1 (MQ=255)
                                                            tgtAACCCAGACCGCTGGCTTGCCACAGCTTCGGATACATACTGATATTAGTGAAGCCAGGCAGGGTGTTGATCTCGTTGATCACCACTTCGTTCTCTGGGGTTAAAAACACGTCTACACGCGCCATGCCTGCGCATCCCAACGTTTGa                                                                                      <  7:98404/149‑1 (MQ=255)
                                                            tgtAACCCAGACCGCTGGCTTGCCACAGCTTCGGATACATACTGATATTAGTGAAGCCAGGCAGGGTGTTGATCTCGTTGATCACCACTTCGTTCTCTGGGGTTAAAAACACGTCTACACGCGCCATGCCTGCGCATCCCAACGTTTGa                                                                                      <  1:372158/149‑1 (MQ=255)
                                                            tgtAACCCAGACCGCTGGCTTGCCACAGCTTCGGATACATACTGATATTAGTGAAGCCAGGCAGGGTGTTGATCTCGTTGATCACCACTTCGTTCTCTGGGGTTAAAAACACGTCTACACGCGCCATGCCTGCGCATCCCAACGTTTGa                                                                                      <  1:142902/149‑1 (MQ=255)
                                                            tgtAACCCAGACCGCTGGCTTGCCACAGCTTCGGATACATACTGATATTAGTGAAGCCAGGCAGGGTGTAGATCTCGTTGATCACCAAATCATTAAAAAAAAAAAAAAACACGTCTACACGCGCCATGCCTGCGCATCCCAACGTTTGa                                                                                      >  4:68159/1‑149 (MQ=255)
                                                                aCCTAGACCGCTGGCTTGCCACAGCTTAGGATACTTACTGATATTAGTGAAGCCAGGCAGGGTGTTGATCTCGTTGATCACCACTTCGTTCTCTGGGGGTAAAAACACGTATACACGCGCCATGCCTGCGCATCCCAACGTTTGATaa                                                                                   <  5:317628/148‑1 (MQ=255)
                                                                  ccAGACCGCTGGCTTGCCACAGCTTCGGATACATACTGATATTAGTGAAGCCAGGCAGGGTGTTGATCTCGTTGATCACCACTTCGTTCTCTGGGGTTAAAAACACGTCTACACGCGCCATGCCTGCGCATCCCAACGTTTGATAAGcc                                                                                >  2:330120/1‑149 (MQ=255)
                                                                     gACCGCTGGCTTGCCACAGCTTCGGATACATACTGATATTAGTGAAGCCAGGCAGGGTGTTGATCTCGTTGATCACCACTTCGTTCTCTGGGGTTAAAAACACGTCTACACGCGCCATGCCTGCGCATCCCAACGTTTGATAAGCCTGa                                                                             <  6:50400/149‑1 (MQ=255)
                                                                                 gCCACAGCTTCGGATACATACTGATATTAGTGAAGCCAGGCAGGGTGTTGATCTCGTTGATCACCACTTCGTTCTCTGGGGTTAAAAACACGTCTACACGCGCCATGCCTGCGCATCCCAACGTTTGATAAGCCTGAACGGCAATCGc                                                                  <  4:124577/148‑1 (MQ=255)
                                                                                   cacaGCTTCGGATACATACTGATATTAGTGAAGCCAGGCAGGGTGTAGATCTCGTTGATCACCAAATCATTAAAAAAAAATAAAAACACGTCAACACGCGCAAGGCCTGCGCATCCCAACGATTGATAAGCCTGAACGGCAATCGccca                                                               >  2:291188/1‑148 (MQ=255)
                                                                                       gCTTCGGATACATACTGATATTAGTGAAGCCAGGCAGGGTGTTGATCTCGTTGATCACCACTTCGTTCTCTGGGGTTAAAAACACGTCTACACGCGCCATGCCTGCGCATCCCAACGTTTGATAAGCCTGAACGGCAATCGCCCGGATc                                                           <  4:241399/149‑1 (MQ=255)
                                                                                       gCTTCGGATACATACTGATATTAGTGAAGCCAGGCAGGGTGTTGATCTCGTTGATCACCACTTCGTTCTCTGGGGTTAAAAACACGTCTACACGCGCCATGCCTGCGCATCCCAACGTTTGATAAGCCTGAACGGCAATCGCCCGGATc                                                           <  8:21223/149‑1 (MQ=255)
                                                                                         ttCGGATACATACTTATATTAGTGAAGCCAGGCAGGGTGTTGAGCTCGTTGATAACCAATTCGTTCTCTGGGGTTAAAAAAACGTCTACACGCGCCATGCCTGCGCATACCAACGTTTGATAAGCCTGAACGGCAATCGCACGGATCtt                                                         <  3:328701/149‑1 (MQ=255)
                                                                                           cGGATACATACTGATATTAGTGAAGCCAGGCAGGGTGTTGATCTCGTTGATCACCACTTCGTTCTCTGGGGTTAAAAACACGTCTACACGCGCCATGCCTGCGCATCCCAACGTTTGATAAGCCTGAACGGCAATCGCCCGGATCtt                                                         <  1:330120/147‑1 (MQ=255)
                                                                                            ggATACATACTGATATTAGTGAAGCCAGGCAGGGTGTTGATCTCGTTGATCACCACTTCGTTCTCTGGGGTTAAAAACACGTCTACACGCGCCATGCCTGCGCATCCCAACGTTTGATAAGCCTGAACGGCAATCGCCCGGATCTTAt                                                       <  2:162066/148‑1 (MQ=255)
                                                                                              atacatacTGATATTAGTGAAGCCAGGCAGGGTGTAGATCTCTTTGATCACCAAATCATTAAAAAAAAAAAAAAACAAGTCTACACGCGCCATGCCTGCGCATCCCAACGTTTGATAAGCCTCCACGGCAATCGCCCGGATCTTCTCGt                                                    >  6:218158/1‑149 (MQ=255)
                                                                                              atacatacTGATATTAGTGAAGCCAGGAAGTGTGTTGATCGCGTTGATCACCACTTCGTTCTCTGGGGTTAAAAACACGTCTACACGCGCCATGCCTGCGCATCCCAACGTTTGATAAGCCTGAACGGCAATCGCCCGGATCTTATCGt                                                    <  1:64013/149‑1 (MQ=255)
                                                                                                 catacTGATATTAGTGAAGCCAGGCAGGGTGTTGATCTCGTTGATCACCACTTCGTTCTCTGGGATTAAAAACACGTCTACACGCGCCATGCCTGCGCATCCCAACGATTGATAAGCCTGAACGGCAATCGCCCGGATCTTATCGttta                                                 >  7:360718/1‑147 (MQ=255)
                                                                                                  atacTGATATTAGTGAAGCCAGGCAGGGTGTTGATCTCGTTGATCACCACTTCGTTCTCTGGGGTTAAAAACACGTCTACACGCGCCATGc                                                                                                          <  1:216309/91‑1 (MQ=255)
                                                                                                  atacTGATATTAGTGAAGCCAGGCAGGGTGTAGATCTCGTTGATCACCAAATCATTAAAAAAAAAAAAAAACACGTCTACACGCGCCATGc                                                                                                          >  2:216309/1‑91 (MQ=255)
                                                                                                   tacTGATCTTCGTGAAGCCAGGCAGGGTGTAGATCTCGTTGATCACCAAATCATTAAAAAAAAAAAAAAACACGTCTACACGCGCCATGCAAGCGCATACCAACGTTTGATAAGCCTGAACGGCAATCTCCCGGATCTGATCGTTGAtt                                               >  6:226039/1‑149 (MQ=255)
                                                                                                    acTGATATTAGTGAAGCCAGGCAGGGTGTTGATCTCGTTGATCACCACTTCGTTCTCTGGGGTTAAAAACACGTCTACACGCGCCATGCCTGCGCATCCCAACGTTTGATAAGCCTGAACGGCAATCGCCCGGATCTTATCGTTGAttt                                              >  1:202853/1‑149 (MQ=255)
                                                                                                      tGATATTAGTGAAGCCAGGCAGGGTGTAGATCTCGTTGATCACCAAATCATTAAAAAAAAATAAAAACACGTCTACACGCGCCATGCCTGCGCATCCCAACGTTTGATAAGCCTGAACGGCAATCGCCCGGATCTTATCGTTGATTTCt                                            >  8:352054/1‑149 (MQ=255)
                                                                                                        atatTAGTGAAGCCAGGCAGGGTGTTGATCTCGTTGATCACCACTTCGTTCTCTGGGGTTAAAAACACGTCTACACGCGCCATGCCTGCGCATCCCAACGTTTGATAAGCCTGAACGGCAATCGCCCGGATCTTATCGTTGATTTCTg                                           <  3:313598/148‑1 (MQ=255)
                                                                                                        atatTAGTGAAGCCAGGCAGGGTGTTGATCTCGTTGATCACCACTTCGTTCTCTGGGGTTAAAAACACGTCTACACGCGCCATGCCTGCGCATCCCAACGTTTGATAAGCCTGAACGGCAATCGCCCGGATCTTATCGTTGATTTCTg                                           <  4:62913/148‑1 (MQ=255)
                                                                                                         tatTAGTGAAGCCAGGCAGGGTGTTGATCTCGTTGATCACCACTTCGTTCTCTGGGGTTAAAAACACGTCTACACGCGCCATGCCTGCGCATCCCAACGTTTGATAAGCCTGAACGGCAATCGCCCGGATCTTATCGTTGATTTCTGgc                                         <  7:122876/149‑1 (MQ=255)
                                                                                                         tatTAGTGAAGCCAGGCAGGGTGTTGATCTCGTTGATCACCACTTCGTTCTCTGGGGTTAAAAACACGTCTACACGCGCCATGCCTGCGCATCCCAACGTTTGATAAGCCTGAACGGCAATCGCCCGGATCTTATCGTTGATTTCTGgc                                         >  7:159808/1‑149 (MQ=255)
                                                                                                              gTGAAGCCAGGCAGGGTGTTGATCTCGTTGATCACCACTTCGTTCTCTGGGGTTAAAAACACGTCTACACGCGCCATGCCTGCGCATCCCAACGTTTGATAAGCCTGAACGGCAATCGCCCGGATCTTATCGTTGATTTCTGGCGCAAt                                    <  8:272912/149‑1 (MQ=255)
                                                                                                               tGAAGCCAGGCAGGGTGTTGATCTCGTTGATCACCACTTCGTTCTCTGGGGTTAAAAACACGTCTACACGCGCCATGCCTGCGCATCCCAACGTTTGATAAGCCTGAACGGCAATCGCCCGGATCTTATCGTTGATTTCTGGCGCAATg                                   <  8:94558/149‑1 (MQ=255)
                                                                                                                             gtgtTGATCTCGTTGATCACCACTTCGTTCTCTGGGGTTAAAAACACGTCTACACGCGCCATGCCTGCGCATCCCAACGTTTGATAAGCCTGAACGGCAATCGCCCGGATCTTATCGTTGATTTCTGGCGCAATGGCTGCCGGAACTAc                     <  8:338881/149‑1 (MQ=255)
                                                                                                                                ttGATCTCGTTGATCACCACTTCGTTCTCTGGGGTTAAAAACACGTCTACACGCGCCATGCCTGCGCATCCCAACGTTTGATAAGCCTGAACGGCAATCGCCCGGATCTTATCGTTGATTTCTGGCGCAATGGCTGCCGGAACTACCa                   <  3:68159/148‑1 (MQ=255)
                                                                                                                                ttGATCTCGTTGATCACCACTTCGTTCTCTGGGGTTAAAAACACGTCTACACGCGCCATGCCTGCGCATCCCAACGTTTGATAAGCCTGAACGGCAATCGCCCGGATCTTATCGTTGATTTCTGGCGCAATGGCTGCCGGAACTACCAc                  >  4:36136/1‑149 (MQ=255)
                                                                                                                                 cgaactcGTTGATCACCACTTCGTTCTCTGGGGTTAAAAACACGTCTACACGCGCCATGCCTGCGCATCCCAACGTTTGATAAGCCTGAACGGCAATCGCCCGGATCTTATCGTTGATTTCTGGCGCAATGGCTGCCGGAACTACCAc                  <  7:270928/144‑1 (MQ=255)
                                                                                                                                   aTCTCGTTGATCACCACTTCGTTCTCTGGGGTTAAAAACACGTCTACACGCGCCATGCCTGCGCATCCCAACGTTTGATAAGCCTGAACGGCAATCGCCCGGATCTTa                                                        <  2:173847/108‑1 (MQ=255)
                                                                                                                                   aTCTCGTTGATCACCACTTCGTTCTCTGGGGTTAAAAACACGTCTACACGCGCCATGCCTGCGCATCCCAACGTTTGATAAGCCTGAACGGCAATCGCCCGGATCTTa                                                        >  1:173847/1‑108 (MQ=255)
                                                                                                                                   aTCGCGTTTATCAACACTTCGTTCTCTGGGGTTAAAAACACGTCTACACGCGCCATGCCTGCGCATCCCAACGTTTGATAAGCCTGAACGGCAATCGCCCGGATCTTATCGTTGATTTCTGGCGCAATGGCTGCCGGAACTACCACttt               <  1:102659/149‑1 (MQ=255)
                                                                                                                                        gTTGATCACCACTTCGTTCTCTGGGGTTAAAAACACGTCTACACGCGCCATGCCTGCGCATCCCAACGTTTGATAAGCCTGAACGGCAATCGCCCGGATCTTATCGTTGATTTCTGGCGCAATGGCTGCCGGAACTACCAc                  >  8:270928/1‑141 (MQ=255)
                                                                                                                                         ttGATCACCACTTCGTTCTCTGGGGTTAAAAACACGTCTACACGCGCCATGCCTGCGCATCCCAACGTTTGATAAGCCTGAACGGCAATCGCCCGGATCTTATCGTTGATTTCTGGCGCAATGGCTGCCGGAACTACCACTTTcgcg           >  7:272060/1‑147 (MQ=255)
                                                                                                                                         ttGATCACCAATTCGGTCTCTTGGGTTAAAAACACGTCTACACGCGCCATGCCTGCGCATCCCAACGTTTGATAAGCCTGAACGGCAATCGCCCGGATCTTATCGTTGATTTCTGGCGCAATGGCTGCCGGAACTACCACTTTcgcg           <  8:272060/147‑1 (MQ=255)
                                                                                                                                           gATCACCACTTCGTTCTCTGGGGTTAAAAACACGTCTACACGCGCCATGCCTGCGCATCCCAACGTTTGATAAGCCTGAACGGCAATCGCCCGGATCTTATCGTTGATTTCTGGCGCAATGGCTGCCGGAACTACCACTTTCGCGCCAt       >  1:332558/1‑149 (MQ=255)
                                                                                                                                             tCACCACTTCGTTCTCTGGGGTTAAAAACACGTCTACACGCGCCATGCCTGCGCATCCCAACGTTTGATAAGCCTGAACGGCAATCGCCCGGATCTTATCGTTGATTTCTGGCGCAATGGCTGCCGGAACTACCACTTTCGCGCCATCt     >  7:10579/1‑149 (MQ=255)
                                                                                                                                             tCACCACTTCGTTCTCTGGGGTTAAAAACACGTCTACACGCGCCATGCCTGCGCATCCCAACGTTTGATAAGCCTGAACGGCAATCGCCCGGATCTTATCGTTGATTTCTGGCGCAATGGCTGCCGGAACTACCACTTTCGCGCCATCt     >  6:311905/1‑149 (MQ=255)
                                                                                                                                                ccacTTCGTTCTCTGGGGTTAAAAACACGTCTACACGCGCCATGCCTGCGCATCCCAACGTTTGATAAGCCTGAACGGCAATCGCCCGGATCTTATCGTTGATTTCTGGCGCAATGGCTGCCGGAACTACCACTTTCGCGCCATCTTCg  >  1:170443/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
ATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGGTGTAACCCAGACCGCTGGCTTGCCACAGCTTCGGATACATACTGATATTAGTGAAGCCAGGCAGGGTGTTGATCTCGTTGATCACCACTTCGTTCTCTGGGGTTAAAAACACGTCTACACGCGCCATGCCTGCGCATCCCAACGTTTGATAAGCCTGAACGGCAATCGCCCGGATCTTATCGTTGATTTCTGGCGCAATGGCTGCCGGAACTACCACTTTCGCGCCATCTTCG  >  NC_000913/399846‑400138

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: