Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,537,490 A→G D35D (GAT→GAC 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,4900AG98.7% 253.8 / ‑5.1 77D35D (GAT→GACnanSN‑acetyl‑9‑O‑acetylneuraminate esterase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (40/36);  total (41/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

GGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCATCTTCCCTGTCCGGTAATGGCAGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAAT  >  NC_000913/4537357‑4537623
                                                                                                                                     |                                                                                                                                     
ggTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCGTCTTcc                                                                                                                                 <  1:368029/140‑1 (MQ=255)
     tAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCGTCTTCCCTGTc                                                                                                                            >  2:729267/1‑140 (MQ=255)
      aaCCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCGTCTTcc                                                                                                                                 >  1:617444/1‑134 (MQ=255)
      aaCCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCGTCTTcc                                                                                                                                 <  2:617444/134‑1 (MQ=255)
       aCCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCGTCTTCCCTGTCCg                                                                                                                          >  1:335584/1‑140 (MQ=255)
         ccTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCGTCTTCCCTGTCCGGt                                                                                                                        <  1:923404/140‑1 (MQ=255)
         ccTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCGTCTTCCCTGTCCGGt                                                                                                                        >  1:967697/1‑140 (MQ=255)
           tGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCGTCTTCCCTGTCCGGTaa                                                                                                                      >  2:672725/1‑140 (MQ=255)
            gCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCGTCTTCCCTGTCCGGTAAt                                                                                                                     >  2:575868/1‑140 (MQ=255)
              atatCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCGTCTTCCCTGTCCGGTAATgg                                                                                                                   >  2:926823/1‑140 (MQ=255)
                  ccTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCGTCTTCCCTGTCCGGTAATGGCAGt                                                                                                               >  1:519341/1‑140 (MQ=255)
                   cTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCGTCTTCCCTGTCCGGTAATGGCAGTc                                                                                                              >  1:872936/1‑140 (MQ=255)
                      aaCATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCGTCTTCCCTGTCCGGTAATGGCAGTCCtt                                                                                                           <  2:90133/140‑1 (MQ=255)
                       aCATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCGTCTTCCCTGTCCGGTAATGGCAGTCCTTc                                                                                                          >  2:667142/1‑140 (MQ=255)
                          tCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCGTCTTCCCTGTCCGGTAATGGCAGTCCTTCGcc                                                                                                       >  2:541314/1‑140 (MQ=255)
                          tCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTAAGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCGTCTTCCCTGTCCGGTAATGGCAGTCCTTCGcc                                                                                                       >  2:291172/1‑140 (MQ=255)
                           cGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCGTCTTCCCTGTCCGGTAATGGCAGTCCtt                                                                                                           <  2:158297/135‑1 (MQ=255)
                           cGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCGTCTTCCCTGTCCGGTAATGGCAGTCCtt                                                                                                           >  1:158297/1‑135 (MQ=255)
                           cGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCGTCTTCCCTGTCCGGTAATGGCAGTCCTTCGCCa                                                                                                      >  2:48118/1‑140 (MQ=255)
                           cGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCGTCTTCCCTGTCCGGTAATGGCAGTCCTTCGCCa                                                                                                      >  1:890453/1‑140 (MQ=255)
                            gtgtGGGCAGTGAGTCAGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCGTCTTCCCTGTCCGGTAATGGCAGTCCTTCGCCAt                                                                                                     <  2:335584/140‑1 (MQ=255)
                            gtgtGGGCAGTGAGTCAGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCGTCTTCCCTGTCCGGTAATGGCAGTCCTTCGCCAt                                                                                                     >  1:154949/1‑140 (MQ=255)
                             tgtgGGCAGTGAGTCAGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCGTCTTCCCTGTCCGGTAATGGCAGTCCTTCGCCATa                                                                                                    >  1:512392/1‑140 (MQ=255)
                                 ggCAGTGAGTCAGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCGTCTTCCCTGTCCGGTAATGGCAGTCCTTCGCCATACGcc                                                                                                >  2:80776/1‑140 (MQ=255)
                                  gCAGTGAGTCAGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCGTCTTCCCTGTCCGGTAATGGCAGTCCTTCGCCATACGCCa                                                                                               <  1:575868/140‑1 (MQ=255)
                                     gTGAGTCAGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCGTCTTCCCTGTCCGGTAATGGCAGTCCTTCGCCATACGCCa                                                                                               <  1:679666/137‑1 (MQ=255)
                                     gTGAGTCAGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCGTCTTCCCTGTCCGGTAATGGCAGTCCTTCGCCATACGCCa                                                                                               >  2:679666/1‑137 (MQ=255)
                                      tGAGTCAGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCGTCTTCCCTGTCCGGTAATGGCAGTCCTTCGCCATACGCCATGGc                                                                                           >  1:836391/1‑140 (MQ=255)
                                      tGAGTCAGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCGTCTTCCCTGTCCGGTAATGGCAGTCCTTCGCCATACGCCATGGc                                                                                           >  1:394953/1‑140 (MQ=255)
                                       gAGTCAGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCATCTTCCCTGTCCGGTAATGGCAGTCCTTCGCCATACGCCATGGCa                                                                                          >  1:373939/1‑140 (MQ=255)
                                            aGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCGTCTTCCCTGTCCGGTAATGGCAGTCCTTCGCCATACGCCATGGCATTAGa                                                                                     >  1:807286/1‑140 (MQ=255)
                                            aGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCGTCTTCCCTGTCCGGTAATGGCAGTCCTTCGCCATACGCCATGGCATTAGa                                                                                     <  2:646739/140‑1 (MQ=255)
                                            aGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCGTCTTCCCTGTCCGGTAATGGCAGTCCTTCGCCATACGCCATGGCATTAGa                                                                                     <  1:353320/140‑1 (MQ=255)
                                                gAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCGTCTTCCCTGTCCGGTAATGGCAGTCCTTCGCCATACGCCATGGCATTAGACTGa                                                                                 >  1:887033/1‑140 (MQ=255)
                                                gAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCGTCTTCCCTGTCCGGTAATGGCAGTCCTTCGCCATACGCCATGGCATTAGACTGa                                                                                 <  2:154949/140‑1 (MQ=255)
                                                  ataatGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCGTCTTCCCTGTCCGGTAATGGCAGTCCTTCGCCATACGCCATGGCATTAGACTGAcc                                                                               >  1:826752/1‑140 (MQ=255)
                                                   taatGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCGTCTTCCCTGTCCGGTAATGGCAGTCCTTCGCCATACGCCATGGCATTAGACTGACCa                                                                              <  1:48118/140‑1 (MQ=255)
                                                      tGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCGTCTTCCCTGTCCGGTAATGGCAGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCa                                                                           <  1:491606/140‑1 (MQ=255)
                                                      tGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCGTCTTCCCTGTCCGGTAATGGCAGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCa                                                                           <  2:967697/140‑1 (MQ=255)
                                                        tCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCGTCTTCCCTGTCCGGTAATGGCAGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCa                                                                           >  2:677995/1‑138 (MQ=255)
                                                        tCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCGTCTTCCCTGTCCGGTAATGGCAGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCa                                                                           <  1:677995/138‑1 (MQ=255)
                                                        tCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCGTCTTCCCTGTCCGGTAATGGCAGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAAc                                                                         >  2:864789/1‑140 (MQ=255)
                                                        tCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCGTCTTCCCTGTCCGGTAATGGCAGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAAc                                                                         <  2:519341/140‑1 (MQ=255)
                                                        tCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCGTCTTCCCTGTCCGGTAATGGCAGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAAc                                                                         >  2:545875/1‑140 (MQ=255)
                                                        tCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCGTCTTCCCTGTCCGGTAATGGCAGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAAc                                                                         >  1:567063/1‑140 (MQ=255)
                                                           ttAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCGTCTTCCCTGTCCGGTAATGGCAGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACggg                                                                      >  2:710887/1‑139 (MQ=255)
                                                           ttAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCGTCTTCCCTGTCCGGTAATGGCAGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGt                                                                      >  2:307355/1‑140 (MQ=255)
                                                           ttAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCGTCTTCCCTGTCCGGTAATGGCAGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGt                                                                      <  1:307355/140‑1 (MQ=255)
                                                             aaaGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCGTCTTCCCTGTCCGGTAATGGCAGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTaa                                                                    >  2:373181/1‑140 (MQ=255)
                                                             aaaGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCGTCTTCCCTGTCCGGTAATGGCAGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTaa                                                                    <  2:807286/140‑1 (MQ=255)
                                                              aaGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCGTCTTCCCTGTCCGGTAATGGCAGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAg                                                                   <  1:541314/140‑1 (MQ=255)
                                                                gTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCGTCTTCCCTGTCCGGTAATGGCAGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGaa                                                                 <  2:909062/140‑1 (MQ=255)
                                                                 tGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCGTCTTCCCTGTCCGGTAATGGCAGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAAc                                                                <  1:672725/140‑1 (MQ=255)
                                                                 tGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCGTCTTCCCTGTCCGGTAATGGCAGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAAc                                                                <  2:198871/140‑1 (MQ=255)
                                                                 tGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCGTCTTCCCTGTCCGGTAATGGCAGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAAc                                                                <  1:864789/140‑1 (MQ=255)
                                                                 tGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCGTCTTCCCTGTCCGGTAATGGCAGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAAc                                                                <  2:394953/140‑1 (MQ=255)
                                                                    cATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCGTCTTCCCTGTCCGGTAATGGCAGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCa                                                                           <  2:127068/126‑1 (MQ=255)
                                                                    cATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCGTCTTCCCTGTCCGGTAATGGCAGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCa                                                                           >  1:127068/1‑126 (MQ=255)
                                                                    cATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCGTCTTCCCTGTCCGGTAATGGCAGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGt                                                                      <  1:165146/131‑1 (MQ=255)
                                                                    cATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGAAGAGGCGCGTCTTCCCTGTCCGGTAATGGCAGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGt                                                                      >  2:165146/1‑131 (MQ=255)
                                                                        gCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCGTCTTCCCTGTCCGGTAATGGCAGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGtaa                                                         >  2:760978/1‑140 (MQ=255)
                                                                                ccGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCGTCTTCCCTGTCCGGTAATGGCAGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCggg                                                 <  1:926823/140‑1 (MQ=255)
                                                                                ccGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCGTCTTCCCTGTCCGGTAATGGCAGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCggg                                                 <  1:729267/140‑1 (MQ=255)
                                                                                ccGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCGTCTTCCCTGTCCGGTAATGGCAGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCggg                                                 <  2:402928/140‑1 (MQ=255)
                                                                                   ggATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCGTCTTCCCTGTCCGGTAATGGCAGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGa                                              >  2:466543/1‑140 (MQ=255)
                                                                                   ggATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCGTCCTCCCTGTCCGGTAATGGCAGTCCTTCGCCATACACCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGa                                              >  2:73761/1‑140 (MQ=255)
                                                                                      tGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCGTCTTCCCTGTCCGGTAATGGCAGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCg                                                   >  1:565740/1‑132 (MQ=255)
                                                                                      tGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCGTCTTCCCTGTCCGGTAATGGCAGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCg                                                   <  2:565740/132‑1 (MQ=255)
                                                                                            tGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCGTCTTCCCTGTCCGGTAATGGCAGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCggg                                                 <  1:493616/128‑1 (MQ=255)
                                                                                            tGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCGTCTTCCCTGTCCGGTAATGGCAGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCggg                                                 >  2:493616/1‑128 (MQ=255)
                                                                                             gcgcAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCGTCTTCCCTGTCCGGTAATGGCAGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCg                                    <  2:744862/140‑1 (MQ=255)
                                                                                             gcgcAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCGTCTTCCCTGTCCGGTAATGGCAGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCg                                    <  2:872936/140‑1 (MQ=255)
                                                                                             gcgcAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCGTCTTCCCTGTCCGGTAATGGCAGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCg                                    <  2:512392/140‑1 (MQ=255)
                                                                                             gcgcAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCGTCTTCCCTGTCCGGTAATGGCAGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCg                                    <  1:466543/140‑1 (MQ=255)
                                                                                                                               aGGCGCGTCTTCCCTGTCCGGTAATGGCAGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAAt  >  2:883360/1‑140 (MQ=255)
                                                                                                                                 gcgcGTCTTCCCTGTCCGGTAATGGCAGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGaa           >  2:572920/1‑129 (MQ=255)
                                                                                                                                 gcgcGTCTTCCCTGTCCGGTAATGGCAGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGaa           <  1:572920/129‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                     
GGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCATCTTCCCTGTCCGGTAATGGCAGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAAT  >  NC_000913/4537357‑4537623

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: