Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,304,796 T→C 9.3% 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,7960TC 9.3% 112.5 / 4.7 56intergenic (+405/+310)eco/mqoecotin, a serine protease inhibitor/malate dehydrogenase, FAD/NAD(P)‑binding domain
Reads supporting (aligned to +/- strand):  ref base T (27/21);  new base C (4/2);  total (32/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.87e-01

TGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGA‑GCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTG‑ATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGT‑AGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATA  >  NC_000913/2304512‑2305083
                                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                                               
tgAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGA‑GCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCACCCCCGTAGGTCGGATACGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTACCTG‑ATGCGACGCTGGCGCGTCGTATCAGGCCTACACCGCTGTGAAGTACTCGACCCCGT‑AGGTTGGATAAGGAGGTTACGCCGCATCCGACATCCAACGTCCGa                                                                                                                                                                                                                                                                                             >  1:68859/1‑289 (MQ=255)
       cGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGA‑GCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTG‑ATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGT‑AGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGtt                                                                                                                                                                                                                                                                                      >  1:187559/1‑289 (MQ=40)
               cGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGA‑GCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGGAGTTCTCCACCCCGTAGGTCGGATAAGGCGGCTACGCCGCACCCGACCTCTAACGCCACAGGCGGTTGCCTG‑AAGCGAAGCTGGCGCGTCTTATCAGCCCCACACCGCTGTGAATATCTCCACCCCGT‑CGGTCTGTGTGGGCGCGTACGCCGCAACCGAGCTCCAACGACCAAACCCGTCGGCTGAt                                                                                                                                                                                                                                                                               >  2:201385/1‑288 (MQ=11)
               cGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGA‑GCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTG‑ATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGT‑AGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGAt                                                                                                                                                                                                                                                                               >  1:776110/1‑288 (MQ=39)
               cGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGA‑GCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTG‑ATGCGACGCTGGCGCGGCTTATCAGGCCTACACCGCTGGGAAGTGCTCCACCCCGT‑AGCTTGGATAAGGCGGTTACGCCGCATCCGACATACAACCTCAGCGCCGGTGGCCTAAt                                                                                                                                                                                                                                                                               >  1:370947/1‑288 (MQ=255)
                 tAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGA‑GCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACCCCGCTGTGACGTGCTCCACCCCGTCGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTTCCTG‑CTGCGACTCTGGCCCGTCTTATCAGGCCTACACCGCTGTGAAGTGCCCCACCGCGG‑CCGTTGGCTAAGGGGGTTACGCCTACCCTGAACTCCAACGTCCAAGCCGGtcacgccac                                                                                                                                                                                                                                                                               >  2:850254/1‑278 (MQ=11)
                         aCAAGGGGGTGACGCCCCGTGCGACGTCCGACGCCCGG‑GCCGGTGGCCGGATGGGAGGTTGGCGGGTCTGATCAGGCCTAAACCGCTGTGAAGTGCCCCACCCGGTAGGTCGGATAAGGGGGTTACGCCGAAGCGAACATCTAACGCCCAAGCCGGTTGCCTG‑ATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGT‑AGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTgg                                                                                                                                                                                                                                                                     <  2:444054/288‑1 (MQ=11)
                         aCAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGA‑GCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTCGGTCGGATAAGGCGGTTACGCCGCATCCGACCTCTAACGCCCAAGCCGGTTGCCTG‑ATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGT‑AGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTgcc                                                                                                                                                                                                                                                                    >  1:106886/1‑287 (MQ=37)
                              gCGGTCACGCCGCATCCGACATCAAACGCCCGA‑GCCGGTTGCCTGATGAGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCTGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTG‑ATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCCGAGAAGTGCTCCACCCCGT‑AGGTTGGATAAGGCGGTTACGCCGCATCCGACATCAAAGGTCCGAGCCGGTTGCGTGATGCGACGCCGGcgcg                                                                                                                                                                                                                                                                 >  1:320179/1‑287 (MQ=255)
                                   gaCGCCGCACCCGGTAGCCAACGCCCGA‑GCGGGTGGCGGAAGGAGACGCTGGCGGGTCTTATCAGGCCTACACCGCTGGGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCGAGCCGGTTGCCTG‑ATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGT‑AGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTAt                                                                                                                                                                                                                                                           <  2:725418/287‑1 (MQ=11)
                                   cACGCCGCGTCCGACATCCAACTCCCGG‑GCCGGTTGCCTGATGCGACGGTGGGGCGTCTTATCAGGCCTACGCCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTG‑ATGCGACGCTGGGGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGT‑AGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTAt                                                                                                                                                                                                                                                           <  2:582260/288‑1 (MQ=255)
                                   cACGCCGCAGCCAAC‑GCGAACGCCCCAGGCCGGTTGCGGGAGGGGAAGATGGCGGGACGTAGCAGGCCTAAACCGCGGTGAAGAGCGCCACCCCGTAGGGCGGAGAAGGCGGGTACGCCGCATCGGACATCTAACGCCCGAGCCGGTGGCCTG‑ATGGGACGCTGGCGCGGCTTATCAGGCCTACACCGCTGTGAAGTGCGCCACCCCGT‑AGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTAt                                                                                                                                                                                                                                                           <  2:771721/288‑1 (MQ=11)
                                   acggcccgcgcctagctccctcgcgcta‑GCCGGTGGCGCGGTGGGTCGCTGTCGCTTCTGCTCAGGCCGTCACCGGTGTGATGTGCTCCCCCCCTTATGTCAGATCAGGCGTTTGCGCCGCGTCCGACATCTACCGCCCAAGCCGTGGGCCAG‑ATGCGAGGCTGCCGCTTCTTATCATACTCACCCCGCTGTGAAGTGCTCCACCCCTG‑AGGTCGGATAAGGCGCTTACGCCGCCTCCGATATCCAACGTCCTGGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATc                                                                                                                                                                                                                                                          <  2:320179/262‑1 (MQ=2)
                                          cATCCGACATCCAACGTCCGA‑GCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTG‑ATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGT‑AGGTTGGATCAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTCCCTGATGCCACGCTGGCGCGTCTTATCAGCCCt                                                                                                                                                                                                                                                    >  1:560694/1‑288 (MQ=12)
                                            cccGCCACCCAACTCCCGA‑GCCGGTTGCCTGAGGCGAGGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTG‑ATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGT‑AGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTaca                                                                                                                                                                                                                                                 <  1:452982/288‑1 (MQ=255)
                                                        aggccgagtgcggtggctggATGGGACTCTGGCGCTTCGTAACGGGCCTAAACCGCGGTGAAGGGCTCCACCCCATAGGTCGGATAAGGAGGTCACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTG‑ATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGCCACCCCGT‑AGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCtgtg                                                                                                                                                                                                                                         <  2:584417/267‑1 (MQ=14)
                                                            cGA‑GCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTG‑ATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGT‑AGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTc                                                                                                                                                                                                                                 <  1:788973/289‑1 (MQ=40)
                                                                gCCGGTGCCCCGAGGCGACGCTGGTGCGGCTTATCAGGCCTAAACCGCTGTGCAGTGCTCCAACCCGTAGGTCGGAAAAGGCGGTTACGCCGCATCCGACATCAAACGCCCAAGCCGGTTGCCTG‑CTGCGACGCTGGTGGGTCTTATCAGGCCTGCACAGCTGGGAAGTGCTCCACCCCGT‑AGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCa                                                                                                                                                                                                                               <  2:776110/288‑1 (MQ=11)
                                                                         cggacgcggcgCTGGCGTGTGTTCTCAGGCCTACACGGGTGGGAGGGGCTCCCCCCCGGAGGTCGGATAAGGCGTTGACGCCGCATCCGACATCTAAGGCCCAAGCCGGTTGCCTG‑ATGCGCCGCTGGCGCGTCGTATCAGGCCTGCACCGCTGTGAAGTGCTCCACCCCGT‑AGGTCGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAgg                                                                                                                                                                                                                      <  2:385258/283‑1 (MQ=11)
                                                                                             cTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCCAGCCGGTTGCCTG‑ATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGT‑AGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTAcgccg                                                                                                                                                                                                 <  1:829665/289‑1 (MQ=39)
                                                                                                      ctgacacCGCTGGGATGCGCTCCACCCCGGAGGTCGGATAAGGCGTTTACGCTGCATCCGACATCTAACGCCCAAGCCGGTTGCCGG‑GTGCGTCGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGT‑AGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGGGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCg                                                                                                                                                                                           <  1:261449/283‑1 (MQ=14)
                                                                                                               cTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTG‑ATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGT‑AGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGCAGTGCTCCACCCCGTAGGTCGGATaa                                                                                                                                                                                                              >  2:541673/1‑258 (MQ=255)
                                                                                                               cTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTG‑ATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGT‑AGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATaa                                                                                                                                                                                                              <  1:541673/258‑1 (MQ=255)
                                                                                                                 gtgAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTG‑ATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGCAGTTCTCCACCCCGT‑AGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTAACAGGCCTAACACGCTGTGAAGTTCTACACCCCGTACGTCGGATAAGGCTGttcacgcgcatcgcccctcacaccccaa                                                                                                                                                                             >  2:633759/1‑263 (MQ=14)
                                                                                                                            cACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTG‑ATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGT‑AGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTCCACCCCTGCGGAGTGCTCCACCGCGTAGGTCGGATAAGGCTGTTACGCCCCACCGCACATCTAACGccgccggcgagtggc                                                                                                                                                                  >  2:179806/1‑276 (MQ=11)
                                                                                                                                   tAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTG‑CTGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGT‑AGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGCGCCGGTTGGCTGATGCGACGCTGGCGCGTCTTATCAGGTCTCACCCGATGTCAAGTTCCCTCCCCCGTCGGGCGGATAAGGCACTTACCCCACACCCGACTTCTAAAGCACCAGCCAGtggacggaggca                                                                                                                                                           >  2:582701/1‑278 (MQ=14)
                                                                                                                                   tAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTG‑ATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGT‑AGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGa                                                                                                                                                                                          >  1:637805/1‑258 (MQ=255)
                                                                                                                                                gCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTG‑ATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGT‑AGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCCCGCCTTATCAGGCCTACCCCGCTGCGAAGTGCTCCACCCCGTCGGTCGGATAAGGCCGTTACGCCGCCCCCGACACCCACCGCCCAAGACGGTTGCCTGATCCGACGCTGGCGCGt                                                                                                                                               >  2:808614/1‑288 (MQ=14)
                                                                                                                                                            cATCCGACATCTAACGCCCAAGCCGGTTGCCTG‑ATGCGACGCTGGCGCGTCTTATCAGCCCTACACCGCTGTGAAATGCTCCACCCCGT‑AGGTTGGATAAGGCGGTTACGCCGCATCCGACCTCCAACGTCCGAGCCGGTTGACTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTCTGACTTCCTCCCCCCCCTTGGTCCGATACGGCGCTTCCGCCGCCTCCCACATCCACACCCCAAGCCGGTTGCCTGCTCCGCCTCTTGCGCCTccttccccacca                                                                                                                                   >  2:27309/1‑277 (MQ=14)
                                                                                                                                                              tCCGACATCTAACGCCCAAGCCGGTTGCCTG‑ATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGT‑AGGTTGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCGAGCCGGTTGCCTGCTGCGACGCTGGCGCGTCTTATCCGGCCTCCACCGCTGTGCAGTGCCCCACCCCGTAGGGCGGATACGCCGGTTGCGCCGCACCCCACCTCTCACGCCCAAGCCGGTTGCCTGGTCCGACGcgcgcgcgcccctgcatccctcc                                                                                                                                 >  2:727600/1‑273 (MQ=11)
                                                                                                                                                              tCCGACACGTAATGCACCAGCCGGTT‑CCTGCATGCTACGCTGGCGCGTCGCATGAGGGCGTCACAGGTGCGACGGGCACCACCCGGT‑AGGCTGGATAAGGCGGTTACGCCGCATCGGACATCCGACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTaca                                                                                                                                <  2:370947/289‑1 (MQ=11)
                                                                                                                                                                      cTAACGCCCAAGCCGGTTGCCTG‑ATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGT‑AGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTCTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGCTGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCtgtg                                                                                                                        >  1:292388/1‑289 (MQ=255)
                                                                                                                                                                         aCGCCCAAGCCGGTTGCCTG‑ATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGT‑AGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGCTGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGCGGAGTGCTCCACCCCGTAGGTCGGATACGGCGGTTACGCCTCACCCCACATCTAACGCCCAAGCCGGTTGCCTGATGCGCCGCCGGCGCGTCTTATCAGGCCTACACCGCAGTGa                                                                                                                       >  2:121705/1‑287 (MQ=14)
                                                                                                                                                                              ccgggcggtgggccggtgggccggCTGGCGCTGCTTAGCAGGCCGAGACGGCTGGGAAGTGCGCCACAGCGC‑AGGTCGGATAAGGCGGTTACGCCGCAGCCGACATCGAACGCCCGAGCCGTTTGCCGGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCt                                                                                                                 <  2:255648/266‑1 (MQ=11)
                                                                                                                                                                                                              tCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGT‑AGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGATGTGCAGTGCCCCACCCCGTAGGTCGGATAAGGCGGTTACGCTGCATCCGACATCTCACGCCCAAGCCGGTTGCCTGATCCGACGCTCGCGCGTCTTATCAGTCCCCCACCGCTGTGCCGGGCTCCACCCCGTAGGTGACAGAAGGTGGTTAccc                                                                                  >  2:186755/1‑286 (MQ=11)
                                                                                                                                                                                                                       gCCTACACCGCGGTGAAGTGCTCCACCCCGT‑ATGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCaa                                                                         <  1:323040/288‑1 (MQ=255)
                                                                                                                                                                                                                                 cTGTGAAGTGCTCCACCCCGT‑AGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGc                                                               >  1:641849/1‑288 (MQ=40)
                                                                                                                                                                                                                                     ggctagtcgtcacccGTAAGGTGGGATAAGCGGGTAACCCCGATGCCGCAATCCAACGTCCGAGCCGTGTTCCCGATGCGACTCTGGCGCGTCGTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATACGGCGGTTACGCCGCTTCCGACATCCAACGCCGAAGCCGGTTGCCTGATGGTACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCg                                                            <  2:505205/276‑1 (MQ=11)
                                                                                                                                                                                                                                     gAAGTGCTCCCCCCCGC‑AGTTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCACGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAg                                                          <  1:730506/289‑1 (MQ=255)
                                                                                                                                                                                                                                           cTCCACCCCGT‑AGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACACCCAAGCCGGTTGCCTGATGCGACGCCCGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCCCCACCCCGTAGGTTGCATAAGGCGGTTACGCCGCAACCATCATCCAACTcccccgaccgtt                                                    >  2:734004/1‑280 (MQ=11)
                                                                                                                                                                                                                                           cTCCACCCCGT‑AGGTGGGATAAGGCGGTTACGCCGCATCCGACTTACAACGGCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGt                                                     <  1:850254/288‑1 (MQ=255)
                                                                                                                                                                                                                                              cACCCCGT‑AGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTg                                                   >  1:368926/1‑287 (MQ=39)
                                                                                                                                                                                                                                                ccccGT‑AGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTg                                                   >  1:514328/1‑285 (MQ=39)
                                                                                                                                                                                                                                                   cGT‑AGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGAt                                             >  1:870843/1‑288 (MQ=39)
                                                                                                                                                                                                                                                   cGT‑AGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATg                                            <  1:631136/289‑1 (MQ=39)
                                                                                                                                                                                                                                                           tGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGAc                                        >  1:867578/1‑286 (MQ=38)
                                                                                                                                                                                                                                                             gATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTAt                                                                                                                                          <  1:812593/186‑1 (MQ=35)
                                                                                                                                                                                                                                                              ataatgcttgtACTCCCCATCCTACATCCAACTTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTgg                                   <  1:282700/278‑1 (MQ=255)
                                                                                                                                                                                                                                                              aTAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTCGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACTCCCAAGCCGGTTGCCTGATGCCACGCTGGCGCGCCTTATCAGGCCTCCACCGCCGTGAAGCTCCTCACCCCTTAGGTTGGCTAAGGCGGTTCAGACGCCTCCCACCTCTGTCGcccacgctgtgtgcttgcgacgaccttggt                                  >  2:237307/1‑262 (MQ=11)
                                                                                                                                                                                                                                                                  ggCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGcgcg                               >  1:520823/1‑288 (MQ=35)
                                                                                                                                                                                                                                                                       ttACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATCAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGGGAAGTGCTCACCCCCGGTGGTTGGATAAGGCGGTTACGCCGCCTCCAACCTCTAACTCCCACTCCGGTTGCCTGGTGCGACGCTGCCGCCTcgca                          >  2:713057/1‑285 (MQ=11)
                                                                                                                                                                                                                                                                       ttACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTaca                                                                                                                                >  2:812593/1‑186 (MQ=32)
                                                                                                                                                                                                                                                                            ccgcATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTaca                                                                                                                                                                                                                                                 <  1:313128/68‑1 (MQ=35)
                                                                                                                                                                                                                                                                             cgcATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCAt                   >  1:30843/1‑289 (MQ=35)
                                                                                                                                                                                                                                                                                aTCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACg                <  1:808614/289‑1 (MQ=37)
                                                                                                                                                                                                                                                                                     aCATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCCTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCAGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAACTGCTTCACCCCGTAGGTTGGATAAGGCGGCTACGCCGCATCCAACATCTAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCGTCTCCGCATAcaccac            >  1:239797/1‑288 (MQ=11)
                                                                                                                                                                                                                                                                                              ctCCGTGCCGGTTGGCGGAGGGGAAGCGGGGTCGGCGAATCAGGCCGAAACCGCTGTGAAGTGCGCCAGCCCGTAGTGCGGAGAAGGAGGTGACGCCGCATCCGACAGCTAACGCCCGTGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATa  <  2:718972/288‑1 (MQ=14)
                                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                                               
TGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGA‑GCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTG‑ATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGT‑AGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATA  >  NC_000913/2304512‑2305083

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