Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,304,796 T→C 17.6% intergenic (+405/+310) eco → / ← mqo ecotin, a serine protease inhibitor/malate dehydrogenase, FAD/NAD(P)‑binding domain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,304,7960TC17.6% 132.1 / 14.0 65intergenic (+405/+310)eco/mqoecotin, a serine protease inhibitor/malate dehydrogenase, FAD/NAD(P)‑binding domain
Reads supporting (aligned to +/- strand):  ref base T (31/21);  new base C (6/6);  total (37/28)
Fisher's exact test for biased strand distribution p-value = 7.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

CTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCG  >  NC_000913/2304509‑2305078
                                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                                          
cTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGTCTCATGCGACGCTGGCCCGTCTTATCAGGCCTACACCGCTGAGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGGTACCCCGGATTCGACATCGAACGt                                                                                                                                                                                                                                                                                            >  1:72807/1‑288 (MQ=255)
cTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGt                                                                                                                                                                                                                                                                                            >  1:279219/1‑288 (MQ=40)
    gacgggCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGCCATCCACCGCCCGCGCCGGTTGCCTGATGCGACGCTGGCCGGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGa                                                                                                                                                                                                                                                                                        <  2:1043044/283‑1 (MQ=255)
     aaGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAg                                                                                                                                                                                                                                                                                       >  1:959788/1‑288 (MQ=255)
     aaGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAg                                                                                                                                                                                                                                                                                       >  1:350481/1‑288 (MQ=21)
      tgtggcggccgccggacgTCATACAAGCGTGTCACGTCGAATCCTAGATCCAACGCCCGTGCTGGGTGCCTGTTGCGTGCCTGGGTCGTCTCATCAGGCCTACACCGTTGCGAAGTGATCCACCCCGTAGGTCGGTTAAGGCGGTTACGCCGCATCCGTCATCTACCGCCCAAGCCGGTTGACTGATGAGCCGCTGGCGCGTCTTATCCGGCCGAGACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGc                                                                                                                                                                                                                                                                                      <  2:132602/271‑1 (MQ=11)
        tGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGCCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCAGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGCAGTTCTCCACCCCGCAGGTTGGATAAGTCGGTTACGCCGCATCCGACCTCCAACGTCCGAGCCgg                                                                                                                                                                                                                                                                                   >  1:283291/1‑289 (MQ=255)
           ggCACCCGGGAGTGCGGACAGGGCGGTAAGGCCGCAGCGGAGGGCAAGCGCCCTAGCCGGTTGCCGGAGGCGAAGCTGGCGGGGCATAACAGGCCTAAACCGCTGTGAAGTGCCGCACCCCGGAGGTCTGATAAGGCGGTTACGCCGCATCCGACATCTATCGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGGCTTATCAGGCCTACACCGCTGTGAAGTGCGCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGtt                                                                                                                                                                                                                                                                                 <  2:508939/288‑1 (MQ=11)
               ccccGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCCACACCGCTGTGAAGTGCCCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGCCATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGCCGCGTCCTATCAGGCGCACACCGCTGCGAAGTGCCCCCCCACGTCGGGCGGATAATGCAGTCACGCCGCAATCAACATCCTACGCCGCAGCCGGTGGCCTg                                                                                                                                                                                                                                                                            >  2:128325/1‑289 (MQ=11)
                    tAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGcc                                                                                                                                                                                                                                                                              >  1:1043044/1‑282 (MQ=39)
                            aCAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGgc                                                                                                                                                                                                                                                               >  1:1082636/1‑289 (MQ=37)
                            aCAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGgc                                                                                                                                                                                                                                                               >  1:850365/1‑289 (MQ=37)
                                  cGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGt                                                                                                                                                                                                                                                           <  1:1126565/287‑1 (MQ=35)
                                      cacaCCGCATCGGACATCAAACGCTCTGGCCGGTGGCCGGAGGCGACGCTGGCGCGTCTTATCAGGCGTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTAt                                                                                                                                                                                                                                                      <  2:850365/288‑1 (MQ=11)
                                      cACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCCGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTAt                                                                                                                                                                                                                                                      <  1:505318/288‑1 (MQ=21)
                                      cACGCCGCATCCGACATCCAACGCCCGAACCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTAt                                                                                                                                                                                                                                                      <  1:266459/288‑1 (MQ=21)
                                          ccgcATCCTAAACCTAGCACCCGGCCCGGTGGCCGAACGTGAGGCTGGGGCTGCTTTTAAGGCTTACAACGCGGGGAAGTGCTCCACCCCGTAGGTCGTATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCTAGCCGCTTGCCTGATGCGAGGCTGGTGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAgca                                                                                                                                                                                                                                                 <  2:279219/289‑3 (MQ=11)
                                                         accGCCCGAGCCGTTGGCCTGGGGCGCGGCTGGGGCGGCTTAACAGGCCTAGACCGCTGTGAAGTGCGCGACCGCGTAGGTCAGAAAAGGCGGTGACGCCGCATCCGACATCTAACGCCCCAGCCGGTGGCCGGATGCGGCGCTGGCGCGTCTTATGAGGCCTACACCGCTGTGAAGTGCGCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGaa                                                                                                                                                                                                                                  <  2:460140/287‑1 (MQ=11)
                                                          aCGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGAGGTTACGTCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCCCGTCTTATCAGGACTACACCCCTGTGAAg                                                                                                                                                                                                                                 >  1:1259567/1‑289 (MQ=255)
                                                              acgggCCGGTTGCGTGAGGCGGAGCTGGCGCGTCTGATCAGGCGTACACCGCTGTGAAGTGATCCACCCCGTAGGTCAGATAAGGCGGTTACGCCGCATCCAGGAGCAAACGCCCAAGCCGGTTGCCTGATGCGACGCGGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCt                                                                                                                                                                                                                             <  2:686541/285‑1 (MQ=11)
                                                                 aGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGCGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCCACTCCGACATCCCAAGCCCCAGCCGGTTGCCTGAAGCGACGGCGGCCCGCATTATCCGGACTACACCCCTTCGAAGTGCTcc                                                                                                                                                                                                                           >  2:1141120/1‑288 (MQ=11)
                                                                  gCCGGTTGCCGGGGCCGACGCTGGGGGGGCTTATCAGGCCGAAACCGCTGTGAAGTGCCCCACCCCGTAGGTCGGATAAGGCGGTTAGGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCGGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCa                                                                                                                                                                                                                          <  2:959788/288‑1 (MQ=14)
                                                                            tGATGCGACGCCGGCGCGTCGTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTc                                                                                                                                                                                                               <  1:341439/289‑1 (MQ=255)
                                                                                  gACGCTGGCGCGTCTGATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATaa                                                                                                                                                                                                         <  1:939292/289‑1 (MQ=255)
                                                                                                                 cTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGACGTGCTCCCCCCCGTAGGTTGGATAAGGCGGTTACGCCGCCTCCTACATCCCACGCCCCAGCCGGTTGCCTGCTGCGCCGCTGCCGCGTCTTATCAGCCCTCCCCCCCTTCGAAGTGCTCCGCCCCGTAGCTCGGCAACCGCGGCTCcgccgcccccccaacataacgc                                                                                                                                                                           >  2:89236/1‑272 (MQ=11)
                                                                                                                 cTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGACGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCCGGCGCGTCTTATCACGCCTACACCGCTGTGACATGCTCCACCCCGTAGGTTGGATAAGGCGCTTACGCCGCACCCGACATCTAACGCCCGAGCCGGTTGCCTGCTGCGAACCTGCCGCGTCTTATCAGCCCTCTAACGCAGTGACGCACTACACCCCCGCTGTCCGACAAGGCCGTCCCGCCGCACCCTATATCCATCGc                                                                                                                                                                           >  2:187789/1‑288 (MQ=11)
                                                                                                                  tgtgAAGCGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCACGCCTACACCTCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCCTCCGACTTCTAACGct                                                                                                                                                                          >  1:2527/1‑287 (MQ=255)
                                                                                                                    tgagagcggcaaCCCCGAGGATCGGATAGGGCGGCAAGGCCCGACCCAAAATCAACGGCCCGAGGCGGTGGCGGGAGGCGAAGCGGGCGCGTCTGAGCGGGCCTACACCGCGGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACAGCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCa                                                                                                                                                                        <  2:818262/277‑1 (MQ=14)
                                                                                                                    tgAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATACGCCGGTTACGCCGAACCCGAGCTCTAACGcccca                                                                                                                                                                       >  2:79495/1‑287 (MQ=255)
                                                                                                                              cACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTg                                                                                                                                                               >  1:537906/1‑287 (MQ=255)
                                                                                                                                ccccGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGTCGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACCCCGCCGTGAAGTGCGCCACCCCGTAGGTAGGATCAGGTGGTCACGCCGCATGCGACATCTATCGCACAAGCCGGTAGcccg                                                                                                                                                           >  1:521605/1‑287 (MQ=14)
                                                                                                                                ccccGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCAACCCCGTAGGGCGGATAAGGCGGTTACGCCTCATCCGACATCTAACGCCCAAGCCGGTTGCCTg                                                                                                                                                           >  1:915334/1‑289 (MQ=255)
                                                                                                                                ccccGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCCCTCCACCGCTGTGAAGTGCCACACCCCGTAGGTCGGATAAGGCGGTTACTCCCCCTTCGACCTCTTCCGCCCCAGCCGGTTGcttc                                                                                                                                                           >  2:629891/1‑286 (MQ=14)
                                                                                                                                   cGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACg                                                                                                                                                                                                                                                                    >  1:1232530/1‑181 (MQ=35)
                                                                                                                                                     gtcaCGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGCAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCAcc                                                                                                                                                                                                                        >  2:1217962/4‑207 (MQ=9)
                                                                                                                                                     gtcaCGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCAcc                                                                                                                                                                                                                        <  1:1217962/204‑1 (MQ=14)
                                                                                                                                                                tccatcttcaaccgccctgcccGGTTGCCTGGGGCTACGCTGGCCAGTCTTGTCAGGCCTACACCTCGGGGAAGTGCTCCACCCCGGAGGTCGGAGAAGGCGGTTACGCCGCATCCGACATCCAACGCCCGGGCCGGTTGCCGGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTaca                                                                                                                           <  2:10743/269‑1 (MQ=11)
                                                                                                                                                                    ccgacaaccggccgggccgcgggCCTGATGCGACGCTGGCGCGGCTTATCAGGCCGACACCGCTGTGAAGTGCGCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCGTACACCGCGGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCg                                                                                                                        <  2:563177/266‑1 (MQ=11)
                                                                                                                                                                      atagaaCGCCCGAGTCGGTGGCCTGATGTGACGCGGGCGTGTCTTTGCAGCCGTACGCCGCTGTGAGTTGCTACAACCCGTAGGTGGGATGAGGCGGTTCCGCCACATCGGACATCCAGCTGCCGAGCCGGCTGCCTGATGCGACGCTGGCGCGTCTGATCAGGCCTACACCGCTGTGAAGGGCTCCACCCCGTAGGTCGGATAAGGCGGTCACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCt                                                                                                                      <  2:706386/284‑1 (MQ=11)
                                                                                                                                                                                  aGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCCCCCCCCCGTAGGTCGGATAAGGCGGTTACGCCGCCTCCCACCTCTAACGCCCAAGCCGGTTGCCTGACGCAACGCTGGCGCGTATTATCACGCCTACCCCCCTGTTAAGTGCTCCa                                                                                                         >  2:1069041/1‑289 (MQ=11)
                                                                                                                                                                                            cagATGCGGCGCTCGCGGGTCTTATCAGGCCTCCACCGCTGTGAATGGCTCCCCCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGCCGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAgg                                                                                                <  2:1082636/286‑1 (MQ=255)
                                                                                                                                                                                                     cGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACCCCGCTGTGAAGTGCTCCACCCCGCAGGTTGGATaca                                                                                       >  1:1269623/1‑286 (MQ=255)
                                                                                                                                                                                                               tCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTCCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTCCACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTAcgc                                                                             >  1:1260712/1‑288 (MQ=16)
                                                                                                                                                                                                                       ggCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTAt                                                                                                                                     >  2:473409/1‑224 (MQ=40)
                                                                                                                                                                                                                       ggCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTAt                                                                                                                                     <  1:473409/224‑1 (MQ=40)
                                                                                                                                                                                                                        gCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGa                                                                                                                                                                                     <  1:1241027/175‑1 (MQ=39)
                                                                                                                                                                                                                        gCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGa                                                                                                                                                                                     >  2:1241027/1‑175 (MQ=39)
                                                                                                                                                                                                                             gacCGCTCTGAAGGGTGCGACCCCGCGGGTGGGGTAGGGCGGTTACGCGGCGTCCGCGATCCACCGCGCGGGCCGGTTGCCGTAGAGGACGCTGGCGAGGCTTATCAGGGCTAAACCGCTGTGAAGTGCTCCACCCCGTAGGTGGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTa                                                              <  2:1225336/288‑1 (MQ=11)
                                                                                                                                                                                                                                  cTGTGCAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTCTCAGGCCTACACTCCTGTGAAGTGCTCCACCCCGTGGGTCGGCTCAGGCGGCTACGCCGCACCCTACCTCTAACGCCCCAGCCGGTTGCCTGATGCGAGGCTGGCGAGTCTTATCAGCCATACGCATCTGTGCCTTGCCCTATCCCGTCGGTGCGGTTAGTCGGTCAcgccgctcccaccacctcactc                                                          >  2:1114597/1‑272 (MQ=11)
                                                                                                                                                                                                                                          tGCTCCATCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTCGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGTCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGCTAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCgg                                                 >  1:26831/1‑289 (MQ=255)
                                                                                                                                                                                                                                                 ccccGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTg                                          >  1:1143505/1‑289 (MQ=39)
                                                                                                                                                                                                                                                    cGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAgg                                                                                      >  1:996455/1‑242 (MQ=37)
                                                                                                                                                                                                                                                      tAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTATGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCg                                     >  1:604769/1‑289 (MQ=39)
                                                                                                                                                                                                                                                      tAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCg                                     >  1:889636/1‑289 (MQ=39)
                                                                                                                                                                                                                                                      tAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCg                                     >  1:288883/1‑289 (MQ=39)
                                                                                                                                                                                                                                                              aTAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATcc                                                                                                                                                                                       <  1:56726/135‑1 (MQ=37)
                                                                                                                                                                                                                                                                  ggCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTaca                                                                                                                                                                                                                                            <  2:585176/78‑1 (MQ=39)
                                                                                                                                                                                                                                                                  ggCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCa                                                                                                                                                                        <  2:1192426/146‑1 (MQ=34)
                                                                                                                                                                                                                                                                  ggCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCa                                                                                                                                                                        >  1:1192426/1‑146 (MQ=34)
                                                                                                                                                                                                                                                                  ggCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGcc                                                                                                                                                             >  2:167203/1‑157 (MQ=33)
                                                                                                                                                                                                                                                                  ggCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGcc                                                                                                                                                             <  1:167203/157‑1 (MQ=33)
                                                                                                                                                                                                                                                                       ttACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTTCTCCACCCCGTCGGTTGGATAAGGCCGTTACGCCCCATCCAACATCCCACCGCCGAGCACGTTGCCTGATGCCAACCCGGCCCGTCTTa                     >  2:823111/1‑288 (MQ=11)
                                                                                                                                                                                                                                                                       ttACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGCCACCCCGTAGGTTGGATAAGGCGGTTACGCTGCATCCAACATCTAACGCCCGAGCCCGTTGCCTGAGGCGACCCTGGCGCGCCTTa                     >  2:800798/1‑288 (MQ=255)
                                                                                                                                                                                                                                                                             cgcATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACCCCGCTGCGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGATACGCCGCATCCAACCTCTAAACCCCGAGCCGGTTGCCCGATGCGCCGCTGGCGGCTCTTATCAGCa               >  2:803362/1‑288 (MQ=11)
                                                                                                                                                                                                                                                                                 tCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACg           >  1:735582/1‑288 (MQ=37)
                                                                                                                                                                                                                                                                                      cATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCcaca      <  1:79495/288‑1 (MQ=39)
                                                                                                                                                                                                                                                                                         caaaCGTCCGGGCGGGTTGCCTGATGCGACGCTGGCGCGTCTTACCAGGCCTACATCGCTGTGAAGGGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCg  <  1:89236/287‑1 (MQ=255)
                                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                                          
CTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCG  >  NC_000913/2304509‑2305078

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