Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 855,411 Δ1 bp 6.9% coding (334/921 nt) ldtB ← L,D‑transpeptidase linking Lpp to murein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913855,4110A. 6.9% 179.6 / 6.8 58coding (334/921 nt)ldtBL,D‑transpeptidase linking Lpp to murein
Reads supporting (aligned to +/- strand):  ref base A (26/27);  new base . (2/2);  total (28/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.07e-01

TGCGTGCATTTTGGCGGTCGGCGTCCAGGTCGGGCCTGCTTTTTTACGCTCAACTTTGGTGGTCCAGTTGATAGGCGTATCTTTGCCTAACTGACCAATGCCGATCGGCAGCACGATAACGGTGTTGGTCCCTTTCGGATAGTAATAAAGACGCATCTCAGCACTGTTAATGACGATGCCTTCATGAACGGTATCCGGCAGGATCAGCTGCTGCGGAATGTTCAGTACAGTACCGCCTTTCGGCAGGAAGGTATCCACACCCGGGTTCGCTTCCATCATATTGGAAAGCC  >  NC_000913/855268‑855557
                                                                                                                                               |                                                                                                                                                  
tGCGTGCATTTTGGCGGTCGGCGTCCAGGTCGGGCCTGCTTTTTTACGCTCAACTTTGGTGGTCCAGTTGATAGGCGTATCTTTGCCTATCTGACCAATGCCGATCGGCAGCACGATAACGGTGTTGGTCCCTTTCGGATAGtaata                                                                                                                                                 >  5:85366/1‑147 (MQ=255)
  cGTGCATTTTGGCGGTCGGCGTCCAGGTCGGGCCTGCTTTATTACGCTCAACTTTGGTGGTCCAGTTGATAGGCGTATCTTTGCCTAACTGACCAATGCCGATCGGCAGCACGATAACGGTGTTGGTCCCTTTCGGATAGTAATaaa                                                                                                                                               <  8:164771/147‑1 (MQ=255)
   gTGCATTTTGGCGGTCGGCGTCCAGGTCGGGCCTGCTTTTTTACGCTCAACTTTGGTGGTCCAGTTGATAGGCGTATCTTTGCCTAACTGACCAATGCCGATCGGCAGCACGATAACGGTGTTGGTCCCTTTCGGATAGTAATAAAGAc                                                                                                                                            <  5:182264/149‑1 (MQ=255)
   gTGCAGTTTGGCGGTCGGCGTCCAGGTCGGGCCTGCTTTTTTAAGCTCAACATAGGAGGTCCAGTTGATAGCCGTATCTTTGCCTAACTGACCAATCCAGATCGGCAGCACGATAACCGTGTTGCTCCCTAACTGATAGTaagaaagag                                                                                                                                            >  8:62004/1‑148 (MQ=255)
        ttttGGCGGTCGGCGTCCAGGTCGGGCCGGCTGTTTTACGCGCAACTTTGGTGGTCCAGTTGATAGGCGTATCTTTGCCTAACTGACCAATGCCGATCGGCAGCACGATAACGGTGTTGGTCCCTTTCGGATAGTAATAAAGACGCAtc                                                                                                                                       <  8:315028/149‑1 (MQ=255)
         tttGGCGGTCGGCGTCCAGGTCGGGCCTGCTTTTTTACGCTCAACTTTGGTGGTCCAGTTGATAGGCGTATCTTTGCCTAACTGACCAATGCCGATCGGCAGCACGATAACGGTGTTGGTCCCTTTCGGATAGTAATAAAGACGCAtct                                                                                                                                      >  7:229605/1‑149 (MQ=255)
                      gTCCAGGTCGGGCCTGCTTTTTTACGCTCAACTTTGGTGGTCCAGTTGATAGGCGTATCTTTGCCTAACTGACCAATGCCGATCGGCAGCACGATAACGGTGTTGGTCCCTTTCGGATAGTAATAAAGACGAATCTCAGCACTGTTAAt                                                                                                                         >  5:155787/1‑149 (MQ=255)
                        ccAGGTCGGGCCTGCTTTTTTACGCTCAACTTTGGTGGTCCAGTTGATAGGCGTATCTTTGCCTAACTGACCAATGCCGATCGGCAGCACGATAACGGTGTTGGTCCCTTTCGGATAGTAATAAAGACGCATCTCAGCACTGTTAATGa                                                                                                                       >  3:235874/1‑149 (MQ=255)
                          aGGTCGGGCCTGCTTTTTTACGCTCAACTTTGGTGGTCCAGTTGATAGGCGTATCTTTGCCTAACTGACCAATGCCGATCGGCAGCACGATAACGGTGTTGGTCCCTTTCGGATAGTAATAAAGACGCATCTCAGCACTGTTAATGACg                                                                                                                     <  4:11123/149‑1 (MQ=255)
                           ggTCGGGCCTGCTTTTTTACGCTCAACTTTGGTGGTCCAGTTGATAGGCGTATCTTTGCCTAACTGACCAATGCCGATCGGCAGCACGATAACGGTGTTGGTCCCTTTCGGATAGTAATAAAGACGCATCTCAGCACTGTTAATGACGa                                                                                                                    <  1:167383/149‑1 (MQ=255)
                           ggTCGGGCCTGCTTTTTTACGCTCAACTTTGGTGGTCCAGTTGATAGGCGTATCTTTGCCTAACTGACCAATGCCGATCGGCAGCACGATAACGGTGTTGGTCCCTTTCGGATAGTAATAAAGACGCATCTCAGCACTGTTAATGACGa                                                                                                                    <  7:262521/149‑1 (MQ=255)
                                   cTGCTTTTTTACGCTCAACTTTGGTGGTCCAGTTGATAGGCGTATCTTTGCCTAACTGACCAATGCCGATCGGCAGCACGATAACGGTGTTGGTCCCTTTCGGATAGTAATAAAGACGCATCTCAGCACTGTTAATGACGATGCCTTCa                                                                                                            >  1:206262/1‑149 (MQ=255)
                                   cTGCTTTTTTACGCTCAACTTTGGTGGTCCAGTTGATAGGCGTATCTTTGCCTAACTGACCAATGCCGATCGGCAGCACGATAACGGTGTTGGTCCCTTTCGGATAGTATAAAGACGCATCTCAGCACTGTTAATGACGATGCCTTCAt                                                                                                           >  5:228533/1‑149 (MQ=255)
                                   cTGCTTTTTTACGCTCAACTTTGGTGGTCCAGTTGATAGGCGTATCTTTGCCTAACTGACCAATGCCGATCGGCAGCACGATAACGGTGTTGGTCCCTTTCGGATAGTATAAAGACGCATCTCAGCACTGTTAATGACGATGCCTTCAt                                                                                                           >  5:233268/1‑149 (MQ=255)
                                             aCGCTCAACTTTGGTGGTCCAGTTGATAGGCGTATCTTTGCCTAACTGACCAATGCCGATCGGCAGCACGATAACGGTGTTGGTCCCTTTCGGATAGTAATAAAGACGCATCTCAGCACTGTTAATGACGATGCCTTCAt                                                                                                           >  4:192101/1‑140 (MQ=255)
                                             aCGCTCAACTTTGGTGGTCCAGTTGATAGGCGTATCTTTGCCTAACTGACCAATGCCGATCGGCAGCACGATAACGGTGTTGGTCCCTTTCGGATAGTAATAAAGACGCATCTCAGCACTGTTAATGACGATGCCTTCAt                                                                                                           <  3:192101/140‑1 (MQ=255)
                                                   aaCTTTGGTGGTCCAGTTGATAGGCGTATCTTTGCCTAACTGACCAATGCCGATCGGCAGCACGATAACGGTGTTGGTCCCTTTCGGATAGTAATAAAGACGCATCTCAGCACTGTTAATGACGATGCCTTCATGAACGGTATCCgg                                                                                              <  2:206262/147‑1 (MQ=255)
                                                   aaCTTTGGTGGTCCAGTTGATAGGCGTATCTTTGCCTAACTGACCAATGCCGATCGGCAGCACGATAACGGTGTTGGTCCCTTTCGGATAGTAATAAAGACGCATCTCAGCACTGTTAATGACGATGCCTTCATGAACGGTATCCGGc                                                                                             <  5:177597/148‑1 (MQ=255)
                                                   aaCTTTGGTGGTCCAGTTGATAGGCGTATCTTTGCCTAACTGACCAATGCCGATCGGCAGCACGATAACGGTGTTGGTCCCTTTCGGATAGTAATAAAGACGCATCTCAGCACTGTTAATGACGATGCCTTCATGAACGGTATCCGGCa                                                                                            <  4:44765/149‑1 (MQ=255)
                                                   aaCTTTGGTGGTCCAGTTGATAGGCGTATCTTTGCCTAACTGACCAATGCCGATCGGCAGCACGATAACGGTGTTGGTCCCTTTCGGATAGTAATAAAGACGCATCTCAGCACTGTTAATGACGATGCCTTCATGAACGGTATCCGGCa                                                                                            <  6:132235/149‑1 (MQ=255)
                                                     cTTTGGTGGTCCAGTTGATAGGCGTATCTTTGCCTAACTGACCAATGCCGATCGGCAGCACGATAACGGTGTTGGTCCCTTTCGGATAGTAATAAAGACGCATCTCAGCACTGTTAATGACGATGCCTTCATGAACGGTATc                                                                                                 <  3:138542/142‑1 (MQ=255)
                                                     cTTTGGTGGTCCAGTTGATAGGCGTATCTTTGCCTAACTGACCAATGCCGATCGGCAGCACGATAACGGTGTTGGTCCCTTTCGGATAGTAATAAAGACGCATCTCAGCACTGTTAATGACGATGCCTTCATGAACGGTATc                                                                                                 >  4:138542/1‑142 (MQ=255)
                                                        tggtggTCCAGTTGATAGGCGTATCTTTGCCTAACTGACCAATGCCGATCGGCAGCACGATAACGGTGTTGGTCCCTTTCGGATAGTATAAAGACGCATCTCAGCACTGTTAATGACGATGCCTTCATGAACGGTATCCGGCAGGAt                                                                                        <  6:228533/147‑1 (MQ=255)
                                                        tggtggTCCAGTTGATAGGCGTATCTTTGCCTAACTGACCAATGCCGATCGGCAGCACGATAACGGTGTTGGTCCCTTTCGGATAGTATAAAGACGCATCTCAGCACTGTTAATGACGATGCCTTCATGAACGGTATCCGGCAGGAt                                                                                        <  6:233268/147‑1 (MQ=255)
                                                         ggaggtCCAGTTGATAGGCGTATCTTTGCCTAACTGACCAATGCCGATCGGCAGCACGATAACGGTGTTGGTCCCTTTCGGATAGTAATAAAGACGCATCTCAGCACTGTTAATGACGATGCCTTCATGAACGGTATCCGGCAGGATc                                                                                       <  2:120047/145‑1 (MQ=255)
                                                             gtCCAGTTGATAGGCGTATCTTTGCCTAACTGACCAATGCCGATCGGCAGCACGATAACGGTGTTGGTCCCTTTCGGATAGTAATAAAGACGCATCTCAGCACTGTTAATGACGATGCCTTCATGAACGGTATCCGGCAGGATCAgct                                                                                   <  2:175344/148‑1 (MQ=255)
                                                              tCCAGTTGATAGGCGTATCTTTGCCTAACTGACCAATGCCGATCGGCAGCACGATAACGGTGTTGGTCCCTTTCGGATAGTAATAAAGACGCATCTCAGCACTGTTAATGACGATGCCTTCATGAACGGTATCCGGCAGGATCAgctgc                                                                                 >  3:293259/1‑149 (MQ=255)
                                                               ccAGTTGATAGGCGTATCTTTGCCTAACTGACCAATGCCGATCGGCAGCACGATAACGGTGTTGGTCCCTTTCGGATAGTAATAAAGACGCATCTCAGCACTGTTAATGACGATGCCTTCATGAACGGTATCCGGCAGGATCAgctgct                                                                                >  5:28922/1‑149 (MQ=255)
                                                                cAGTTGATAGGCGTATCTTTGCCTAACTGACCAATGCCGATCGGCAGCACGATAACGGTGTTGGTCCCTTTCGGATAGTAATAAAGACGCATCTCAGCACTGTTAATGACGATGCCTTCATGAACGGTATCCGGCAGGATCAGctgctg                                                                               >  5:148050/1‑149 (MQ=255)
                                                                cAGTTGAAAGGCCTATCATTGCCTAACTGACCAATGCCGCTCGGCAGCACGATAACGGTCTAGGTCCCTTTCGCATAGTAATAAAGACGCATCTCAGAACTGTTAATGACGATTCCTTCATGAACCCTATCCGGCAGGAACAGcagctg                                                                               >  8:87787/1‑149 (MQ=255)
                                                                         ggCGTATCTTTGCCTAACTGACCAATGCCGATCGGCAGCACGATAACGGTGTTGGTCCCTTTCGGATAGTAATAAAGACGCATCTCAGCACTGTTAATGACGATGCCTTCATGAACGGTATCCGGCAGGATCAGCTGCTGCGGAATGtt                                                                      <  5:67427/149‑1 (MQ=255)
                                                                          gCGTATCTTTGCCTATCTGACCAATGCCGATCGGCAGCACGATAACGGTGTTGGTCCCTTTCGGATAGTAATAAAGACGCATCTCAGCACTGTTAATGACGATGCCTTCATGAACGGTATCCGGCAGGATCAGCTGCTGCGGAATGTTc                                                                     <  6:85366/149‑1 (MQ=255)
                                                                            gTATCTTTGCCTAACTGACCAATGCCGATCGGCAGCACGATAACGGTGGTGGTCCCTTTCGGATAGTAATAAAGACGCATCTCAGCACTGTTAATGACGATGCCTTCATGAACGGTATCCGGCAGGATCAGCTGCTGCGGAATGTTcag                                                                   <  1:89037/149‑1 (MQ=255)
                                                                              aTCTTTGCCTAACTGACCAATGCCGATCGGCAGCACGATAACGGTGTTGGTCCCTTTCGGATAGTAATAAAGACGCATCTCAGCACTGTTAATGACGATGCCTTCATGAACGGTATCCGGCAGGATCAGCTGCTGCGGAATGTTcagta                                                                 >  4:215700/1‑149 (MQ=255)
                                                                               tCTTTTCCTAACTGACCAATGCCGATCGGCAGCACGATAACGGTGTTGGTCCCTTTCGGATAGTAATAAAGACGCATCTCAGCACTGTTAATGACGATGCCTTCATGAACGGTATCAGGCAGGATCAGCTGCTGCGGAATGTTcagta                                                                 <  8:115051/148‑1 (MQ=255)
                                                                                cTTTGCCTAACTGACCAATGCCGATCGGCAGCACGATAACGGTGTTGGTCCCTTTCGGATAGTAATAAAGACGCATCTCAGCACTGTTAATGACGATGCCTTCATGAACGGTATCCGGCAGGATCAGCTGCTGCGGAATGTTcagta                                                                 >  8:29634/1‑147 (MQ=255)
                                                                                     ccTAACTGACCAATGCCGATCGGCAGCACGATAACGGTGTTGGTCCCTTTCGGATAGTAATAAAGACGCATCTCAGCACTGTTAATGACGATGCCTTCATGAACGGTATCCGGCAGGATCAGCTGCTGCGGAATGTTCAGTACAGTAcc                                                          >  7:194702/1‑149 (MQ=255)
                                                                                     ccTAACTGACCAATGCCGATCGGCAGCACGATAACGGTGTTGGTCCCTTTCGGATAGTAATAAAGACGCATCTCAGCACTGTTAATGACGATGCCTTCATGAACGGTATCCGGCAGGATCAGCTGCTGCGGAATGTTCAGTACAGTAcc                                                          <  6:148050/149‑1 (MQ=255)
                                                                                        aaCTGACCAATGCCGATCGGCAGCACGATAACGGTGTTGGTCCCTTTCGGATAGTAATAAAGACGCATCTCAGCACTGTTAATGACGATGCGTTCATGAACGGTATCCGGCAGGATCAGCTGCTGCGGAATGTTCAGTACAGTACCGcc                                                       >  2:186123/1‑149 (MQ=255)
                                                                                             tacAATGGCGATCGGCTGCACGATAACTGTGTTGGTGCCTTTCTGATAGTTATAAAGACGCATCTCAGCACTGTTAATGAAGATTCCTTCATGAACGGTATCCGTCAGGATGAGCTGCTGCGGAATGTTCAGTACAGTACCGCCTTTc                                                   <  5:316664/146‑1 (MQ=255)
                                                                                             aCCAATGCCGATCGGCAGCACGATAACGGTGTTGGTCCCTTTCGGATAGTAATAAAGACGCATCTCAGCACTGTTAATGACGATGCCTTCATGAACGGTATCCGGCAGGATCAGCTGCTGCGGAATGTTCAGTACAGTACCGCCTTTCg                                                  <  7:29634/149‑1 (MQ=255)
                                                                                               cAATGCCGATCGGCAGCACGATAACGGTGTTGGTCCCTTTCGGATAGTAATAAAGACGCATCTCAGCACTGTTAATGACGATGCCTTCATGAACGGTATCCGGCAGGATCAGCTGCTGCGGAATGTTCAGTACAGTACCGCCTTTCgg                                                 >  1:254032/1‑148 (MQ=255)
                                                                                               cAATGCCGATCGGCAGCACGATAACGGTGTTGGTCCCTTTCGGATAGTAATAAAGACGCATCTCAGCACTGTTAATGACGATGCCTTCATGAACGGTATCCGGCAGGATCAGCTGCTGCGGAATGTTCAGTACAGTACCGCCTTTCgg                                                 >  5:324474/1‑148 (MQ=255)
                                                                                                     cGATCGGCAGCACGATAACGGTGTTGGTCCCTTTCGGATAGTAATAAAGACGCATCTCAGCACTGTTAATGACGATGCCTTCATGAACGGTATCCGGCAGGATCAGCTGCTGCGGAATGTTCAGTACAGTACCGCCTTTCGGCAGGAAg                                          >  1:46142/1‑149 (MQ=255)
                                                                                                     cGATCGGCAGCACGATAACGGTGTTGGTCCCTTTAGGATAGTAATAAAGACGCATCTCAGCACTGTTAATGACGATGCCTTCATGAACGGTATCCGGCAGGATCAGCTGCTGCGGAATGTTCAGTACAGTACCGCCTTTCGGCAGGAAg                                          >  5:103066/1‑149 (MQ=255)
                                                                                                       aTCGGCAGCACGATAACGGTGTTGGTCCCTTTCGGATAGTAATAAAGACGCATCTCAGCACTGTTAATGACGATGCCTTCATGAACGGTATCCGGCAGGATCAGCTGCTGCGGAATGTTCAGTACAGTACCACCTTTCGGCAGGAAGGt                                        <  4:293259/149‑1 (MQ=255)
                                                                                                          ggCAGCACGATAACGGTGTTGGTCCCTTTCGGATAGTAATAAAGACGCATCTCAGCACTGTTAATGACGATGCCTTCATGAACGGTATCCGGCAGGATCAGCTGCTGCGGAATGTTCAGTACAGTACCGCCTTTCGGCAGGAAGGTATc                                     <  3:87646/149‑1 (MQ=255)
                                                                                                             agcaCGATAACGGTGTTGGTCCCTTTCGGATAGTAATAAAGACGCATCTCAGCACTGTTAATGACGATGCCTTCATGAACGGTATCCGGCAGGATCAGCTGCTGCGGAATGTTCAGTACAGTACCGCCTTTCGGCAGGAAGGTATCca                                   >  3:264678/1‑148 (MQ=255)
                                                                                                              gcaCGATAACGGTGTTGGTCCCTTTCGGATAGTAATAAAGACGCATCTCAGCACTGTTAATGACGATGCCTTCATGAACGGTATCCGGCAGGATCAGCTGCTGCGGAATGTTCAGTACAGTACCGCCTTTCGGCAGGAAGGTATCcaca                                 <  3:215700/149‑1 (MQ=255)
                                                                                                                  gATAACGGTGTTGGTCCCTTTCGGATAGTAATAAAGACGCATCTCAGCACTGTTAATGACGATGCCTTCATGAACGGTATCCGGCAGGATCAGCTGCTGCGGAATGTTCAGTACAGTACCGCCTTTCGGCAGGAAGGTATCCACACCCg                             >  5:206965/1‑149 (MQ=255)
                                                                                                                         gtgtTGGTCCCTTTCGGATAGTAATAAAGACGCATCTCAGCACTGTTAATGACGATGCCTTCATGAACGGTATCCGGCAGGATCAGCTGCTGCGGAATGTTCAGTACAGTACCGCCTTTCGGCAGGAAGGTATCCACACCCGGGTTCGc                      >  3:32911/1‑149 (MQ=255)
                                                                                                                          tgtTGGTCCCTTTCGGATAGTAATAAAGACGCATCTCAGCACTGTTAATGACGATGCCTTCATGAACGGTATCCGGCAGGATCAGCTGCTGCGGAATGTTCAGTACAGTACCGCCTTTCGGCAGGAAGGTATCCACACCCGGGTTCGCt                     <  4:32911/149‑1 (MQ=255)
                                                                                                                            ttGGTCCCTTTCGGATAGTAATAAAGACGCATCTCAGCACTGTTAATGACGATGCCTTCATGAACGGTTTCCGGCAGGATCAGCTGCTGCGGAATGTTCAGTACAGTACCGCCTTTCGGCAGGAAGGTATCCACACCCGGGTTCGCTTc                   <  8:229605/149‑1 (MQ=255)
                                                                                                                                tCCCTTTCGGATAGTAATAAAGACGCATCTCAGCACTGTTAATGACGATGCCTTCATGAACGGTATCCGGCAGGATCAGCTGCTGCGGAATGTTCAGTACAGTACCGCCTTTCGGCAGGAAGGTATCCACACCCGGGTTCGCTTCcat                >  2:338172/1‑148 (MQ=255)
                                                                                                                                tCCCTTTCGGATAGTAATAAAGACGCATCTCAGCACTGTTAATGACGATGCCTTCATGAACGGTATCCGGCAGGATCAGCTGCTGCGGAATGTTAAGTACAGTACCGCCTTTCGGCAGGAAGGTATCCACACCCGGGTTCGCTTCcat                >  1:102773/1‑148 (MQ=255)
                                                                                                                                  ccTTTCGGATAGTAATAAAGACGCATCTCAGCACTGTTAATGACGATGCCTTCATGAACGGTATCCGGCAGGATCAGCTGCTGCGGAATGTTCAGTACAGTACCGCCTTTAGGCAGGAAGGTATCCACACCCGGTTTCGCTTCcatca              >  3:9187/1‑148 (MQ=255)
                                                                                                                                   cTTTCGGATAGTAATAAAGACGCATCTCAGCACTGTTAATGACGATGCCTTCATGAACGGTATCCGGCAGGATCAGCTGCTGCGGAATGTTCAGTACAGTACCGCCTTTCGGCAGGAAGGTATCCACACCCGGGTTCGCTTCCATCata            <  4:264678/149‑1 (MQ=255)
                                                                                                                                    tttCGGATAGTAATAAAGACGCATCTAAGCACTGTTAATGACGATGCCTTCATGAACGGTATCCGGCAGGATCAGCTGCTGCGGAATGTTCAGTACAGTACCGCCTTTCTGAAGGAAGGTATCCACACCCGGGTTCGCTTCCATCatat           >  3:268061/1‑149 (MQ=255)
                                                                                                                                          agaGTAATAAAGACGCCTCTCAGCACTGTTAATGACGATGCCTTCATGCACGGTATCCGGCAGGATCAGCTGCTGCGGAATGTTCAGTACAGTACCGCGTTTCGGCAGGAAGCGATCCCCACCCGGGTTCGCTTCCATCATATTGGa       >  5:182124/3‑147 (MQ=255)
                                                                                                                                             gTAATAAAGACGCATCTCAGCACTGTTAATGACGATGCCTTCATGAACGGTATCCGGCAGGATCAGCTGCTGCGGAATGTTCAGTACAGTACCGCCTTTCGGCAGGAAGGTATCCACACCCGGGTTCGCTTCCATCATATTGGAAAGcc  <  1:338172/149‑1 (MQ=255)
                                                                                                                                              taataaAGACGCATCTCAGCACTGTTAATGACGATGCCTTCATGAACGGTATCCGGCAGGATCAGCTGCTGCGGAATGTTCAGTACAGTACCGCCTTTCGGCAGGAAGGTATc                                     >  1:91605/1‑113 (MQ=255)
                                                                                                                                              taataaAGACGCATCTCAGCACTGTTAATGACGATGCCTTCATGAACGGTATCCGGCAGGATCAGCTGCTGCGGAATGTTCAGTACAGTACCGCCTTTCGGCAGGAAGGTATc                                     <  2:91605/113‑1 (MQ=255)
                                                                                                                                               aatacAGACGCATCTCAGCACTGTTAATGACGATGCCTTCATGAACGGTATCCGGCAGGATCAGCTGCTGCGGAATGTTCAGTACAGTACCGCCTTTCGGCAGGAAGGTATCCACACCCGGGTTCGCTTCCATCATATTg         <  1:45902/140‑1 (MQ=255)
                                                                                                                                               aataaAGACGCATCTCAGCACTGTTAATGACGATGCCTTCATGAACGGTATCCGGCAGGATCAGCTGCTGCGGAATGTTCAGTACAGTACCGCCTTTCGGCAGGAAGGTATCCACACCCGGGTTCGCTTCCATCATATTg         >  2:45902/1‑140 (MQ=255)
                                                                                                                                               |                                                                                                                                                  
TGCGTGCATTTTGGCGGTCGGCGTCCAGGTCGGGCCTGCTTTTTTACGCTCAACTTTGGTGGTCCAGTTGATAGGCGTATCTTTGCCTAACTGACCAATGCCGATCGGCAGCACGATAACGGTGTTGGTCCCTTTCGGATAGTAATAAAGACGCATCTCAGCACTGTTAATGACGATGCCTTCATGAACGGTATCCGGCAGGATCAGCTGCTGCGGAATGTTCAGTACAGTACCGCCTTTCGGCAGGAAGGTATCCACACCCGGGTTCGCTTCCATCATATTGGAAAGCC  >  NC_000913/855268‑855557

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