Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 11,481 G→A S102S (AGC→AGT yaaI ← DUF2541 domain‑containing protein YaaI

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091311,4810GA92.9% 113.2 / ‑1.5 44S102S (AGC→AGTyaaIDUF2541 domain‑containing protein YaaI
Reads supporting (aligned to +/- strand):  ref base G (2/1);  new base A (19/20);  total (23/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

AGAGGCGAGTTTTGCCTTTTTCATTATGGGTGAGCAATCGGGCGAAATTTGCCAACTGTTCCTCACTACAATGCTGAAGAAAATCCAGATCTGAATCATTCAGGTAATTAACATTCATTTTTTGTGGCTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTATGACCCGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCGCTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTCTGAAGGAATATTCAGGCTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGATCGTTGGTTTCGTTACGCGGCATTGCAATGGCGCCGAGGAGTTTATGGTCGTTTGCCTGCG  >  NC_000913/11239‑11725
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                    
atagtagAGTTTTGCCTTATTCATTATGGTTGAGCAATCTGGCGAAATTTGCCATCTGTTCTTCACTACAATTCTGAAGAAAAGCCAGATCTGAATCATTCAGGTAATTAACATACATTTTTTGTGGTTTCTATATTCTGGCGTTAGTCGTCGCCTATAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTATGTCCGGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTTTCACTGTTAAt                                                                                                                                                                                                                                              <  2:429904/245‑1 (MQ=255)
  agGCGAGTTTTGCCTTTTTCATTATGGGTGAGCAATCGGGCGAAATTTGCCAACTGTTCCTCACTACAATGCTGCAGAAAATCCAGATCTGAATCATTCAGGTACTTAACATTCATTTTTTGTGGTTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTGGCCATATCCGATGCGTTCACCGTATGACCGGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCGCTGTTACTGt                                                                                                                                                                                                                                            >  1:499484/1‑251 (MQ=255)
      gaattttGCCTTTTTCATTACGGGTGAGCAATCGGGCGAAATTTGCCAACTGTTACTCACTACAATGCTGAAGAAAAGCCAGATCTGAATCATTCAGGTAATTAACAGTCATTTTTGGTGGTTTCTATATGCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTATGACCGGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCACTGTTAATGTTGAt                                                                                                                                                                                                                                        <  1:5061/248‑1 (MQ=255)
        gTTTTGCCTTTTTCATGATGGGTGAGCAATCGGTCGAAATTTGCCAACTGTTCCTCACTACAATTCTGAAGAAAATCCAGATCTGAATCATTCAGGTAATTAACATTCATTTTTTGTGGTTTCTATAGTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTATGACCGGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCACTGTTAATGTTGATc                                                                                                                                                                                                                                       <  2:638642/250‑1 (MQ=255)
                 ttttCATTATGGGTGTGCAATCGGGCGAAATTTGCCAACTGTTCCTCACTACAATGCTTAAGAAAATCCAGATCTGAATCATTCAGGTAATTAACATTCATTTTTTGTGGTTTCTATATTCTGTCGTTAGTCGTCGCCGATAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTATTACCGGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCACTGTTAATGTTGATCcagtcagt                                                                                                                                                                                                                               <  2:38436/249‑1 (MQ=255)
                   ttCATTATGGGTGAGCAATCGGGCGAAATATGCCACCTGTTCCTCCCTACAAAGCTGAAGAATATCCAGATCTGAATCATTCAGGTAATTAACATTCATTTTTTGAGGTTTCTATATTCTGGCATTAGTCTTCGCCGATCATTTTAAGCGTGTCAATATCCGCTGAGTTCACCGTATGCACGGTAAAGGTGATTTTTGATACGCAGCGTTTATTGACATTATATCTGTTAATTTTTATCCAGTAAGTgatt                                                                                                                                                                                                                           >  2:473471/1‑248 (MQ=255)
                     cATTATGGGCGAGCAATCGGGCGAAATTTGCCAACTGTGCCTCACTACAATGCTTAAGAAAATCCAGAGCTGAATTATTCAGGTAATTAACATGCATTTTTTGTGGTTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTATGACCGGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCACTGTTAATGTTGATCCAGTCAGTGGttt                                                                                                                                                                                                                          <  2:381524/250‑1 (MQ=255)
                     cATTATGGGAGAGCAATCAGGCGAAATATGCCAACTGTTCCTCACTCCAATGCTTAAGAAAACCCAGTTCTGAAACATTCAGGTAATTAACATTCATTTTTTGTGGGTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTGGCCATATCCTACGAGTTCACCGTATGACCGGAAAAGGTGATTTTCTAGACGCAGCGTTTATTTCCGTTATCACTGTTAATTTTGATCCAGTCAGTGGttt                                                                                                                                                                                                                          <  2:348685/250‑1 (MQ=255)
                               gAGCAATCGGGCGAAATTTGCCAACTGTTCCTCACTACAATGCTGAAGAAAATCCAGATCTGAATCATTCAGGTAATTAACATTCATTTTTTGTGGTTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTATGACCGGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCACTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCtttt                                                                                                                                                                                                               >  1:104789/1‑251 (MQ=255)
                                    aTCGGGCGAAATTAGCCAACTGTTCCTCACTACAATGCTGAAGAAAATCCAGATCTGACTCATTCAGGTAATTAACATTCATTTTTTGTGATTTCTATATTCTGGAATTAGTCGTCGCCGATAATTTTCAGCGTGGCCATATCCGTTGAGTTCACCGTATGACCGGAAAAGGTGATTTTTTAGACGAAGCGTTTCTTGTCGTTATCACTGTTAATGTTTATCCAGTCAGTGGTTTGCCCTTCTTTTAtttt                                                                                                                                                                                                          >  2:537259/1‑250 (MQ=255)
                                       gggCGAAATTTGCCAACTGTTCCTCACTACAATGCTGAAGAAAATCCAGATCTGAATCATTCAGGTAATTAACATTCATTTTTTGTGGTTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTATGACCGGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCACTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTc                                                                                                                                                                                                          <  1:179599/248‑1 (MQ=255)
                                       gggCGAAATTTGCCAACTGTTCCTCACTACAATGCTGAAGAAAATCCAGATCTGAATCATTCAGGTAATTAACATTCATTTTTTGTGGTTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTATGACCGGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCACTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTc                                                                                                                                                                                                          >  2:179599/1‑248 (MQ=255)
                                       gggCGAAATTTGCCAACTGTTCCTCACTACAATGCTGAAGAAAATCCAGATCTGAATCATTCAGGTAATTAACATTCATTTTTTGTGGTTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTATGACCGGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCACTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTCTg                                                                                                                                                                                                        >  1:600325/1‑250 (MQ=255)
                                        ggCGAAATTTGCCAACTGTTCCTCACTACAATGCTGAAGAAAATCCAGATCTGAATCATTCAGGTAATTAACATTCATTTTTTGTGGTTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTATGACCGGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCACTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTCTGaa                                                                                                                                                                                                      <  1:195586/251‑1 (MQ=255)
                                                        tGTTCCTCACTACACTGCTGCCGAAACTCCAGATCTGAATCATTCAGGTCATTAACATTCATTTTTTGTGGTTTCTATATTCTTGCGATAGTCGTCGCCGATAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTCTGACCGGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCACTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTCTGTAGGATTATTCAGActct                                                                                                                                                                                      >  2:94184/1‑251 (MQ=255)
                                                          ttCCTCACTACAATGCTGAAGAAAATCCAGATCTGAATCATTCAGGTAATTAACATTCATTTTTTGTGGTTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTATGACCGGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCACTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTCTGAAGGAATATTCaaactctg                                                                                                                                                                                     >  2:264661/1‑243 (MQ=255)
                                                             cTCACTACAATGCTGAAGAAAATCCAGATCTGAATCATTCAGGTAATTAACATTCATTTTTTGTGGTTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTATGACCGGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCACTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTCTGAAGGAATATTCAGACTCTGACTg                                                                                                                                                                                 >  1:676304/1‑251 (MQ=255)
                                                                     aaTGCTGAAGAAAATCCAGATCTGAATCATTCAGGTAATTAACATTCATTTTTTGTGGTTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTATGACCGGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCACTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTCTGAAGGAATATTCAGACTCTGACTGGCGCTACg                                                                                                                                                                         <  1:388952/251‑1 (MQ=255)
                                                                                     cAGATCGGAATCATTCAAGTAATTAACATTCAGTTTTTGTGGTTTCTATATTCTGGCGTTAGGCGTCGCCGAGAATTTTCAGCGTGTCCATATCCGATGAGTTCACCGTATGACCGGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCACTGTTAATGTTGATCCAGTTAGTGGTTTGCTCTGCTTTTATTTCTGAAGGAATATTCAGACTCTGACTGGCGTTACGGGCGGCTTTGAAATaa                                                                                                                                                         <  2:618039/251‑1 (MQ=255)
                                                                                                   tCAGGTAATTAACATTCATTTTTTGTGGTTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTATGACCGGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCACTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTCTGAAGGAATATTCAGACTCTGACTGGCGCTACGGGCGGCTTTGAATTAAACCGATgacccgct                                                                                                                                           >  1:435430/1‑244 (MQ=255)
                                                                                                        tAATTAACATTAATTTTGTGTGGTTTCTATATTCTGGCGTTAGTCGTAGCCGATAATTTTCAGCGTGTCCATATCCGATGAGTTCACCGTATGAACGGAAAAGGTGTTTTTTGAGACGCAGCGTTTATTGTCGTTATCACTGTTAATGATGCTCCAGTCAGTGGTTTTACCTTCTTTTATTTCTGATGGAATATACAGACTCTGACTGGCGCTACGGGCGGCTTTGACATACACCGATTCACCGTTTAACt                                                                                                                                      >  2:152350/1‑251 (MQ=255)
                                                                                                           ttAACATTCATTTTTTGTGGCTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTATGACCCGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCGCTGTTAATGtt                                                                                                                                                                                                                                           <  1:121627/147‑1 (MQ=255)
                                                                                                               cattcattTTTTGTGGTTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTGGCCAAATCCGATGAGTTCACCGTATGACCGGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCACTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTCTGAAGGAATATTCAGACTCTGACTGGCGCTACGGGAGGCTTTGAAATAAACCGATGCACCGCTTAACTGTAAATc                                                                                                                               <  1:433074/251‑1 (MQ=255)
                                                                                                                atacaatctttGTGGTTGCTATATTATGGCGTTAGTCGACGCAGATATTTTTCATAGTGGCCATATTCGATGAGTTCACCGTATGCCAGGAATAGGTGAGTTTTGAGACGCAGCGCTTATTGTCGCTATCACTGTGAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTAGTTCTGCAGGAATATTCAGCCTCTGATTGGCGCTACGAGCGGCTTTGTAATAAACCGATGCACCGATTAACTGTatctcg                                                                                                                              >  2:282005/9‑246 (MQ=255)
                                                                                                                   cattTTTTGTGGTTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTATGACCGGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCACTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTCTGAAGGAATATTCAGACTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGTAAATCGcc                                                                                                                            >  1:317993/1‑250 (MQ=255)
                                                                                                                      tttttGTGGTTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTATGACCGGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCACTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTCTGAAGGAATATTCAGACTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGTAAATCGCCATg                                                                                                                         <  1:59572/250‑1 (MQ=255)
                                                                                                                          tgtgGTTTCTATATTCTGGAGTTAGTCGTCGCCGATAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTATGACCGGAAAAGGTGATTTTTGAGACGCAGCTTTTATTGTCGTTATCACTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTCTGAAGGAATATTCAGACTCTGACTGGCGCTACGGGCGGTTTTGAAATAAACCGATGCACCGCTTAACTGTAAATCGCCATGGTCg                                                                                                                     <  2:54947/250‑1 (MQ=255)
                                                                                                                                                tAGTCGTCGCCGATAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTATGACCGGAAAAGGTGATTTTTGAGACGCCGCGTTTATTGTCGTTATCCCTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTCTGACGGAATATTCAGACTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACAGATGCACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTcac                                                                                              >  1:416217/1‑251 (MQ=255)
                                                                                                                                                       cgCCGATAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTATGACCGGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCACTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTCTGAAGGAATATTCAGACTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCg                                                                                        >  1:48544/1‑250 (MQ=255)
                                                                                                                                                            aTAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTATGCCCGGTAAAGGTGATTTTTGAGACGCAGCGATTATTGTCGTTATCACTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTCTGAAGGAATATTCAGACTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGCCAAAc                                                                                  >  2:108446/1‑251 (MQ=255)
                                                                                                                                                            aTAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTATGACCGGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCACTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTCTGAAGGAATATTCAGACTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAAc                                                                                  <  1:108446/251‑1 (MQ=255)
                                                                                                                                                                       ggTGGCCATATCCGATGAGTTCAACGTATGACCGGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCACTGGTAATGTTGATCCAGTCAGTGGGTTGCCCTTCGTTTATTTCTGAAGGAATATTCAGACTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCa                                                                       <  2:676304/250‑1 (MQ=255)
                                                                                                                                                                          ggCCATATCCGATGAGTTCACCGTATGACCGGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCACTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTCTGAAGGAATATTCAGACTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAgcg                                                                    <  2:104789/251‑1 (MQ=255)
                                                                                                                                                                              atatCCGATGAGTTCACCGTATGACCGGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCACTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTCTGAAGGAATATTCAGACTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCTATGCACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTTTATGCGTTTCACAATGCGACCAACAGGAAGTTTCAGCGcc                                                                  >  2:436019/1‑249 (MQ=255)
                                                                                                                                                                                  ccGATGAGTTCACCGTATGACCGGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCACTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTCTGAAGGAATATTCAGACTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGATCg                                                            <  1:309258/251‑1 (MQ=255)
                                                                                                                                                                                               cGTATGACCGGCAATGGTTATTTTTGAACCTCAGCGTTTATTGTCGTTATCAATGTTAATGTTGATCCAGTCAGTGGTTTGCACTTCTTTTATTTATGAAGGAATATTCAGCCTCTGACTGGCGCTACGGGCGGCATTGCAATAAACCGATGCACCGCTTAACTGTAAATCGCCATGGTAGGCAGAGAGATGTGTGCGTTTCACAATGCGCCAAACAGGAGGTTTCAgcgaacgctggtgggtttcttta                                                >  2:86248/1‑230 (MQ=255)
                                                                                                                                                                                                        ggAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCACTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTCTGAAGGAATATTCAGACTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGATCGTTGGTTTCGTTACGCGGCATTg                                      <  1:185307/250‑1 (MQ=255)
                                                                                                                                                                                                                 gATTTTTGAGACGCAGCGTTTATTGTCGTTATCACTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTCTGAAGGAATATTCAGACTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGATCGTTGGTTTCGTTACGCGGCATTGCAATgacg                              <  1:364139/250‑4 (MQ=255)
                                                                                                                                                                                                                     tttGAGACGCAGCGTTTATTGTCGTTATCACTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTCTGAAGGAATATTCAGACTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGATCGTTGGTTTCGTTACGCGGCATTGCAATGACGCCg                           >  1:171420/1‑249 (MQ=255)
                                                                                                                                                                                                                              cAGCGTTTATTGTCGTTATCGCTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTCTGAAGGAATATTCAGGCTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCGGCTCGTTGGTTTCGTTACGCGGCATTGCAATGGCGCCGGGGGGTTTATg                >  2:89729/1‑251 (MQ=255)
                                                                                                                                                                                                                                  gTTTATTGTCGTTATCACTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTCTGAAGGAATATTCAGACTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGATCGTTGGTTTCGTTACGCGGCATTGCAATGACGCCGAGGAGTTTATGGTCg            <  1:436019/251‑1 (MQ=255)
                                                                                                                                                                                                                                    ttATTGTCGTTATCACTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTCTGAAGGAATATTCAGACTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGATCGTTGGTTTCGTTACGCGGCATTGCAATGACGCCGAGGAGTTTATGGACGtt          >  2:672777/1‑251 (MQ=255)
                                                                                                                                                                                                                                    ttATTGACGTTATCACTGTTAATGTTGATCCAGTCAGTGGTTTGCCGTTCTTTTATTTCTGATGGAATATTCAGACTCTGACTGGCGCTACGGGAGGCTTTGATATAAACCGATGCACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCACCGCCAGATCGTTGGTTTCGTTACGCGGCATTGCAAGGAAGCCGAGGAGTTTATGGTCGtt          >  2:553861/1‑251 (MQ=255)
                                                                                                                                                                                                                                            gTTATCACTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTCTGAAGGAATATTCAGACTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGATCGTTGGTTTCGTTACGCGGCATTGCAATGACGCCGAGGAGTTTATGGTCGTTTGCCTgcg  >  2:578647/1‑251 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                    
AGAGGCGAGTTTTGCCTTTTTCATTATGGGTGAGCAATCGGGCGAAATTTGCCAACTGTTCCTCACTACAATGCTGAAGAAAATCCAGATCTGAATCATTCAGGTAATTAACATTCATTTTTTGTGGCTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTATGACCCGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCGCTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTCTGAAGGAATATTCAGGCTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGATCGTTGGTTTCGTTACGCGGCATTGCAATGGCGCCGAGGAGTTTATGGTCGTTTGCCTGCG  >  NC_000913/11239‑11725

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