Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,529,039 C→T G402D (GGT→GAT)  sgcC ← putative PTS enzyme IIC component SgcC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,529,0390CT97.6% 126.8 / ‑6.2 41G402D (GGT→GAT) sgcCputative PTS enzyme IIC component SgcC
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (28/12);  minor base G (0/1);  total (28/13)
Fisher's exact test for biased strand distribution p-value = 3.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.39e-01

CTTTGTCGATAACCCAGTTCATCCCCAACTGCGCATCAAAACTCGGGTCGCCGGGCAATGCACGTAGATGGCACATGGCGATGACCGCTTTTTCCGTCCCAATAACCTCTTTCAGCCAACTCATTGAAAGTACTCCTTATTCCGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAACCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACGTCCTGACCAGGTCTCCGCGGTGGATAACCGTAGCCATACAGATAAAGAACGGTGCTACCGGCAGGTCTGCGAGGGGTAACACTTTGTTGC  >  NC_000913/4528794‑4529279
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                
cTTTGTCGATAACCCAGTTCATCCCCAACTGCGCATCAAAACTCGGGTCGCCGGGCAATGCACGTAGATGGCACATGGCGATGACCGCTTTTTCCGTCCCAATAACCTCTTTCAGCCAACTCATTGAAAGTACTCCTTATTCCGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGaa                                                                                                                                                                                                                                              >  1:212389/1‑250 (MQ=255)
  ttGTCGATAACCCAGTTCATCCCCAACTGCGCATCAAAACTCGGGTCGCCGGGCAATGCACGTAGATGGCACATGGCGATGACCGCTTTTTCCGTCCCAATAACCTCTTTCAGCCAACTCATTGAAAGTACTACTTATTCCGATAATTCCCGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCATGACAACGCCAATAATCACAAGCGACATCAATTCAGAAATCGACCAATCGAACa                                                                                                                                                                                                                                            >  1:372164/1‑250 (MQ=255)
           ccccaattcccccccAACTGCCCATCAAAACCCGGGTCGCCGTGAAATGCACGTAGAGGGCACAGGGCGATGACCGCTTTTTCCATCCCAATAACCTCTTCCAGCCAATCTATTGAAAGTACCCCTTATACCGATATTCCCCGTTTTCGAAAGAAAAGCACGCCGCTGGCCAAAATCCCAACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTTCCCAc                                                                                                                                                                                                                                  <  2:185368/240‑1 (MQ=255)
             ccAGTCCACCCCCAACTGCGCAACAAAACTCGGGCCGCCGGCCAATGCCCGTAGTGGGAACAGGGGGATGCCCGCTTTTTCCGCCCCAATACCCCCTTTCAGCCAACCCATTGAAATTCCTCTTTATTCCGATAATTCCCTTTGCCGAAAGAACAGCACGCCGCTGGAGACAATACAGACCGCCACGACAACGCCAATAACCCCAAGCGACATCAATTCAGAAATCGACCAAGCGAACATATGGCCCACgg                                                                                                                                                                                                                                <  2:70992/251‑1 (MQ=255)
                     cccccAGCTGCGCATCAAAACTCGGGTCGGCGGGCAATGCGCGTAGAGGGTACATGGCGATGACCTTTTTTTCCGTCCCAATAACCTCTTTCAGCCACCTCATTGAAAGTCCCCCTTATTCCGAGAATTCTCGTTTACGTAAGAACAGCAGGACGCGGGCGACAAGACCAACCGCCACGACAACGCCAATAATCCCAAGCGACAGCAATTCAGAAATCGACCAATCGAACATAGGGCCCACTGACAACGCt                                                                                                                                                                                                                        <  2:383532/250‑1 (MQ=255)
                                    cAAAACTCGGGTCGCCGGGCAATGCACGTAGATGGCACATGGCGATGACCGCTTTTTCCGTCCCAATAACCTCTTTCAGCCAACTCATTGAAAGAACTCCTTATTCCGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGAGCACtt                                                                                                                                                                                                         >  1:248191/1‑251 (MQ=255)
                                      aaactcggtgttgtgggaaaTGCACGTAGATTGCTAATGGTGTTCACCGTTTTTTCCGCCCCAATAACCTCTTCCATCCAACCCATGAAAAGTACTCCTTATCCCGATAATTACCGTTTAAGTAAGACCTGCACGCCGCTGGGGCAAATCCCTACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTc                                                                                                                                                                                                       <  1:99468/233‑1 (MQ=255)
                                       aaCTCGGGTCGCCGGGCAATGCACGTAGATGGCACATGGCGATGACCGCTTTTTCCGTCCCAATAACCTCTTTCAGCCAACTCATTGAAAGTACTCCTTATTCCGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCa                                                                                                                                                                                                      <  1:180675/251‑1 (MQ=255)
                                                       cAATGCACGTAGATGGCACATGGCGATGACCGCTTTTTCCGTCCCAATAACCTCTTTCAGCCAACTCATTGAAAGTACTCCTTATTCCGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAAc                                                                                                                                                                                      >  1:400248/1‑251 (MQ=255)
                                                                            gccGATGACCGTTTTTTCCGTCCAAATAACCTCTTTCAGCAAACTCATTAAAAGTACTCCTTATCCCGAAAATCCCCGTTACCGAAAGAACACCCCGCCGCTGGCAAAAAACCCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAg                                                                                                                                                                 <  2:121909/249‑1 (MQ=255)
                                                                                aTGACCGCTTTTTCCGTCCCAATAACCTCTTTCAGCCAACTCATTGAAAGTACTCCTTATTCCGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTgc                                                                                                                                                              >  1:322002/1‑249 (MQ=255)
                                                                                   aCCGCTTTTTCGGTCCCAATAACCTCTTTAAGCCAACTCATTGAAAGTACTCCTTCTTCCGATAATTCCAGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGacgt                                                                                                                                                           >  2:21560/1‑247 (MQ=255)
                                                                                      cgtttttCCCGCCCAAAAACCCTCTTTCAGCCAACACATTGAAAGTACTCCTTATTCCGATAATTCACGTTAACGTAAGAACAGCACGACGCTGGCGACAATCCCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAg                                                                                                                                                        <  1:320289/249‑1 (MQ=255)
                                                                                       cTTTTTCCGTCCCAATAACCTCTTTCAGCCAACTCATTGAAAGTACTCCTTATTCCGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGgc                                                                                                                                                      >  2:360802/1‑251 (MQ=255)
                                                                                         ttttCCGTCCCAATAACCTCTTTCAGCCAACTCATTGAAAGTACTCCTTATTCCGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACtt                                                                                                                                                                                                         >  2:325206/1‑198 (MQ=255)
                                                                                         ttttCCGTCCCAATAACCTCTTTCAGCCAACTCATTGAAAGTACTCCTTATTCCGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACtt                                                                                                                                                                                                         <  1:325206/198‑1 (MQ=255)
                                                                                                   cAATAACCTCTTTCAGCCAACTCATTGAAAGTACTCCTTATTCCGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTg                                                                                                                                          >  1:292383/1‑251 (MQ=255)
                                                                                                            tctTTCAGCCAACTCATTGAAAGTACTCCTTATTCCGATAATACACGTTTACGTAAGAACAGCACGACGCTGACGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGAGACCTCAATTCAAACATCGACCCATCGAACATATTGCCCACGAACAACGCTGAGCTTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAACGCGCGAACTGGGTGGCGATc                                                                                                                                    >  1:228733/1‑249 (MQ=255)
                                                                                                              tttaagacAACTCATTGAAAGTACTCCTTATTACGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACACTCCCGACCGCCGCGACAACGCCAATAATCCCAAGCGCCATCAATTAAGAAATCGACCAATCGACCCTATTGCCCACGGACATCGCTGTGCTTTGCGCCCTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGCAGTAAGGCGCGAACTGCGTGGCGATcagc                                                                                                                                 >  2:256108/8‑250 (MQ=255)
                                                                                                               ttCAGCCAACTCATGAAAAGTACTCCTTATCCCAATAATTCCCGTTTCCGAAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATcagc                                                                                                                                 <  1:48376/249‑1 (MQ=255)
                                                                                                                     cAACTCATTGAAAGTACTCCTTATTCCGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAAAGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGt                                                                                                                         >  2:152956/1‑251 (MQ=255)
                                                                                                                        cTCATTGAAAGTACTCCTTATTCCGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAat                                                                                                                      >  2:57042/1‑251 (MQ=255)
                                                                                                                              agAAAGTACTCTTTTCCCGAATAATTCACTTTCACGAAAGAACAGCCCGCGGTGGGCAAAAATCCCGCCCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGATATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAatcatc                                                                                                                  <  2:138354/249‑1 (MQ=255)
                                                                                                                                             ttCCGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGCAATCCACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACAGCCTTTTAGCGCCCTTTCCGTGAAGTAAGGCGCGAACTGGGTGGCGCTCAGCAGTACGGTAATCATCACGATGACCCCGCtt                                                                                                   >  1:49790/1‑250 (MQ=255)
                                                                                                                                                cGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCAACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGAGAACTGGGAGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATc                                                                                                >  1:80876/1‑250 (MQ=255)
                                                                                                                                                 gATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGAAATCATCACGATGACCCCGCTTATCa                                                                                               >  1:346747/1‑250 (MQ=255)
                                                                                                                                                           cGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACGTCCTGAc                                                                                     >  1:183088/1‑251 (MQ=255)
                                                                                                                                                                           gCACGACGCTGTCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATAGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCAGCAAAACTTAAACCGCCTTTGAGAGCCATTTCAGTGACGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGAAATCATCACGATGACCTCGCTGATCACCGTCATACCCAGGTCTCCGCggggg                                                                     >  2:271957/1‑248 (MQ=255)
                                                                                                                                                                               gacgCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCACTTCAGACATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGGGAAGTAAGGCGCGAACTGGGCGGCGATCAGCAGGACGATAATCATCACGATGACCACGATTATCAACGTCCTGACCAGGTCTCCGCGCTGGATa                                                                  >  2:278594/1‑250 (MQ=255)
                                                                                                                                                                                    tagCGAAAATACCGCCCGCCACGACACCGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTCTGAGCGCCATTTCAGTGAAGTAAGTCGCGAACTGGGGGGCGATCAACAGGACGGTAATCATCACGATGACCCCGCTTATCAACGTCATGACCAGGTCTCCGCGGTGGAGAACCGtt                                                            >  2:369516/3‑250 (MQ=255)
                                                                                                                                                                                                             cAACGCCAATCATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACGTCCTGACCAGGTCTCCGCGGTGGATAACCGTAGCCATACAGATAAAGAAAGGTGCt                                    >  2:156272/1‑251 (MQ=255)
                                                                                                                                                                                                                  ccAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCCGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACGTCCTGACCAGGTCTCCGCGGTGGATAACCGTAGCCATACAGATAAAGAACGGTGCTACCgg                               >  1:35551/1‑251 (MQ=255)
                                                                                                                                                                                                                    aataatCCCAAGCGACTTCCATGCAGAAATCGCCCAATCAAACATATTGCCCACGTACAACGCTGAGATTTGCGCACTTCCCCCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACGTCCTGACCAGGTCTCCGCGGTGGATAACCGTAGCCATACAGATAAAGAACGGTGCCACCGGCa                             <  2:33417/251‑1 (MQ=255)
                                                                                                                                                                                                                          cccAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACGTCCTGACCAGGTCTCCGCGGTGGATAACCGTAGCCATACAGATAAAGAACGGTGCTACCGGCAGGTCTg                       <  1:264097/251‑1 (MQ=255)
                                                                                                                                                                                                                                           cAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAAACGACTTTGAGCGCCATTTCAGTGAAGAAAGGCGCGAACTGGGAGGCGATCAGCAGGACGGTAATCATCACGATGCCCCCGCTTATCAACGTCCTGACCAGGTCTCCGCGGTGGATAACCGTAGCCATACAGATAAAGAACGGTGCTACCGGCAGGTCTGCGAGGGGTAACACTttg      >  1:243353/1‑251 (MQ=255)
                                                                                                                                                                                                                                           cAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGAATTGCGCACTTTCACCGGCAAACCTGAACCCGCCTTTGAGCGCCATTTCAGTGACGTAAGGCGCGAACTGGGTGGCGATCAGCAGGCCGGTAATCATCACGATGACACCGCTTATCAACGTCCTGACCAGGTCTCCGCGGTGGATAACCGTAGCCATACAGCTAAAGAACGGCGCTACCGGCAGGACTCCGAGGGGAAACACtttt      >  1:170801/1‑250 (MQ=255)
                                                                                                                                                                                                                                           cAGAAATCGACCAATCGAACATATTGCCCACGAACAACGCTGAGATTTGAGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACGTCCTGACCAGGTCTCCGCGGTGGATAACCGTAGCCATACAGATAAAGAACGGTGCTACCGGCAGGTCTGCGAGGGGTaaaaacttg      >  1:242493/1‑244 (MQ=255)
                                                                                                                                                                                                                                            aGAAATCGACCAATCGAACATATTGCCCACGGACAACGCGGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACGTCCTGACCAGGTCTCCGCGGTGGATAACCGTAGCCATACAGATAAAGAACGGTGCTACCGGCAGGTCTGCGAGGGGTAACACTttgt     <  1:136855/251‑1 (MQ=255)
                                                                                                                                                                                                                                             gAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACGTCCTGACCAGGTCTCCGCGGTGGATAACCGTAGCCATACAGATAAAGAACGGTGCTACCGGCAGGTCTGCGAGGGGTAACACTttgtt    >  1:70400/1‑251 (MQ=255)
                                                                                                                                                                                                                                              aaaTCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACGTCCTGACCAGGTCTCCGCGGTGGATAACCGTAGCCATACAGATAAAGAACGGTGCTACCGGCAGGTCTGCGAGGGGTAACACTttgttg   >  1:166952/1‑251 (MQ=255)
                                                                                                                                                                                                                                               aaTCGACCAATCGAACATATTGCCCACGGCCAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGCACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACGTCCTGACCAGGTCTCCGCGGTGGATAACCGTAGCCATACAGATAAAGAACGGTGCTACCGGCAGGGCTGCGAGGGGTAAAACTTTTTTGc  >  2:182086/1‑251 (MQ=255)
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                
CTTTGTCGATAACCCAGTTCATCCCCAACTGCGCATCAAAACTCGGGTCGCCGGGCAATGCACGTAGATGGCACATGGCGATGACCGCTTTTTCCGTCCCAATAACCTCTTTCAGCCAACTCATTGAAAGTACTCCTTATTCCGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAACCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACGTCCTGACCAGGTCTCCGCGGTGGATAACCGTAGCCATACAGATAAAGAACGGTGCTACCGGCAGGTCTGCGAGGGGTAACACTTTGTTGC  >  NC_000913/4528794‑4529279

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