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

ATCGCCTACATCGTCAACCACGCGTTCGCTAAAAGCCTGTTTTTCCTTGTAGCAGGTGCGCTGAGTTACAGCTGCGGCACGCGCTTGTTGCCGCGTCTGCGTGGCGTATTGCACACCCTGCCGTTGCCAGGCGTGGGTTTCTGCGTGGCAGCGCTGGCGATTACCGGCGTGCCGCCGTTCAACGGCTTCTTCAGTAAATTCCCGCTGTTTGCTGCCGGTTTTGCGTTGTCAGTGGAGTACTGGATCCTGCTGCCCGCCATGATTCTTCTGATGATTGAATCGGTCGCCAGTTTCGCCTGGTTTATTCGCTGGTTTGGTCGCGTTGTGCCTGGCAAACCGAGCGAGGCCGTCGCCGATGCCGCACCGCTGCCAGGATCAATGCGCCTGGTGTTGATTGTACTGATTGTGATGTCGCTGATTTCCAGCGTAATCGCCGCGACCTGGTTGCAGTAAGGAGATGATGAATGAC  >  NC_000913/2605796‑2606264
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                          
ctcgcttactatggcacccctgtgTTCGCTAAAGGCCTGTGTGTCCTTGTAGCAGGTGCTCTGAGTTACGGCTGCGGAACGCGGTTGTTGCCGCGTCTGCGTGGCGTATTGCACACCCTGCCGTTGCCAGGCGTGGGTTTCTGCGTGGCAGCGCTGGCGATTACCGGCGTGCCGCCGTTCAACGGCTTCTTCAGTAAATTCCCGCTGTTTGCTGCCGGTTTTGCGTTGTCAGTGGAGTa                                                                                                                                                                                                                                         <  2:437279/216‑1 (MQ=255)
atccgcTACATCCTCAACCACGCGTTCGCTAAAAGCCTGTTTTTCCTTGTAGCAGGTGCGCTGAGTTACAGCTGCGGCACGCGCTTGTTGCCGCGTCTGCGTGGCGTATTGCACACCCTGCCGTTGCCAGGCGTGGGTTTCTGCGTGGCAGCGCTGGCGATTACCGGCGTGCCGCCGTTCAACGGCTTCTTCAGTAAATTCCCGCTGTTTGCTGCCGGTTTTGCGTTGTCAGTGGAGTa                                                                                                                                                                                                                                         <  2:481300/234‑1 (MQ=255)
  ctgtctactcttCAACCCCGCTGTCGCTAAAAGCCGTTTTTTCCTTGTAGCAGGGTGGCTGAGTTACAGCTGCGGCCCGGGCTTGTTTCCGCGTCTGCGTGGCGTATTGCAACCCCTGCCGTTGCCAGGCGTGGGTTTCTGCGTGGCAGCGCTGGCGATTACCGGCGTGCCGCCGTTCAACGTCTTCTCCAGTAAATTCCCGCTGTTTGCTGCCGGTTTTGCGTTGTCAGTGGAGTACt                                                                                                                                                                                                                                       <  1:863988/228‑1 (MQ=255)
               aaCCACGCGTTCGCTAAAAGCCTGTTTTTCCTTGTAGCAGGTGCGCTGAGTTACAGCTGCGGCACGCGCTTGTTGCCGCGTCTGCGTGGCGTCTTGCACACCCTGCCGTTGCCAGGCGTGGGTTTCTGCGTGGCAGCGCTGGCGATTACCGGCGTGCCGCCGTTCAACGGCTTTTTCCGTAAATTTCCCCCTTTTGTTGCCGGGTTTTCGTGTTCAGGGGGGGACTGGGGCCggctgc                                                                                                                                                                                                                           >  1:758874/1‑238 (MQ=255)
                  ggtgtaaaggagacaggCCTGTTTTTCCTTGTAGCGGGTGCGTTGAGTTACAGCTGCGGCCCGCGCTGGTTGCCGCGTCGGCGTGGCGTATTGCACCCCCTGCCGTTGCCAGGCGTGGGTTTCTGCGTGGCAGCGCTGGCGATTACCGGCGTGCCGCCGTTCAACGGCTTCTTCAGTAAATTCCCGCTGTTTGCTGCCGGTTTTGCGTTGTCAGTGGAGTACTGGATCCTGCTGCCCGc                                                                                                                                                                                                                       <  2:317570/223‑1 (MQ=255)
                   aCGCGTTCGCTAAAAGCCTGTTTTTCCTTGTAGCAGGTGCGCTGAGTTACAGCTGCGGCACGCGCTTGTTGCCGCGTCTGCGTGGCGTATTGCACACCCTGCCGTTGCCAGGCGTGGGTTTCTGCGTGGCAGCGCTGGCGATTACCGGCGTGCCGCCGTTCAACGGCTTCTTCAGTAAATTCCCCCTGTTTGCTGCCGGTTTTGCGTTGTCAGTGGAGTACTGGATCCTCCTTCCCGc                                                                                                                                                                                                                       >  2:405938/1‑238 (MQ=255)
                    cgcgTTCGCTAAAAGCCTGTTTTTCCTTGTAGCAGGTGCGCTGAGTTACAGCTGCGGCACGCGCTTGTTGCCGCGTCTGCGTGGCGTATTGCACACCCTGCCGTTGCCAGGCGTGGGTTTCTGCGTGGCAGCGCTGGCGATTACCGGCGTGCCGCCGGTCAACGGCTTCTTCAGTAAATTCCCGCTGTTTGCTGCCGGTTTTGCGTTTTCAGTGGGGTAATGGATCCTGCTGCCCGccc                                                                                                                                                                                                                     >  1:170295/1‑238 (MQ=255)
                    cgcgTTCGCTAAAAGCCTGTTTTTCCTTGTAGCAGGCGCGCTGAGTTACAGCTGCGGCACGCGCTTGTTGCCGCGTCTGCGTGGCGTATTGCACAACCTGCCGTTGCCAGGCGTGGGTTTCTGCGTGGCATGGCTGGCGCTTACCAGCGTGCCGCCCATTACAGGTGTCTTCAGTAAGTTCCCCCTGTTTGATTCCGCGTCTCGCTGGTCAGTGGAGtcatcgctccagtacccccccc                                                                                                                                                                                                                     >  1:579169/1‑218 (MQ=255)
                     ggtctCGCTAAACGCTTTTTTTTCCTTGGGGCAGGGGGCTTGAGTTCCCGCGGGGGCCGGCGCTTGTTCCCGCTTCTGGGTGGCTTATTGCCCCCCCTGCCGTTGCCGGGCTTGGGTTTCTGCGTGGCAGCGCTGGCGATTCCCGGGGTGCCGCCGTTCAACGGTTTCTTCAGTAAATTCCCGCTGTTTGCTGCCGGTTTTGCGTTGTCAGTGGAGTACGGGATCCTGCTGCCCGCCAt                                                                                                                                                                                                                    <  1:668852/235‑1 (MQ=255)
                     gcgTTCGCTAAAAGCCTGTTGTTCTTTTGAGCCGGTGCGCTGAGTTACAGCTGCGGCACGCGCTTGTTGCCGTGTCTGCGTGGCGTATTGCACACCCTGCCGTTGCCAGGCGTGGTTTTCTGCGTGGCAGCGCTGGCGATTACCGGCGTGCCGCCGTTCAACGGCTTCTTCAGTAAATTCCCGCTGTTTGCTGCCGGTTTTGCGTTGTCAGTGGAGTACTGGATCCTGCTGCCCGCCAt                                                                                                                                                                                                                    <  2:186824/239‑1 (MQ=255)
                       gTTCGCTAAAAGCCTGTTTTTCCTTGTAGCAGGTGCGCTGAGTTACAGCTGCGGCACGCGCTTGTTGCCGCGTCTGCGTGGCGTATTGCACACCCTGCCGTTGCCAGGCGTGGGGTTCTGCGTGGCAGCGCTGGCGATTACCGGCGTGACGCCGTTCAACGGCTTTTTAAGGAAATTTCCGCTTGTTGCTGCCGGGTTTGGGTTTTCAGGGGGGGACTGGGTCcttgggccccccagct                                                                                                                                                                                                                  >  1:73467/1‑225 (MQ=255)
                        ttCGCTAAAAGCCTGTTTTTCCTTGTAGCAGGTGCGCTGAGTTACAGCTGCGGCACGCGCTTGTTGCCGCGTCTGCGTGGCGTATTGCACACCCTGCCGTTGCCAGGCGTGGGTTTCTGCGTGGCAGCGCTGGCGATTACCGGCGTGCCGCCGTTCAACGGCTTCTTCAGTAAATTCCCGCTGTTTGCTGCCGGTTTTGCGTTTTCAGTGGGGTACTGGATCCTGCTGCCCGCCATGAt                                                                                                                                                                                                                 >  2:47756/1‑239 (MQ=255)
                            cTAAAAGCCTGTTTTTCCTTGTAGCAGGTGCGCTGAGTTACAGCTGCGGCACGCGCTTGTTGCCGCGTCTGCGTGGCGTATTGCACACCCTGCCGTTGCCAGGCGTGGGTTTCTGCGTGGCAGCGCTGGCGATTACCGGCGTGCCGCCGTTCAACGGCTTCTTCCGTAAATTCCCGCTGTTTGCTGCCGGTTTTCCGTTTTACGTGGGGTACTGGATCCTGCTGCCCGCCATGCttct                                                                                                                                                                                                              >  2:51169/1‑238 (MQ=255)
                                             cTTGTAGCAGGTGCGCTGAGTTACAGCTGCGGCACGCGCTTGTTGCCGCGTCTGCGTGGCGTATTGCACACCCTGCCGTTGCCAGGCGTGGGTTTCTGCGTGGCAGCGCTGGCGATTACCGGCGTGCCGCCGTTCAACGGCTTCTTCAGTAAATTCCCGCTGTTTGCTGCCCGTTTTGCGTTTTCAGTGGGGTACAGGATCCTGCTGGCCCCCCTGGttcttccggtctcttataaaca                                                                                                                                                                                            >  1:317570/1‑223 (MQ=255)
                                                 tAGCAGGTGCGCTGAGTTACAGCTGCGGCACGCGCTTGTTGCCGCGTCTGCGTGGCGTATTGCACACCCTGCCGTTGCCAGGCGTGGGTTTCTGCGTGGCAGCGCTGGCGATTACCGGCGTGCCGCCGTTCAACGGCTTCTTCAGTAAATTCCCGCTGTTTGCTGGCGGGTTTGCGTTGTCAGTGGAGGACTGGATCCTGCCGCCCGCCCAGATTCTTCTGTTGATGGaaaggtggggc                                                                                                                                                                                        >  1:390995/1‑230 (MQ=255)
                                                 gaGCAGGTGCGCTGAGTTACATCGGCGGCACGCGCTTGTTGCCGCGTCTGCGTGGCGTATTGCACACCCTGCCGTCGCCAGGCGTGGGTTTCTGCGTGGCAGCGCTGGCGATTACCGGCGTGCCGCCGTTCAACGGCTTCTTCAGTACATTCCCGCGGTTTTCTGCCGGGTTTTGGTTGTCCGTGGGGGACCGGATCCTCGTGCCCGCCAGCATTCTTCTGATTGTTGAAACGGTCGc                                                                                                                                                                                         >  1:166032/2‑238 (MQ=255)
                                                       gTGCGCTGAGTTACAGCTGCGGCACGCGCTTGTTGCCGCGTCTGCGTGGCGTATTGCACACCCTGCCGTTGCCAGGCGTGGGTTTCTGCGTGGCAGCGCTGGCGATTACCGGCGTGCCGCCGTTCAACGGCTTCTTCAGTAAATTCCCGCTTGTTGCTGCCGGTTTTGCGTTGTCACTGGAGTACTGGGTCCTGCTGCCCGCCATGATTCTCTTGATGATTGAATCGGTCGCCAGttt                                                                                                                                                                                   >  2:119109/1‑238 (MQ=255)
                                                          cgcTGAGTTCCGGGTGCGGCACGCGTTTGTTGCCGCGTCGGCGTGGCGTTTTGCACACCCTGCCGTTGCCGGGCGGGGGTTTCTGCGTGGCAGCGCTGGCGATTACCGGCGTGCCGCCGTTCAACGGCTTCTTCAGTAAATTCCCGCTGTTTGCTGCCGGTTTTGCGTTGTCAGTGGAGTACTGGATCCTGCTGCCCGCCATGATTCTTCTGATGATTGAATCGGTCGCCAGTTTCGc                                                                                                                                                                                <  2:703320/238‑1 (MQ=255)
                                                            cTGAGTTACAGCTGCGGCACGCGCTTGTTGCCGCGTCTGCGTGGCGTATTGCACACCCTGCCGTTGCCAGGCGTGGGTTTCTGCGTGGCAGCGCTGGCGATTACCGGCGTGCCGCCGTTCAACGGCTTCTTCATTAAATTCCCCCTGTTTGTCGCCGGTTTTGCGTTGTCCGGGGCGGACTGGGTCCTTGCGCCCCCCCTGCTTTTTTTTGTGGTTGAAACGGtgcgcaggttccccgg                                                                                                                                                                             >  1:881743/1‑224 (MQ=255)
                                                                 ttACAGCTGCGGCACGCGCTTGTTGCCGCGTCTGCGTGGCGTATTGCACACCCTGCCGTTGCCAGGCGTGGGTTTCTGCGTGGCAGCGATGGCGATTACCGGCGTGCCGCCGGTCAACTGCTTCTTCAGTAAATTACCGCTGGTTGCTGCACGTTTTGAGTTTTCAGTGGGGTACAGGGTCCGCGCGCCCGCACGCATTCTTCTGATGATTTAATCGGCGGCCGATTTCCCCCGGGTTa                                                                                                                                                                        >  1:598534/1‑239 (MQ=255)
                                                                      tcctggggcctgccTTGTTGGCCGGTTCGGCGTGGGGTTTTCCCCCCCCGTCCGTTGCCGGGGGTGGTTTTCGGCGTGGCAGCGCTGGCTATTCCCGGCGTGCCGCCGTTCAACGGTTTTTTCAGTAAATTCCCGCTGTTTGCTGCCGGTTTTGCGTTGTCAGTGGAGTACTGGATCCTGCTGCCCGCCATGATTCTTCTGATGATTGAATCGGTCGCCAGTTTCGCCTGGTTTATTCg                                                                                                                                                                    <  1:629396/226‑1 (MQ=255)
                                                                         gCGGCACGCGCTTGTTGCCGCGTCTGCGTGGCGTATTGCACACCCTGCCGTTGCCAGGCGTGGGTTTCTGCGTGGCAGCGCTGGCGATTACCGGCGTGCCGCCGTTCAACGGCTTCTTCAGTAAATTCCCGCTGTTTGCTGCCGGTTTTGCGTTGTCAGTGGAGTACTGGATCCTGCTGCCCGCCATGATTCTTCTGTTGATTGAATCGGTCGCCAGTTTCGCCTGGTTTATTCGCTgg                                                                                                                                                                >  2:480449/1‑239 (MQ=255)
                                                                         gCGGCACGCGCTTGTTGCCGCGTCTGCGTGGCGTATTGCACACCCTGCCGTTGCCAGGCGTGGGTTTCTGCGTGGCAGCGCTGGCGATTACCGGCGTGCCGCCGTTCAACGGCTTCTTCAGTAAATTCCCGCTGTTTGCTGCCGGTTTTGCGTTGTCAGTGGAGTACTGGATCCTGCTGCCCGCCACGATTCTTCTGATGATTGAATCGGTCGCCAGTTTCGCCTGGTTTATTCGCTgg                                                                                                                                                                >  2:257730/1‑239 (MQ=255)
                                                                                    ttgttgCCGCGTCTGCGTGGCGTATTGCACGCCCGGCCGTTGCCAGGCGTGGGTTTCTGCGTGGCAGCGCTGGCGATTACCGGCGTGCCGCCGTTCAACGGCTTCTTCAGTAAATTCCCGCTGTTTGCTGCCGGTTTTGCGTTGTCAGTGGAGTACTGGATCCTGCTGCCCGCCATGATTCTTCTGATGATTGAATCGGTCGCCAGTTTCGCCTGGTTTATTCGCTGGTTTGGTcgcg                                                                                                                                                      <  2:237337/238‑1 (MQ=255)
                                                                                    ttgttgCCGCGTCTGCGTGGCGTATTGCACACCCTGCCGTTGCCAGGCGTGGGTTTCTGCGTGGCAGCGCTGGCGATTACCGGCGTGCCGCGGTTCGACGGCTTCTTCAGTAAATTCCCGCTGTTTGCTGCCGGTTTTGCGTTGTCAGTGGAGTACTGGATCCTGCTGCCCGCCATGATTCTTCTGATGATTGAATCGGTCGCCAGTTTCGCCTGGTTTATTCGCTGGTTTGGTcgcg                                                                                                                                                      <  2:174546/238‑1 (MQ=255)
                                                                                               gTCTGCGTGGCGTATTGCACACCCTGCCGTTGCCAGGCGTGGGTTTCTGCGTGGCAGCGCTGGCGATTACCGGCGTGCCGCCGTTCAACGGCTTCTTCAGTAAATTCCCGCTGTTTGCTGCCGGTTTTGCGTTGTCAGTGGAGTACTGGATCCTGCTGCCCGCCATGATTCTTCTGGTGATTGAATCGGTCGCCAGTTTCGCCTGGGTTTTTCCCTGGTTTGGTTGCGTTGGCCCTggt                                                                                                                                           >  2:676441/1‑238 (MQ=255)
                                                                                                         cGTATTGCACACCCTGCCGTTGCCAGGCGTGGGTTTCTGCGTGGCAGCGCTGGCGATTACCGGCGTGCCGCCGTTCAACGGCTTCTTCAGTAAATTCCCGCTGTTTGCTGCCGGTTTTTCGTTGTCAGTGGAGTACTGGATCCTGCTGCCCGCCATGATTTTTCTGATGATTGAATCGGTCGCCAGTTTTGCCTGGTTTATTTCCTGGTTTGGTCCCGTTGTGCCTGGCAAACcggccg                                                                                                                                 >  2:256354/1‑235 (MQ=255)
                                                                                                              tGCACACCCTGCCGTTGCCAGGCGTGGGTTTCTGCGTGGCAGCGCTGGCGATTACCGGCGTGCCGCCGTTCAACGGCTTCTTCAGTAAATTCCCGCTGTTTGCTTCCGGTTTTGCGTTTTCAGTGGGGTACTGGGTTCTCCCGCCCGCCATGATTCTTCTGCTGATTGGATCGGTCGCCCGTTTCGCCTCGTTTTTTCGCGGGGTTGGTCCCGTGGTGGCCGGCGAACCGCGCGAGGcc                                                                                                                            >  1:547228/1‑239 (MQ=255)
                                                                                                                   ccccTGCCGTTGCCAGGCTTGGGTTTCTGCGTGGCAGCGCTGGCGATTGCCGGCGTGCCGCCGTTCAACGGCTTCTTCAGTAAATTCCCGCTGTTTGCTGCCGGTTTTGCGTTGTCAGTGGAGTACTGGATCCTGCTGCCCGCCATGATTCTTCTGATGATTGAATCGGTCGCCAGTTTCGCCTGGTTTATTCGCTGGTTTGGTCGCGTTGTGCCTGGCAAACCGAGCGAGGCCGTCGc                                                                                                                       <  2:166032/238‑1 (MQ=255)
                                                                                                                   acCCTGCCGTTGCCAGGCGTGGGTTTCTGCGTGGCAGCGCTGGCGATTACCGGCGTGCCGCCGTTCAACGGCTTCTTCAGTAACTTCCCGCTGTTTGCTGCCCGTTTTGCGTTGTCAGTGGAGTACTGGATCCTGCTGCCCGCCATGATTCTTCTGATTATTTAATCGGTTGCCCGTTTCGCCTGGTTTTTTGACTGGTTTGGTTGCGGTTTGCCTGGCAAACCGGGCGGGGcactccc                                                                                                                       >  1:402889/1‑233 (MQ=255)
                                                                                                                     ccTGCCGTTGCCAGGCGTGGGTTTCTGCGTGGCAGCGCTGGCGATTACCGGCGTGCCGCCGTTCAACGGCTTCTTCAGTAAATTCCCGCTGTTTGCTGCCGGTTTTGCGTTGTCAGTGGAGTACTGGATCCTGCTGCCCGCCCTGATTCTTCTGATGATTGAATCGGTCGCCAGTTTCGCCTGGTTTATTCGCTGGTTTGGTCGCGTTTGGCCTGGCAAACCCGCCGAGGCCGTCGCCg                                                                                                                     >  2:479324/1‑239 (MQ=255)
                                                                                                                       tGCCGTTGCCAGGCGTGGGTTTCTGCGTGGCAGCGCTGGCGATTACCGGCGTGCCGCCGTTCAACGGCTTCTTCAGTAAATTCCCGCTGTTTGCTGCCGGTTTTGCGTTGTCCGTGGAGTACTGGATCCTGCTGCCCGCCAGGCTTCTTCTGATGATTGAATCGGTGGCCAGTTTCGCCTGGTTTATTTGCTGGGTTGGTCGCGGTGGGCCTGGCGAACcgcgggggggcggcgcggca                                                                                                                   >  1:563753/1‑221 (MQ=255)
                                                                                                                        gCCGTTGCCAGGCGTGGGTTTCTGCGTGGCAGCGCTGGCGATTACCGGCGTGCCGCCGTTCAACGGCTTCTTCAGTAAATTCCCGCTGTTTGCTGCCGGTTTTGGGTTGTCAGTGGAGGACTGGATCCTGCTTCCCGCCATGATTCTTTTGGTTATTGAATCGGTCGCCAGGTTTGCCCGGGTTTTTTGCTGGTTTGGGTGCGTTTGGCCGTGCAAACCGGGCGGGGCCGTcccggagg                                                                                                                  >  1:675630/1‑232 (MQ=255)
                                                                                                                        gCCGTTGCCAGGCGTGGGTTTCTGCGTGGCAGCGCTGGCGATTACCGGCGTGCCGCCGTTCAACGGCTTCTTCAGTAAATTCCCGCTGTTTGCTGCCGGTTTTGCGTTGTCAGTGGAGTACTGGATCCTGCTGCCCGCCATGATTCTTCTGATGATTGAATCGGTCGCCAGTTTCGCCTGGTTTATTCGCTGATTTGGTCGCGTTGTGCCTGGCAAACCGCGCGAGGCCGTCGCCCAt                                                                                                                   >  2:737942/1‑238 (MQ=255)
                                                                                                                                         gTTTCTGCGTGGCAGCGCTGGCGATTACCGGCGTGCCGCCGTTCAACGGCTTCTTCAGTAAATTCCCGCTGTTTTCTGCCGGGTTTTCGTTTTCAGTGGAGTAATGGATCCTGCTGCCCGCCATGATTCTTCTTGTTATTGAATCGTTCGCCATTTTTTCCTGGTTTATTCGCTGGGTTTGTCGTGTTGTTCCTGGCAAACCGAGCGAGGCCGTCGCCGATTCCCCACCGCCGCCCgg                                                                                                  >  1:501868/1‑238 (MQ=255)
                                                                                                                                                       gcgcgggggtttcccgggttcccgcctTTCAACGGCTTTTTCAGTAAATCCCTGGTGTTTGCGGCCGGTTTGGCGTTGGCAGTGGAGTACGGGATCCTGTTGCCCGCCATGATTCTTCTGATGATTGAATCGGTCGCCAGTTTCGCCTGGTTTATTCGCTGGTTTGGTCGCGTTGTGCCTGGCAAACCGAGCGAGGCCGTCGCCGATGCCGCACCGCTGCCAGGATCAATGCGCCTGgt                                                                                  <  1:480449/218‑1 (MQ=255)
                                                                                                                                                        gtctggtgtgtcgcgggggggcgcctTTCAGCGGCTTCTTAAGTAACTTCCGCCTTTTGTCTCCCGGTTTTGCGTTTTCAGTGGAGTACTGGGTCCTGCTGCCCGCCATGATTCTTCTGATGATTGAATCGGTCGCCAGTTTCGCCTGGTTTATTCGCTGGTTTGGTCGCGTTGTGCCTGGCAAACCGAGCGAGGCCGTCGCCGATGCCGCCCCGCTGCCAGGATCAATGCGCCTGgt                                                                                  <  1:229924/217‑1 (MQ=255)
                                                                                                                                                        gggtGGGGATTCCCGGGGTGCCCGCGGTGCAGGGTTTTTTCGGGAAATTCCCGCGTGTGTCGGCCGGTTTTGGGTTGTCGGTGGGTTCGTGGTCCCGGCGGCCCGCCATGATTCTTCGGATGATTGAATCGGTCGCCAGTTTCGCCTGGTTTATTCGCTGGTTTGGTCGCGTTGTGCCTGGCAAACCGAGCGAGGCCGTCGCCGATGCCGCACCGCTGCCAGGATCAATGCGCCTGgt                                                                                  <  1:637993/235‑1 (MQ=255)
                                                                                                                                                        cgcTGGCGATTCCCGGCGTGCCGCCGTTCAACGGCTTCTTCAGTAAATTCCCGCTGTTTGCTGCCGGTTTTGCGTTGTCAGTGGAGTACTGGATCCTGCTGCCCGCCATGATTCTTCTGATGATTGAATCGGTCGCCAGTTTCGCCTGGTTTATTCGCTGGTTTGGTCGCGTTGTGCCTGGCAAACCGAGCGAGGCCGTCGCCGATGCCGCACCGCTGCCAGGATCAATGCGCCTGgt                                                                                  <  2:170295/238‑1 (MQ=255)
                                                                                                                                                        cgcTGGCGATTACCGGCGTGCCGCCGTTCAACGGCTTCTCCAGGAAATTCCCGCTGTTTGCTGCGGGTTTGGCGTTGTCAGTGGAGTACTGGATCCTGCTGCCCGCCATGATTCTTCTGATGATTGAATCGGTCGCCAGTTTCGCCTGGTTTATTCGCTGGTTTGGTCGCGTTGTGCCTGGCAAACCGAGCGAGGCCGTCGCCGATGCCGCACCGCTGCCAGGATCAATGCGCCTGgt                                                                                  <  2:758874/238‑1 (MQ=255)
                                                                                                                                                                       gCGTGCCGCCGTTCAACGGCTTCTTCAGGAAATTCCCGCTGTTTGCTGCCGGTTTTGCGTTGTCAGGGGAGTCATGGCTCCTGCTGCCCGCCATGATTCTTCTGATGATTGAATCGGTCGCCCGATTTGCCTGGTTTTTTCGCTGGGTTGGTCGCGTTGTGCCTTGCAAACCGGGCGAGGCCCTCGCCGGTTCCGCCCCCGTGCCGGGATCAATACGCCTCGGTGTGATTGTACTTAtt                                                                  >  1:763795/1‑239 (MQ=255)
                                                                                                                                                                                        ggtttttcaattAAATTCCCGTTGTTTGCTGCCGGTTTTGCGTTGTCAGTGGAGTACTGGTTCCTGCTGCCCGCCATGATTCTTCTGATGATTGAATCGGTCGCCAGTTTCGCCTGGTTTATTCGCTGGTTTGGTCGCGTTGTGCCTGGCAAACCGAGCGAGGCCGTCGCCGATGCCGCACCGCTGCCAGGATCAATGCGCCTGGTGTTGATTGTACTGATTGTGATGTCGCTGATTTc                                                 <  2:392428/228‑1 (MQ=255)
                                                                                                                                                                                                          ccGGTGTTTGCTGCCGGTTTTGCGTTGTCAGTGGGGTACTGGATCCTTCTGCCCGCCATGATTCTTCTGATGATTGAATCGGTCGCCAGTTTCGCCTGGTTTATTCGCTGGTTTGGTCGCGTTGTGCCTGGCAAACCGAGCGAGGCCGTCGCCGATGCCGCACCGCTGCCAGGATCAATGCGCCTGGTGTTGATTGTACTGATTGTGATGTCGCTGATTTCCAGCGTAATCGCCGCGAc                               <  2:881743/239‑1 (MQ=255)
                                                                                                                                                                                                             cTGTTTGCTGCCGGTTTTGCGTTGTCAGTGGAGTACTGGATCCTGCTGCCCGCCATGATTCTTCTGATGATTGAATCGGTCGCCAGTTTCGCCTGGTTTATTCGCTGGTTTGGTCGCGTTGTGCCTGGCAAACCGAGCGAGGCCGTCGCCGATGCCGCACCGCTGCCAGGATCAATGCGCCTGGTGTTGATTGTACTGATTGTGATGTCGCTGATTTCCAGCGAAATCGCCGCGACCTg                            >  2:769612/1‑239 (MQ=255)
                                                                                                                                                                                                              tttttGCTCGCGGTTTTGCTTTGTCAGTGGATTACTGGTTCCTGCTGCCCGCCATGATTCTTCGGATGATTGAATCGGTCGCCAGTTTCGCCTGGTTTATTCGCTTGTTTGGTCGCGTTGTGCCTGGCAAACCGAGCGAGGCCGTCGCCGATGCCGCACCGCTGCCAGGATCAATGCGCCTGGTGTTGATTGTACTGATTGTGATGTCGCTGATTTCCAGCGTAATCGCCGCGACCTgg                           <  2:192399/237‑1 (MQ=255)
                                                                                                                                                                                                              tGTTTGCTGCCGGTTTTGCGTTGTCAGGGGAGTACTGGATCCTGCTGCCCGCCATGATTCTTCTGATGATTGAATCGGTCGCCAGTTTCGCCTGGTTTATTCGCTGGTTTGGTCGCGTTGTGCCTGGCAAACCGAGCGAGGCCGTCGCCGATGCCGCACCGCTGCCAGGATCAATGCGCCTGGTGTTGATTGTACTGATTGTGATGTCGCTGATTTCCAGCGTAATCGCCGCGACCTgg                           <  2:675630/239‑1 (MQ=255)
                                                                                                                                                                                                              ggTTTGATGCCGGTTTTGCGTTTTCAGTGGAGGACTGGATCCGGCTGCCCGCCATGATTCTTCTGATGATTGAATCGGTCGGCAGTTTCGCCTGGTTTATTCGCTGGTTTGGTCGCGTTGTGCCTGGCAAACCGAGCGAGGCCGTCGCCGATGCCGCACCGCTGCCAGGATCAATGCGCCTGGTGTTGATTGTACTGATTGTGATGTCGCTGATTTCCAGCGTAATCGCCGCGACCTgg                           <  2:658329/238‑1 (MQ=255)
                                                                                                                                                                                                               ggttGCGGCCGGTTTTGGGTTGTCAGGGGAGTATTGGATCCTGCTGCCCGCCATGATTCTTCTGATGATTGAATCGGTCGCCAGTTTCGCCTGGTTTATTCGCTGGTTTGGTCGCGTTGTGCCTGGCAAACCGAGCGAGGCCGTCGCCGATGCCGCACCGCTGCCAGGATCAATGCGCCTGGTGTTGATTGTACTGATTGTGATGTCGCTGATTTCCAGCGTAATCGCCGCGACCTgg                           <  2:553989/236‑1 (MQ=255)
                                                                                                                                                                                                                ttttggtgtgggTGTTGCGTTGTCAGTCGTGTACTGGATCCTGTGGCCCGCCATGATTCTTCTGATGATTGAATCGGTCGCCAGTTTCGCCTGGTTTATTCGCTGGTTTGGTCGCGTTGTGCCTGGCAAACCGAGCGAGGCCGTCGCCGATGCCGCACCGCTGCCAGGATCAATGCGCCTGGTGTTGATTGTACTGATTGTGATGTCGCTGATTTCCAGCGTAATCGCCGCGACCTgg                           <  2:833660/228‑1 (MQ=255)
                                                                                                                                                                                                                tgttggcctgggtttttGCTTTGGCAGGGGAGTGCGGGTGCCTGTCGCCAGCCATGATTCTTCTGATGATTGAATCGGTCGCAGGTTTCCCCTGGTTTATTCGCTGGTTTGGTCGCGTTGTGCCTGGCAAACCGAGCGAGGCCGTCGCCGATGCCGCACCGCTGCCAGGATCAATGCGCCTGGTGTTGATTGTACTGATTGTGATGTCGCTGATTTCCAGCGTAATCGCCGCGACCTgg                           <  1:566240/226‑1 (MQ=255)
                                                                                                                                                                                                                  ttgtcgcCTGTTTTGCCTTGTCAGTGGAGTCCTGGATCCTTGTGCCCGCCATGATTCTTCTGATGATTGAATCGGTCGCCAGTTTCGCCTGGTTTATTCGCTGGTTTGGTCGCGTTGTGCCTGGCAAACCGAGCGAGGCCGTCGCCGATGCCGCACCGCTGCCAGGATCAATGCGCCTGGTGTTGATTGTACTGATTGTGATGTCGCTGATTTCCAGCGTAATCGCCGCGACCTGGTTg                        <  2:563753/234‑1 (MQ=255)
                                                                                                                                                                                                                  ttgcgtccGGTTTGGGTTTGTCGTGGGAGTCCTGGATCCTGCTGCCCGCCATGATTCTTCTGATGATTGAATCGGTCGCCAGTTTCGCCTGGTTTATTCGCTGGTTTGGTCGCGTTGTGCCTGGCAAACCGAGCGAGGCCGTCGCCGATGCCGCACCGCTGCCAGGATCAATGCGCCTGGTGTTGATTGTACTGATTGTGATGTCGCTGATTTCCAGCGTAATCGCCGCGACCTGGTTg                        <  2:383570/233‑1 (MQ=255)
                                                                                                                                                                                                                  tgctgcCGGTTTTGCGTTGTCAGTGGGGTACTGGATCCTGCTGCCCGCCATGATTCTTCTGATGATTGAATCGGTCGCCAGTTTCGCCTGGTTTATTCGCTGGTTTGGTCGCGTTGTGCCTGGCAAACCGAGCGAGGCCGTCGCCGATGCCGCACCGCTGCCAGGATCAATGCGACTGGTGTTGATTGTACTGATTGTGATGTCGCTGCTTTCCCGCGTAATCACCGCGACTGGGTTGc                       >  2:829948/1‑239 (MQ=255)
                                                                                                                                                                                                                  tgctgcCGGTTTTGCGTTGTCAGTGGAGTACTGGATCCTGCTGCCCGCCATGATTCTTCTGATGATTGAATCGGTCGCCAGTTTCGCCTGGTTTATTCGCTGGTTTGGTCGCGTTGTGCCTGTCAAACCGAGCGAGGCCGTCGCCGATGCCGCACCGCTGCCAGGATCAATGCGCCTGGTGTTGATTGTACTGTTTGTGATGTCCCTGATTTCCAGCGGAATCGCCCCGACCGGGTTGc                       >  1:553989/1‑239 (MQ=255)
                                                                                                                                                                                                                       ggcGGTTTTGCGTTGTAAGTGGGGGATTGGAGCCTGCTGCCCGCCATGATTCTTCGTATGATTGAGTCGGTCGCCAGTTTCGCCTGGTTTATTCGCTGGTTTGGTCGCGTTGTGCCTGGCAAACCGAGCGAGGCCGTCGCCGATGCCGCACCGCTGCCAGGATCAATGCGCCTGGTGTTGATTGTACTGATTGTGATGTCGCTGATTTCCAGCGTAATCGCCGCGACCTGGTTGCAGt                    <  2:501868/236‑1 (MQ=255)
                                                                                                                                                                                                                       gcgggTTTTGCGTTTCCATGGGATTCCTGGATCCTGCTGCCCGCCATGATTCTTCTGATGATTGAATCGGTCGCCAGTTTCGCCTGGTTTATTCGCTGGTTTGGTCGCGTTGTGCCTGGCAAACCGAGCGAGGCCGTCGCCGATGCCGCACCGCTGCCAGGATCAATGCGCCTGGTGTTGATTGTACTGATTGTGATGTCGCTGATTTCCAGCGTAATCGCCGCGACCTGGTTGCAGTa                   <  2:390995/236‑1 (MQ=255)
                                                                                                                                                                                                                              tttGCGTTGTCAGGGGAGTACTGGATCCTGCTGCCCGCCATGATTCTTCTGATGATTGAATCGGTCGCCAGTTTCGCCTGGTTTATTCGCTGGTTTGGTCGCGTTGTGCCTGGCAAACCGAGCGAGGCCGTCGCCGATGCCGCACCGCTGCCAGGATCAATGCGCCTGGTGTTGATTGTACTGATTGTGATGTCGCTGATTTCCAGCGTAATCGCCGCGACCTGGTTGCAGTAAGgaga             <  2:505570/239‑1 (MQ=255)
                                                                                                                                                                                                                                       tCAGTGGAGTACTGGATCCTGCTGCCCGCCATGTTTCTTCTGATGATTGAATCGGTCGCCAGTTTCGCCTGGTTTATTCGCTGGTTTGGTCGCGTTGTGCCTGGCAAACCGAGCGAGGCCGTCGCCGATGCCGCACCGCTGCCAGGATCAATGCGCCTGGGGTTGATTGTACTGATTGTGATGTCGCTTATTTCCAGCGTAATCGCCGCGCCCCGGTTGCAGTATGGAgattgtgcat     >  2:257618/1‑231 (MQ=255)
                                                                                                                                                                                                                                         aGTGGAGTTCTGGATCCTGCTTCCCGCCATTATTCTTCTGATGATTGAATCGGTCGCCAGTTTCGCCTGGTTTATTCGCTGGTTTGGTCGCGTTGTGCCTGGCAAACCGAGCGAGGCCGTCGCCGATGCCGCACCGCTGCCAGGATCAATGCGCCTGGTGTTGATTGTACTGATTGTGATGTCGCTGATTTCCAGCGTAATCGCCGCGACCTGGTTGCAGTAAGGAGATGATGAATGAc  <  2:73467/239‑1 (MQ=255)
                                                                                                                                                                                                                                         aGTGGAGTATGGTATCCTGCGGCCCGCCTTGATTCTTCTGGTGATTGAATCGGTCGCCAGTTTCGGCTGGTTTATTCGCTGGTTTGGTCGCGTTGTGCCTGGCAAACCGAGCGAGGCCGTCGCCGATGCCGCACCGCTGCCAGGATCAATGCGCCTGGTGTTGATTGTACTGATTGTGATGTCGCTGATTTCCAGCGTAATCGCCGCGACCTGGTTGCAGTAAGGAGATGATGAATGAc  <  2:388889/239‑1 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                          
ATCGCCTACATCGTCAACCACGCGTTCGCTAAAAGCCTGTTTTTCCTTGTAGCAGGTGCGCTGAGTTACAGCTGCGGCACGCGCTTGTTGCCGCGTCTGCGTGGCGTATTGCACACCCTGCCGTTGCCAGGCGTGGGTTTCTGCGTGGCAGCGCTGGCGATTACCGGCGTGCCGCCGTTCAACGGCTTCTTCAGTAAATTCCCGCTGTTTGCTGCCGGTTTTGCGTTGTCAGTGGAGTACTGGATCCTGCTGCCCGCCATGATTCTTCTGATGATTGAATCGGTCGCCAGTTTCGCCTGGTTTATTCGCTGGTTTGGTCGCGTTGTGCCTGGCAAACCGAGCGAGGCCGTCGCCGATGCCGCACCGCTGCCAGGATCAATGCGCCTGGTGTTGATTGTACTGATTGTGATGTCGCTGATTTCCAGCGTAATCGCCGCGACCTGGTTGCAGTAAGGAGATGATGAATGAC  >  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: