Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,133,273 G→A G234S (GGC→AGC)  flgE → flagellar hook protein FlgE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,133,2730GA92.5% 100.6 / 1.1 40G234S (GGC→AGC) flgEflagellar hook protein FlgE
Reads supporting (aligned to +/- strand):  ref base G (2/1);  new base A (20/17);  total (22/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.04e-01

AGATCAACCTGAATTCCAGTGATCCGCTTCCTACTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTG‑TGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAG‑ACAGCGACAACACTGGAATTTAATGCTAATGGCACATTAGTGGATGGTGCGATGGCGAATAATATCGCAACCGGCGCAATTAACGGTGCAGAACCCGCCACGTTTAGTCTGAGCTTCCTCAACTCCATGCAGCAAAATACCGGCGCTAACAATATTGTGGCAACCACCCAGAACGGCTACAAACCGGGCGATCTGGTGAGTTATCAAATCAATGATGACGGTACGGTTGTCGGCAACTATTCCAACGAACAAACCCAACT  >  NC_000913/1133037‑1133502
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                     
aGATCAACCTGAATTCCAGTGATCCGCTTCCTACTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTG‑TGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAG‑ACAGCGACAACACTGGAATTTAATGCTAATAGCACATTAGTGGAt                                                                                                                                                                                                                         <  1:152414/251‑1 (MQ=255)
                cAGTGATCCGCTTCCTACTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTG‑TGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAG‑ACAGCGACAACACTGGAATTTAATGCTAATAGCACATTAGTGGATGGTGCGATGGCGaat                                                                                                                                                                                                          >  1:8414/1‑250 (MQ=255)
                       ccGCTTCCTACTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTG‑TGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAG‑ACAGCGACAACACTGGAATTTAATGCTAATAGCACATTAGTGGATGGTGCGATGGCGAATAATATCGc                                                                                                                                                                                                  >  1:251882/1‑251 (MQ=255)
                            tCCTACTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTG‑TGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAG‑ACAGCGACAACACTGGAATTTAATGCTAATAGCACATTAGTGGATGGTGCGATGGCGAATAATATCGCAAcc                                                                                                                                                                                              <  2:291431/250‑1 (MQ=255)
                                  tGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTG‑TGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAG‑ACAGCGACAACACTGGAATTTAATGCTAATAGCACATTAGTGGATGGTGCGATGGCGAATAATATCGCAACCGGCGCa                                                                                                                                                                                        >  1:364354/1‑250 (MQ=255)
                                              cAGCGCCAGCCATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTG‑TGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAG‑ACAGCGAAAACACTGGAATTTAATGCTAATAGCACATTAGTGGATGGTGCGATGGCGAATAATATAGCAACCGGCGCAATTAACGGTga                                                                                                                                                                             >  2:325073/1‑248 (MQ=255)
                                                 cgcCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTG‑TGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAG‑ACAGCGACAACACTGGAATTTAATGCTAATAGCACATTAGTGGATggggc                                                                                                                                                                                                                    >  2:392522/1‑204 (MQ=255)
                                                 cgcCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTG‑TGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAG‑ACAGCGACAACACTGGAATTTAATGCTAATAGCACATTAGTGGATGGTGc                                                                                                                                                                                                                    <  1:392522/207‑1 (MQ=255)
                                                       cAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGATATGAGCGTCTACTTTG‑TGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAG‑ACAGCGACAACACTGGAATTTAATGCTAATAGCACATTAGTGGATGGTGCGATGGCGAATAATATCGCAACCGGCGCAATTAACGGTGCAGAACCCGccc                                                                                                                                                                  >  1:387015/1‑250 (MQ=255)
                                                       aaaTGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTG‑TGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAG‑ACAGCGACAACACTGGAATTTAATGCTAATAGCACATTAGTGGATGGTGCGATGGCGAATAATATCGCAACCGGCGCAATTAACGGTGCAGAACCCGCCa                                                                                                                                                                  <  1:164015/250‑1 (MQ=255)
                                                                  gttaaaaCAAAAAAGGTTGGGGGACTGTTTCAGACAGTCAGGGTAATGCTCAGGACATGAGCGT‑TTCTTTGGTGAAGACCGGGAATAATAACGGGCAGGTCTACCCCCAGGATAGCAGTGATCCAAACAGCATTGGGAAG‑ACAGCGACAACACGGAAATTTAATGCTAATAGCACATTAGGGGATGGTGCGAGGGCGAATAATATCGCAACCGGCGCAATTAACGGTGCAGAACCCGCCACGTTTAGTCTg                                                                                                                                                       <  2:453108/246‑1 (MQ=255)
                                                                        aCAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTG‑TGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAG‑ACAGCGACAACACTGGAATTTAATGCTAATAGCACATTAGTGGATGGTGCGATGGCGAATAATATCGCAACCGGCGCAATTAACGGTGCAGAACCCGCCACGTTTAGTCTGAGCTTc                                                                                                                                                 >  1:67482/1‑251 (MQ=255)
                                                                              aaGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTG‑TGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAG‑ACAGCGACAACACTGGAATTTAATGCTAATAGCACATTAGTGGATGGTGCGATGGCGAATAATATCGCAACCGGCGCAATTAACGGTGCAGAACCCGCCACGTTTAGTCTGAGCTTCCTCAAc                                                                                                                                           >  1:461397/1‑251 (MQ=255)
                                                                               aGGTCCGGTGACTGTTTTCGAAAGACAGGGTAGGGCTCATGACAGGAGCGTCTCCTTTG‑TGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGAATTGCGAAG‑ACAGCGACAACACTGGAATTTAATGCTAATAGCACATTAGTGGATGGTGCGATGGCGAATAATATCGCAACCGGCGCAATTAACGTTGCAGAACCCGCCACGTTTAGTCTGAGCTTCCTCAACt                                                                                                                                          <  1:96723/251‑1 (MQ=255)
                                                                                        gccTGTTTTCGACAGTCAGGGAAATGCTCATGACATGAGCGTCTACTTTG‑TGAAGACCGGGGATAATAACTGGCAGGTCTACCCCCAGGATAGCAGTGATCCAAACAGCATTGCGAAG‑ACAGCGACAACACTGGAATTTAATGCTAATAGCACATTAGTGGATGGTGCGATGGCGAATAATATCGCAACCGGCGCAATTAACGGTGCAGAACCCGCCACGTTTAGTCTGAGCTTCCTCAACTCCATgcagc                                                                                                                                 <  1:346039/249‑1 (MQ=255)
                                                                                              tttCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTAATTTG‑TGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAG‑ACAGCGACAACACTGGAATTTAATGCTAATAGCACATTAGTGGATGGTGCGATGGCGAATAATATCGCAACCGGCGCAATTAACGGTGCAGAACCCGCCACGTTTAGTCTGAGCTTCCTCAACTCCATGCAGCAAAAt                                                                                                                            >  1:18388/1‑250 (MQ=255)
                                                                                                      gTCAGGGTAATGCTCATGACATGAGCGTCTACTTTG‑TGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAG‑ACAGCGACAACACTGGAATTTAATGCTAATAGCACATTAGTGGATGGTGCGATGGCGAATAATATCGCAACCGGCGCAATTAACGGGGCAGAACCCGCCCCGTTTAATCTGAGCTTCCTCAACTCCATGCAGCCAAATACCCGCGCt                                                                                                                   >  2:97378/1‑251 (MQ=255)
                                                                                                         aGGGTAATGCTCATGACATGAGCGTCTACTTTG‑TGAAGACCGAGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAG‑ACAGCGACAACACTGGAATTTAATGCTAATGGCACATTAGTGGATGGTGCGATGGCGAATAATATCGCAACCGGCGCAATTAACGGGTCAGAAACCGCCACGTTTAGTCTGAGCTTCCTCAACTCCATGCAGCAAAATACCGGCGCTAAc                                                                                                                >  2:422485/1‑251 (MQ=255)
                                                                                                                tGCTCATGACATGCGCGTCTACTTTG‑TGAAGACCGGGGATAATAACTGGCAGGTATACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAG‑ACAGCGACAACACTGGAATTTAATGCTAATAGAACATTAGTGGATGGTGCGATGGCGAATAATATCGCAACAGGCGCAATTAACGGTGCAGAACCCGCCACGTTTAGTCTGAGCTTCCTCAACTCCATGCAGCAAAATACCGGCGCTAACAATATtg                                                                                                         >  1:105365/1‑251 (MQ=255)
                                                                                                                 gCTCATGACATGAGCGTCTACTTTG‑TGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAG‑ACAGCGACAACACTGGAATTTAATGCTAATAGCACATTAGTGGATGGTGCGATGGCGAATAATATCGCAACCGGCGCAATTAACGGTGCAGAACCCGCCACGTTTAGTCTGAGCTTCCTCAACTCCATGCAGCAAAATACCGGCGCTAACAATATtgt                                                                                                        >  1:199004/1‑251 (MQ=255)
                                                                                                                  cTCATGACATGAGCGTCTACTTTG‑TGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAG‑ACAGCGACAACACTGGAATTTAATGCTAATAGCACATTAGTGGATGGTGCGATGGCGAATAATATCGCAACCGGCGCAATTAACGGTGCAGAAACCGCCACGTTTAGTCTGAGCTTCCTCAACTCCATGCAGCAAAATACCGGCGCTAACAATATtgtg                                                                                                       >  2:464874/1‑251 (MQ=255)
                                                                                                                         catgaGCGTCTACTTTG‑TGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGACCCAAACAGCATTGCGAAG‑ACAGCGACAACACTGGAATTTAATGCTAATAGCACATTAGTGGATGGTGCGATGGCGAATAAAATCGCAACCGGCGCAATTAACGGTGCAGAACCCGCCACGTTTAGTCTGAGCTTCCTCAACTCCATGCAGCAAAATACCGGCGCTAACAATATTGTGGCAacca                                                                                                <  1:97378/251‑1 (MQ=255)
                                                                                                                               cGTCTACTCTG‑TGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCAAAG‑ACAGCGACAACACTGGAATTTAATGCTAATGGCACATTAGTGGATGGTGCTATGGCGAATAATATCGCAACCGGCGCAATTAACGGTGCAGAACCCGCCACGTTTAGTCTGAGCTTCCTCa                                                                                                                                             >  1:135426/1‑200 (MQ=255)
                                                                                                                               cGTCTACTCTG‑TGAAGACCGGGGATAATAACGGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGAAAAG‑ACAGCGACAACACTGGAATTAAATGCTAATGGCACATTAGTGGATGGTGCTAGGGCGAATAATATCGCACCCGGCCCAATTAACGGTGCAGAACCCGCCACGTTTAGTCTGAGCTTCCTCa                                                                                                                                             <  2:135426/200‑1 (MQ=255)
                                                                                                                                  cTACTTTG‑TGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAG‑ACAGCGACAACACTGGAATTTAATGCTAATAGCACATTAGTGGATGGTGCGATGGCGAATAATATCGCAACCGGCGCAATTAACGGTGCAGAACCCGCCACGTTTAGTCTGAGCTTCCTCAACTCCATGCAGCAAAATACCGGCGCTAACAATATTGTGGCAACCACCCAGAACg                                                                                       >  1:320981/1‑251 (MQ=255)
                                                                                                                                   tACTTTG‑TGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAG‑ACAGCGACAACACTGGAATTTAATGCTAATAGCACATTAGTGGATGGTGCGATGGCGAATAATATCGCAACCGGCGCAATTAACGGTGCAGAACCCGCCACGTTTAGTCTGAGCTTCCTCAACTCCATGCAGCAAAATACCGGCGCTAACAATATTGTGGCAACCACCCAGAACgg                                                                                      >  1:77446/1‑251 (MQ=255)
                                                                                                                                        gtgggAAGACCGGGGAAAAAAACGGGCAGTTCTACCCCCGGGATAGTAGTTATCCAAACAGCATTG‑GAAGGACAGGGAAAACACTGGAATTAAATGCAAAAAGCACATTAGTGGATGGTGCGATGGCGAATAATATCGCACCCGGCGCAATTAACGGTGCATACCCCGCCACGTTTAGTCTGAGCTTCCTCAACTCCATGCAGCAAAATACCGGCGCTAACAATATTGTGGCAACCACCCAGAACGGCTa                                                                                   <  2:18388/246‑1 (MQ=255)
                                                                                                                                              aGACCGGGGATAAAAACGGGCAGTTCTCCACCCAGGATAGAAGTGATCCAAACAGCATTGCGAAG‑ACAGCGACAACACTGGAATTTAATGCTAATAGCACATTAGTGGATGGTGCGATGGCGAATAATATCGCACCCGGCGCAATTAACGGTGCAGAACCCGCCACGTTTAGTCTGAGCTTCCTCAACTCCATGCAGCAAAATACCGGCGCTAACAATATTGTGGCAACCACCCAGAACGGCTACAAACCg                                                                            <  1:343590/251‑1 (MQ=255)
                                                                                                                                                 ccGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAG‑ACAGCGACAACACTGGAATTTAATGCTAATAGCACATTAGTGGATGGTGCGATGGCGAATAATATCGCAACCGGCGCAATTAACGGTGCAGAACCCGCCACGTTTAGTCTGAGCTTCCTCAACTCCATGCAGCAAAATACCGGCGCTAACAATATTGTGGCAACCACCCAGAACGGCTACAAACCgggg                                                                         >  2:327564/1‑250 (MQ=255)
                                                                                                                                                      aataataACTGCCAGGTCTACACCAAGGATAGCAGTTATCCAAACAGCATTGCTAAT‑ACAGCAACAACACTGGAATTTAATGCTAATACCACATTAGTGGATGGGGCGATGGCGAATAATATCGCAACCGGCGCAATAAACGGTGCAGAACCCGCCACGTTTAGTCTGAGCTTCCTCAACTCCATGCAGCAAAATACCGGCGCTAACAATATTGTGGCAACCACCCAGAACGGCTACAAACCGGGCGATCt                                                                    <  1:126006/250‑1 (MQ=255)
                                                                                                                                                                ggCAGGTCTACACCAAGGATAGCAGTGATCCAAACAGCATTGCGAAG‑ACAGCGACAACACTGGAATTTAATGCTAATAGCACATTAGTGGATGGTGCGATGGCGAATAATATCGCAACCGGCGCAATTAACGGTGCAGAACCCGCCACGTTTAGTCTGAGCTTCCTCAACTCCATGCAGCAAAATACCGGCGCTAACAATATTGTGGCAACCACCCAGAACGGCTACAAACCGGGCGATCTGGTGAGTTat                                                          <  1:52340/251‑1 (MQ=255)
                                                                                                                                                                  cAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAG‑ACAGCGACAACACTGGAATTTAATGCTAATAGCACATTAGTGGATGGTGCGATGGCGAATAATATCGCAACCGGCGCAATTAACGGTGCAGAACCCGCCACGTTTAGTCTGAGCTTCCTCAACTCCATGCAGCAAAATACCGGCGCTAACAATATTGTGGCAACCACCCAGAAAGGCTACAAACCAGGCGATCTGGTGAGTTatc                                                         >  1:15813/1‑250 (MQ=255)
                                                                                                                                                                         acacCCAGGATAGAATTGATCCAAACAGCATTGCGAAG‑ACAGCGAAAACACGGGAATTTAATGCTAATAGCACATTAGTGGATGGTGCGATGGCGAATAATATCGCAACCGGCGCAATTAACGGGGCAGAACCCGCCACGTTTAGTCTGAGCTGCCTCAACTCCATGCAGCAAAATACCGGCGCTAACAATATTGTGGCAACCACCCAGAACGGCTACAAACCGGGCGATCTGGTGAGTTatcaaatc                                                    <  1:325073/248‑1 (MQ=255)
                                                                                                                                                                                          gacccAAACAGCATTGCAAAG‑ACAGCGACAACACTGGAATTTAAGGCTAAAAGCACATTAGTGGAGGGGGGGATGGCGAATAATATCGCACCCGGCGCAATTAACGGTGCAGAACCCGCCACGTTTAGTCTGAGCTTCCTCAACTCCATGCAGCAAAATACCGGCGCTAACAATATTGTGGCAACCACCCAGAGCGGCTACAAACCGGGCGATCTGGTGAGTTATCAAATCAATGATGACGGTACGGTTg                                 <  1:23548/247‑1 (MQ=255)
                                                                                                                                                                                                          cGAAG‑ACAGCGACAACACTGGAATTTAATGCTAATAGCACATTAGTGGATGGTGCGATGGCGAATAATGTCGCAACCGGCGCAATTAACGGTGCAGAACCCGCCACGTTTAGTCTGAGCTTCCTCAACTCCATGCAGCAAAATACCGGAGCTAACAATATTGTGGCAACCACCCAGAACGGCTACAAACCGGGCGATCTGGTGAGTTATCAAATCAATGCTGACGGTCCGATTGTCGGCAACTATTCCAAc                >  1:297391/1‑251 (MQ=255)
                                                                                                                                                                                                             aG‑ACAGCGAAAACACTGGAATTTAATGCTAATAGCACATTAGTGGATGGTGCGATGGCAAATAATATCGCAACCGGCGCAATTAACGGTGCAGAACCCGCCACGTTTAGTCTGAGCTTCCTCAACTCCATGCAGCAAAATACCGGCGCTAACAATATTGTGGCAACCACCCAGAACGGCTACAAACCGGGCGATCTGGTGAGTTATCAAATCAATGATGACGGTACGGTTGTCGGCAACTATTCCAACGaa             <  1:208789/251‑1 (MQ=255)
                                                                                                                                                                                                                aCAGCGACAACACGGAAATTTAATGCTAATAGCACATTAGTGGATGGTGCGATGGCGAATAATATCGCAACCGGCGCAATTAACGGTGCAGAACCCGCCACGTTTAGTCTGAGCTTCCTCAACTCCATGCAGCAAAATACCGGCGCTAACAATATTGTGGCAACCACCCAGAACGGCTACAAACCGGGCGATCTGGTGAGTTATCAAATCAATGATGACGGTACGGTTGTCGGCAACTATTCCAACGAACa           <  1:327564/251‑1 (MQ=255)
                                                                                                                                                                                                                 cAGCGACAACACTGGAATTTAATGCTAATAGCACATTAGTGGATGGTGCGATGGCGAATAATATCGCAACCGGCGCAATTAACGGTGCAGAACCCGCCACGTTTAGTCTGAGCTTCCTCAACTCCATGCAGCAAAATACCGGCGCTAACAATATTGTGGCAacc                                                                                                 >  1:100196/1‑164 (MQ=255)
                                                                                                                                                                                                                 cAGCGACAACACTGAAATTTAATGCTAATAGCACATTAGTGGATGGTGCGATGGCGAATAATATCGCAACCGGCGCAATTAACGGTGCAGAACCCGCCACGTTTAGTCTGAGCTTCCTCAACTCCATGCAGCAAAATACCGGCGCTAACAATATTGTGGCAacc                                                                                                 <  2:100196/164‑1 (MQ=255)
                                                                                                                                                                                                                         acacTGGAATTTAATGCTAATAGCACATTAGTGGATGGTGCGATGGCGAATAATATCGCAACCGGCGCAATTAACGGTGCAGAACCCGCCACGTTTAGTCTGAGCTTCCTCAACTCCATGCAGCAAAATACCGGCGCTAACAATATTGTGGCAACCACCCAGAACGGCTACAACCCGGGCGATCTGGTGAGTTATCAAATCAATGATGACGGTACGGTTGACGGCAACTATTACAACGAACAAACCCAAct  >  1:80759/1‑251 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                     
AGATCAACCTGAATTCCAGTGATCCGCTTCCTACTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTG‑TGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAG‑ACAGCGACAACACTGGAATTTAATGCTAATGGCACATTAGTGGATGGTGCGATGGCGAATAATATCGCAACCGGCGCAATTAACGGTGCAGAACCCGCCACGTTTAGTCTGAGCTTCCTCAACTCCATGCAGCAAAATACCGGCGCTAACAATATTGTGGCAACCACCCAGAACGGCTACAAACCGGGCGATCTGGTGAGTTATCAAATCAATGATGACGGTACGGTTGTCGGCAACTATTCCAACGAACAAACCCAACT  >  NC_000913/1133037‑1133502

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