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

AGATCAACCTGAATTCCAGTGATCCGCTTCCTACTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAGACAGCGACAACACTGGAATTTAATGCTAATGGCACATTAGTGGATGGTGCGATGGCGAATAATATCGCAACCGGCGCAATTAACGGTGCAGAACCCGCCACGTTTAGTCTGAGCTTCCTCAACTCCATGCAGCAAAATACCGGCGCTAACAATATTGTGGCAACCACCCAGAACGGCTACAAACCGGGCGATCTGGTGAGTTATCAAATCAATGATGACGGTACGGTTGTCGGCAACTATTCCAACGAACAAACCCAACT  >  NC_000913/1133037‑1133502
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                     
aGATCAACCTGAATTCCAGTGATCCGCTTCCTACTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAGACAGCGACAACACTGGAATTTAATGCTAATAGCACATTAGTGGAt                                                                                                                                                                                                                         <  1:152414/251‑1 (MQ=255)
                cAGTGATCCGCTTCCTACTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAGACAGCGACAACACTGGAATTTAATGCTAATAGCACATTAGTGGATGGTGCGATGGCGaat                                                                                                                                                                                                          >  1:8414/1‑250 (MQ=255)
                       ccGCTTCCTACTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAGACAGCGACAACACTGGAATTTAATGCTAATAGCACATTAGTGGATGGTGCGATGGCGAATAATATCGc                                                                                                                                                                                                  >  1:251882/1‑251 (MQ=255)
                            tCCTACTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAGACAGCGACAACACTGGAATTTAATGCTAATAGCACATTAGTGGATGGTGCGATGGCGAATAATATCGCAAcc                                                                                                                                                                                              <  2:291431/250‑1 (MQ=255)
                                  tGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAGACAGCGACAACACTGGAATTTAATGCTAATAGCACATTAGTGGATGGTGCGATGGCGAATAATATCGCAACCGGCGCa                                                                                                                                                                                        >  1:364354/1‑250 (MQ=255)
                                              cAGCGCCAGCCATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAGACAGCGAAAACACTGGAATTTAATGCTAATAGCACATTAGTGGATGGTGCGATGGCGAATAATATAGCAACCGGCGCAATTAACGGTga                                                                                                                                                                             >  2:325073/1‑248 (MQ=255)
                                                 cgcCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAGACAGCGACAACACTGGAATTTAATGCTAATAGCACATTAGTGGATggggc                                                                                                                                                                                                                    >  2:392522/1‑204 (MQ=255)
                                                 cgcCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAGACAGCGACAACACTGGAATTTAATGCTAATAGCACATTAGTGGATGGTGc                                                                                                                                                                                                                    <  1:392522/207‑1 (MQ=255)
                                                       cAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGATATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAGACAGCGACAACACTGGAATTTAATGCTAATAGCACATTAGTGGATGGTGCGATGGCGAATAATATCGCAACCGGCGCAATTAACGGTGCAGAACCCGccc                                                                                                                                                                  >  1:387015/1‑250 (MQ=255)
                                                       aaaTGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAGACAGCGACAACACTGGAATTTAATGCTAATAGCACATTAGTGGATGGTGCGATGGCGAATAATATCGCAACCGGCGCAATTAACGGTGCAGAACCCGCCa                                                                                                                                                                  <  1:164015/250‑1 (MQ=255)
                                                                  gttaaaaCAAAAAAGGTTGGGGGACTGTTTCAGACAGTCAGGGTAATGCTCAGGACATGAGCGTTTCTTTGGTGAAGACCGGGAATAATAACGGGCAGGTCTACCCCCAGGATAGCAGTGATCCAAACAGCATTGGGAAGACAGCGACAACACGGAAATTTAATGCTAATAGCACATTAGGGGATGGTGCGAGGGCGAATAATATCGCAACCGGCGCAATTAACGGTGCAGAACCCGCCACGTTTAGTCTg                                                                                                                                                       <  2:453108/246‑1 (MQ=255)
                                                                        aCAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAGACAGCGACAACACTGGAATTTAATGCTAATAGCACATTAGTGGATGGTGCGATGGCGAATAATATCGCAACCGGCGCAATTAACGGTGCAGAACCCGCCACGTTTAGTCTGAGCTTc                                                                                                                                                 >  1:67482/1‑251 (MQ=255)
                                                                              aaGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAGACAGCGACAACACTGGAATTTAATGCTAATAGCACATTAGTGGATGGTGCGATGGCGAATAATATCGCAACCGGCGCAATTAACGGTGCAGAACCCGCCACGTTTAGTCTGAGCTTCCTCAAc                                                                                                                                           >  1:461397/1‑251 (MQ=255)
                                                                               aGGTCCGGTGACTGTTTTCGAAAGACAGGGTAGGGCTCATGACAGGAGCGTCTCCTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGAATTGCGAAGACAGCGACAACACTGGAATTTAATGCTAATAGCACATTAGTGGATGGTGCGATGGCGAATAATATCGCAACCGGCGCAATTAACGTTGCAGAACCCGCCACGTTTAGTCTGAGCTTCCTCAACt                                                                                                                                          <  1:96723/251‑1 (MQ=255)
                                                                                        gccTGTTTTCGACAGTCAGGGAAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACCCCCAGGATAGCAGTGATCCAAACAGCATTGCGAAGACAGCGACAACACTGGAATTTAATGCTAATAGCACATTAGTGGATGGTGCGATGGCGAATAATATCGCAACCGGCGCAATTAACGGTGCAGAACCCGCCACGTTTAGTCTGAGCTTCCTCAACTCCATgcagc                                                                                                                                 <  1:346039/249‑1 (MQ=255)
                                                                                              tttCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTAATTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAGACAGCGACAACACTGGAATTTAATGCTAATAGCACATTAGTGGATGGTGCGATGGCGAATAATATCGCAACCGGCGCAATTAACGGTGCAGAACCCGCCACGTTTAGTCTGAGCTTCCTCAACTCCATGCAGCAAAAt                                                                                                                            >  1:18388/1‑250 (MQ=255)
                                                                                                      gTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAGACAGCGACAACACTGGAATTTAATGCTAATAGCACATTAGTGGATGGTGCGATGGCGAATAATATCGCAACCGGCGCAATTAACGGGGCAGAACCCGCCCCGTTTAATCTGAGCTTCCTCAACTCCATGCAGCCAAATACCCGCGCt                                                                                                                   >  2:97378/1‑251 (MQ=255)
                                                                                                         aGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGAGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAGACAGCGACAACACTGGAATTTAATGCTAATGGCACATTAGTGGATGGTGCGATGGCGAATAATATCGCAACCGGCGCAATTAACGGGTCAGAAACCGCCACGTTTAGTCTGAGCTTCCTCAACTCCATGCAGCAAAATACCGGCGCTAAc                                                                                                                >  2:422485/1‑251 (MQ=255)
                                                                                                                tGCTCATGACATGCGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTATACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAGACAGCGACAACACTGGAATTTAATGCTAATAGAACATTAGTGGATGGTGCGATGGCGAATAATATCGCAACAGGCGCAATTAACGGTGCAGAACCCGCCACGTTTAGTCTGAGCTTCCTCAACTCCATGCAGCAAAATACCGGCGCTAACAATATtg                                                                                                         >  1:105365/1‑251 (MQ=255)
                                                                                                                 gCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAGACAGCGACAACACTGGAATTTAATGCTAATAGCACATTAGTGGATGGTGCGATGGCGAATAATATCGCAACCGGCGCAATTAACGGTGCAGAACCCGCCACGTTTAGTCTGAGCTTCCTCAACTCCATGCAGCAAAATACCGGCGCTAACAATATtgt                                                                                                        >  1:199004/1‑251 (MQ=255)
                                                                                                                  cTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAGACAGCGACAACACTGGAATTTAATGCTAATAGCACATTAGTGGATGGTGCGATGGCGAATAATATCGCAACCGGCGCAATTAACGGTGCAGAAACCGCCACGTTTAGTCTGAGCTTCCTCAACTCCATGCAGCAAAATACCGGCGCTAACAATATtgtg                                                                                                       >  2:464874/1‑251 (MQ=255)
                                                                                                                         catgaGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGACCCAAACAGCATTGCGAAGACAGCGACAACACTGGAATTTAATGCTAATAGCACATTAGTGGATGGTGCGATGGCGAATAAAATCGCAACCGGCGCAATTAACGGTGCAGAACCCGCCACGTTTAGTCTGAGCTTCCTCAACTCCATGCAGCAAAATACCGGCGCTAACAATATTGTGGCAacca                                                                                                <  1:97378/251‑1 (MQ=255)
                                                                                                                               cGTCTACTCTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCAAAGACAGCGACAACACTGGAATTTAATGCTAATGGCACATTAGTGGATGGTGCTATGGCGAATAATATCGCAACCGGCGCAATTAACGGTGCAGAACCCGCCACGTTTAGTCTGAGCTTCCTCa                                                                                                                                             >  1:135426/1‑200 (MQ=255)
                                                                                                                               cGTCTACTCTGTGAAGACCGGGGATAATAACGGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGAAAAGACAGCGACAACACTGGAATTAAATGCTAATGGCACATTAGTGGATGGTGCTAGGGCGAATAATATCGCACCCGGCCCAATTAACGGTGCAGAACCCGCCACGTTTAGTCTGAGCTTCCTCa                                                                                                                                             <  2:135426/200‑1 (MQ=255)
                                                                                                                                  cTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAGACAGCGACAACACTGGAATTTAATGCTAATAGCACATTAGTGGATGGTGCGATGGCGAATAATATCGCAACCGGCGCAATTAACGGTGCAGAACCCGCCACGTTTAGTCTGAGCTTCCTCAACTCCATGCAGCAAAATACCGGCGCTAACAATATTGTGGCAACCACCCAGAACg                                                                                       >  1:320981/1‑251 (MQ=255)
                                                                                                                                   tACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAGACAGCGACAACACTGGAATTTAATGCTAATAGCACATTAGTGGATGGTGCGATGGCGAATAATATCGCAACCGGCGCAATTAACGGTGCAGAACCCGCCACGTTTAGTCTGAGCTTCCTCAACTCCATGCAGCAAAATACCGGCGCTAACAATATTGTGGCAACCACCCAGAACgg                                                                                      >  1:77446/1‑251 (MQ=255)
                                                                                                                                        gtgggAAGACCGGGGAAAAAAACGGGCAGTTCTACCCCCGGGATAGTAGTTATCCAAACAGCATTGGAAGGACAGGGAAAACACTGGAATTAAATGCAAAAAGCACATTAGTGGATGGTGCGATGGCGAATAATATCGCACCCGGCGCAATTAACGGTGCATACCCCGCCACGTTTAGTCTGAGCTTCCTCAACTCCATGCAGCAAAATACCGGCGCTAACAATATTGTGGCAACCACCCAGAACGGCTa                                                                                   <  2:18388/246‑1 (MQ=255)
                                                                                                                                              aGACCGGGGATAAAAACGGGCAGTTCTCCACCCAGGATAGAAGTGATCCAAACAGCATTGCGAAGACAGCGACAACACTGGAATTTAATGCTAATAGCACATTAGTGGATGGTGCGATGGCGAATAATATCGCACCCGGCGCAATTAACGGTGCAGAACCCGCCACGTTTAGTCTGAGCTTCCTCAACTCCATGCAGCAAAATACCGGCGCTAACAATATTGTGGCAACCACCCAGAACGGCTACAAACCg                                                                            <  1:343590/251‑1 (MQ=255)
                                                                                                                                                 ccGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAGACAGCGACAACACTGGAATTTAATGCTAATAGCACATTAGTGGATGGTGCGATGGCGAATAATATCGCAACCGGCGCAATTAACGGTGCAGAACCCGCCACGTTTAGTCTGAGCTTCCTCAACTCCATGCAGCAAAATACCGGCGCTAACAATATTGTGGCAACCACCCAGAACGGCTACAAACCgggg                                                                         >  2:327564/1‑250 (MQ=255)
                                                                                                                                                      aataataACTGCCAGGTCTACACCAAGGATAGCAGTTATCCAAACAGCATTGCTAATACAGCAACAACACTGGAATTTAATGCTAATACCACATTAGTGGATGGGGCGATGGCGAATAATATCGCAACCGGCGCAATAAACGGTGCAGAACCCGCCACGTTTAGTCTGAGCTTCCTCAACTCCATGCAGCAAAATACCGGCGCTAACAATATTGTGGCAACCACCCAGAACGGCTACAAACCGGGCGATCt                                                                    <  1:126006/250‑1 (MQ=255)
                                                                                                                                                                ggCAGGTCTACACCAAGGATAGCAGTGATCCAAACAGCATTGCGAAGACAGCGACAACACTGGAATTTAATGCTAATAGCACATTAGTGGATGGTGCGATGGCGAATAATATCGCAACCGGCGCAATTAACGGTGCAGAACCCGCCACGTTTAGTCTGAGCTTCCTCAACTCCATGCAGCAAAATACCGGCGCTAACAATATTGTGGCAACCACCCAGAACGGCTACAAACCGGGCGATCTGGTGAGTTat                                                          <  1:52340/251‑1 (MQ=255)
                                                                                                                                                                  cAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAGACAGCGACAACACTGGAATTTAATGCTAATAGCACATTAGTGGATGGTGCGATGGCGAATAATATCGCAACCGGCGCAATTAACGGTGCAGAACCCGCCACGTTTAGTCTGAGCTTCCTCAACTCCATGCAGCAAAATACCGGCGCTAACAATATTGTGGCAACCACCCAGAAAGGCTACAAACCAGGCGATCTGGTGAGTTatc                                                         >  1:15813/1‑250 (MQ=255)
                                                                                                                                                                         acacCCAGGATAGAATTGATCCAAACAGCATTGCGAAGACAGCGAAAACACGGGAATTTAATGCTAATAGCACATTAGTGGATGGTGCGATGGCGAATAATATCGCAACCGGCGCAATTAACGGGGCAGAACCCGCCACGTTTAGTCTGAGCTGCCTCAACTCCATGCAGCAAAATACCGGCGCTAACAATATTGTGGCAACCACCCAGAACGGCTACAAACCGGGCGATCTGGTGAGTTatcaaatc                                                    <  1:325073/248‑1 (MQ=255)
                                                                                                                                                                                          gacccAAACAGCATTGCAAAGACAGCGACAACACTGGAATTTAAGGCTAAAAGCACATTAGTGGAGGGGGGGATGGCGAATAATATCGCACCCGGCGCAATTAACGGTGCAGAACCCGCCACGTTTAGTCTGAGCTTCCTCAACTCCATGCAGCAAAATACCGGCGCTAACAATATTGTGGCAACCACCCAGAGCGGCTACAAACCGGGCGATCTGGTGAGTTATCAAATCAATGATGACGGTACGGTTg                                 <  1:23548/247‑1 (MQ=255)
                                                                                                                                                                                                          cGAAGACAGCGACAACACTGGAATTTAATGCTAATAGCACATTAGTGGATGGTGCGATGGCGAATAATGTCGCAACCGGCGCAATTAACGGTGCAGAACCCGCCACGTTTAGTCTGAGCTTCCTCAACTCCATGCAGCAAAATACCGGAGCTAACAATATTGTGGCAACCACCCAGAACGGCTACAAACCGGGCGATCTGGTGAGTTATCAAATCAATGCTGACGGTCCGATTGTCGGCAACTATTCCAAc                >  1:297391/1‑251 (MQ=255)
                                                                                                                                                                                                             aGACAGCGAAAACACTGGAATTTAATGCTAATAGCACATTAGTGGATGGTGCGATGGCAAATAATATCGCAACCGGCGCAATTAACGGTGCAGAACCCGCCACGTTTAGTCTGAGCTTCCTCAACTCCATGCAGCAAAATACCGGCGCTAACAATATTGTGGCAACCACCCAGAACGGCTACAAACCGGGCGATCTGGTGAGTTATCAAATCAATGATGACGGTACGGTTGTCGGCAACTATTCCAACGaa             <  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)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                     
AGATCAACCTGAATTCCAGTGATCCGCTTCCTACTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAGACAGCGACAACACTGGAATTTAATGCTAATGGCACATTAGTGGATGGTGCGATGGCGAATAATATCGCAACCGGCGCAATTAACGGTGCAGAACCCGCCACGTTTAGTCTGAGCTTCCTCAACTCCATGCAGCAAAATACCGGCGCTAACAATATTGTGGCAACCACCCAGAACGGCTACAAACCGGGCGATCTGGTGAGTTATCAAATCAATGATGACGGTACGGTTGTCGGCAACTATTCCAACGAACAAACCCAACT  >  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: