Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,537,366 C→T G77S (GGT→AGT)  nanS ← N‑acetyl‑9‑O‑acetylneuraminate esterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,537,3660CT97.6% 136.5 / ‑3.5 42G77S (GGT→AGT) nanSN‑acetyl‑9‑O‑acetylneuraminate esterase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (20/21);  total (20/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 = 3.77e-01

GCTGACTAAATCCTGGTACAGCGGAGTATCCGTTCCCCAACGACAAGCATCATGGCTGGCTCCGTGCCGTTCTGAATATGTCCCTTCGCTGCCCGCGGTAAAAGCCGATCCGCCACGGCAACACGGAACGATAAGAATCCCTGCATTATCAGGAATGAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGG‑TGCCGTACTGTGTTTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCATCTTCCCTGTCCGGTAATGGCAGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCG  >  NC_000913/4537121‑4537574
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                
gctgACTAAATCCTGGTACAGCGGAGTATCCGTTCCCCAACGACAAGCATCATGGCTGGCTCCGTGCCGTTCTGAATATGTCCCTTCGCTGCCCGCGGTAAAAGCCGATCCGCCACGGCAACACGGAACGATAAGAATCCCTGCATTATCAGGAATGAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGG‑TGCCGTACTGTGTTTGATGATTCGTTGCCAGAGGATGGTGATAACTCTGc                                                                                                                                                                                                              >  1:308691/1‑250 (MQ=255)
  tgACTAAATCCTGGTACAGCGGAGTATCCGTTCCCCAACGACAAGCATCATGGCTGGCTCCGTGCCGTTCTGAATATGTCCCTTCGCTGCCCGCGGTAAAAGCCGATCCGCCACGGCAACACGGAACGATAAGAATCCCTGCATTATCAGGAATGAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGG‑TGCCGTACTGTGTTTGATGATTCGTTGCCAGAGGATGGTGATAACTCTGCat                                                                                                                                                                                                            >  1:4646/1‑250 (MQ=255)
           cgtGGTACAGAGGACTATCCGTTCCCCAGCGACCAGAATCCTGGCTGACTCCGAGCCGTTCTGAATATGTCCCTTCGCTGCCCGCGGTAAAAGCAGATCCCCCACGGCAACACGGAACGATAAGAATCCCTGAATTATCAGGAATGAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGG‑TGCCGTACTGTGTTTGATGATTCGTTGCCAGAGGATGGTGATAACTCTGCATATCCTGAAc                                                                                                                                                                                                   >  2:241789/3‑250 (MQ=255)
           ccTGGTACAGCGGAGTATCCGTTCCCCAACGACAAGCATCATGGCTGGCTCCGTGCCGTTCTGAATATGTCCCTTCGCTGCCCGCGGTAAAAGCCGATCCGCCACGGCAACACGGAACGATAAGAATCCCTGCATTATCAGGAATGAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGG‑TGCCGTACTGTGTTTGATGATTCGTTGCCAGAGGATGGTGATAACTCTGCATATCCTGAAc                                                                                                                                                                                                   <  2:228297/250‑1 (MQ=255)
             tGGTAAAGCGGAGTATACGTTCCCCAACGACAAGCATCATGGCTGGCTCCGTGCAGTTCTGAATATGTCCCTTCGCTGCCCGCGGTAAAAGCCGATCCGCCACGGCAACACGGAACGATACGAATCCCTGCATTATCAGGAATGACGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGG‑TGCCGTACTGTGTTTGATGATTCGTTGCCAGAGGATGGTGATAACTCTGCATATCCTGAACATc                                                                                                                                                                                                >  2:62270/1‑251 (MQ=255)
                      ggAGTATCCGTTCCCCAACGACAAGCATCATGGCTGGCTCCGTGCCGTTCTGAATATGTCCCTTCGCTGCCCGCGGTAAAAGCCGATCCGCCACGGCAACACGGAACGATAAGAATCCCTGCATTATCAGGAATGAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGG‑TGCCGTACTGTGTTTGATGATTCGTTGCCAGAGGATGGTGATAACTCTGCATATCCTGAACATCGTGTGGGc                                                                                                                                                                                        >  2:386341/1‑250 (MQ=255)
                           aTCCGTTCCCCAACGACAAGCATCATGGCTGGCTCCGTGCCGTTCTGAATATGTCCCTTCGCTGCCCGCGGTAAAAGCCGATCCGCCACGGCAACACGGAACGATAAGAATCCCTGCATTATCAGGAATGAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGG‑TGCCGTACTGTGTTTGATGATTCTTTGCCAGAGGATGGTGATAACTCTGCATATCCTGAACATCGTGTGGGCAGTGa                                                                                                                                                                                   <  1:71993/250‑1 (MQ=255)
                                                atcatGGCTGGCTCCGTGCCGTTCTGAATATGTCCCTTCGCTGCCCGCGGTAAAAGCCGATCCGCCACGGCAACACGGAACGATAAGAATCCCTGCATTATCAGGAATGAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGG‑TGCCGTACTGTGTTTGATGATTCGTTGCCAGAGGATGGTGATAACTCTGCATATCCTGAACATCGTGTGGGCAGGGAGTCAGTGGAATAATGTCGTTa                                                                                                                                                              >  2:15847/1‑250 (MQ=255)
                                                   atGGCTGGCTCCGTGCCGTTCTGAATATGTCCCTTCGCTGCCCGCGGTAAAAGCCGATCCGCCACGGCAACACGGAACGATAAGAATCCCTGCATTATCAGGAATGAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGG‑TGCCGTACTGTGTTTGATGATTCGTTGCCAGAGGATGGTGATAACTCTGCATATCCTGAACATCGTGTGGGCCGGGAGTCAGTGGAATAATGTCGTTAAAGt                                                                                                                                                          >  1:368488/1‑251 (MQ=255)
                                                         ggctCCGTGCCGTTCTGAATATGTCCCTTCGCTGCCCGCGGTAAAAGCCGATCCGCCACGGCAACACGGAACGATAAGAATCCCTGCATTATCAGGAATGAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGG‑TGCCGTACTGTGTTTGATGATTCGTTGCCAGAGGATGGTGATAACTCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTAAAGTGACATg                                                                                                                                                    <  1:365902/251‑1 (MQ=255)
                                                                   cGTTCTGAATATGTCCCTTCGCTGCCCGCGGTAAAAGCCGATCCGCCACGGCAACACGGAACGATAAGAATCCCTGCATTATCAGGAATGAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGG‑TGCCGTACTGTGTTTGATGATTCGTTGCCAGAGGATGGTGATAACTCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTcc                                                                                                                                          <  2:143595/251‑1 (MQ=255)
                                                                   cGTTCTGAATATGTCCCTTCGCTGCCCGCGGTAAAAGCCGATCCGCCACGGCAACACGGAACGATAAGAATCCCTGCATTATCAGGAATGAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGG‑TGCCGTACTGGTTTTGATGATTCGTTGCCAGAGGATGGTGATAACTCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTAAAGTGACATGTCGGGCCTcc                                                                                                                                          <  2:173216/251‑1 (MQ=255)
                                                                     ttCTGAATATGTCCCTTCGCTGCCCGCGGTAAAAGCCGATCCGCCACGGCAACACGGAACGATAAGAATCCCTGCATTATCAGGAATGAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGG‑TGCCGTACTGTGTTTGATGATTCGTTGCCAGAGGATGGTGATAACTCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTCCgg                                                                                                                                        >  1:120333/1‑251 (MQ=255)
                                                                         gAATATGTCCCTTCGCTGCCCGCGGTAAAAGCCGATCCGCCACGGCAACACGGAACGATAAGAATCCCTGCATTATCAGGAATGAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGG‑TGCCGTACTGTGTTTGATGATTCGTTGCCAGAGGATGGTGATAACTCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGATg                                                                                                                                    <  2:141067/251‑1 (MQ=255)
                                                                              gggcccTTCGCTGCCCGCGGGAAAAGCCTCTCCGCCACGGCAACACGGAACGATAAGAATCCCTGCATTATCAGGAATGAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGG‑TGCCGTACTGTGTTTGATGATTCGTTGCCAGAGGATGGTGATAACTCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTAAAGTGACATGGCGTGCCTCCGGGATGCGTAt                                                                                                                               <  2:212095/248‑1 (MQ=255)
                                                                                                    aaaGCCGATCCGCCACGGCAACACGGAACGATAAGAATCCCTGCATTATCAGGAATGAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGG‑TGCCGTACTGTGTTTGATGATTCGTTGCCAGAGGATGGTGATAACTCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAAGGTCGTTAAAGTTACAGGGCTGGCCTCCGGGATGCGCATGCGCAAATCTCGCTTATTGttt                                                                                                         <  2:16528/251‑1 (MQ=255)
                                                                                                     aaGCCGATCCGTCACCGCAACACGGAACGATAAGAATCCCTGCATTATCAGGAATGAAGGGCAGTAATTTCCGTGCAATACGCAGTGCCTTGCCAACGG‑TGCCGTACTGTGTTTGATGATTCGTTGCCAGAGGATGGTGATAACTCTGCATATCCTGAACATCGTGTGTTCAGTGAGTCAGTGGAATAATGTCGTTAAAGTGACATGTCGGGCCTCCTGGATGCCTACGCGCAAATCGCGTTAATTGTTTa                                                                                                        <  2:253909/251‑1 (MQ=255)
                                                                                                        ccGATCCGCCACGGCAACACGGAACGATAAGAATCCCTGCATTATCAGGAATGAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGG‑TGCCGTACTGTGTTTGATGATTCGTTGCCAGAGGATGGTGATAACTCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAAtt                                                                                                     <  1:270744/251‑1 (MQ=255)
                                                                                                            tCCGCCACGGCAACACGGAACGATAAGAATCCCTGCATTATCAGGAATGAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGG‑TGCCGTACTGTGTTTGATGATTCGTTGCCAGAGGATGGTGATAACTCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTgg                                                                                                 >  1:207676/1‑251 (MQ=255)
                                                                                                                   cGGCAACACGGAACGATAAGAATCCCTGCATTATCAGGAATGAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGG‑TGCCGTACTGTGTTTGATGATTCGTTGCCAGAGGATGGTGATAACTCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAgg                                                                                          <  2:120333/251‑1 (MQ=255)
                                                                                                                      cAACACGGAACGATAAGAATCACTGCATTATCAGGAATGAAGGGCAGTAATTTCCGTGCCATATGCAGTGCCTGGCCAACGG‑TGCCGTACTGTGTTTGATGATTCGTTGCCAGAGGATGGTGATAACTCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCCTCGGGATGCGTATGCGCAAATCTCGCTCATTGTTTAATTCTGGGATGAGgcgc                                                                                       >  2:20651/1‑251 (MQ=255)
                                                                                                                        acacGGAACGATAAGAATCCCTGCATTATCAGGAATGAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGG‑TGCCGTACTGTGTTTGATGATTAGTTGCCAGAGGACGGAGATAACTCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGCATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGAGCAt                                                                                     >  1:255582/1‑251 (MQ=255)
                                                                                                                          acGGAACGATAAGAATCCCTGCATTATCAGGAATGAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGG‑TGCCGTACTGTGTTTGATGATTCGTTGCCAGAGGATGGTGATAACTCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCATCt                                                                                   >  1:330079/1‑251 (MQ=255)
                                                                                                                            ggAACGATAAGAATCCCTGCATTATCAGGAATGAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGG‑TGCCGTACTGTGTTTGATGATTCGTTGCCAGAGGATGGTGATAACTCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCATCTTc                                                                                 >  1:32650/1‑251 (MQ=255)
                                                                                                                               aCGATAAGAATCCCTGCATTATCAGGAATGAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGG‑TGCCGTACTGTGTTTGATGATTCGTTGCCAGAGGATGGTGATAACTCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATTAGGCGCATCTTCCCt                                                                              >  1:110324/1‑251 (MQ=255)
                                                                                                                                            ccgCATTATCAGGAATGAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGG‑TGCCGTACTGTGTTTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCGGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCATCTTCCCTGTCCGGTAATGGc                                                                 <  1:412480/249‑1 (MQ=255)
                                                                                                                                              gCATTATCAGGAATGAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGG‑TGCCGTACTGTGTTTGATGATTCGTTGCCAGAGGATGGTGATAACTCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCATCTTCCCTTTCCGGTAATGGCAg                                                               <  1:108495/251‑1 (MQ=255)
                                                                                                                                                   aTCAGGAATGAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGG‑TGCCGTACTGTGTTTGATGATTCGTTGCCAGAGGATGGTGATAACTCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCATCTTCCCTGTCCGGTAATGGCAGTCCt                                                           <  2:165529/250‑1 (MQ=255)
                                                                                                                                                        gAATGAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGG‑TGCCGTACTGTGTTTGATGATTCGTTGCCAGAGGATGGTGATAACTCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCATCTTCCCTGTCCGGTAATGGCCATCCTTCGccc                                                     >  2:142440/1‑250 (MQ=255)
                                                                                                                                                          aTGAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGG‑TGCCGTACTGTGTTTGATGATTCGTTGCCAGAGGATGGTGATAACTCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCATCTTCCCTGTCCGGTAATGGCAGTCCTTCGCCAt                                                    >  1:69532/1‑250 (MQ=255)
                                                                                                                                                          aTGAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGG‑TGCCGTACTGTGTTTGATGATTCGTTGCCAGAGGATGGAGATAACTCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGAAAATCTCGCTAATTGTTTAATTCTGGGAGGAGGCGCATCTTCCCTGTCCGGTAATGGCAGTCCTTCGCCAt                                                    >  2:15757/1‑250 (MQ=255)
                                                                                                                                                            gAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGG‑TGCCGTACTGTGTTTGATGATTCGTTGCCAGAGGATGGTGATAACTCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCATCTTCCCTGTCCGGTAATGGCAGTCCTTCGCCATACg                                                 <  1:39534/251‑1 (MQ=255)
                                                                                                                                                                  cAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGG‑TGCCGTACTGTGTTTGATGATTCGTTGCCAGAGGATGGTGATAACTCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCATCTTCCCTGTCCGGTAATGGCAGTCCTTCGCCCTACGCCATgg                                           >  1:408060/1‑251 (MQ=255)
                                                                                                                                                                          tCNGTGCAATATGCAGTGCCTGGCCAACGG‑TGCCGTACTGTGTTTGATGATTCGTTGCCAGAGGATGGTGATAACTCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCATCTTCCCTGTCCGGTAATGGCAGGCCTCCGCCGTACGCCATGGCATTAGa                                    <  2:368488/250‑1 (MQ=255)
                                                                                                                                                                                  atatTCAGTTCCTGGCCAACGGGTGC‑GTACTGTGTTTGATGATTCGTTGCCAGATGATGGTGATAACTCTGCATATCCTGAACATCGTGTGGGCACTGAGTCAGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTTATTCTGGGATGAGGCGCATCTTCCCTGTCCGTTAATGGCAGGCCTTCGCCAGACGCCATGGCATTAGACTGTCCAGc                           <  2:268565/251‑1 (MQ=255)
                                                                                                                                                                                           gCCTGGCCAACGG‑TGCCGTACTGTGTTTGATTATTCGTTGCCAGAGGATGGTGATAACTCTGCATATCCTGAACATCGTGTGTGCAGTGAGTCAGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATTCGCAAATCTCGCTAATGGTTTAAGTCTGGGATGAGGCGCATCTTCCCTGTCCGGTAATGGCAGGCCTTCGCCATACGCCAGGGCATTAGACTGTCTAGCAACGGtgag                  <  2:232974/251‑4 (MQ=255)
                                                                                                                                                                                            ccTGGCCAACGG‑TGCCGTACTGTGTTTGATGATTCGTTGCCAGAGGATGGTGATAACTCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCATCTTCCCTGTCCGGTAATGGCAGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGa                 <  2:207676/251‑1 (MQ=255)
                                                                                                                                                                                              tGGCCAACGG‑TGCCGTACTGTGTTTGATGATTCGTTGCCAGAGGATGGTGATAACTCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCATCTTCCCTGTCCGGTAATGGCAGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAAc               >  1:217735/1‑251 (MQ=255)
                                                                                                                                                                                                    aCGG‑TGCCGTACTGTGTTTGATGATTCGTTGCCAGAGGATGGTGATAACTCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCATCTTCCCTGTCCGGTa                                                                      >  1:5632/1‑190 (MQ=255)
                                                                                                                                                                                                    aCGG‑TGCCGTACTGTGTTTGATGATTCGTTGCCAGAGGATGGTGATAACTCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCATCTTCCCTGTCCGGTa                                                                      <  2:5632/190‑1 (MQ=255)
                                                                                                                                                                                                           ccGTACTGTGTTTGATGATTCGTTGCCAGAGGATGGTGATAACTCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCATCTTCCCTGTCCGGTAATGGCAGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATc   <  1:386341/251‑1 (MQ=255)
                                                                                                                                                                                                            cGTACTGTGTTTGATGATTCGTTGCCAGAGGATGGTGATAACTCTGCATATCCTGAACATCGTGTGGGCAGTTAGTCAGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCATCTTCCCTGTCCGGTAATGGCAGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCg  <  2:269477/251‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                
GCTGACTAAATCCTGGTACAGCGGAGTATCCGTTCCCCAACGACAAGCATCATGGCTGGCTCCGTGCCGTTCTGAATATGTCCCTTCGCTGCCCGCGGTAAAAGCCGATCCGCCACGGCAACACGGAACGATAAGAATCCCTGCATTATCAGGAATGAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGG‑TGCCGTACTGTGTTTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCATCTTCCCTGTCCGGTAATGGCAGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCG  >  NC_000913/4537121‑4537574

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