Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,537,030 A→G 7.6% Y189H (TAC→CAC)  nanS ← 9‑O‑acetyl N‑acetylneuraminic acid esterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,537,0300AG 7.6% 142.8 / 2.8 53Y189H (TAC→CAC) nanS9‑O‑acetyl N‑acetylneuraminic acid esterase
Reads supporting (aligned to +/- strand):  ref base A (28/20);  new base G (2/2);  total (31/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 = 9.23e-01

ACCGTAATATCCCGTGCTTAAATCGTCCGGGTCTTCATCAGGCGCGTTCGTCAGTCCTCTTTCACCTTGTTGCTGGAAGTCGACGAAAATAATATTGGCTAAAACATTATTTTGATAATTGCCATAAATAGCTTCATACGAATGCGGGAAATTTTCTTTCCAGTACCAGGTGGTATCGCCGCAAAACCACGGCGCGTCAGTAATATTATTAAGCTGAGAATGGTATTGTTTTAGATCCCTACGAAAGGCTTCAACCATATGATTAAAGTGTTGAGGGTGTGACGCGTAGTCACTGGTCATTAAGTCAAATTCGCCTTGCATCCAGCATGCGCCGAGGAATTTGTTCTGCGGATTTTTTGCCAGTGCGGCTCGCGTTCTGCTGACTAAATCCTGGTACAGCGGAGTATCCGTTCCCCAACGACAAGCATCATGGCTGGCTCCGTGCCGTTCTGAATATGTCCCTTCGCTGCCCGCGGTAAAAGCCGATCCGCCACGGCAACACGGAACGATAAGAATCCCTGCATTATCAGGAATGAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTG  >  NC_000913/4536743‑4537330
                                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                                                            
aCCGTAATATCCCGTGCTTAAATCGTCCGGGTCTTCATCAGGCGCGTTCGTCAGTCCTCTTTCACCTTGTTGCTGGAAGTCGACGAAAATAATATTGGCTAAAACATTATTTTGATAATTGCCATAAATAGCTTCATACGAATGCGGGAAATTTTCTTTCCAGTACCAGGTGGTATCGCCGCAAAACCACGGCGCGTCAGTAATATTATTAAGCTGAGAATGGTATTGTTTTAGATCCCTACGAAAGGCTTCAACCATATGATTAAAGTGTTGAGGGTGTGACGCGTAGTCACTGGTCAtt                                                                                                                                                                                                                                                                                                 <  2:631878/301‑1 (MQ=255)
 ccGTAATATCCCGTGCTTAAATCGTCCGGGTCTTCATCCGGCGCGTTCGTCAGTCCTCTTTCACCTTGTTGCTGGAAGTCGACGAAAATAATATTTGCTAAAACATTATTTTGATAATTGCCATAAATAGCTTCATACGAATGCGGGAAATTTTCTTTCCAGTACCAGGTGGTATCGCCGCAAAACCACGGCGCGTCAGTAATATTATTAAGCTGAGAATGGTATTGTTTTAGATCCCTACGAAAGGCTTCAACCATATGATTAAAGTGTTGAGGGTGTGACGCGTAGTCACTGGTCATTa                                                                                                                                                                                                                                                                                                <  2:231827/301‑1 (MQ=255)
   gTAATATGCCGTGTTTTATTTGTCCGGGTCTTCATCAGGCGCGTTCGTCAGTCCTCTTTCACCTTGTTGCTGGAAGTCGACGAAAATAATATTGGCTAAAACATTATTTTGATAATTGCCATAAATAGCTTCATACGAATGCGGGAAA‑TTTCTTTCCAGTACCAGGTGGTATCGCCGCAAAACCACGGCGCGTCAGTAATATTATTAAGCTGAGAATGGTATTGTTTTAGATCCCTACGAAAGGCTTCAACCATATGATTAAAGTGTTGAGGGTGTGACGCGTAGTCACTGGTCATTAAGt                                                                                                                                                                                                                                                                                             <  1:277257/301‑1 (MQ=255)
   gTAATATCCCGTGCTTAAATCGTCCGGGTCTTCATCAGGCGCGTTCGTCAGTCCTCTTTCACCTTGTTGCTGGAAGTCGACGAAAATAATATTGGCTAAAACATTATTTTGATAATTGCCATAAATAGCTTCATACGAATGCGGGAAATTTTCTTTCCAGTACCAGGTGGTATCGCCGCAAAACCACGGCGCGTCAGTAATACTATTAAGCTGAGAATGGTATTGTTTTAGATCCCTACGAAAGGCTTCAACCATATGATTAAAGTGTTGAGGGTGTGACGCGTAGTCACTGGTCATTaaa                                                                                                                                                                                                                                                                                              >  2:549878/1‑300 (MQ=255)
            cGTGCTTAAATCGTCCGGGTCTTCATCAGGCGCGTTCGTCAGTCCTCTTTCACCTTGTTGCTGGAAGTCGACGAAAATAATATTGGCTAAAACATTATTTTGATAATTGCCATAAATAGCTTCATACGAATGCGGGAAATTTTCTTTCCAGTACCAGGTGGTATCGCCGCAAAACCACGGCGCGTCAGTACTATTATTAAGCTGAGAATGGTATTGTTTTAGATCCCTACGACAGGCTTCAACCATATGATTAAAGTGTTCAGGATATGACGCGAAGTCACTGGTCATTCAGTCACAATCg                                                                                                                                                                                                                                                                                     >  1:608698/1‑301 (MQ=255)
               ttctAAATCGTCCGGGTCTTCATCAGGCGCGTTCGTCAGTCCTCTTTCACCTTGTTGCTGGAAGTCGACGAAAATAATATTGGCTAAAACATTATTTTGATAATTGCCATAAATAGCTTCATACGAATGCGGGAAATTTTCTTTCCAGTACCAGGTGGTATCGCCGCAAAACCACGGCGCGTCAGTAATATTATTAAGCTGAGAATGGTATTGTTTTAGATCCCTACGAAAGGCTTCAACCATATGATTAAAGTGTTGAGGGTGTGACGCGTAGTCACTGGTCATTAAGTCAAATTCGCCt                                                                                                                                                                                                                                                                                  <  2:332937/298‑1 (MQ=255)
                cTTAAATCGTCCGGGTCTTCATCAGGCGCGTTCGTCAGTCCTCTTTCACCTTGTTGCTGGAAGTCGACGAAAATAATATTGGCTAAAACATTATTTTGATAATTGCCATAAATAGCTTCATACGAATGCGGGACATTTTCTTTCCAGTACCAGGTGGTATCGCCGCAAAACCACGGCGCGTCAGTAATATTATTAAGCTGAGAATGGTCTTGTTTTAGATCCCTACGAAAGGCTTCAACCATATGATTAAAGTGTTGAGGGTGTGACACGTAGTCACTGGTCATTAAGTAAAATACGCCtt                                                                                                                                                                                                                                                                                 >  1:188143/1‑301 (MQ=255)
                   taaTCGTCCGGGTCTTCATCAGGCGCGTTCGTCAGTCCTCTTTCACCTTGTTGCTGGAAGTCGACGAAAATAATATTGGCTAAAACATTATTTTGATAATTGCCATAAATAGCTTCATACGAATGCGGGAAATTTTCTTTCCAGTACCAGGTGGTATCGCCGCAAAACCACGGCGCGTCAGTAATATTATTAAGCTGAGAATGGTATTGTTTTAGATCCCTACGAAAGGCTTCAACCATATGATTAAAGTGTTGAGGGTGTGACGCGTAGTCACTGGTCATTAAGTCAAATTCGCCTTGCa                                                                                                                                                                                                                                                                              <  2:623044/300‑1 (MQ=255)
                   aaaTCGTCCGGGTCTTCATCAGGCGCGTTCGTCAGTCCTCTTTCACCTTGTTGCTGGAAGTCGACGAAAATAATATTGGCTAAAACATTATTTTGATAATTGCCATAAATAGCTTCATACGAATGCGGGAAATTTTCTTTCCAGTACCAGGTGGTATCGCCGCAAAACCACGGCGCGTCAGTAATATTATTAAGCTGAGAATGGTATTGTTTTAGATCCCTACGAAAGGCTTCAACCATATGATTAAAGTGTTGAGGGAGTGAAGCATAGGCAAGGGACATACAGTAAAATTCGCCTTCCa                                                                                                                                                                                                                                                                              >  1:337797/1‑301 (MQ=255)
                   aaaTCGTCCGGGTCTTCATCAGGCGCGTTCGTCAGTCCTCTTTCACCTTGTTGCTGGAAGTCGACGAAAATAATATTGGCTAAAACATTATTTTGATAATTGCCATAAATAGCTTAATACGAATGCGGGAAATTTTCTTTCCAGTACCAGGTGGTATCGCCGCAAACCCACGGCGCGTCAGTAATATTATTAAGCTGCGAATGGTATTGTTTTAGATCCCTACGAAAGGCTTCAACCATATGATTAACGTGTTGAGGGTTTGACGCGTATTCACTGGTCTTTAAGTCACATTCGCCtatat                                                                                                                                                                                                                                                                              >  1:332532/1‑297 (MQ=255)
                    aaTCGTCCGGGTCTTCATCAGGCGCGTTCGTCAGTCCTCTTTCACCTTGTTGCTGGAAGTCGACGAAAATAATATTGGCTAAAACATTATTTTGATAATTGCCATAAATAGCTTCATACGAATGCGGGAAATTTTCTTTCCAGTACCAGGTGGTATCGCCGCAAAACCACGGCGCGTCAGTAATATTATTAAGCTGAGAATGGTATTGTTTTAGATCCCTACGAAAGGCTTCAACCATATGCTTAAAGTGTTGAGGGTGTGACGCGTAGTCACTGGTCATTAAGTCAAATTCGCCTTGCaa                                                                                                                                                                                                                                                                             >  1:622185/1‑300 (MQ=255)
                                              ttcttcAGTCTGCTTTGTCCTTGTTGCTGGTAGTTGATGTAAATAATTTTGGCTAAAACATTATTTTGATAATTGCCATAAATAGCTTCATACGAATGCGTGAAAGTGTCTTTCCAGTACCAGGTTGTATCGCCGCAAAACCACGGCGCGTCAGTAATATTATTAAGCTGAGAATGGTAGTGTGTTAGATCCCTACGAGAGGCTTCAACCATATGATTAAAGTGTTGAGGGTGTGACGCGTAGTCACTGGTCATTAAGTCAAATTCGCCTTGCAGCCAGCATGCGCCGAGGAATTTGGTCt                                                                                                                                                                                                                                                   <  1:519921/301‑1 (MQ=255)
                                                   cAGTCCTCTTTCACCTTGTTGCTGGAAGTCGACGAAAATAATATTGGCTAAAACATTATTTTGATAATTGCCATAAATAGCTTCATACGAATGCGGGAAATTTTCTTTCCAGTACCAGGTGGTATCGCCGCAAAACCACGGCGCGTCAGTAATATTATTAAGCTGAGAATGGTATTGTTTTAGATCCCTACGAAAGGCTTCAACCATATGATTAAAGTGTTGAGGGTGTGACGCGTAGTCACTGGTCATTAAGTCAAATTCGCCTTGCATCCAGCATGCGCCGAGGAATTTGTTCTGCGGa                                                                                                                                                                                                                                              >  2:126047/1‑301 (MQ=255)
                                                              taCCTTGTTGCTGGAAGTCGACGAAAATAATATTGGCTAAAACATTATTTTGATAATTGCCATAAATAGCTTCATACGAATGCGGGAAATTTTCTTTCCAGTACCAGGTGGTATCGCCGCAAAACCACGGCGCGTCAGTAATATTATTAAGCTGAGAATGGTATTGTTTTAGATCCCTACGAAAGGCTTCAACCATATGATTAAAGTGTTGAGGGTGTGACGCGTAGTCACTGGTCATTAAGTCAAATTCGCCTTGCATCCAGCATGCGCCGAGGAATTTGTTCTGCGGATTTTTTGCCAg                                                                                                                                                                                                                                   <  2:54229/300‑1 (MQ=255)
                                                                ccTTGTTGCTGGAAGTCGACGAAAATAATATTGGCTAAAACATTATTTTGATAATTGCCATAAATAGCTTCATACGAATGCGGGAAATTTTCTTTCCAGTACCAGGTGGTATCGCCGCAAAACCACGGCGCGTCAGTAATATTATTAAGCTGAGAATGGTATTGTTTTAGATCCCTACGAAAGGCTTCAACCATATGATTAAAGTGTTGAGGGTGTGACGCGTAGTCACTGGTCATTAAGTCAAATTCGCCTTGCATCCAGCATGCGCCGAGGAATTTGTTCTGCGGATTTTTTGCCAGTg                                                                                                                                                                                                                                 >  2:151208/1‑301 (MQ=255)
                                                                              gTCGACGAAAATAATATTGGCTAAAACATTATTTTGATAATTGCCATAAATAGCTTCATACGAATGCGGGAAATTTTCTTTCCAGTACCAGGTGGTATCGCCGCAAAACCACGGCGCGTCAGTAATATTATTAAGCTGAGAATGGTATTGTTTTAGATCCCTACGAAAGGCTTCAACCATATGATTAAAGTGTTGAGGGTGTGACGCGTAGTCACTGGTCATTAAGTCAAATTCGCCTTGCATCCAGCATGCGCCGAGGAATTTGTTCTGCGGATTTTTTGCCAGTGCGGCCCGCGTTctg                                                                                                                                                                                                                   >  1:573242/1‑301 (MQ=255)
                                                                                 gacgaAAATAATATTGGCTAAAACATTATTTTGATAATTGCCATAAATAGCTTCATACGAATGCGGGAAATTTTCTTTCCAGTACCAGGTGGTATCGCCGCAAAACCACGGCGCGTCAGTAATATTATTAAGCTGAGAATGGTATTGTTTTAGATCCCTACGAAAGGCTTCAACCATATGATTAAAGTGTTGAGGGTGTGACGCGTAGTCACTGGTCATTAAGTCAAATTCGCCTTGCATCCAGCATGCGCCGAGGAATTTGTTCTGCGGATTTTTTGCCAGTGCGGCTCGCGTTctgctg                                                                                                                                                                                                                <  2:573957/301‑1 (MQ=255)
                                                                                 gacgaAAATAATATTGGCTAAAACATTATTTTGATAATTGCCATAAATAGCTTCATACGAATGCGGGAAATTTTCTTTCCAGTACCAGGTGGTATCGCCGCAAAACCACGGCGCGTCAGTAATATTATTAAGCTGAGAATGGTATTGTTTTAGATCCCTACGAAAGGCNNNNNNCATATGATTAAAGTGTTGAGGGTGTGACCNGTAGTCACTGGTCATNAAGTCAAATTCGCCTTGCATCCAGCATGCGCCGAGGAAATTGTACTTCGGATTTTTTGACCGTGCGGCTCGCGTTatgctt                                                                                                                                                                                                                >  1:589010/1‑300 (MQ=255)
                                                                                         taataTTGGCTAAAACATTATTTTGATAATTGCCATAAATAGCTTCATACGAATGCGGGAAATTTTCTTTCCAGTACCAGGTGGTATCGCCGCAAAACCACGGCGCGTCAGTAATATTATTAAGCTGAGAATGGTATTGTTTTAGATCCCTACGAAAGGCTTCAACCATATGATTAAAGTGTTGAGGGTGTGACGCGTAGTCACTGGTCATTAAGTCAAATTCGCCTTGCATCCAGCATGCGCCGAGGAATTTGTTCTGCGGATTTTTTGCCAGTGCGGCTCGCGTTCTGCTGACTAAATc                                                                                                                                                                                                        >  2:480453/1‑301 (MQ=255)
                                                                                             tttGGCTAAAACATTATTTTGATAATTGCCATAAATAGCTTCATACGAATGCGGGAAATTTTCTTTCCAGTACCAGGTGGTATCGCCGCAAAACCACGGCGCGTCAGTAATATTATTAAGCTGAGAATGGTATTGTTTTAGATCCCTACGAAAGGCTTCAACCATATGATTAAAGTGTTGAGGGTGTGACGCGTAGTCACTGGTCATTAAGTCAAATTCGCCTTGCATCCAGCATGCGCCGAGGAATTTGTTCTGCGGATTTTTTGCCAGTGCGGCTCGCGTTCTGCTGACTAAATCCTgg                                                                                                                                                                                                    <  2:589734/300‑1 (MQ=255)
                                                                                               tGGCTAAAACATTATTTTGATAATTGCCATAAATAGCTTCATACGAATGCGGGAAATTTTCTTTCCAGTACCAGGTGGTATCGCCGCAAAACCACGGCGCGTCAGTAATATTATTAAGCTGAGAATGGTATTGTTTTAGATCCCTACGAAAGCCTTCAACCATATGATTAAAGTGTTGAGGGTGTGACGCGTCGTCACTGGTCATTAAGTCAAATTCGCCTTGCATCCAGCATCCGCCGAGGAATTTGTTCTCAGGATCTTTTGCCAGTGCGGCTCGCGTTCCGCTCACTaaccccccgtc                                                                                                                                                                                                  >  1:293658/1‑292 (MQ=255)
                                                                                                    aaaaCATTATTTTGATAATTGCCATAAATAGCTTCATACGAATGCGGGAAATTTTCTTTCCAGTACCAGGTGGTATCGCCGCAAAACCACGGCGCGTCAGTAATATTATTAAGCTGAGAATGGTATTGTTTTAGATCCCTACGAAAGGCTTCAACCATATGATTAAAGTGTTGAGGGTGTGACGCGTGGTCACTGGTCATTAAGTCAAATTCGCCTTGCATCCAGCATGCGCCGAGGAATTTGTTCTGCGGATTTTTTGCCAGTGCGGCTCGCGTTCTGCTGACTAAATCCTGGTACAGCg                                                                                                                                                                                             >  2:204191/1‑301 (MQ=255)
                                                                                                       aCATTATTTTGATAATTGCCATAAATAGCTTCATACGAATGCGGGAAATTTTCTTTCCAGTAACAGGTGGTATCGCCGCAAAACCACGGCGCGTCAGTAATATTATTAAGCTGAGAATGGTATTGTTTTAGATCCCTACGAAAGGCTTCAACCATATGATTAAAGTGTTGAGGGTGTGACGCGTAGTCACTGGTCATTAAGTCAAATTCGCCTTGCATCCAGCATGCGCCGAGGAATTTGTTCTGCGGATTTTTTGCCAGTGCGGCTCGCGTTCTGCTGACTAAATCCTGGTACAGCGGAg                                                                                                                                                                                          <  1:529880/301‑1 (MQ=255)
                                                                                                       aCATTAGTTTGATAATTGCCATATATAGCTTCATACGAATGCGGGAAATTTTCTTTCCAGTACCAGGTGGTATCGCCGCAAAACCACGGCGCGTCAGTAATACTATTAAGCTGAGAATGGTATTGTTTTAGATCCCTACGAAAGGCTTCAACCATATGATTAAAGTGTTGAGGGTGTGACGCGTAGTCACTGGTCATTAAGTCAAATTCGCCTTGCATCCAGCATGCGCCGAGGAATTTGTTCTGCGGATTTTTTGCCAGTGCGGCTCGCGTTCTGCTGACTAAATCCTGGTACAGCGGAg                                                                                                                                                                                          <  1:549205/301‑1 (MQ=255)
                                                                                                        cATTATTTTGATAATTGCCATAAATAGCTTCATACGAATGCGGGAAATTTTCTTTCCAGTACCAGGTGGTATCGCCGCAAAACCACGGCGCGTCAGTAATATTATTAAGCTGAGAATGGTATTGTTTTAGATCCCTACGAAAGGCTTCAACCATATGATTAAAGTGTTGAGGGTGTGACGCGTAGTCACTGGTCATTAAGTCCAATTCGCCTTGCATCCAGCATGCGCCGAGGAATTTGTTCTGCGGATTTTTTGCCAGTGCGGCTCGCGTTCTGCTGACTAAATCCTGGTACAGCGGAGt                                                                                                                                                                                         >  1:77710/1‑301 (MQ=255)
                                                                                                                                  gCTTCATACGAATGCGGGAAATTTTCTTTCCAGTACCAGGTGGTATCGCCGCAAAACCACGGCGCGTCAGTAATATTATTAAGCTGAGAATGGTATTGTTTTAGATCCCTACGAAAGGCTTCAACCATATGATTAAAGTGTTGAGGGTGTGACGCGTAGTCACTGGTCATTAAGTCAAATTCGCCTTGCATCCAGCATGCGCCGAGGAATTTGTTCTGCGGATTTTTTGCCAGTGCGGCTCGCGTTCTGCTGACTAAATCCTGCTACAGCGGAGTATCCGTTCCCCAACGACAAGcatcaa                                                                                                                                                               >  2:98618/1‑300 (MQ=255)
                                                                                                                                        tACGAATGCGGGAAATTTTCTATCCAGTACCAGGTGGTATCGCCGCAAAACCACGGCGCGTCAGTAATATTATTAAGCTGAGAATGGTATTGTTTTAGATCCCTACGAAAGGCTTCAACCATATGATTAAAGTGTTGAGGGTGTGACGCGTAGTCACTGGTCATTAAGTCAAATTCGCCTTGCATCCAGCATGCGCCGAGGAATTTGTTCTGCGGATTTTTTGCCAGTGCGGCTCGCGTTCTGCTGACTAAATCCTGATACAGCGGAGTATCCGTTCCCCAACGACAAGCATCAtggctgg                                                                                                                                                         >  2:39703/1‑301 (MQ=255)
                                                                                                                                                 ttgAAATTTTCTTTCCAGTACCAGGTGGTATCGCCGCAAAACCACGGCGCGTCAGTAATATTATTAAGCTGAGAATGGTATTGTTTTAGATCCCTACGAAAGGCTTCAACCATATGATTAAAGTGTTGAGGGTGTGACGCGTAGTCACTGGTCATTCAGTCAAATTCGCCTTGCATCCAGCATGCGCCGAGGAATTTGTTCTGCGGATTTTTTGCCAGTGCGGCTCTCTTTCTGCTGACTAAATCCTGGTACAGCGTAGTATCCGTTCCCCAACGACAAGCATCATTGCTGGCTCCGTGcc                                                                                                                                                <  2:203028/299‑1 (MQ=255)
                                                                                                                                                   gAAATTTTCTTTCCAGTACCAGGTGGTATCGCCGCAAAACCACGGCGCGTCAGTAATATTATTAAGCTGAGAATGGTATTGTTTTAGATCCCTACGAAAGGCTTCAACCATATGATTAAAGTGTTGAGGGTGTGACGCGTAGTCACTGGTCATTAAGTCAAATTCGCCTTGCATCCAGCATGCGCCGAGGAATTTGTTCTGCGGATTTTTTGCCAGTGCGGCTCGCGTTCTGCTGACTAAATCCTGGTACAGCGGAGTATCCGTTCCCCAACGACAAGCACAATCGCTGGCCCCGTGccag                                                                                                                                              >  1:202785/1‑299 (MQ=255)
                                                                                                                                                           ttttcCAGTACTAGGTGGTATCGCCGCAAAACCACGGCGCGTCAGTAATATTATTAAGCTGAGAATGGTATTGTTTTAGATCCCTACGAAAGGCTTCAACCATATGATTAAAGTGTTGAGGGTGTGACGCGTAGTCACTGGTCATTAAGTCAAATTCGCCTTGCATCCAGCATGCGCCGAGGAATTTGTTCTGCGGATTTTTTGCCAGTGCGGCTCGCGTTCTGCTGACTAAATCCTGGTACAGCGGAGTATCCGTTCCCCAACGACAAGCATCATGGCTGGCTCCGTGCCGTTCTGAata                                                                                                                                      <  2:338260/300‑1 (MQ=255)
                                                                                                                                                              tcCAGTACCAGGTGGTATCGCCGCAAAACCACGGCGCGTCAGTAATATTATTAAGCTGAGAATGGTATTGTTTTAGATCCCTACGAAAGGCTTCAACCATATGATTAAAGTGTTGAGGGTGTGACGCGTAGTCACTGGTCATTAAGTCAAATTCGCCTTGCATCCAGCATGCGCCGAGGAATTTGTTCTGCGGATTTTTTGCCAGTGCGGCTCGCGTTCTGCTGACTAAATCCTGGTACAGCGGAGTATCCGTTCCCCAACGACAAGCATCAAGGCTGGCTCCGTGCCGTTCTGAATATGt                                                                                                                                   >  2:447473/1‑301 (MQ=255)
                                                                                                                                                               ccAGTACCAGTTGGTATCGCCGCAAAACCACGGCGCGTCAGTAATATTATTAAGCTTAGAATGGTATTGTTTTAGATCCCTACGAAAGGCTTCAACCATATGATTAAAGTGTTGAGGGTGTGACGCGTAGTCACTGGTCATTAAGTCAAATTCGCCTTGCATCCAGCATGCGCCGAGGAATTTGTTCTGCGGATTTTTTGCCAGTGCGGCTCGCGTTCTGCTGACTAAATCCTGGTACAGCGGAGTATCCGTTCCCCAACGACAAGCATCATGGCTGGCTCCGTGCCGTTCTGAATATGTc                                                                                                                                  <  1:151057/301‑1 (MQ=255)
                                                                                                                                                                       aGGTGGTATCGCCGCAAAACCACGGCGCGTCAGTAATATTATTAAGCTGAGAATGGTATTGTTTTAGATCCCTACGAAAGGCTTCAACCATATGATTAAAGTGTTGAGGGTGTGACGCGTAGTCACTGGTCATTAAGTCAAATTCGCCTTGCATCCAGCATGCGCCGAGGAATTTGTTCTGCGGATTTTTTGCCAGTGCGGCTCGCGTTCTGCTGACTAAATCCTGGTACAGCGGAGTATCCGTTCCCCAACGACAAGCATCATGGCTGGCTCCGTGCCGTTCTGAATATGTCCCTTCGCt                                                                                                                          >  2:226498/1‑301 (MQ=255)
                                                                                                                                                                               tCGCCGCAAAACCACGGCGCGTCAGTAATATTATTAAGCTGAGAATGGGATTGTTTTAGATCCCTACGAAAGGCTTCAACCATATGATTAAAGTGTTGAGGGTGTGACGCGTAGTCACTTGTCATTAAGTCAAATTCGCCTTGCATCCAGCATGCGCCGAGGAATTTGTTCTGCGGATTTTTGGCCAGTGCGGCTCGCGTTCTGCTGACTAAATCCTGGTACAGCGGAGTATCCGTTCCCCAACGACAAGCATCATGGCTGGCTCCGTGCCGTTCTGAATATGTCCCTTCGCTGCCCGCgg                                                                                                                  <  2:609556/301‑1 (MQ=255)
                                                                                                                                                                                     aatatccgcggcggggtaGTAATTTTATTAAGCTGATAATGGTATTGTTTTAGATCCCTACGAAAGGCTTCAACCAGATGATTAAAGTGTTGAGGGTGTGACGCGTAGTCACTGGTCATTAAGTCAAATTCGCCTTGCATCCAGCATGCGCCGAGGAATTTGTTCTGCGGATTTTTTGCCAGTGCGGCTCGCGTTCTGCTGACTAAATCCTGGTACAGCGGAGTATCCGTTCCCCAACGACAAGCATCATGGCTGGCTCCGTGCCGTTCTGAATATGTCCCTTCGCTGCCCGCGGTAAAAg                                                                                                            <  1:140173/284‑1 (MQ=255)
                                                                                                                                                                                        aaCCACGGCGCGTCAGTAATATTATTAAGCTGAGAATGGTATTGTTTTAGATCCCTACGAAAGGCTTCAACCATATGATTAAAGTGTTGAGGGTGTGACGCGTAGTCACTGGTCATTAAGTCAAATTCGCCTTGCATCCAGCATGCGCCGAGGAATTTGTTCTGCGGATTTTTTGCCAGTGCGGCTCGCGTTCTGCTGACTAAATCCTGGTACAGCGGAGTATCCGTTCCCCAACGACAAGCATCATGGCTGGCTCCGTGCCGTTCTGAATATGTCACTTCGCTGCCCGCGGTAAAAgaca                                                                                                         >  2:594688/1‑298 (MQ=255)
                                                                                                                                                                                           cACGGCGCGTCAGTAATATTATTAAGCTGAGAATGGTATTGTTTTAGATCCCTACGAAAGGCTTCAACCATATGATTAAAGTGTTGAGGGTGTGACGCGTAGTCACTGGTCATTAAGTCAAATTCGCCTTGCATCCAGCATGCGCCGAGGAATTTGTTCTGCGGATTTTTTGCCAGTGCGGCTCGCGTTCTGCTGACTAAATCCTGGTACAGCGGAGTATCCGTTCCCCAACGACAAGCATCATGGCTGGCTCCGTGCCGTTCTGAATATGTCCCTTCGCTGCCCGCGGTAAAAGCCGATc                                                                                                      <  1:479852/301‑1 (MQ=255)
                                                                                                                                                                                                     cAGTAATATTATTAAGCTGAGAATGGTATTGTTTTAGATCCCTACGAAAGGCTTCAACCATATGATTAAAGTGTTGAGGGTGTGACGCGTAGTCACTGGTCATTAAGTCAAATTCGCCTTGCATCCAGCATGCGCCGAGGAATTTGTTCTGCGGATTTTTTGCCAGTGCGGCTCGCGTTCTGCTGACTAAATCCTGGTACAGCGGAGTATCCGTTCCCCAACGACAAGCATCATGGCTGGCTCCGTGCCGTTCTGAATATGTCCCTTCGCAGCCCGCGGTAAAAGCCGATCCGCCACgccc                                                                                            >  1:633587/1‑298 (MQ=255)
                                                                                                                                                                                                       atAATATTATTAAGCTGAGATTGGTATTGTTTTAGATCCCTACGAAAGGCTTCAACCATATGATTAAAGTGTTGAGGGTGTGACGCGTGGTCACTGGTCATTAAGTCAAATTCGCCTTGCATCCAGCATGCGCCGAGGAATTTGTTCTGCGGATTTTTTGCCAGTGCGGCTCGCGTTCTGCTGACTAAATCCGGGTACAGCGGAGTATCCGTTCCCCAACGACAAGCATCATGGCTGGCTCCGTGCCGTTCTGAATATGTCCCTTCGCTGCCCGCGGTAAAAGCCGATCCGCCACGGCAac                                                                                          <  1:203948/300‑1 (MQ=255)
                                                                                                                                                                                                                      tGAGAATGGTATTGTTTTAGATCCCTACGAAAGGCTTCAACCATATGATTAAAGTGTTGAGGGTGTGACGCGTAGTCACTGGTCATTAAGTCAAATTCGCCTTGCATCCAGCATGCGCCGAGGAATTTGTTCTGCGGATTTTTTGCCAGTGCGGCTCGCGTTCTGCTGACTAAATCCTGGTACAGCGGAGTATCCGTTCCCCAACGACAAGCATCATGGCTGGCTCCGTGCCGTTCTGAATATGTCCCTTCGCTGCCCGCGGTAAAAGCCGATCCGCCACGGCAACACgacaacacaagaa                                                                           >  1:307549/1‑289 (MQ=255)
                                                                                                                                                                                                                              gTATTGTTTTAGATTCCTACGAAAGGCTTCAACCATATGATTAAAGTGTTGAGGGTGTGACGCGTAGTCACTGGTCATTAAGTCAAATTCGCCTTGCATCCAGCATGCGCCGAGGAATTTGTTCTGCGGATTTTTTGCCAGTGCGGCTCGCGTTCTGCTGACTAAATCCTGGTACAGCGGAGTATCCGTTCCCCAACGACAAGCATCATGGCTGGCTCCGTGCCGTTCTGAATATGTCCCTTCGCTGCCCGCCGTAAAAGCCGATCCGCCACGACAACCCGGaaacaaaaaaaacccctcc                                                                   >  1:618165/1‑284 (MQ=255)
                                                                                                                                                                                                                                   gTTTTAGATCCCTACGAAAGGCTTCAACCATATGATTAAAGTGTTGAGGGTGTGACGCGTAGTCACTGGTCATTAAGTCAAATTCGCCTTGCATCCAGCATGCGCCGAGGAATTTGTTCTGCGGATTTTTTGCCAGTGCGGCTCGCGTTCTGCTGACTAAATCCTGGTACAGCGGAGTATCCGTTCCCCAACGACAAGCATCATGGCTGGCTCCGTGCCGTTCTGAATATGTCCCTTCGCTGCCCGCGGTAAAAGCCGATCCGCCACGGCAACACGGAACGATAAGAATCCCTGCATTAta                                                              >  2:385448/1‑300 (MQ=255)
                                                                                                                                                                                                                                             ccTTCGTAAGGCTTCAACCATATGATTAAAGTGTTGAGGGTGTGGCGCGTAGTCACTGGTCATTAAGTCAAATTCGCCTTGCATCCAGCATGCGCCGAGGAATTTGTTCTGCGGATTTTTTGCCAGTGCGGCTCGCGTTCTGCTGACTAAATCCTGGTACAGCGGAGTATCCGTTCCCCAACGACAAGCATCATGGCTGGCTCCGTGCCGTTCTGAATATGTCCCTTCGCTGCCCGCGGTAAAAGCCGATCCGCCACGGCAACACGGAACGATAAGAATCCCTGCATTATCAGGAATGAAg                                                    <  2:243504/301‑1 (MQ=255)
                                                                                                                                                                                                                                                aCGAAAGGCTTCAACCATATGATTAAAGTGTTGAGGGTGTGACGCGTAGTCACTGGTCATTAAGTCAAATTCGCCTTGCATCCAGCATGCGCCGAGGAATTTGTTCTGCGGATTTTTTGCCAGTGCGGCTCGCGTTCTGCTGACTAAATCCTGGTACAGCGGAGTATCCGTTCCCCAACGACAAGCATCATGGCTGGCTCCGTGCCGTTCTGAATATGTCCCTTCGCTGCCCGCGGTAAAAGCCGATCCGCCACGGCAACACGGAACGATAAGAATCCCTGCATTATCAGGAATGAAGGGc                                                 <  2:307953/301‑1 (MQ=255)
                                                                                                                                                                                                                                                      ggCTTCAACCATATGATTAAAGTGTTGAGGGTGTGACGCGTGGTCACTGGTCATTAAGTCAAATTCGCCTTGCATCCAGCATGCGCCGAGGAATTTGTTCTGCGGATTTTTTGCCAGTGCGGCTCGCGTTCTGCTGACTAAATCCTGGTACAGCGGAGTATCCGTTCCCCAACGACAAGCATCATGGCTGGCTCCGTGCCGTTCTGAATATGTCCCTTCGCTGCCCGCGGTAAAAGCCGATCCGCCACGGCCACACGGAACGATAAGAAACCATGCATTAACAGGAATGAAGGGCAGaaaa                                           >  2:402958/1‑300 (MQ=255)
                                                                                                                                                                                                                                                            aaCCATATGATTAAAGTGTTGAGGGTGTGACGCGTAGTCACTGGTCATTAAGTCAAATTCGCCTTGCATCCAGCATGCGCCGAGGAATTTGTTCTGCGGATTTTTTGCCAGTGCGGCTCGCGTTCTGCTGACTAAATCCTGGTACAGCGGAGTATCCGTTCCCCAACGACAAGCATCATGGCTGGCTCCGTGCCGTTCTGAATATGTCCCTTCGCTGCCCGCGGTAAAAGCCGATCCGCCACGGCAACACGGAACGATAAGAATCCCTGCATTATCAGGAATGAAGGGCAGTAATTTCCGt                                     <  2:77800/301‑1 (MQ=255)
                                                                                                                                                                                                                                                            aaCCATATGATTAAAGTGTTGAGGGTGTGACGCGTAGTCACTGGTCATTAAGTCAAATTCGCCTTGCATCCAGCATGCGCCGAGGAATTTGTTCTGCGGATTTTTTGCCAGTGCGGCTCGCGTTCTGCTGACTAAATCCTGGTACAGCGGAGTATCCGTTCCCCAACGACAAGCATCATGGCTGGCTCCGTGCCGTTCTGAATATGTCCCTTCGCTGCCCGCGGTAAAAGCCGATCCGCCACGGCAACACGGAACGATAAGAATCCCTGCATTATCAGGAATGAAGGGAAGTAATTTCCGt                                     >  1:291958/1‑301 (MQ=255)
                                                                                                                                                                                                                                                             tccATATGTTTATGTTTTTGTGGGTGTTACGTGTGTTCACTGGTCATTAAGTTAATGTCGCCGGTCATCCAGCAGGCGCCGAGGAATTGGTTCTGCGGATTTTTTGCCAGTGCGGCTCGCGTTCTGCTGACTAAATCCTGGTACAGCGGAGTATCCGTTCCCCAACGACACGCATCATGGCTGGCTCCGTGCCGTTCTGAATATGTCCCTTCGCTGCCCGCGGTAAAAGCCGATCCGCCACGGCAACACGGAACGATAAGAATCCCTGCATTATCAGGAATGAAGGGCAGTAATTTCCGTg                                    <  1:402395/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                    gATTAAAGTGTTGAGGGTGTGACGCGTAGTCACTGGTCATTAAGTCAAATTCGCCTTGCATCCAGCATGCGCCGAGGAATTTGTTCTGCGGATTTTTTGCCAGTGCGGCTCGCGTTCTGCTGACTAAATCCTGGTACAGCGGAGTATCCGTTCCCCAACGACAAGCATCATGGCTGGCTCCGTGCCGTTCTGAATATGTCCCTTCGCTGCCCGCGGTAAAAGCCGATCCGCCACGGCAACACGGAACGATAAGAATCCCTGCATTACCAGGAATGAAGGGCAGTAATTACCGTGCAatatt                             >  1:299024/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                             gtTGAGGGTGTGACGCGTAGTCACTGGTCATTAAGTCAAATTCGCCTTGCATCCAGCATGCGCCGAGGAATTTGTTCTGCGGATTTTTTGCCAGTGCGGCTCGCGTTCTGCTGACTAAATCCTGGTACAGCGGAGTATCCGTTCCCCAACGACAAGCATCATGGCTGGCTCCGTGCCGTTCTGAATATGTCCCTTCGCTGCCCGCGGTAAAAGCCGATCCGCCACGGCAACACGGAACGATAAGAATCCCTGCATTATCAGGAATGAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTg                    >  2:583327/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                gAGGGTGTGACGCGTAGTCACTGGTCATTAAGTCAAATTCGCCTTGCATCCAGCATGCGCCGAGGAATTTGTTCTGCGGATTTTTTGCCAGTGCGGCTCGCGTTCTGCTGACTAAATCCTGGTACAGCGGAGTATCCGTTCCCCAACGACAAGCATCATGGCTGGCTCCGTGCCGTTCTGAATATGTCCCTTCGCTGCCCGCGGTAAAAGCCGATCCGCCACGGCAACACGGAACGATAAGAATCCCTGCATTATCAGGAATGAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGcc                 >  2:220446/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                          cgcgTAGTCACTGGTCATTAAGTCAAATTCGCCTTGCATCCAGCATGCGCCGAGGAATTTGTTCTGCGGATTTTTTGCCAGTGCGGCTCGCGTTCTGCTGACTAAATCCTGGTACAGCGGAGTATCCGTTCCCCAACGACAAGCATCATGGCTGGCTCCGTGCCGTTCTGAATATGTCCCTTCGCTGCCCGCGGTAAAAGCCGATCCGCCACGGCAACACGGAACGATAAGAATCCCTGCATTATCAGGAATGAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCg       >  1:171195/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                            cgTAGTCACTGGTCATTAAGTCAAATTCGCCTTGCATCCAGCATGCGCCGAGGAATTTGTTCTGCGGATTTTTTGCCAGTGCGGCTCGCGTTCTGCTGACTAAATCCTGGTACAGCGGAGTATCCGTTCCCCAACGACAAGCATCATGGCTGGCTCCGTGCCGTTCTGAATATGTCCCTTCGCTGCCCGCGGTAAAAGCCGATCCGCCACGGCAACACGGAACGATAAGAATCCCTGCATTATCAGGAATGAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTa     >  2:176494/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                               aGTCACTGGTCATTAAGTCAAATTCGCCTTGCATCCAGCATGCGCCGAGGAATTTGTTCTGCGGATTTTTTGCCAGTGCGGCTCGCGTTCTGCTGACTAAATCCTGGTACAGCGGAGTATCCGTTCCCCAACGACAAGCATCATGGCTGGCTCCGTGCCGTTCTGAATATGTCCCTTCGCTGCCCGCGGTAAAAGCCGATCCGCCACGGCAACACGGAACGATAAGAATCCCTGCATTATCAGGAATGAAGGGCAGTAATTTCCGTGAAATATGCAGTGCCTGGCCAACGGTGCCGTAAtg  >  2:30663/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                                                            
ACCGTAATATCCCGTGCTTAAATCGTCCGGGTCTTCATCAGGCGCGTTCGTCAGTCCTCTTTCACCTTGTTGCTGGAAGTCGACGAAAATAATATTGGCTAAAACATTATTTTGATAATTGCCATAAATAGCTTCATACGAATGCGGGAAATTTTCTTTCCAGTACCAGGTGGTATCGCCGCAAAACCACGGCGCGTCAGTAATATTATTAAGCTGAGAATGGTATTGTTTTAGATCCCTACGAAAGGCTTCAACCATATGATTAAAGTGTTGAGGGTGTGACGCGTAGTCACTGGTCATTAAGTCAAATTCGCCTTGCATCCAGCATGCGCCGAGGAATTTGTTCTGCGGATTTTTTGCCAGTGCGGCTCGCGTTCTGCTGACTAAATCCTGGTACAGCGGAGTATCCGTTCCCCAACGACAAGCATCATGGCTGGCTCCGTGCCGTTCTGAATATGTCCCTTCGCTGCCCGCGGTAAAAGCCGATCCGCCACGGCAACACGGAACGATAAGAATCCCTGCATTATCAGGAATGAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTG  >  NC_000913/4536743‑4537330

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