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,4810GA91.7% 121.7 / 0.7 49S102S (AGC→AGTyaaIDUF2541 domain‑containing protein YaaI
Reads supporting (aligned to +/- strand):  ref base G (2/2);  new base A (23/21);  total (25/24)
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.73e-01

CTGGAGAGGCGAGTTTTGCCTTTTTCATTATGGGTGAGCAATCGGGCGAAATTTGCCAACTGTTCCTCACTACAATGCTGAAGAAAATCCAGATCTGAATCATTCAGGTAATTAACATTCATTTTTTGTGGCTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTATGACCCGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCGCTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTCTGAAGGAATATTCAGGCTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGATCGTTGGTTTCGTTACGCGGCATTGCAATGGCGCCGAGGAGTTTATGGTCGTTTGCCTGC  >  NC_000913/11235‑11724
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                   
cTGGAGAGGCGAGTTTTGCCTTTTTCATTATGGGTGAGCAATCGGGCGAAATTTGCCAACTGTTCCTCACTACAATGCTGAAGAAAATCCAGATCTGAATCATTCAGGTAATTAACATTCATTTTTTGTGGTTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTATGACCGGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCACTGt                                                                                                                                                                                                                                                 >  1:305758/1‑251 (MQ=255)
           aGTTTTGCCTTTTTCATTATGGGTGAGCAATCGGGCGAAATTTGCCAACTGTTCCTCACTACAATGCTGAAGAAAATCCAGATCTGAATCATTCAGGTAATTAACATTCATTTTTTGTGGTTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTATGACCGGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCACTGTTAATGTTGATc                                                                                                                                                                                                                                      <  1:302540/251‑1 (MQ=255)
           aGTTTTGCCTTTTTCATTATGGGGGAGCAATCGGGCGAAATTTGCTAACTGTTCCTCACTACAATTCTGAAGAAAATCCAGATCTGAATCATTCAGGTAATTAACATTCATTTTTTGTGGTTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTATGACCGGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCACTGTTAATGTTGATc                                                                                                                                                                                                                                      <  2:406614/251‑1 (MQ=255)
              tttGCCTTTTTCATTATGGGTGAGCAATCGGGCGAAATTTGCCAACTGTTCCTCACTACAATGCTGAAGAAAATCCAGATCTGAATCATTCAGGTAATTAACATTCATTTTTTGTGGCTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTATGACCCGAAAAGGGGGTTTTTGAGACGCAGCGTTTATTTTCGTTATCGCTGTTAATGTTTATcc                                                                                                                                                                                                                                     >  2:235346/1‑249 (MQ=255)
                tGCCTTTTTCATTATGGGTGAGCAATCGGGCGAAATTTGCCAACTGTTCCTCACTACAATGCTGAAGAAAATCCAGATCTGAATCATTCAGGTAATTAACATTCATTTTTTGTGGTTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTATGACCGGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCACTGTTAATGTTGATCcagt                                                                                                                                                                                                                                  <  1:465603/250‑1 (MQ=255)
                 gCCTTTTTCATTATGGGTGAGCAATCGGGCGAAATTTGCCAACTGTTCCTCACTACAATGCTGAAGAAAATCCAGATCTGAATCATTCAGGTAATTAACATTCATTTTTTGTGGTTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTATGACCGGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCACTGTTAATGTTGATCcagtc                                                                                                                                                                                                                                 <  1:1825/250‑1 (MQ=255)
                    tttttCATTATGGGTGAGCAATCGGGCGAAATTTGCCAACTGTTCCTCACTACAATGCTGAAGAAAATCCAGATCTGAATCATTCAGGTAATTAACATTCATTTTTTGTGGTTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTATGACCGGAAACGGTGCTTTTTGAGACGCAGCGTTTCTTGTCGTTATCATTGTTATTGTTGTTCCGGTCAGTg                                                                                                                                                                                                                             >  2:118942/1‑251 (MQ=255)
                     ttttCATTATGGGTGAGCAATCGGGCGAAATTTGCCAACTGTTCCTCACTACAATGCTGAAGAAAATCCAGATCTGAATCATTCAGGTAATTAACATTCATTTTTTGTGGTTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTATGACCGGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCACTGTTAATGTTGATCcagtcagt                                                                                                                                                                                                                              <  1:72271/249‑1 (MQ=255)
                     ttttCATTATGGGTGAGAAATTGGGCTAAATTGTCCAACTTTGCCACACTACAAGGCTTAAGAAAATCCAGATCTGAATCATTCAGGTAATTAACATTCATTTTTTGTGGGTTCTATATTCGGGCGTTATTCGTCGCTGATAATTTTCAGCGTGGTCATATCCGATTAGTTCACCGTATGACCGGAAAAGGTGATTTTTGAGACGCAGCGGTTATGGGCGTTATCACTGTTAATGTTGACCCAGTCAGTg                                                                                                                                                                                                                             <  2:311120/250‑1 (MQ=255)
                        tCATTATGGGTGAGCAATCGGGCGAAATTTGCCAACTGTTCCTCACTACAATGCTGAAGAAAATCCAGATCTGAATCATTCAGGTAATTAACATTCATTTTTTGTGGTTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTATGACCGGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCTCTGTTAATGTTGATCCAGTCAGTGGtt                                                                                                                                                                                                                          <  1:200077/250‑1 (MQ=255)
                               tggggAGCAATCGGGCGAAATTTGCCAACTGTTCCTCACTACAATGCTGAATAAAATCCAGATCTGAATCATTCAGGTAATTAACATTCATTTTTTGTGGTTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTATGACCGGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCACTGTTAATGTTGATCCAGTCAGTGGTTTGCCctt                                                                                                                                                                                                                   <  2:220465/246‑1 (MQ=255)
                                        aTCGGGCGAAATTTGCCAACTGTTCCTCACTACAATGCTGAAGAAAATCCAGATCTGAATCATTCAGGTAATTAACATTCATTTTTTGTGGTTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTATGACCGGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCACTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTc                                                                                                                                                                                                         >  1:231937/1‑251 (MQ=255)
                                          cGGGCGAAATTTGCCAACTGTTCCTCACTACAATGTGGAAGAAAATCCAGATCTGAATCATTCAGGTAATTAACATTCATTTTTTGTGGTTTCTATATTCTGGCGTTAGTCGTCGCCGATAAGTTTCAGCTTTGCCATTTCCGATGAGTTCACCGTATTGCCGGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATTACTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTCt                                                                                                                                                                                                        <  1:207401/250‑1 (MQ=255)
                                                            tGTTCCTCACTACAATGCTGAAGAAAATCCAGATCTGAATCATTCAGGTAATTAACATTCATTTTTTGTGGTTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTATGACCGGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCACTg                                                                                                                                                                                                                                                  <  2:364850/190‑1 (MQ=255)
                                                            tGTTCCTCACTACAATGCTGAAGAAAATCCAGATCTGAATCATTCAGGTAATTAACATTCATTTTTTGTGGTTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTATGACCGGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCACTg                                                                                                                                                                                                                                                  >  1:364850/1‑190 (MQ=255)
                                                                             cTGAAGAAAATCCAGATCTGAATCATTCAGGTAATTAACATTCATTTTTTGTGGTTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTATGACCGGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCACTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTCTGAAGGAATATTCAGACTCTGACTGGCGCTACGggcg                                                                                                                                                                    >  1:238904/1‑251 (MQ=255)
                                                                                   aaaaTCCAGATCTGAATCATTCAGGTAATTAACATTCATTTTTTGTGGTTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTATGACCGGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCACTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTCTGAAGGAATATTCAGACTCTGACTGGCGCTACGGGCGGCTTTg                                                                                                                                                              >  1:301696/1‑251 (MQ=255)
                                                                                    aaaTCCAGATCTGAATCATTCAGGTAAATAACATTCATTTTTTGTGGTTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTATGACCGGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCACTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTCTGAAGGAATATTCAGACTCTGACTGGCGCTACGGGCGGCTTTGa                                                                                                                                                             <  2:338026/251‑1 (MQ=255)
                                                                                         cAGATCTGAATCATTCAGGTAATTAACATTCATTTTTTGTGGCTTCTATATTCTGGCGTTAGTCGTCGCCGCTAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTATGACCCGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTTTCGCTGTTAATGTTTATCCAGTCAGGGGTTTGCCCTTTTTTTATTTCTGAAGGAATATTCAGGCTCTGACTGGCGCTACGGGCGGCTTTGAAATa                                                                                                                                                         >  2:190788/1‑250 (MQ=255)
                                                                                                    cATTCAGGTAATTAACATTCATTTTTTGTGGTTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTATGACCGTAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCACTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTCTGAAGGAATATTCAGACTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCAcc                                                                                                                                             <  1:62587/251‑1 (MQ=255)
                                                                                                       tCAGGTAATTAACATTCATTTTTTGTGGTTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTATGACCGGCAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCACTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTCTGACGGCATATTCAGACTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCt                                                                                                                                          >  1:464684/1‑251 (MQ=255)
                                                                                                             aaTTAACATTCATTTTTTGTGGTTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTATGACCGGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCACTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTCTGAAGGAATATTCAGACTCTGACTGACGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTg                                                                                                                                    >  2:198659/1‑251 (MQ=255)
                                                                                                              aTTAACATTCATTTTTTGTGGTTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTATGACCGGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCACTGTTAATGTTGATCCAGTCAGTg                                                                                                                                                                                                                             >  1:217313/1‑161 (MQ=255)
                                                                                                              aTTAACATTCATTTTTTGTGGTTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTATGACCGGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCACTGTTAATGTTGATCCAGTCAGTg                                                                                                                                                                                                                             <  2:217313/161‑1 (MQ=255)
                                                                                                                   cattcattTTTTGTGGCTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTATGACCCGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCGCTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTCTGAAGGAATATTCAAGCTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGTAAATc                                                                                                                              <  1:235346/251‑1 (MQ=255)
                                                                                                                    attcattTTTTGTGGTTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTATGACCGGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCACTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTCTGAAGGAATATTCAGACTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGCTGCACCGCTTAACTGTAAATCg                                                                                                                             >  1:434008/1‑251 (MQ=255)
                                                                                                                          tttttGTGGTTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTATGACCGGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCACTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTCTGAAGGAATATTCAGACTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGTAAATCGCCATg                                                                                                                        <  1:108081/250‑1 (MQ=255)
                                                                                                                                     tCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTATGACCGGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCACTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTCTGAAGGAATATTCAGACTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGTAAATCGCCATGGTCGGCagaga                                                                                                             <  1:156886/250‑1 (MQ=255)
                                                                                                                                                     aGTCGTCGCCGATAATTTTCAGCGTGGCCATATCCTATGAGTTCACCGTATGACCGGAAAAGGTGATTGTTGAGACGCAGCGTTTATTGTCGTTATCACTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTCTGAAGGAATATTCAGACTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTcac                                                                                             <  2:227001/250‑1 (MQ=255)
                                                                                                                                                        cgtcgCCGATAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTATGACCGGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCACTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTCTGAAGGAATATTCAGACTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATg                                                                                         <  1:106425/251‑1 (MQ=255)
                                                                                                                                                         gtcgCCGATAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTATGACCGGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCACTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTCTGAAGGAATATTCAGACTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGc                                                                                        >  1:375659/1‑251 (MQ=255)
                                                                                                                                                               gATAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTATGACCGGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCACTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTCTGATGGAATATTAAGACTCTGACTGGCGCAACGGGCGGCTTTGAAATAAACAGATGCACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACaa                                                                                   >  2:405988/1‑250 (MQ=255)
                                                                                                                                                               gATAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTATGACCGGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCACTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTCTGAAGGAATATTCAGACTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACaaa                                                                                  <  2:305758/251‑1 (MQ=255)
                                                                                                                                                                aTAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTATGACCGGAAAAGGGGATTTTTGAGACGCAGCGTTTATTGTCGTTATCACTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTCTGAAGGAATATTCAGACTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAAc                                                                                 >  2:221978/1‑251 (MQ=255)
                                                                                                                                                                                           gAGTTCACCGTATGACCGGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCACTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTCTGAAGGAATATTCAGACTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGATCGTTGGt                                                      >  1:30140/1‑251 (MQ=255)
                                                                                                                                                                                                   cGTATGACCGGACAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCACTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTCTGAAGGAATATTCAGACTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGATCGTTGGTTTCGTTAc                                              >  2:106119/1‑251 (MQ=255)
                                                                                                                                                                                                    gtaggaacggAAAAGGTGATTTTTGAGACGCAGAGTTTATTGTCGTTAACACTATTAATGTTGATCGAAACAGTGGTATGCCCTTATTTTATTTCTGAAGGAATATTCAGACTCTGACTGTCGCTACGGGCGGCATTGAAATAAACCGCTGCACCGCTTAACTTTAAATCGACATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCAACAAACAGGAAGGTTCTGCACCAGATCGTTGGATTCGTTa                                               >  2:394057/10‑249 (MQ=255)
                                                                                                                                                                                                      aTGACCCGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCGCTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTCTGAAGGAATATTCAGGCTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGATCGTTGGTTTCGTTAcgcg                                           <  1:190788/251‑1 (MQ=255)
                                                                                                                                                                                                               aaaGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCACTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTCTGAAGGAATATTCAGACTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGATCGTTGGTTTCGTTACGCGGCATTGCa                                   >  1:177024/1‑250 (MQ=255)
                                                                                                                                                                                                               aaaGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCACTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTCTGAAGGAATATTCAGACTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGATCGTTGGTTTCGTTACGCGGCATTGCa                                   >  1:427913/1‑250 (MQ=255)
                                                                                                                                                                                                                 aGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCACTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTCTGAAGGAATATTCAGACTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGTAAATCGCCATGGGCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTGCAGCGCCAGATCGTTGGTTTCGTTACGCGGCATTGCAATg                                <  2:327282/251‑1 (MQ=255)
                                                                                                                                                                                                                       tttttGAGACGCAGCGTTTATTGTCGTTATCACTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTCTGCAGGAATATTCAGACTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACAGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGATCGTTGGTTTCGTTACGCGGCATTGCaagaacgccg                          >  2:174847/1‑243 (MQ=255)
                                                                                                                                                                                                                           tGAGACGCAGCGTTTATTGTCGTTATCACTGATAATGTTGCTCAAGTCAGTGGTTTGCCCTTCTTTTATTTCTGAAGGAATATTCAGACTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGGTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGATCGTTGGTTTCGTTACGCGGCATTGCAATgaagacgagg                       >  1:416853/1‑241 (MQ=255)
                                                                                                                                                                                                                                   ggCGTTTATGGTGGTATCACTGTTAATGTTTATCCCGTCAGTGGGTTGCCCTTCTTTTATTTCTGAAGGAATATTCAGGCTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGGAAATCGCCCTGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGATCGTTGGTTTCGTTACGCGGCATTGCAATGACGCCGAGGAGTTTATGGt             <  1:118942/250‑1 (MQ=255)
                                                                                                                                                                                                                                             gTCGTTATCACTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTCTGAAGGAATATTCAGACTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGATCGTTGGTTTCGTTACGCGGCATTGCAATGACGCCGAGGAGTTTATGGTCGTTTGCCt    <  1:106119/251‑1 (MQ=255)
                                                                                                                                                                                                                                               cGTTATCACTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTCTGAAGGAATATTCAGACTCTGACTGGCGCTACGGTCGGCTTTGAAATAAACCGATGCACCGCTTAACTGTAACTCGCCATGGTCGGCAGCGAGTTGTATGCTTTTCACAATGCGACAACCAGGAAGTTTCAGCGCCAGATCGTTGGTTGCGTAATACGGCATTGCATTGACGCCGAGGAGTTTCTGATCGTTTGCCTgc  >  2:371128/1‑251 (MQ=255)
                                                                                                                                                                                                                                               aggtATCACTGTTAATGTTGATCCAGTAAGTGGTTTGCCATTCTTTTATTTCTGAAGGAATATTCAGACTCTGGCCGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGTAAATAGCCATGGTCGGCAGGGAGATGTATGAGTTTCACAATGCACCACACAGGAAGTTTCAGCGCCCGCTTGTTGTTTTCGTTACGCGGCATTGAAATGACGCCGAGGAGATTATGGGAGTTTCCCTgc  >  2:253303/4‑251 (MQ=255)
                                                                                                                                                                                                                                                  tATCACTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTCTGAAGGAATATTCAGACTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGATCGTTGGTTTCGTTACGCGGCATTGCAATgacg                             >  1:367994/1‑218 (MQ=255)
                                                                                                                                                                                                                                                  tATCACTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTCTGAAGGAATATTCAGACGCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTGTCAGCGCCAGATCGTTGGTTTCGTTACGCGGCATTGCAATgacg                             <  2:367994/221‑4 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                   
CTGGAGAGGCGAGTTTTGCCTTTTTCATTATGGGTGAGCAATCGGGCGAAATTTGCCAACTGTTCCTCACTACAATGCTGAAGAAAATCCAGATCTGAATCATTCAGGTAATTAACATTCATTTTTTGTGGCTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTATGACCCGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCGCTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTCTGAAGGAATATTCAGGCTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGATCGTTGGTTTCGTTACGCGGCATTGCAATGGCGCCGAGGAGTTTATGGTCGTTTGCCTGC  >  NC_000913/11235‑11724

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