Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,372,125 T→A 11.9% I150F (ATT→TTT)  nanT ← sialic acid transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,372,1250TA11.9% 162.3 / 4.4 66I150F (ATT→TTT) nanTsialic acid transporter
Reads supporting (aligned to +/- strand):  ref base T (29/27);  new base A (2/6);  total (33/33)
Fisher's exact test for biased strand distribution p-value = 2.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGCTTCCGGGATGTTTTTACGCAGCCAGAGAGCAAAGATGATTGGCAAAATGCCGATAAAGAACAGCGCACGCCAGCCCCAGACCGGAACCACCAGGCTATAGACCT‑GAGCGGCAACGACGGCCCCCACAGAGAAGCCTGAAATCAAAAAACCACTGGCTTTGTTACGCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAGCCTGGCGCAAAGCCGCAGGCC  >  NC_000913/3371981‑3372273
                                                                                                                                                 |                                                                                                                                                    
ccgcatccgGGATGTTTTTACGCAGCCAGAGAGCAAAGATGATTGGCAAAATGCCGATAAAGAACAGCGCACGCCAGCCCCAGACCGGAACCACCAGGCTATAGACCT‑GAGCGACAACGACGGCCCCCACAGAGAAGCCTGCAAACaac                                                                                                                                                  >  2:434250/1‑148 (MQ=255)
    ttccgGGATGTTTTTACGCAGCCAGAGAGCAAAGATGATTGGCAAAATGCCGATAAAGAACAGCGCACGCCAGCCCCAGACCGGAACCAGCAGGCTATAGACCT‑GAGCGGCAACGACGGCCCCCACAGAGAAGCCTGAAATCAAAAAAc                                                                                                                                              <  2:555052/149‑1 (MQ=255)
      ccgGGATGTTTTTACGCAGCCAGAGAGCAAAGATGATTGGCAAAATGCCGATAAAGAACAGCGCACGCCAGCCCCAGACCGGAACCACCAGGTTATAGACCT‑GAGCGGCAACGACGGCCCCCACAGAGAAGCCTGAAATCAAAAAACCa                                                                                                                                            <  1:24336/149‑1 (MQ=255)
      ccgGGATGCTTTTAGACAGCCAAAGAGCAGAGATGATTGGCCAAATGCCGATAAAGAACGGCGCACGCCAGCCCCAGACCGGAACCACCAGGCTATAGACCT‑GAGCGGCAACGCCGGCCCCCACAGAGACGCCCGAAACCAAAAAAccc                                                                                                                                            >  2:69750/1‑148 (MQ=255)
        gggATGTTTTTACGCAGCCAGAGAGCAAAGATGATTGGCAAAATGCCGATAAAGAACAGCGCACGCCAGCCCCAGACCGGAACCACCAGGCTATAGACCT‑GAGCGGCAACGACGGCCCCCACAGAGAAGCCTGAAATCAAAAAACCACt                                                                                                                                          <  1:602134/149‑1 (MQ=255)
            tGTTTTTACGCAGCCAGAGAGCAAAGATGATTGGCAAAATGCCGATAAAGAACAGCGCACGCCAGCCCCAGACCGGAACCACCAGGCTATAGACCT‑GAGCGGCACCGACGGCCCCCACCGCGACGATCTAAATCAAACAACCACTGGCt                                                                                                                                      >  1:15395/1‑149 (MQ=255)
               ttttACGCTGCCAGGGAGCAAAGATGATTGGCAAAATGCCGATAAAGAACAGCGCACGCCAGCCCCAGATCGGAACTACCTGTCTATAGACTT‑GAGCTTCAACGACGGCCCCCACAGATAAGCCTGAAATCAAAAAACCACTGGCTTTg                                                                                                                                   <  2:339034/149‑1 (MQ=255)
                gttACGCAGCCACAGAGCAAAGATAATGAGCAAAAAGCCGATAAAGAACAGCGCACGCCAGCCCCAGACCGGAACCACCAGGCTATACACCT‑GAGCGGCAACGACGGCCCCCACAGAGAAGAAAAAAAAAAAAAAACCACTGGCTTTGt                                                                                                                                  <  2:384060/148‑1 (MQ=255)
                      cAGCCAGAGAGCAAAGATGATTGGCAAAATACCGATAAAGAACAGCGCGCGCCAGCCCCAGACCGGAACCACCAGGCTATAGACCT‑GAGCGGCAACGACGGCCCCCCCAGAGAAGCCTGAAATCAAAAAACCACTGGCTTTGTTACGCa                                                                                                                            >  1:383786/1‑149 (MQ=255)
                           agagaaGAAAGAAGATGGGCAAGAGGCCGATAAAGAAGAGTGCACGCCAGCACCAGACCGGAACCACCCCGCTATACACCT‑GCGCGGCTACGACGGCCCCCACAGAGAAGAAAAAAAAAAAAAAACCACTGGCTTTGTTACGCAGATGt                                                                                                                       <  2:666124/149‑1 (MQ=255)
                             acctcAAAGAGGAGTGGCAAAACGCCGATAAAGAGCAGCGCGCGCCAGCCCCAGACCGGAACCACCAGGCTATAGACCT‑GAGCGGCAACGACGGCCCCCACAGCGAAGCCTGAAATCAAAAAACCACTGGCTTTGTTACGCAGATGttt                                                                                                                     <  2:561115/145‑1 (MQ=255)
                              gagCAAAGATGATTGGCAAAATGCCGATAAAGAACAGCGCACGCCAGCCCCAGACCGGAACCACCAGGCTATAGACCT‑GAGCGGCAACGACGGCCCCCACAGAGAAGCCTGAAATCAAAAAACCACTGGCTTTGTTACGCAGATGtttt                                                                                                                    <  1:379731/149‑1 (MQ=255)
                                  aaaGATGATTGGCAAAATGCCGATAAAGAACAGCGCACGCCAGCCCCAGACCGGAACCACCAGGCTATAGACCT‑GAGCGGCAACGACGGCCCCCACAGAGTGGCCTGAAATCAAAAAACCCCTGGCTTTGTTACGCAGATGTTTTGGcc                                                                                                                >  1:217657/1‑149 (MQ=255)
                                    aGATGATTGGCAAAATGCCGATAAAGAACAGCGCACGCCAGCCCCAGACCGGAACCACCAGGCTATAGACCT‑GAGCGGCAACGACGGCCCCCACAGAGAAGCCTGAAATCAAAAAACCACTGGCTTTGTTACGCAGATGTTTTGACCAg                                                                                                              >  1:272503/1‑149 (MQ=255)
                                       tgatTGGCAAAATGCCGATAAAGAACAGCGCACGCCAGCCCCAGACCGGAACCACCAGGCTATAGACCT‑GAGCGGCAACGACGGCCCCCACAGAGAAGAAAAAAAAAAAAAAACCACTGGCTTTGTTACGCAGATGTTTTGGCCAGCtt                                                                                                           <  2:477890/149‑1 (MQ=255)
                                         atTGGCAATATGCCGCTAAAGAACAGCGCACTCCAGCCCCAGCCCGGAACCACCAGGCTATAGACCT‑GAGCTGAAACGACGGCCCCCACATATAAGCCTGAAATCAAAAAACCCCTGGCTTTGTTACGTAGATGTATTGGCCAGCTTTc                                                                                                         <  2:162719/149‑1 (MQ=255)
                                          ttGGCAAAATGCCGATAAAGAACAGCGCACGCCAGCCCCAGACCGGAACCACCAGGCTATAGACCT‑GAGCGGCAACGACGGCCCCCACAGAGAAGCCTGAAATCAAAAAACCACTGGCTTTGTTACGCAGATGTTTTGGCCAGCTTTCa                                                                                                        <  1:349206/149‑1 (MQ=255)
                                           tGGCAAAATGCCGATAAAGAACAGCGCACGCCAGCCCCAGACCGGAACCACCAGGCTATAGACCT‑GAGCGGCAACGACGGCCCCCACAGAGAAGAATGAAATCAAAAAACCACTGGCTTTGTTACGCAGATGTTTTGGCCAGCTTTCaa                                                                                                       <  2:170429/149‑1 (MQ=255)
                                            ggCAAAATGCCGATAAAGAACAGCGAACGACAGCCCAAGACCGGAACCACCAGGCTATAGAGCTGGAGC‑GCAACGACGGCCCCCACCGAGAAGCCTGAAATCAAAAAACCACTGGCATTGTTACGCAGAAGTTCGGGCGAGCTTCCagt                                                                                                      >  2:629697/1‑147 (MQ=255)
                                                aagtGCCGAAAAAGCACGGCGAACGCCGGCCGCAGACCGTAACACCCCGCCTATAGACTT‑GAGTTGCAACGACGGCCCCCTCCGAGAAGAATGAAAACAAAAAACCACTGGCTTTGTTACGCAGATGTTTTGGCCAGCTTTCAATGACa                                                                                                  <  2:214276/146‑1 (MQ=255)
                                                 aaTGCCGATAAAGAACAGCGCACGCCAGCCCCAGACCGGAACCACCAGGCTATAGACCT‑GAGCGGCAAGGACGGCCCCCACAGAGAAGCCTGAAATCAAAAAACCACTGGCTTTGTTACGCAGATGTTTTGGCCAGCTTGCAATGGCag                                                                                                 >  1:626458/1‑148 (MQ=255)
                                                   tGCCGATAAAGAACAGCGCACGCCAGCCCCAGACCGGAACCACCAGGCTATAGACCT‑GAGCGGCAACGACGGCCCCCACAGAGAAGCCTGAAATCAAAAAACCACTGGCTTTGTTACGCAGATGTTTTGGCCAGCTTTCAATGACATAg                                                                                               <  1:578721/149‑1 (MQ=255)
                                                       gATAAAGAACAGCGCACGCCAGCCCCAGACCGGAACCACCAGGCTATAGACCT‑GAGCGGCAACGACGGCCCCCACAGAGAAGCCTGAAATCAAAAAACCACTGGCTTTGTTACGCAGATGTTTTGGCCAGCTTTCAATGACATAGGTgg                                                                                           >  2:703309/1‑149 (MQ=255)
                                                             gAACAGCGCACGCCAGCCCCAGACCGGAACCACCAGGCTATAGACCT‑GAGCGGCAACGACGGCCCCCACAGAGAAGCCTGAAATCAAAAAACCACTGGCTTTGTTACGCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTgg                                                                                     >  1:371600/1‑149 (MQ=255)
                                                              aaCAGCGCACGCCAGCCCCAGACCGGAACCACCAGGCTATAGACCT‑GAGCGGCAACGACGGCCCCCGCAAAGAAGCCTGAACACAAAAAACCACTGGCTTTGTTACGCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCTCTGGa                                                                                    >  1:175956/1‑149 (MQ=255)
                                                              aaCAGCGCACGCCAGCCCCAGACCGGAACCACCAGGCTATAGACCT‑GAGCGGCAACGACGGCCCCCACAGAGAAGCCTGAAATCAAAAAACCACTGGCATTGTTACGCAGATGTATTGGACAGCTTTCCATGACATAGGTGGCGCGGGa                                                                                    >  2:478788/1‑149 (MQ=255)
                                                                  gcgcACGCCAGCCCCAGACCGGAACCACCAGGCTATAGCCCT‑GAGCGGCAACGACGGCCCCCCCAGAGAAGCCTGAAATCAAAAAACCACTGGCGTTGATAGGCAGATGGGTAGGCCAGCTTTCAAAGCCCCAGGAGGCGCTGGGACCg                                                                                >  2:609848/1‑149 (MQ=255)
                                                                     cACGCCAGCCCCAGACCGGCACCACCAGGCCATAGCCCT‑GAGCGGCAACGACGGCCCCCACAGAGAAAAAAAAAAAAAAAAAACCACTGGCTTTGTTACGCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTAt                                                                             <  2:454753/149‑1 (MQ=255)
                                                                        gCCAGCCCCAGACCGGAACCACCAGGCTATAGACCT‑GAGCGGCAACGACAGCCCCCACAGAGAAGCCTGAAATCAAAAAACCACTGGCTTTGTTACGCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCAGGAACCGCATTcc                                                                          >  1:539057/1‑148 (MQ=255)
                                                                         ccAGCCCCAGACCGGAACCACCAGGCTATAGACCT‑GAGCGGCAACGACGGCCCCCACAGAGAAGCCTGAAATCAAAAAACCACTGGCTTTGTTACGCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCAc                                                                         <  1:652663/149‑1 (MQ=255)
                                                                         ccAGACCCCGACCGGCACCACCAGGCTATAGACCT‑GAGCGGCAACGACGGCCCCCACGGAGAAGAAAAAAATAAAAAAACCACTGGCTTTGTTACGCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCAc                                                                         <  2:692585/149‑1 (MQ=255)
                                                                          cAGCCCCAGACCGGAACCACCAGGCTATAGACCT‑GAGCGGCAACGACGGCCCCCACAGAGAAGCCTGAAATCAAAAAACCACTGGCTTTGTTACGCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTcccc                                                                        >  1:353336/1‑149 (MQ=255)
                                                                           aGCCCCAGACCGGAACCACCAGGCTATCGACCT‑GAGCGGCAACGACGGCCCCCACAGAGAAGCCTGAAATCAAAAAACCACTGGCTTTGTTACGCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCAccc                                                                       <  1:524593/149‑1 (MQ=255)
                                                                           aGCCCCAGACCGGAACCACCAGGCTATAGACCT‑GAGCGGCAACGACGGCCCCCACAGAGACGCCGGAAATCAAAAAACCACTGGCCTCGTTACCCAGATGCTTTGGCCAGCTTCCAATGACAGCGGTGGGGGGGGAACCGCATTCAccc                                                                       >  2:551709/1‑149 (MQ=255)
                                                                                cAGACCGGAACCACCAGGCTATAGACCT‑GAGCGGCAACGACGGCCCCCACAGAGAAGCCTGAAATCAAAAAACCACTGGCTTCGTCACGCAGACGCTTTGGCCAGCTGTCAATGACATAGGTGGCGCTGGAACCGTAATCACGCGCCAt                                                                  >  2:611731/1‑149 (MQ=255)
                                                                                 aGACCGTAACCACCAGGCTATAGACCT‑GAGCGGCAACGACGGCCCCCACAGAGAAGCCTGAAATCAAAAAACCACTGGCTTTCCTACGCACCTGTTTTGGCCAGCTTTCAATGACATACGTGGCGCTGGAACCGTATTCACCCGCCATc                                                                 <  1:388891/149‑1 (MQ=255)
                                                                                 aGACCGGAACCACCAGGCTATAGACCT‑TAGCGGCAACGACGGCCCCCACAGAGAAGCCTGAAATCAAAAAACCACTGGCTTTGTTACGCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATc                                                                 <  1:699169/149‑1 (MQ=255)
                                                                                     cGGAACCACCAGGCTATAGACCT‑GAGCGGCAACGACGGCCCCCACAGAGAAGCCTGAAATCAAAAAACCACTGGCTTTGTTACGCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCa                                                             <  2:15395/149‑1 (MQ=255)
                                                                                       gAAGCACCCGGCTATAGACCT‑GAGCGGCAACGACGGCCCCCACAGAGAAGCCTGAAATCAAAAAACCACTGGCTTTGTTACGCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATg                                                           <  1:418527/149‑1 (MQ=255)
                                                                                       gAACCACCAGGCTATAGACCT‑GAGTGGCAACGACGGCCCCCACAGAGAAGCCTGAAATCAAAAAACCACTGGCTTTGTTACGCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCCTg                                                           >  1:221624/1‑149 (MQ=255)
                                                                                           caccaGGCTATAGACCT‑GAGCGGCAACGACGGCCCCCACAGAGAAGCCTGAAATCAAAAAACCACTGGCTTTGTTACGCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTGTTCACCCGCCATCCCCATGCCGa                                                       >  1:396328/1‑149 (MQ=255)
                                                                                           caccaGGCTATAGACCT‑GAGCGGCAACGACGGCCCCCACAGAGAAGCCTGAAATCAAAAAACCACTGGCTTTGTTACGCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGGCGa                                                       >  2:695430/1‑149 (MQ=255)
                                                                                              caGGCTATAGACCT‑GAGCGGCAACGACGGCCCCCACAGCGAAGCCTGAAATCAAAAAACCCCTGGCTCTGGTACGCCGACGTACTCGCCAGCTTTCAGTGACATAGGTGGCGCTGGAAGCGTATGCACCCGCGCTCCCCATGCCGaggc                                                    >  2:817/1‑146 (MQ=255)
                                                                                              caGGCTATAGACCT‑GAGCGGCAACGACGGCCCCCACAGAGGAGCCTGAAATCAAAAAACCACCGGCCACGGTACGCAGACGCCTTGGCCAGCCTTCAATGACATCGGTGGCGCTGGAACCGGATTCACCCGCCGTCCCCATCCCGATgc                                                    >  2:704878/1‑148 (MQ=255)
                                                                                              caGGCTATAGACCT‑GAGCGGCAACGACGGCCCCCACAGAGAAGCCTGAAATCAAAAAACCACTGGCTTTGTTACGCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGa                                                    <  2:394916/149‑1 (MQ=255)
                                                                                              caGGCTATAGACCT‑GAGCGGCAACGACGGCCCCCACAGAGAAGCCTGAAATCAAAAAACCACTGGCTTTGTTACGCAGATGTTTTGGCCAGCCTTTAAGCACATATGTTGCGCTGGCACCGTATTCACCCGCCATCCCCATGCCGATGa                                                    >  1:328395/1‑149 (MQ=255)
                                                                                               aGGCTATAGACCT‑GAGCGGCAACGACGGCCCCCACAGAGAAGCAAAAAATCAAAAAACCACTGGCTTTGTTACGCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGAc                                                   <  2:577845/149‑1 (MQ=255)
                                                                                                   tataGACCT‑GAGCGGCAACGACGGCCCCCACAGAGAAGCCTGAAATCAAAAAACCACTGGCTTTGTTACGCAGATGTTTTGGCCAGCTTTCAATGTCATAGGTGGCGCTGGAACCGTATTCACCCTCCATCCCCATGCCGATGACCAGa                                               <  1:412910/149‑1 (MQ=255)
                                                                                                      aGACCT‑GAGCGGCAACGACGGCCCCCACAGAGAAGCCTGAAATCAAAAAACCACTGGCTTTGTTACGCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGa                                            <  1:422530/149‑1 (MQ=255)
                                                                                                              aGCGGCAACGACGGCCCCCACAGAGAAGCCTGAAAACAAAAAACCACTTGCTTTGTTACGCAGATGTTTTGGCCAGCTTTCAATGATATAGGTGGCGCTGGAACCGTATCCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATaa                                     <  2:460041/149‑1 (MQ=255)
                                                                                                                cGGCAACGACGGCCCCCACAGAGAAGCCTGAAATCAAAAAACCACTGGCTTTGTTACGCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAAc                                   <  1:432831/149‑1 (MQ=255)
                                                                                                                  gCAACGACGGCCCCCACAGAGAAGCCAGAAATCAAAAACCCCCTGGCATTGTTACGCGGAGGCTGGGACCGGATTCCAACGACATCCGCGGGGCGGGGACCATACGAGCCAGCCAGCCCCAGGCCGAGGCCCAGACGAGGGaggaaggc                                 >  2:376558/1‑142 (MQ=255)
                                                                                                                   cAACGACGGCCCCCACAGAGAAGCCTGAAATCAAAAAACCACTGGCTTTGTTACGCAGATGTTTTGGCCAGCTTTCAATGCCATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAAccgg                                >  1:204236/1‑146 (MQ=255)
                                                                                                                   cAACGACGGCCCCCACAGAGAAGCCTGAAAGCAAAAAACCACTGGCTTTGTTACGCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGCCGAGCGATAAACATg                                >  1:182354/1‑149 (MQ=255)
                                                                                                                       gacgGCCCCCACAGAGAAGGCTGAAATCAAAAAACCACTGGCTTTGATACGCAGATGTTTTGGCGAGAACTCAATGACCTAGGTGGCGCTGGAACCGGACTCACCCGGCATCCCAATGCCGATGACCAGACGAGCGAAAAAcccggggc                            >  2:165684/1‑142 (MQ=255)
                                                                                                                          ggCCCCCACAGAGAAGCCTGAAATCAAAAAACCACTGGCTTTGTTACGCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGt                         >  1:643467/1‑149 (MQ=255)
                                                                                                                           gCCCCCACAGAGAAGCCTGAAATCAAAAAACCACTGGCTTTGTTACGCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTa                        <  1:645248/149‑1 (MQ=255)
                                                                                                                            cccccACAGAGAAGCCTGAAATCAAAAAACCACTGGCTTTGTTACGCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAg                       <  1:629697/149‑1 (MQ=255)
                                                                                                                                cacaGAGAAGCCTGAAATCAAAAAACCACTGGCTTTGTTACGCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGCGATGTAGCCTg                   >  2:238951/1‑149 (MQ=255)
                                                                                                                                  caGAGAAGCCTGAAATCAAAAAACCACTGGCTTTGTTACGCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAGCCTGgc                 <  2:92296/149‑1 (MQ=255)
                                                                                                                                  caGAGAAGCCTGAAATCAAAAAACCACTGGCTTTGTTACGCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAGCCTGgc                 <  1:71017/149‑1 (MQ=255)
                                                                                                                                            tGAAATCAAAAAACCACTGGCTTTGTTACGCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGGTAGACATGGTGATGTAGCCTGGCGCGAcgccgc       >  2:123150/1‑149 (MQ=255)
                                                                                                                                            tGAAATCAAAAAACCACTGGCTTTGTTACGCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAGCCTGGCGCAAAGCCGc       >  1:576257/1‑149 (MQ=255)
                                                                                                                                            tGAAATCAAAAAACCACTGGCTTTGTTACGCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAGCCTGGCGCAAAGCCGc       >  1:404467/1‑149 (MQ=255)
                                                                                                                                              aaaTCAAAAAACCACTGGCTTTGTTACGCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAGCCTGGCGCAAAGCCGCAg     >  2:477873/1‑149 (MQ=255)
                                                                                                                                               aaTCAAAAAACCACGTGCTTTGTTTCCAGGCTGTTTTGGCCAGCTTTCAAGTGCACAGGTGGCGCTGTAACCGTATTCACCCGCCAACCCCATGCCGATGACCAGACGAGCGATAATCATGGTGATGTAGCCTGGCGCAAAGCCGCAgg    <  2:228693/149‑1 (MQ=255)
                                                                                                                                                 tCAAAAAACCACTGGCTTTGTTACGCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAGCCTGGCGCAAAGCCGCAGGcc  >  1:681027/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
CCGCTTCCGGGATGTTTTTACGCAGCCAGAGAGCAAAGATGATTGGCAAAATGCCGATAAAGAACAGCGCACGCCAGCCCCAGACCGGAACCACCAGGCTATAGACCT‑GAGCGGCAACGACGGCCCCCACAGAGAAGCCTGAAATCAAAAAACCACTGGCTTTGTTACGCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAGCCTGGCGCAAAGCCGCAGGCC  >  NC_000913/3371981‑3372273

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