Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,309,712 C→T W423* (TGG→TGA)  alsA ← D‑allose ABC transporter ATP binding subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,309,7120CT94.7% 99.7 / ‑3.8 38W423* (TGG→TGA) alsAD‑allose ABC transporter ATP binding subunit
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (15/21);  total (16/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCCATGCCATAATCTCCTCTTCGCTCATGTCATCGCGATTCGTCAGGATTTGCGTCAGTCGTCCTTCGCAGAACACGGCGATGCGGTCGCAGACGGTGATAATTTCAGGTAGTTCAGATGACACCATCAGGATGACTTTTCCGTCGTCCGCCAGTTGGCGCATCACTTTGTAAATTTCGGCTTTCGCGCCAACGTCGATGCCGCGGGTAGGTTCATCGAAAATAATCACTTCCGGGCAACAGCACAGCCATTTGGAGATCAGGACTTTCTGCTGATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTAC‑GTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTA  >  NC_000913/4309465‑4309957
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                      
cccATGCCATAATCTCCTCTTCGCTCATGTCATCGCGATTCGTCAGGATTTGCGTCAGTCGTCCTTCGCAGAACACGGCGATGCGGTCGCAGACGGTGATAATTTCAGGTAGTTCAGATGACACCATCAGGATGACTTTTCCGTCGTCCGCCAGTTGGCGCATCACTTTGTAAATTTCGGCTTTCGCGCCAACGTCGATGCCGCGGGCAGGTTCATCGCAAATAATCACTTACGGGCAACAGCACAGTCAt                                                                                                                                                                                                                                                     >  1:64682/1‑251 (MQ=255)
 ccATGCCATAATCGCCGCTTCGCTCATGTCATCGCGATTCGTCAGGATTTGCGGCAGTCGTCCACCGCAGAACACGGCGATGCGGTCGCAGACGGTGATAATTTCAGGTAGTTCAGATGACACCATCAGTAGGACTTTTCCGTCGTCCGCCAGTTGGCGCATCACTTTGTAAATTTCGGCTTTCGCGCCAACGTCGATGCCGCGGGTAGGTTCATCGAAAATATTCACTTCCGGGCAACAGCACAGTCAtt                                                                                                                                                                                                                                                    <  2:361394/251‑1 (MQ=255)
   aTGCCATCATCTCCTCTTCGATCATGTCCTCGCGTTTAACCAGGATTTGCGTCAGACGGCCTTCGCATAACACGGCGATGCGGTCGCAGACGGTGATAACATCAGGCAGTTAAGATGACACCATCAGGATGACTTTTCCGTAGACCGCCTATTGGCGCATAACTTTGTAAATTTCGGCTTTCGCGCCCACGTCGATGCCGCGGTTCGGTTCAACGCAAATAATCACTTACGGGCACAAGCACATTCATTTg                                                                                                                                                                                                                                                  >  2:282573/1‑251 (MQ=255)
                   ttCGCTCATGTCATCGCGATTCGTCAGGATTTGCGTCAGTCGTCCTTCGCAGAACACGGCGATGCGGTCGCAGACGGTGATAATTTCAGGTAGTTCAGATGACACCATCAGGATGACTTTTCCGTCGTCCGCCAGTTGGCGCATCACTTTGTAAATTTCGGCTTTCGCGCCAACGTCGATGCCGCGGGTAGGTTCATCGAAAATAATCACTTCCGGGCAACAGCACAGTCATTTGGAGATCAGGACTTTct                                                                                                                                                                                                                                  >  1:254019/1‑251 (MQ=255)
                     cGCTCATGTCATCGCGATTCGTCAGGATTTGCGTCAGTCGTCCTTCGCAGAACACGGCGATGCGGTCGCAGACGGTGATAATTTCAGGTAGTTCAGATGACACCATCAGGATGACTTTTCCGTCGTCCGCCAGTTGGCGCATCACTTTGTAAATTTCGGCTTTCGCGCCAACGTCGATGCCGCGGGTAGGTTCATCGAAAATAATCACTTCCGGGCAACAGCACAGTCATTTGGAGATCAGGACTTTctgc                                                                                                                                                                                                                                <  1:400882/251‑1 (MQ=255)
                      gCTCATGTCATCGCGATTCGTCAGGATTTGCGTCAGTCGTCCTTCGCAGAACACGGCGATGCGGTCGCAGACGGTGATAATTTCAGGTAGTTCAGATGACACCATCAGGATGACTTTTCCGTCGTCCGCCAGTTGGCGCATCACTTTGTAAATTTCGGCTTTCGCGCCAACGTCGATGCCGCGGGTAGGTTCATCGAAAATAATCACTTCCGGGCAACAGCACAGTCATTTGGAGATCAGGACTTTctgct                                                                                                                                                                                                                               <  1:155934/251‑1 (MQ=255)
                       cTCATGTCATCGCGATTCGTCAGGATTTGCGTCAGTCGTCCTTCGCAGAACACGGCGATGCGGTCGCAGACGGTGATAATTTCAGGTAGTTCAGATGACACCATCAGGATGACTTTTCCGTCGTCCGCCAGTTGGCGCATCACTTTGTAAATTTCGGCTTTCGCGCCCACGTCGATGCCGCGGGTAGGTTCATCGAAAATAATCACTTCCGGGCAACAGCACAGTCATTTGGAGATCAGGACTTTctgctg                                                                                                                                                                                                                              >  1:343339/1‑251 (MQ=255)
                                        cGTCAGGATTTGCGTCAGTCGTCCTTCGCAGAACACGGCGATGCGGTCGCAGACGGTGATAATTTCAGGTAGTTCAGATGACACCATCAGGATGACTTTTCCGTCGTCCGCCAGTTGGCGCATCACTTTGTAAATTTCGGCTTTCGCGCCAACGTCGATGCCGCGGGTAGGTTCATCGAAAATAATCACTTCCGGGCAACAGCACAGTCATTTGGAGATCAGGACTTTCTGCTGATTTCCCCCGGAGAGtt                                                                                                                                                                                                             >  1:372125/1‑251 (MQ=255)
                                            aGGATTTGCGTCAGTCGTCCTTCGCAGAACACGGCGATGCGGTCGCAGACGGTGATAATTTCAGGTAGTTCAGATGACACCATCAGGATGACTTTTCCGTCGTCCGCCAGTTGGCGCATCACTTTGTAAATTTCGGCTTTCGCGCCAACGTCGATGCCGCGGGTAGGTTCATCGAAAATAATCACTTCCGGGCAACAGCACAGTCATTTGGAGATCAGGACTTTCTGCTGATTTCCCCCGGAGAGTTCGGt                                                                                                                                                                                                         >  1:16334/1‑251 (MQ=255)
                                            aGGATTTGCGTCAGTCGTCCTTCGCAGAACACGGCGATGCGGTCGCAGACGGTGATAATTTCAGGTAGTTCAGATGACACCATCAGGATGACTTTTCCGTCGTCCGCCAGTTGGCGCATCACTTTGTAAATTTCGGCTTTCGCGCCAACGTCGATGCCGCGGGTAGGTTCATCGAAAATAATCACTTCCGGGCAACAGCACAGTCATTTGGAGATCAGGACTTTCTGCTGATTTCCCCCGGAGAGTTCGGt                                                                                                                                                                                                         <  1:326872/251‑1 (MQ=255)
                                               aTTTGCGTCAGTCGTCCTTCGCAGAACACGGCGATGCGGTCGCAGACGGTGATAATTTCAGGTAGTTCAGATGACACCATCAGGATGACTTTTCCGTCGTCCGCCAGTTGGCGCATCACTTTGTAAATTTCGGCTTTCGCGCCAACGTCGATGCCGCGGGTAGGTTCATCGAAAATAATCACTTCCGGGCAACAGCACAGTCATTTGGAGATCAGGACTTTCTGCTGATTTCCCCCGGAGAGTTCGGTGat                                                                                                                                                                                                      >  1:371915/1‑251 (MQ=255)
                                                   gCGTCTGTCGTCCTTCGCAGATCACGGCGATGCGGTCTCAGATGTTGATAATTTCATGTAGTTCAGATGACACCATCAGGATGACTTTTCCGTCGCCCGCCAGTTGGCGCATCTCTTTGTAAATTTCGGCTTTCGCGCCAACGTCGATGCCGCGGGTAGGTTCATCGAAAATAATCACTTCCGGGCAACAGCACAGTCATTTGGAGATCAGGACTTTCGGCTGATTTCCCCCGGAGAGTTCGGTGATATTc                                                                                                                                                                                                  <  2:342395/251‑1 (MQ=255)
                                                          ccgtcCTTCGCAGAACACTGCGATGCGGTCGCAGACGGTGATAATTTCAGGTAGTTCAGATGACACCATCAGGATGACTTTTCCGTCGTCCGCCAGTTGGCGCATCACTTTGTAAATTTCGGCTTTCGCGCCAACGTCGATGCCGCGGGTAGGTTCATCGAAAATAATCACTTCCGGGCAACAGCACAGTCATTTGGAGATCAGGACTTTCTGCTGATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTa                                                                                                                                                                                           <  1:14967/250‑1 (MQ=255)
                                                            gtcCTTCGCAGAACACGGCGATGCGGTCGCAGACGGTGATAATTTCAGGTAGTTCAGATGACACCATCAGGATGACTTTTCCGTCGTCCGCCAGTTGGCGCATCACTTTGTAAATTTCGGCTTTCGCGCCAACGTCGATGCCGCGGGTAGGTTCATCGAAAATAATCACTTCCGGGCAACAGCACAGTCATTTGGAGATCAGGACTTTCTGCTGATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTAcc                                                                                                                                                                                         <  1:452035/251‑1 (MQ=255)
                                                               cTTCGCAGAACACGGCGATGCGGTCGCAGACGGTGATAATTTCAGGTAGTTCAGATGACACCATCAGGATGACTTTTCCGTCGTCCGCCAGTTGGCGCATCACTTTGTAAATTTCGGCTTTCGCGCCAACGTCGATGCCGCGGGTAGGTTCATCGAAAATAATCACTTCCGGGCAACAGCACAGTCATTTGGAGATCAGGACTTTCTGCTGATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGa                                                                                                                                                                                       >  1:416377/1‑250 (MQ=255)
                                                                 tCGCAGAACACGGCGATGCGGTCGCAGACGGTGATAATTTCAGGTAGTTCAGATGACACCATCAGGATGACTTTTCCGTCGTCCGCCAGTTGGCGCATCACTTTGTAAATTTCGGCTTTCGCGCCAACGTCGATGCCGCGGGTAGGTTCATCGAAAATAATCACTTCCGGGCAACAGCACAGTCATTTGGAGATCAGGACTTTCTGCTGATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGaa                                                                                                                                                                                      <  1:266135/249‑1 (MQ=255)
                                                                  cGCAGAACTCGGCGATGCTGGCGTAGTGTGTGATAATTTCAGGGAGTTCAGGTGACCCCATCAGGATGACTTTTCCGTCGTCCGCCAGGTGGCGCAGCACTTTGTAAATTTAGGATTTCGCGGCAACTTCGATGCCGCGGGTAGGTTCATCGAAAATAATCACTTCCGGGCAACAGCACAGTCATTTGGAGATCATGACTTTCTGCTGATTTCCTCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGa                                                                                                                                                                                   <  2:27895/251‑1 (MQ=255)
                                                                    cAGAACACGGCGATGCGGTCGCAGACGGTGATAATTTCAGGTAGTTCAGATGACACCATCAGGATGACTTTTCCGTCGTCCGCCAGTTGGCGCATCACTTTGTAAATTTCGGCTTTCGCGCCAACGTCGATGCCGCGGGTAGGTTCATCGAAAATAATCACTTCCGGGCAACAGCACAGTCATTTGGAGATCAGGACTTTCTGCTGATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACa                                                                                                                                                                                 <  1:294020/251‑1 (MQ=255)
                                                                       aaCACGGCGATGCGGTCGCAGACTGTGATAATTTCAGGTAGTTCAGATGACACCATCAGGATGACTTTTCCGTCGTCCGCCAGTTGGCGCATCACTTTGTAAATTTCGGCTTTCGCGCCAACGTCGATGCCGCGGGTAGGTTCATCGAAAATAATCACTTCCGGGCAACAGCACAGTCATTTGGAGATCAGGACTTTCTGCTGATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACAttt                                                                                                                                                                              <  2:372125/251‑1 (MQ=255)
                                                                          acGGCGATGCGGTCGCAGACGGTGATAATTTCAGGTAGTTCAGATGACACCATCAGGATGACTTTTCCGTCGTCCGCCAGTTGGCGCATCACTTTGTAAATTTCGGCTTTCGCGCCAACGTCGATGCCGCGGGTAGGTTCATCGAAAATAATCACTTCCGGGCAACAGCACAGTCATTTGGAGATCAGGACTTTCTGCTGATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTc                                                                                                                                                                             >  1:179615/1‑249 (MQ=255)
                                                                                       cGCAGACGGTGATAATTTCAGGTAGTTCAGATGACACCATCAGGATGGCTTTTCCGTCGTCCGCCAGTTGGCGCATCACTTTGTAAATTTCGGCTTTCGCGCCAACGTCGATGCCGCGGGTAGGTTCATCGAAAATAATCACTTCCGGGCAACAGCACAGTCATTTGGAGATCAGGACTTTCTGCTGATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTc                                                                                                                                                               <  1:284602/250‑1 (MQ=255)
                                                                                               gtgttAATTTCAGGGAGTTCAGATGACACCATCAGGATGACTTTTCCGTCGTCCGCCAGTTGGCGCATCACTTTGTAAATTTCGGCTTTCGCGCCAACGTCGATGCCGCGGGTAGGTTCATCGAAAATAATCACTTCCGGGCAACAGCACAGTCATTTGGAGATCAGGACTTTCTGCTGATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGAtt                                                                                                                                                      <  2:371915/247‑1 (MQ=255)
                                                                                                         cAGGTAGTTCAGATGACACCATCAGGATGGATTTTCCGTCGTCCGCCAGTTGGCGCATCACTTTGTAAATTTCGGCTTTCGCGCCAACGTCGATGCCGCGGGTAGGTTCATCGAAAATAATCACTTCAGGGCAACACAACAGTCATTTGGAGATCAGGACTTTCTGCTGATTTCCCCCGGAGCGTTCGGGGATATTATGGTTTACCGAATGCCATTTCTGCGCCAGCGGTCCGCGTTCATTTTctgaggt                                                                                                                                             >  2:210922/1‑244 (MQ=255)
                                                                                                                                         tttCCGTCGTCCGCCAGTTGGCGCATCACTTTGTAAATTTCGGCTTTCGCGCCAACGTCGATGCCGCGGGTAGGTTCATCGAAAATAATCACTTCCGGGCAACAGCACAGTCATTTGGAGATCAGGACTTTCTGCTGATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTAC‑GTTGCTCGTCAACTTCATGAAACAAGCCCa                                                                                                            <  1:81952/250‑1 (MQ=255)
                                                                                                                                         tttCCGTCGTCCGCCAGGTGGCGCATCCCTTTGTAAATTTCGGCTGTCGCGCCAACGTCGATGCCGCGGGTATGTTCATCGAAAATAATCACTTCCGGGCAACAGCACAGTCATTTGGAGATCAGGACTTTCTGCTGATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTAC‑GTTGCTCGTCAACTTCATGAAACAAGCCCa                                                                                                            <  2:416377/250‑1 (MQ=255)
                                                                                                                                           gccttcgtcCGCCAGGTTGCGCATCACTTTGTAAAGTTCGGCTTTCGCGCTACCGGCGATGCCGCGGGTAGGTTCATCGAAAATAATCACTTCCGGGCAACAGCACAGTCATTTTGAGATCAGGACTTTCTGCTGATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGAGTTTCAGCGGTAC‑GTTGCTCGTCAACTTCATGAATCAAGCCTATc                                                                                                          <  2:425904/246‑1 (MQ=255)
                                                                                                                                                  tcCGCCAGTTGGCGCATCACTTTGTAAATTTCGGCTTTCGCGCCAACGTCGATGCCGCGGGTAGGTTCATCGAAAATAATCACTTCCGGGCAACAGCACAGTCATTTGGAGATCAGGACTTTCTGCTGATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTAC‑GTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCttt                                                                                                  >  1:310681/1‑251 (MQ=255)
                                                                                                                                                    cGCCAGTTGGCGCATCACTTTGTAAATTTCGGCTTTCGCGCCAACGTCGATGCCGCGGGTAGGTTCATCGAAAATAATCACTTCCGGGCAACAGCACAGTCATTTGGAGATCAGGACTTTCTGCTGATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTAC‑GTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTtat                                                                                                >  1:362178/1‑251 (MQ=255)
                                                                                                                                                               gcATCACTTTGTAAATTTCGGCTTTCGCGCCAACGTCGATGCCGCGGGTAGGTTCATCGAAAATAATCACTTCCGGGCAACAGCACAGTCATTTGGAGATCAGGACTTTCTGCTGATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCACCGCCAGCAGTTCGCGTTGATTTTCAGCGGTAC‑GTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCt                                                                                     >  1:391674/1‑251 (MQ=255)
                                                                                                                                                                         gTAAATTTCGGCTTGCGCGCCAGCGTCGATTGCGCGGGTAGGTTCATCGAAAAGAATCACTTCCGGGCAACAGCACAGTCATTTGGAGATCAGGACTTTCTGCTGATTTCCCCCGTAGAGGTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTAC‑GTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGc                                                                           <  2:26516/251‑1 (MQ=255)
                                                                                                                                                                           aaaTTTCGGCTTTCGCGCCAACGTCGATGCCGCGGGTAGGTTCATCGAAAATAATCACTTCCGGGCAACAGCACAGTCATTTGGAGATCAGGACTTTCTGCTGATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTAC‑GTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCgg                                                                         >  1:186400/1‑251 (MQ=255)
                                                                                                                                                                                       tCGCGCCAACGGCGATGCCGCGGGTAGGTTCATCGAAAATAATCACTTCCGGGCAACAGCACAGCCATTTGGAGATCAGGACTTTCTGCTGATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCTGCAGTTCGCGTTGATTTTCAGCGGTAC‑GTGGCTCGTCAACTTCATGAAACAAGCCCATCGCGTCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCAt                                                              <  2:94039/250‑1 (MQ=255)
                                                                                                                                                                                                                      aTCGAAAATAATCACTTCCGGGCAACAGCACAGTCATTTGGAGATCCGGACTTTCTGCTGATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCATCGGTAC‑GTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAAc                              <  2:152017/251‑1 (MQ=255)
                                                                                                                                                                                                                               aatCACTTCCGGGCAACAGCACAGTCATTTGGAGATCAGGACTTTCTGCTGATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTAC‑GTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTTGGGAAAAAACCGTTATccc                     <  2:233200/251‑1 (MQ=255)
                                                                                                                                                                                                                                      tCCGGGCAACAGCACAGTCATTTGGAGATCAGGACTTTCTGCTGATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTAC‑GTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCcggcgg                <  1:257632/249‑1 (MQ=255)
                                                                                                                                                                                                                                         ggacAACAGCACAGTCATTTGGAGATCAGGACTTTCTGCTGGTTTCCCCCGGAGAGTTCGGTGTTATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGAGGTACCGTT‑TTCGTCAACTTACTGAAACAAGCCCATATAGCCTTTATAGTCGCCGACCTTCAGACTGCGTCTGGTCGCCATGATCTGAGCGATGAACCAGTTGGGGACAAACCCTATATCCcggcgggcta            >  2:369740/4‑246 (MQ=255)
                                                                                                                                                                                                                                          ggCAACAGCACAGTCATTTGGATATCAGGACTTTCTGCTGATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTAC‑GTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCa          <  2:16334/251‑1 (MQ=255)
                                                                                                                                                                                                                                          ggCAACAGCACAGCCATTTGGAGATCAGGACTTTCTGCTGATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTAC‑GTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTc           >  1:16619/1‑250 (MQ=255)
                                                                                                                                                                                                                                                   acagcCATTTGGAGATCAGGACTTTCTGCTGATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTAC‑GTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTa  >  1:100660/1‑250 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                      
CCCATGCCATAATCTCCTCTTCGCTCATGTCATCGCGATTCGTCAGGATTTGCGTCAGTCGTCCTTCGCAGAACACGGCGATGCGGTCGCAGACGGTGATAATTTCAGGTAGTTCAGATGACACCATCAGGATGACTTTTCCGTCGTCCGCCAGTTGGCGCATCACTTTGTAAATTTCGGCTTTCGCGCCAACGTCGATGCCGCGGGTAGGTTCATCGAAAATAATCACTTCCGGGCAACAGCACAGCCATTTGGAGATCAGGACTTTCTGCTGATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTAC‑GTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTA  >  NC_000913/4309465‑4309957

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