Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,606,031 A→G 20.0% E408G (GAG→GGG)  hyfD → hydrogenase 4, membrane subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,606,0310AG20.0% 87.3 / 3.2 60E408G (GAG→GGG) hyfDhydrogenase 4, membrane subunit
Reads supporting (aligned to +/- strand):  ref base A (18/27);  new base G (10/3);  total (29/31)
Fisher's exact test for biased strand distribution p-value = 2.75e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCGCCTACATCGTCAACCACGCGTTCGCTAAAAGCCTGTTTTTCCTTGTAGCAGGTGCGCTGAGTTACAGCTGCGGCACGCGCTTGTT‑GCCGCGTCTGCGTGGCGTATTGCACACCCTGCCGTTGCCAGGCGTGGGTTTCTGCGTGGCAGCGCTGGCGATTACCGGCGTGCCGCCGTTCAACGGCTTCTTCAGTAAATTCCCGCTGTTTG‑CTGCCG‑GTTTTGCGTTGTCAGTGGAGTACTGGATCCTGCTGCCCGCCATGATTCTTCTGATGATTGAATCGGTCGCCAGTTTCGCCTGGTTTATTCGCTGGTTTGGTCGCGTTGTGCCTGGCAAACCGAGCGAGGCCGTCGCCGATGCCGCACCGCTGCCAGGATCAATGCGCCT‑GGTGTTGATTGTACTGATTGTGATGTCGCTGATTTCCAGCGTAATCGCCGCGACCTGGTTGCAGTAAGGAGATGATGAATGAC  >  NC_000913/2605796‑2606264
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                          
ctcgcttactatggcacccctgtgTTCGCTAAAGGCCTGTGTGTCCTTGTAGCAGGTGCTCTGAGTTACGGCTGCGGAACGCGGTTGTT‑GCCGCGTCTGCGTGGCGTATTGCACACCCTGCCGTTGCCAGGCGTGGGTTTCTGCGTGGCAGCGCTGGCGATTACCGGCGTGCCGCCGTTCAACGGCTTCTTCAGTAAATTCCCGCTGTTTG‑CTGCCG‑GTTTTGCGTTGTCAGTGGAGTa                                                                                                                                                                                                                                         <  2:437279/216‑1 (MQ=255)
atccgcTACATCCTCAACCACGCGTTCGCTAAAAGCCTGTTTTTCCTTGTAGCAGGTGCGCTGAGTTACAGCTGCGGCACGCGCTTGTT‑GCCGCGTCTGCGTGGCGTATTGCACACCCTGCCGTTGCCAGGCGTGGGTTTCTGCGTGGCAGCGCTGGCGATTACCGGCGTGCCGCCGTTCAACGGCTTCTTCAGTAAATTCCCGCTGTTTG‑CTGCCG‑GTTTTGCGTTGTCAGTGGAGTa                                                                                                                                                                                                                                         <  2:481300/234‑1 (MQ=255)
  ctgtctactcttCAACCCCGCTGTCGCTAAAAGCCGTTTTTTCCTTGTAGCAGGGTGGCTGAGTTACAGCTGCGGCCCGGGCTTGTT‑TCCGCGTCTGCGTGGCGTATTGCAACCCCTGCCGTTGCCAGGCGTGGGTTTCTGCGTGGCAGCGCTGGCGATTACCGGCGTGCCGCCGTTCAACGTCTTCTCCAGTAAATTCCCGCTGTTTG‑CTGCCG‑GTTTTGCGTTGTCAGTGGAGTACt                                                                                                                                                                                                                                       <  1:863988/228‑1 (MQ=255)
               aaCCACGCGTTCGCTAAAAGCCTGTTTTTCCTTGTAGCAGGTGCGCTGAGTTACAGCTGCGGCACGCGCTTGTT‑GCCGCGTCTGCGTGGCGTCTTGCACACCCTGCCGTTGCCAGGCGTGGGTTTCTGCGTGGCAGCGCTGGCGATTACCGGCGTGCCGCCGTTCAACGGCTTTTTCCGTAAATTTCCCCCTTTTG‑TTGCCG‑GGTTTTCGTGTTCAGGGGGGGACTGGGGCCggctgc                                                                                                                                                                                                                           >  1:758874/1‑238 (MQ=255)
                  ggtgtaaaggagacaggCCTGTTTTTCCTTGTAGCGGGTGCGTTGAGTTACAGCTGCGGCCCGCGCTGGTT‑GCCGCGTCGGCGTGGCGTATTGCACCCCCTGCCGTTGCCAGGCGTGGGTTTCTGCGTGGCAGCGCTGGCGATTACCGGCGTGCCGCCGTTCAACGGCTTCTTCAGTAAATTCCCGCTGTTTG‑CTGCCG‑GTTTTGCGTTGTCAGTGGAGTACTGGATCCTGCTGCCCGc                                                                                                                                                                                                                       <  2:317570/223‑1 (MQ=255)
                   aCGCGTTCGCTAAAAGCCTGTTTTTCCTTGTAGCAGGTGCGCTGAGTTACAGCTGCGGCACGCGCTTGTT‑GCCGCGTCTGCGTGGCGTATTGCACACCCTGCCGTTGCCAGGCGTGGGTTTCTGCGTGGCAGCGCTGGCGATTACCGGCGTGCCGCCGTTCAACGGCTTCTTCAGTAAATTCCCCCTGTTTG‑CTGCCG‑GTTTTGCGTTGTCAGTGGAGTACTGGATCCTCCTTCCCGc                                                                                                                                                                                                                       >  2:405938/1‑238 (MQ=255)
                    cgcgTTCGCTAAAAGCCTGTTTTTCCTTGTAGCAGGTGCGCTGAGTTACAGCTGCGGCACGCGCTTGTT‑GCCGCGTCTGCGTGGCGTATTGCACACCCTGCCGTTGCCAGGCGTGGGTTTCTGCGTGGCAGCGCTGGCGATTACCGGCGTGCCGCCGGTCAACGGCTTCTTCAGTAAATTCCCGCTGTTTG‑CTGCCG‑GTTTTGCGTTTTCAGTGGGGTAATGGATCCTGCTGCCCGccc                                                                                                                                                                                                                     >  1:170295/1‑238 (MQ=255)
                    cgcgTTCGCTAAAAGCCTGTTTTTCCTTGTAGCAGGCGCGCTGAGTTACAGCTGCGGCACGCGCTTGTT‑GCCGCGTCTGCGTGGCGTATTGCACAACCTGCCGTTGCCAGGCGTGGGTTTCTGCGTGGCATGGCTGGCGCTTACCAGCGTGCCGCCCATTACAGGTGTCTTCAGTAAGTTCCCCCTGTTTG‑ATTCCGCGTCTCGC‑TGGTCAGTGGAGtcatcgctccagtacccccccc                                                                                                                                                                                                                     >  1:579169/1‑218 (MQ=255)
                     ggtctCGCTAAACGCTTTTTTTTCCTTGGGGCAGGGGGCTTGAGTTCCCGCGGGGGCCGGCGCTTGTT‑CCCGCTTCTGGGTGGCTTATTGCCCCCCCTGCCGTTGCCGGGCTTGGGTTTCTGCGTGGCAGCGCTGGCGATTCCCGGGGTGCCGCCGTTCAACGGTTTCTTCAGTAAATTCCCGCTGTTTG‑CTGCCG‑GTTTTGCGTTGTCAGTGGAGTACGGGATCCTGCTGCCCGCCAt                                                                                                                                                                                                                    <  1:668852/235‑1 (MQ=255)
                     gcgTTCGCTAAAAGCCTGTTGTTCTTTTGAGCCGGTGCGCTGAGTTACAGCTGCGGCACGCGCTTGTT‑GCCGTGTCTGCGTGGCGTATTGCACACCCTGCCGTTGCCAGGCGTGGTTTTCTGCGTGGCAGCGCTGGCGATTACCGGCGTGCCGCCGTTCAACGGCTTCTTCAGTAAATTCCCGCTGTTTG‑CTGCCG‑GTTTTGCGTTGTCAGTGGAGTACTGGATCCTGCTGCCCGCCAt                                                                                                                                                                                                                    <  2:186824/239‑1 (MQ=255)
                       gTTCGCTAAAAGCCTGTTTTTCCTTGTAGCAGGTGCGCTGAGTTACAGCTGCGGCACGCGCTTGTT‑GCCGCGTCTGCGTGGCGTATTGCACACCCTGCCGTTGCCAGGCGTGGGGTTCTGCGTGGCAGCGCTGGCGATTACCGGCGTGACGCCGTTCAACGGCTTTTTAAGGAAATTTCCGCTTGTTG‑CTGCCG‑GGTTTGGGTTTTCAGGGGGGGACTGGGTCcttgggccccccagct                                                                                                                                                                                                                  >  1:73467/1‑225 (MQ=255)
                        ttCGCTAAAAGCCTGTTTTTCCTTGTAGCAGGTGCGCTGAGTTACAGCTGCGGCACGCGCTTGTT‑GCCGCGTCTGCGTGGCGTATTGCACACCCTGCCGTTGCCAGGCGTGGGTTTCTGCGTGGCAGCGCTGGCGATTACCGGCGTGCCGCCGTTCAACGGCTTCTTCAGTAAATTCCCGCTGTTTG‑CTGCCG‑GTTTTGCGTTTTCAGTGGGGTACTGGATCCTGCTGCCCGCCATGAt                                                                                                                                                                                                                 >  2:47756/1‑239 (MQ=255)
                            cTAAAAGCCTGTTTTTCCTTGTAGCAGGTGCGCTGAGTTACAGCTGCGGCACGCGCTTGTT‑GCCGCGTCTGCGTGGCGTATTGCACACCCTGCCGTTGCCAGGCGTGGGTTTCTGCGTGGCAGCGCTGGCGATTACCGGCGTGCCGCCGTTCAACGGCTTCTTCCGTAAATTCCCGCTGTTTG‑CTGCCG‑GTTTTCCGTTTTACGTGGGGTACTGGATCCTGCTGCCCGCCATGCttct                                                                                                                                                                                                              >  2:51169/1‑238 (MQ=255)
                                             cTTGTAGCAGGTGCGCTGAGTTACAGCTGCGGCACGCGCTTGTT‑GCCGCGTCTGCGTGGCGTATTGCACACCCTGCCGTTGCCAGGCGTGGGTTTCTGCGTGGCAGCGCTGGCGATTACCGGCGTGCCGCCGTTCAACGGCTTCTTCAGTAAATTCCCGCTGTTTG‑CTGCCC‑GTTTTGCGTTTTCAGTGGGGTACAGGATCCTGCTGGCCCCCCTGGttcttccggtctcttataaaca                                                                                                                                                                                            >  1:317570/1‑223 (MQ=255)
                                                 tAGCAGGTGCGCTGAGTTACAGCTGCGGCACGCGCTTGTT‑GCCGCGTCTGCGTGGCGTATTGCACACCCTGCCGTTGCCAGGCGTGGGTTTCTGCGTGGCAGCGCTGGCGATTACCGGCGTGCCGCCGTTCAACGGCTTCTTCAGTAAATTCCCGCTGTTTG‑CTGGCG‑GGTTTGCGTTGTCAGTGGAGGACTGGATCCTGCCGCCCGCCCAGATTCTTCTGTTGATGGaaaggtggggc                                                                                                                                                                                        >  1:390995/1‑230 (MQ=255)
                                                 gaGCAGGTGCGCTGAGTTACATCGGCGGCACGCGCTTGTT‑GCCGCGTCTGCGTGGCGTATTGCACACCCTGCCGTCGCCAGGCGTGGGTTTCTGCGTGGCAGCGCTGGCGATTACCGGCGTGCCGCCGTTCAACGGCTTCTTCAGTACATTCCCGCGGTTTT‑CTGCCG‑GGTTTTGGTTGTCCGTGGGGGACCGGATCCTCGTGCCCGCCAGCATTCTTCTGATTGTTGAAACGGTCGc                                                                                                                                                                                         >  1:166032/2‑238 (MQ=255)
                                                       gTGCGCTGAGTTACAGCTGCGGCACGCGCTTGTT‑GCCGCGTCTGCGTGGCGTATTGCACACCCTGCCGTTGCCAGGCGTGGGTTTCTGCGTGGCAGCGCTGGCGATTACCGGCGTGCCGCCGTTCAACGGCTTCTTCAGTAAATTCCCGCTTGTTG‑CTGCCG‑GTTTTGCGTTGTCACTGGAGTACTGGGTCCTGCTGCCCGCCATGATTCTCTTGATGATTGAATCGGTCGCCAGttt                                                                                                                                                                                   >  2:119109/1‑238 (MQ=255)
                                                          cgcTGAGTTCCGGGTGCGGCACGCGTTTGTT‑GCCGCGTCGGCGTGGCGTTTTGCACACCCTGCCGTTGCCGGGCGGGGGTTTCTGCGTGGCAGCGCTGGCGATTACCGGCGTGCCGCCGTTCAACGGCTTCTTCAGTAAATTCCCGCTGTTTG‑CTGCCG‑GTTTTGCGTTGTCAGTGGAGTACTGGATCCTGCTGCCCGCCATGATTCTTCTGATGATTGAATCGGTCGCCAGTTTCGc                                                                                                                                                                                <  2:703320/238‑1 (MQ=255)
                                                            cTGAGTTACAGCTGCGGCACGCGCTTGTT‑GCCGCGTCTGCGTGGCGTATTGCACACCCTGCCGTTGCCAGGCGTGGGTTTCTGCGTGGCAGCGCTGGCGATTACCGGCGTGCCGCCGTTCAACGGCTTCTTCATTAAATTCCCCCTGTTTG‑TCGCCG‑GTTTTGCGTTGTCCGGGGCGGACTGGGTCCTTGCGCCCCCCCTGCTTTTTTTTGTGGTTGAAACGGtgcgcaggttccccgg                                                                                                                                                                             >  1:881743/1‑224 (MQ=255)
                                                                 ttACAGCTGCGGCACGCGCTTGTT‑GCCGCGTCTGCGTGGCGTATTGCACACCCTGCCGTTGCCAGGCGTGGGTTTCTGCGTGGCAGCGATGGCGATTACCGGCGTGCCGCCGGTCAACTGCTTCTTCAGTAAATTACCGCTGGTTG‑CTGCAC‑GTTTTGAGTTTTCAGTGGGGTACAGGGTCCGCGCGCCCGCACGCATTCTTCTGATGATTTAATCGGCGGCCGATTTCCCCCGGGTTa                                                                                                                                                                        >  1:598534/1‑239 (MQ=255)
                                                                      tcctggggcctgccTTGTTGGCCGGTTCGGCGTGGGGTTTTCCCCCCCCGTCCGTTGCCGGGGGTGGTTTTCGGCGTGGCAGCGCTGGCTATTCCCGGCGTGCCGCCGTTCAACGGTTTTTTCAGTAAATTCCCGCTGTTTG‑CTGCCG‑GTTTTGCGTTGTCAGTGGAGTACTGGATCCTGCTGCCCGCCATGATTCTTCTGATGATTGAATCGGTCGCCAGTTTCGCCTGGTTTATTCg                                                                                                                                                                    <  1:629396/226‑1 (MQ=255)
                                                                         gCGGCACGCGCTTGTT‑GCCGCGTCTGCGTGGCGTATTGCACACCCTGCCGTTGCCAGGCGTGGGTTTCTGCGTGGCAGCGCTGGCGATTACCGGCGTGCCGCCGTTCAACGGCTTCTTCAGTAAATTCCCGCTGTTTG‑CTGCCG‑GTTTTGCGTTGTCAGTGGAGTACTGGATCCTGCTGCCCGCCATGATTCTTCTGTTGATTGAATCGGTCGCCAGTTTCGCCTGGTTTATTCGCTgg                                                                                                                                                                >  2:480449/1‑239 (MQ=255)
                                                                         gCGGCACGCGCTTGTT‑GCCGCGTCTGCGTGGCGTATTGCACACCCTGCCGTTGCCAGGCGTGGGTTTCTGCGTGGCAGCGCTGGCGATTACCGGCGTGCCGCCGTTCAACGGCTTCTTCAGTAAATTCCCGCTGTTTG‑CTGCCG‑GTTTTGCGTTGTCAGTGGAGTACTGGATCCTGCTGCCCGCCACGATTCTTCTGATGATTGAATCGGTCGCCAGTTTCGCCTGGTTTATTCGCTgg                                                                                                                                                                >  2:257730/1‑239 (MQ=255)
                                                                                    ttgtt‑gCCGCGTCTGCGTGGCGTATTGCACGCCCGGCCGTTGCCAGGCGTGGGTTTCTGCGTGGCAGCGCTGGCGATTACCGGCGTGCCGCCGTTCAACGGCTTCTTCAGTAAATTCCCGCTGTTTG‑CTGCCG‑GTTTTGCGTTGTCAGTGGAGTACTGGATCCTGCTGCCCGCCATGATTCTTCTGATGATTGAATCGGTCGCCAGTTTCGCCTGGTTTATTCGCTGGTTTGGTcgcg                                                                                                                                                      <  2:237337/238‑1 (MQ=255)
                                                                                    ttgtt‑gCCGCGTCTGCGTGGCGTATTGCACACCCTGCCGTTGCCAGGCGTGGGTTTCTGCGTGGCAGCGCTGGCGATTACCGGCGTGCCGCGGTTCGACGGCTTCTTCAGTAAATTCCCGCTGTTTG‑CTGCCG‑GTTTTGCGTTGTCAGTGGAGTACTGGATCCTGCTGCCCGCCATGATTCTTCTGATGATTGAATCGGTCGCCAGTTTCGCCTGGTTTATTCGCTGGTTTGGTcgcg                                                                                                                                                      <  2:174546/238‑1 (MQ=255)
                                                                                               gTCTGCGTGGCGTATTGCACACCCTGCCGTTGCCAGGCGTGGGTTTCTGCGTGGCAGCGCTGGCGATTACCGGCGTGCCGCCGTTCAACGGCTTCTTCAGTAAATTCCCGCTGTTTG‑CTGCCG‑GTTTTGCGTTGTCAGTGGAGTACTGGATCCTGCTGCCCGCCATGATTCTTCTGGTGATTGAATCGGTCGCCAGTTTCGCCTGGGTTTTTCCCTGGTTTGGTTGCGTTGGCCCTggt                                                                                                                                           >  2:676441/1‑238 (MQ=255)
                                                                                                         cGTATTGCACACCCTGCCGTTGCCAGGCGTGGGTTTCTGCGTGGCAGCGCTGGCGATTACCGGCGTGCCGCCGTTCAACGGCTTCTTCAGTAAATTCCCGCTGTTTG‑CTGCCG‑GTTTTTCGTTGTCAGTGGAGTACTGGATCCTGCTGCCCGCCATGATTTTTCTGATGATTGAATCGGTCGCCAGTTTTGCCTGGTTTATTTCCTGGTTTGGTCCCGTTGTGCCTGGCAAACcggccg                                                                                                                                 >  2:256354/1‑235 (MQ=255)
                                                                                                              tGCACACCCTGCCGTTGCCAGGCGTGGGTTTCTGCGTGGCAGCGCTGGCGATTACCGGCGTGCCGCCGTTCAACGGCTTCTTCAGTAAATTCCCGCTGTTTG‑CTTCCG‑GTTTTGCGTTTTCAGTGGGGTACTGGGTTCTCCCGCCCGCCATGATTCTTCTGCTGATTGGATCGGTCGCCCGTTTCGCCTCGTTTTTTCGCGGGGTTGGTCCCGTGGTGGCCGGCGAACCGCGCGAGGcc                                                                                                                            >  1:547228/1‑239 (MQ=255)
                                                                                                                   ccccTGCCGTTGCCAGGCTTGGGTTTCTGCGTGGCAGCGCTGGCGATTGCCGGCGTGCCGCCGTTCAACGGCTTCTTCAGTAAATTCCCGCTGTTTG‑CTGCCG‑GTTTTGCGTTGTCAGTGGAGTACTGGATCCTGCTGCCCGCCATGATTCTTCTGATGATTGAATCGGTCGCCAGTTTCGCCTGGTTTATTCGCTGGTTTGGTCGCGTTGTGCCTGGCAAACCGAGCGAGGCCGTCGc                                                                                                                       <  2:166032/238‑1 (MQ=255)
                                                                                                                   acCCTGCCGTTGCCAGGCGTGGGTTTCTGCGTGGCAGCGCTGGCGATTACCGGCGTGCCGCCGTTCAACGGCTTCTTCAGTAACTTCCCGCTGTTTG‑CTGCCC‑GTTTTGCGTTGTCAGTGGAGTACTGGATCCTGCTGCCCGCCATGATTCTTCTGATTATTTAATCGGTTGCCCGTTTCGCCTGGTTTTTTGACTGGTTTGGTTGCGGTTTGCCTGGCAAACCGGGCGGGGcactccc                                                                                                                       >  1:402889/1‑233 (MQ=255)
                                                                                                                     ccTGCCGTTGCCAGGCGTGGGTTTCTGCGTGGCAGCGCTGGCGATTACCGGCGTGCCGCCGTTCAACGGCTTCTTCAGTAAATTCCCGCTGTTTG‑CTGCCG‑GTTTTGCGTTGTCAGTGGAGTACTGGATCCTGCTGCCCGCCCTGATTCTTCTGATGATTGAATCGGTCGCCAGTTTCGCCTGGTTTATTCGCTGGTTTGGTCGCGTTTGGCCTGGCAAACCCGCCGAGGCCGTCGCCg                                                                                                                     >  2:479324/1‑239 (MQ=255)
                                                                                                                       tGCCGTTGCCAGGCGTGGGTTTCTGCGTGGCAGCGCTGGCGATTACCGGCGTGCCGCCGTTCAACGGCTTCTTCAGTAAATTCCCGCTGTTTG‑CTGCCG‑GTTTTGCGTTGTCCGTGGAGTACTGGATCCTGCTGCCCGCCAGGCTTCTTCTGATGATTGAATCGGTGGCCAGTTTCGCCTGGTTTATTTGCTGGGTTGGTCGCGGTGGGCCTGGCGAACcgcgggggggcggcgcggca                                                                                                                   >  1:563753/1‑221 (MQ=255)
                                                                                                                        gCCGTTGCCAGGCGTGGGTTTCTGCGTGGCAGCGCTGGCGATTACCGGCGTGCCGCCGTTCAACGGCTTCTTCAGTAAATTCCCGCTGTTTG‑CTGCCG‑GTTTTGGGTTGTCAGTGGAGGACTGGATCCTGCTTCCCGCCATGATTCTTTTGGTTATTGAATCGGTCGCCAGGTTTGCCCGGGTTTTTTGCTGGTTTGGGTGCGTTTGGCCGTGCAAACCGGGCGGGGCCGTcccggagg                                                                                                                  >  1:675630/1‑232 (MQ=255)
                                                                                                                        gCCGTTGCCAGGCGTGGGTTTCTGCGTGGCAGCGCTGGCGATTACCGGCGTGCCGCCGTTCAACGGCTTCTTCAGTAAATTCCCGCTGTTTG‑CTGCCG‑GTTTTGCGTTGTCAGTGGAGTACTGGATCCTGCTGCCCGCCATGATTCTTCTGATGATTGAATCGGTCGCCAGTTTCGCCTGGTTTATTCGCTGATTTGGTCGCGTTGTGCCTGGCAAACCGCGCGAGGCCGTCGCCCAt                                                                                                                   >  2:737942/1‑238 (MQ=255)
                                                                                                                                         gTTTCTGCGTGGCAGCGCTGGCGATTACCGGCGTGCCGCCGTTCAACGGCTTCTTCAGTAAATTCCCGCTGTTTT‑CTGCCG‑GGTTTTCGTTTTCAGTGGAGTAATGGATCCTGCTGCCCGCCATGATTCTTCTTGTTATTGAATCGTTCGCCATTTTTTCCTGGTTTATTCGCTGGGTTTGTCGTGTTGTTCCTGGCAAACCGAGCGAGGCCGTCGCCGATTCCCCACCGCCGCCCgg                                                                                                  >  1:501868/1‑238 (MQ=255)
                                                                                                                                                       gcgcgggggtttcccgggttcccgcctTTCAACGGCTTTTTCAGTAAATCCCTGGTGTTTG‑CGGCCG‑GTTTGGCGTTGGCAGTGGAGTACGGGATCCTGTTGCCCGCCATGATTCTTCTGATGATTGAATCGGTCGCCAGTTTCGCCTGGTTTATTCGCTGGTTTGGTCGCGTTGTGCCTGGCAAACCGAGCGAGGCCGTCGCCGATGCCGCACCGCTGCCAGGATCAATGCGCCT‑Ggt                                                                                  <  1:480449/218‑1 (MQ=255)
                                                                                                                                                        gtctggtgtgtcgcgggggggcgcctTTCAGCGGCTTCTTAAGTAACTTCCGCCTTTTGT‑CTCCCG‑GTTTTGCGTTTTCAGTGGAGTACTGGGTCCTGCTGCCCGCCATGATTCTTCTGATGATTGAATCGGTCGCCAGTTTCGCCTGGTTTATTCGCTGGTTTGGTCGCGTTGTGCCTGGCAAACCGAGCGAGGCCGTCGCCGATGCCGCCCCGCTGCCAGGATCAATGCGCCT‑Ggt                                                                                  <  1:229924/217‑1 (MQ=255)
                                                                                                                                                        gggtGGGGATTCCCGGGGTGCCCGCGGTGCAGGGTTTTTTCGGGAAATTCCCGC‑GTGTGTCGGCCG‑GTTTTGGGTTGTCGGTGGGTTCGTGGTCCCGGCGGCCCGCCATGATTCTTCGGATGATTGAATCGGTCGCCAGTTTCGCCTGGTTTATTCGCTGGTTTGGTCGCGTTGTGCCTGGCAAACCGAGCGAGGCCGTCGCCGATGCCGCACCGCTGCCAGGATCAATGCGCCT‑Ggt                                                                                  <  1:637993/235‑1 (MQ=255)
                                                                                                                                                        cgcTGGCGATTCCCGGCGTGCCGCCGTTCAACGGCTTCTTCAGTAAATTCCCGCTGTTTG‑CTGCCG‑GTTTTGCGTTGTCAGTGGAGTACTGGATCCTGCTGCCCGCCATGATTCTTCTGATGATTGAATCGGTCGCCAGTTTCGCCTGGTTTATTCGCTGGTTTGGTCGCGTTGTGCCTGGCAAACCGAGCGAGGCCGTCGCCGATGCCGCACCGCTGCCAGGATCAATGCGCCT‑Ggt                                                                                  <  2:170295/238‑1 (MQ=255)
                                                                                                                                                        cgcTGGCGATTACCGGCGTGCCGCCGTTCAACGGCTTCTCCAGGAAATTCCCGCTGTTTG‑CTGCGG‑GTTTGGCGTTGTCAGTGGAGTACTGGATCCTGCTGCCCGCCATGATTCTTCTGATGATTGAATCGGTCGCCAGTTTCGCCTGGTTTATTCGCTGGTTTGGTCGCGTTGTGCCTGGCAAACCGAGCGAGGCCGTCGCCGATGCCGCACCGCTGCCAGGATCAATGCGCCT‑Ggt                                                                                  <  2:758874/238‑1 (MQ=255)
                                                                                                                                                                       gCGTGCCGCCGTTCAACGGCTTCTTCAGGAAATTCCCGCTGTTTG‑CTGCCG‑GTTTTGCGTTGTCAGGGGAGTCATGGCTCCTGCTGCCCGCCATGATTCTTCTGATGATTGAATCGGTCGCCCGATTTGCCTGGTTTTTTCGCTGGGTTGGTCGCGTTGTGCCTTGCAAACCGGGCGAGGCCCTCGCCGGTTCCGCCCCCGTGCCGGGATCAATACGCCTCGGTG‑TGATTGTACTTAtt                                                                  >  1:763795/1‑239 (MQ=255)
                                                                                                                                                                                        ggtttttcaattAAATTCCCGTTGTTTG‑CTGCCG‑GTTTTGCGTTGTCAGTGGAGTACTGGTTCCTGCTGCCCGCCATGATTCTTCTGATGATTGAATCGGTCGCCAGTTTCGCCTGGTTTATTCGCTGGTTTGGTCGCGTTGTGCCTGGCAAACCGAGCGAGGCCGTCGCCGATGCCGCACCGCTGCCAGGATCAATGCGCCT‑GGTGTTGATTGTACTGATTGTGATGTCGCTGATTTc                                                 <  2:392428/228‑1 (MQ=255)
                                                                                                                                                                                                          ccGGTGTTTG‑CTGCCG‑GTTTTGCGTTGTCAGTGGGGTACTGGATCCTTCTGCCCGCCATGATTCTTCTGATGATTGAATCGGTCGCCAGTTTCGCCTGGTTTATTCGCTGGTTTGGTCGCGTTGTGCCTGGCAAACCGAGCGAGGCCGTCGCCGATGCCGCACCGCTGCCAGGATCAATGCGCCT‑GGTGTTGATTGTACTGATTGTGATGTCGCTGATTTCCAGCGTAATCGCCGCGAc                               <  2:881743/239‑1 (MQ=255)
                                                                                                                                                                                                             cTGTTTG‑CTGCCG‑GTTTTGCGTTGTCAGTGGAGTACTGGATCCTGCTGCCCGCCATGATTCTTCTGATGATTGAATCGGTCGCCAGTTTCGCCTGGTTTATTCGCTGGTTTGGTCGCGTTGTGCCTGGCAAACCGAGCGAGGCCGTCGCCGATGCCGCACCGCTGCCAGGATCAATGCGCCT‑GGTGTTGATTGTACTGATTGTGATGTCGCTGATTTCCAGCGAAATCGCCGCGACCTg                            >  2:769612/1‑239 (MQ=255)
                                                                                                                                                                                                              tttttG‑CTCGCG‑GTTTTGCTTTGTCAGTGGATTACTGGTTCCTGCTGCCCGCCATGATTCTTCGGATGATTGAATCGGTCGCCAGTTTCGCCTGGTTTATTCGCTTGTTTGGTCGCGTTGTGCCTGGCAAACCGAGCGAGGCCGTCGCCGATGCCGCACCGCTGCCAGGATCAATGCGCCT‑GGTGTTGATTGTACTGATTGTGATGTCGCTGATTTCCAGCGTAATCGCCGCGACCTgg                           <  2:192399/237‑1 (MQ=255)
                                                                                                                                                                                                              tGTTTG‑CTGCCG‑GTTTTGCGTTGTCAGGGGAGTACTGGATCCTGCTGCCCGCCATGATTCTTCTGATGATTGAATCGGTCGCCAGTTTCGCCTGGTTTATTCGCTGGTTTGGTCGCGTTGTGCCTGGCAAACCGAGCGAGGCCGTCGCCGATGCCGCACCGCTGCCAGGATCAATGCGCCT‑GGTGTTGATTGTACTGATTGTGATGTCGCTGATTTCCAGCGTAATCGCCGCGACCTgg                           <  2:675630/239‑1 (MQ=255)
                                                                                                                                                                                                              ggTTTG‑ATGCCG‑GTTTTGCGTTTTCAGTGGAGGACTGGATCCGGCTGCCCGCCATGATTCTTCTGATGATTGAATCGGTCGGCAGTTTCGCCTGGTTTATTCGCTGGTTTGGTCGCGTTGTGCCTGGCAAACCGAGCGAGGCCGTCGCCGATGCCGCACCGCTGCCAGGATCAATGCGCCT‑GGTGTTGATTGTACTGATTGTGATGTCGCTGATTTCCAGCGTAATCGCCGCGACCTgg                           <  2:658329/238‑1 (MQ=255)
                                                                                                                                                                                                               ggttG‑CGGCCG‑GTTTTGGGTTGTCAGGGGAGTATTGGATCCTGCTGCCCGCCATGATTCTTCTGATGATTGAATCGGTCGCCAGTTTCGCCTGGTTTATTCGCTGGTTTGGTCGCGTTGTGCCTGGCAAACCGAGCGAGGCCGTCGCCGATGCCGCACCGCTGCCAGGATCAATGCGCCT‑GGTGTTGATTGTACTGATTGTGATGTCGCTGATTTCCAGCGTAATCGCCGCGACCTgg                           <  2:553989/236‑1 (MQ=255)
                                                                                                                                                                                                                ttttggtgtgg‑gTGTTGCGTTGTCAGTCGTGTACTGGATCCTGTGGCCCGCCATGATTCTTCTGATGATTGAATCGGTCGCCAGTTTCGCCTGGTTTATTCGCTGGTTTGGTCGCGTTGTGCCTGGCAAACCGAGCGAGGCCGTCGCCGATGCCGCACCGCTGCCAGGATCAATGCGCCT‑GGTGTTGATTGTACTGATTGTGATGTCGCTGATTTCCAGCGTAATCGCCGCGACCTgg                           <  2:833660/228‑1 (MQ=255)
                                                                                                                                                                                                                tgttggcctgggtttttGCTTTGGCAGGGGAGTGCGGGTGCCTGTCGCCAGCCATGATTCTTCTGATGATTGAATCGGTCGCAGGTTTCCCCTGGTTTATTCGCTGGTTTGGTCGCGTTGTGCCTGGCAAACCGAGCGAGGCCGTCGCCGATGCCGCACCGCTGCCAGGATCAATGCGCCT‑GGTGTTGATTGTACTGATTGTGATGTCGCTGATTTCCAGCGTAATCGCCGCGACCTgg                           <  1:566240/226‑1 (MQ=255)
                                                                                                                                                                                                                  ttgtcgcCT‑GTTTTGCCTTGTCAGTGGAGTCCTGGATCCTTGTGCCCGCCATGATTCTTCTGATGATTGAATCGGTCGCCAGTTTCGCCTGGTTTATTCGCTGGTTTGGTCGCGTTGTGCCTGGCAAACCGAGCGAGGCCGTCGCCGATGCCGCACCGCTGCCAGGATCAATGCGCCT‑GGTGTTGATTGTACTGATTGTGATGTCGCTGATTTCCAGCGTAATCGCCGCGACCTGGTTg                        <  2:563753/234‑1 (MQ=255)
                                                                                                                                                                                                                  ttgcgtccG‑GTTTGGGTTTGTCGTGGGAGTCCTGGATCCTGCTGCCCGCCATGATTCTTCTGATGATTGAATCGGTCGCCAGTTTCGCCTGGTTTATTCGCTGGTTTGGTCGCGTTGTGCCTGGCAAACCGAGCGAGGCCGTCGCCGATGCCGCACCGCTGCCAGGATCAATGCGCCT‑GGTGTTGATTGTACTGATTGTGATGTCGCTGATTTCCAGCGTAATCGCCGCGACCTGGTTg                        <  2:383570/233‑1 (MQ=255)
                                                                                                                                                                                                                  tg‑ctgcCG‑GTTTTGCGTTGTCAGTGGGGTACTGGATCCTGCTGCCCGCCATGATTCTTCTGATGATTGAATCGGTCGCCAGTTTCGCCTGGTTTATTCGCTGGTTTGGTCGCGTTGTGCCTGGCAAACCGAGCGAGGCCGTCGCCGATGCCGCACCGCTGCCAGGATCAATGCGACT‑GGTGTTGATTGTACTGATTGTGATGTCGCTGCTTTCCCGCGTAATCACCGCGACTGGGTTGc                       >  2:829948/1‑239 (MQ=255)
                                                                                                                                                                                                                  tg‑ctgcCG‑GTTTTGCGTTGTCAGTGGAGTACTGGATCCTGCTGCCCGCCATGATTCTTCTGATGATTGAATCGGTCGCCAGTTTCGCCTGGTTTATTCGCTGGTTTGGTCGCGTTGTGCCTGTCAAACCGAGCGAGGCCGTCGCCGATGCCGCACCGCTGCCAGGATCAATGCGCCT‑GGTGTTGATTGTACTGTTTGTGATGTCCCTGATTTCCAGCGGAATCGCCCCGACCGGGTTGc                       >  1:553989/1‑239 (MQ=255)
                                                                                                                                                                                                                       ggcG‑GTTTTGCGTTGTAAGTGGGGGATTGGAGCCTGCTGCCCGCCATGATTCTTCGTATGATTGAGTCGGTCGCCAGTTTCGCCTGGTTTATTCGCTGGTTTGGTCGCGTTGTGCCTGGCAAACCGAGCGAGGCCGTCGCCGATGCCGCACCGCTGCCAGGATCAATGCGCCT‑GGTGTTGATTGTACTGATTGTGATGTCGCTGATTTCCAGCGTAATCGCCGCGACCTGGTTGCAGt                    <  2:501868/236‑1 (MQ=255)
                                                                                                                                                                                                                       gcgg‑gTTTTGCGTTTCCATGGGATTCCTGGATCCTGCTGCCCGCCATGATTCTTCTGATGATTGAATCGGTCGCCAGTTTCGCCTGGTTTATTCGCTGGTTTGGTCGCGTTGTGCCTGGCAAACCGAGCGAGGCCGTCGCCGATGCCGCACCGCTGCCAGGATCAATGCGCCT‑GGTGTTGATTGTACTGATTGTGATGTCGCTGATTTCCAGCGTAATCGCCGCGACCTGGTTGCAGTa                   <  2:390995/236‑1 (MQ=255)
                                                                                                                                                                                                                              tttGCGTTGTCAGGGGAGTACTGGATCCTGCTGCCCGCCATGATTCTTCTGATGATTGAATCGGTCGCCAGTTTCGCCTGGTTTATTCGCTGGTTTGGTCGCGTTGTGCCTGGCAAACCGAGCGAGGCCGTCGCCGATGCCGCACCGCTGCCAGGATCAATGCGCCT‑GGTGTTGATTGTACTGATTGTGATGTCGCTGATTTCCAGCGTAATCGCCGCGACCTGGTTGCAGTAAGgaga             <  2:505570/239‑1 (MQ=255)
                                                                                                                                                                                                                                       tCAGTGGAGTACTGGATCCTGCTGCCCGCCATGTTTCTTCTGATGATTGAATCGGTCGCCAGTTTCGCCTGGTTTATTCGCTGGTTTGGTCGCGTTGTGCCTGGCAAACCGAGCGAGGCCGTCGCCGATGCCGCACCGCTGCCAGGATCAATGCGCCT‑GGGGTTGATTGTACTGATTGTGATGTCGCTTATTTCCAGCGTAATCGCCGCGCCCCGGTTGCAGTATGGAgattgtgcat     >  2:257618/1‑231 (MQ=255)
                                                                                                                                                                                                                                         aGTGGAGTTCTGGATCCTGCTTCCCGCCATTATTCTTCTGATGATTGAATCGGTCGCCAGTTTCGCCTGGTTTATTCGCTGGTTTGGTCGCGTTGTGCCTGGCAAACCGAGCGAGGCCGTCGCCGATGCCGCACCGCTGCCAGGATCAATGCGCCT‑GGTGTTGATTGTACTGATTGTGATGTCGCTGATTTCCAGCGTAATCGCCGCGACCTGGTTGCAGTAAGGAGATGATGAATGAc  <  2:73467/239‑1 (MQ=255)
                                                                                                                                                                                                                                         aGTGGAGTATGGTATCCTGCGGCCCGCCTTGATTCTTCTGGTGATTGAATCGGTCGCCAGTTTCGGCTGGTTTATTCGCTGGTTTGGTCGCGTTGTGCCTGGCAAACCGAGCGAGGCCGTCGCCGATGCCGCACCGCTGCCAGGATCAATGCGCCT‑GGTGTTGATTGTACTGATTGTGATGTCGCTGATTTCCAGCGTAATCGCCGCGACCTGGTTGCAGTAAGGAGATGATGAATGAc  <  2:388889/239‑1 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                          
ATCGCCTACATCGTCAACCACGCGTTCGCTAAAAGCCTGTTTTTCCTTGTAGCAGGTGCGCTGAGTTACAGCTGCGGCACGCGCTTGTT‑GCCGCGTCTGCGTGGCGTATTGCACACCCTGCCGTTGCCAGGCGTGGGTTTCTGCGTGGCAGCGCTGGCGATTACCGGCGTGCCGCCGTTCAACGGCTTCTTCAGTAAATTCCCGCTGTTTG‑CTGCCG‑GTTTTGCGTTGTCAGTGGAGTACTGGATCCTGCTGCCCGCCATGATTCTTCTGATGATTGAATCGGTCGCCAGTTTCGCCTGGTTTATTCGCTGGTTTGGTCGCGTTGTGCCTGGCAAACCGAGCGAGGCCGTCGCCGATGCCGCACCGCTGCCAGGATCAATGCGCCT‑GGTGTTGATTGTACTGATTGTGATGTCGCTGATTTCCAGCGTAATCGCCGCGACCTGGTTGCAGTAAGGAGATGATGAATGAC  >  NC_000913/2605796‑2606264

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