Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,610,107 G→A E541K (GAG→AAG)  hyfG → hydrogenase 4, subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,610,1070GA99.5% 1403.4 / ‑5.0 364E541K (GAG→AAG) hyfGhydrogenase 4, subunit
Reads supporting (aligned to +/- strand):  ref base G (0/1);  major base A (150/211);  minor base T (1/1);  total (151/213)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 377 total aligned reads displayed.

GGCCGGTGTTGCGTTACATGCTGCGCGGCAATACCGTTTCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAGAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACGGTGAAATATCCCTTCGCGCCACTGGAGGTCAGCCCTGGCTTTCGCG  >  NC_000913/2609961‑2610254
                                                                                                                                                  |                                                                                                                                                   
ggCCGGTGTTGCGTTACATGCTGCGCGGCAATACCGTTTCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAaaga                                                                                                                                                  <  6:982538/150‑1 (MQ=255)
 ggcGGTGTTGCGTTACATGCTGCGCGGCAATACCGTTTCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATc                                                                                                                                                <  3:1070116/149‑1 (MQ=255)
   cGGTGTTGCGTTACATGCTGCGCGGCAATACCGTTTCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGa                                                                                                                                              <  3:229554/151‑1 (MQ=255)
    ggTGTTGCGTTACATGCTGCGCGGCAATACCGTTTCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGaa                                                                                                                                             <  8:378205/151‑1 (MQ=255)
           cGTTACATGCTGCGCGGCAATACCGTTTCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACg                                                                                                                                      <  3:1233122/151‑1 (MQ=255)
            gTTACATGCTGCGCGGCAATACCGTTTCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACgg                                                                                                                                     >  7:431828/1‑151 (MQ=255)
              tACATGCTGCGCGGCAATACCGTTTCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGc                                                                                                                                    <  6:816503/150‑1 (MQ=255)
                cATGCTGCGCGGCAATACCGTTTCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCAtt                                                                                                                                 <  6:1507657/151‑1 (MQ=255)
                  tgctgcGCGGCAATACCGTTTCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTg                                                                                                                                >  5:1014776/1‑150 (MQ=255)
                      gcgcGGCAATACCGTTTCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCg                                                                                                                            <  2:332931/150‑1 (MQ=255)
                      gcgcGGCAATACCGATTCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCg                                                                                                                            <  1:1039575/150‑1 (MQ=255)
                        gcgGCAATACCGTTTCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTaa                                                                                                                         <  2:171060/151‑1 (MQ=255)
                         cgGCAATACCGTTTCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAAc                                                                                                                        <  5:863867/151‑1 (MQ=255)
                          ggCAATAGCGTTTCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAAcc                                                                                                                       <  3:75173/151‑1 (MQ=255)
                          ggCAATACCGTTTCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGGGACGCTGGTAGATGTGCGCAAGGGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAAcc                                                                                                                       <  4:109223/151‑1 (MQ=255)
                           gCAATACCGTTTCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCg                                                                                                                      >  5:440360/1‑151 (MQ=255)
                            cAATACCGTTTCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGt                                                                                                                     <  5:680364/151‑1 (MQ=255)
                            cAATACCGTTTCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGt                                                                                                                     <  6:597387/151‑1 (MQ=255)
                            cAATACCGTTTCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGt                                                                                                                     <  3:6292/151‑1 (MQ=255)
                             aaTACCGTTTCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGtt                                                                                                                    >  1:853415/1‑151 (MQ=255)
                              aTACCGTTTCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTc                                                                                                                   <  6:130562/151‑1 (MQ=255)
                               tACCGTTTCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTcg                                                                                                                  >  2:1003172/1‑151 (MQ=255)
                                aCCGTTTCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCCTAACCGTTcgc                                                                                                                 >  4:1226736/1‑151 (MQ=255)
                                      tCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTg                                                                                                            >  8:1013607/1‑150 (MQ=255)
                                         gACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGaa                                                                                                          <  7:1279748/149‑1 (MQ=255)
                                         gACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAg                                                                                                         <  4:133036/150‑1 (MQ=255)
                                         gACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGt                                                                                                        <  7:295029/151‑1 (MQ=255)
                                          aCGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTa                                                                                                       <  7:185084/151‑1 (MQ=255)
                                          aCGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTa                                                                                                       <  8:799752/151‑1 (MQ=255)
                                                cGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGAc                                                                                                  >  8:104418/1‑150 (MQ=255)
                                                cGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACa                                                                                                 <  4:1595679/151‑1 (MQ=255)
                                                cGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACa                                                                                                 <  6:434450/151‑1 (MQ=255)
                                                 gCTGATTCTCGGTAGCCTTGATCCCTGCTACTCGTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACag                                                                                                <  3:1351375/151‑1 (MQ=255)
                                                 gCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGAc                                                                                                  >  6:1200740/1‑149 (MQ=255)
                                                 gCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACag                                                                                                <  7:591901/151‑1 (MQ=255)
                                                  cTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGAc                                                                                                  >  4:610658/1‑148 (MQ=255)
                                                  cTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGAc                                                                                                  >  2:1233606/1‑148 (MQ=255)
                                                    gATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACagaag                                                                                             <  6:936200/151‑1 (MQ=255)
                                                    gATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACagaa                                                                                              <  2:1481324/150‑1 (MQ=255)
                                                      ttATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACaga                                                                                               >  1:471471/1‑147 (MQ=255)
                                                        aTCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAGAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGc                                                                                         <  4:408231/151‑1 (MQ=255)
                                                          cGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTg                                                                                       >  8:1330589/1‑151 (MQ=255)
                                                           ggTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGa                                                                                      <  6:1330851/151‑1 (MQ=255)
                                                           ggTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGa                                                                                      <  2:139843/151‑1 (MQ=255)
                                                            gTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGaa                                                                                     <  4:706611/151‑1 (MQ=255)
                                                            gTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGaa                                                                                     <  5:602823/151‑1 (MQ=255)
                                                              aGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGt                                                                                   <  8:48915/151‑1 (MQ=255)
                                                                ccTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGtt                                                                                  >  8:930025/1‑150 (MQ=255)
                                                                ccTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTa                                                                                 <  1:1047092/151‑1 (MQ=255)
                                                                ccTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTa                                                                                 <  7:944552/151‑1 (MQ=255)
                                                                ccTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTa                                                                                 <  5:512728/151‑1 (MQ=255)
                                                                    gATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTg                                                                              >  7:1308497/1‑150 (MQ=255)
                                                                    gATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTg                                                                              <  5:123773/150‑1 (MQ=255)
                                                                     aTCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGaa                                                                            <  7:237511/151‑1 (MQ=255)
                                                                       cccTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGaaa                                                                           <  2:463053/150‑1 (MQ=255)
                                                                             tACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACtttt                                                                     >  1:763668/1‑150 (MQ=255)
                                                                             tACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACtatt                                                                     >  5:1074783/1‑150 (MQ=255)
                                                                                 ccTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCCTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATg                                                                  <  6:692425/149‑1 (MQ=255)
                                                                                    gTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATgcgc                                                               <  5:945749/149‑1 (MQ=255)
                                                                                    gTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATgcgc                                                               <  7:497570/149‑1 (MQ=255)
                                                                                     tACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGcc                                                            >  3:556557/1‑151 (MQ=255)
                                                                                        cgaccgTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCgg                                                          >  8:1518338/1‑150 (MQ=255)
                                                                                         gaccgTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGa                                                         >  5:1304510/1‑150 (MQ=255)
                                                                                          accgTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAAc                                                       <  1:1346040/151‑1 (MQ=255)
                                                                                             gtgtGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACcgc                                                    >  2:980041/1‑151 (MQ=255)
                                                                                              tgtgACGCTGGTAGATGTGAGCTAGCGCCAGTCAAAATCCGTGCCTTATAAAAAGATAGAACGCTACTTCATTGATCGTAACAGTTCTCCGCTGAAGTAAGTACAGAAGATGCTGTAGTTACTGAAAACTATTATGCGCGACGGAAcct                                                     <  3:1293293/149‑2 (MQ=255)
                                                                                               gtgACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGCTGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACcgcg                                                   <  8:19305/150‑1 (MQ=255)
                                                                                                 gACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACg                                                <  5:1165087/151‑1 (MQ=255)
                                                                                                       ggTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACGGTGaa                                           <  6:310388/150‑1 (MQ=255)
                                                                                                       ggTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACGGTGaa                                           >  8:420888/1‑150 (MQ=255)
                                                                                                        gTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACGGTGAAat                                         <  1:212031/151‑1 (MQ=255)
                                                                                                              gtgCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACGGTGAAATATCCCt                                   <  6:1117352/151‑1 (MQ=255)
                                                                                                                gcgcAAGCGCCATTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACGGTGAAATATCCCTTc                                 <  8:717398/151‑1 (MQ=255)
                                                                                                                gcgcAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACGGTGAAATATCCCTTc                                 <  8:1227940/151‑1 (MQ=255)
                                                                                                                gcgcAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACGGTGAAATATCCCTTc                                 <  5:1442338/151‑1 (MQ=255)
                                                                                                                 cgcAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTAGTGAAAACTATTATGCGCGCCGGAACCGCGACGGTGAAATATCCCTTcg                                >  3:1160638/1‑151 (MQ=255)
                                                                                                                 cgcAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACGGTGAAATATCCCTTcg                                >  7:1134577/1‑151 (MQ=255)
                                                                                                                  gcAAGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACGGTGAAATATCCCTTcgc                               <  4:297823/151‑1 (MQ=255)
                                                                                                                   cAAGCGCCAGTCAAAAACCGTGACGTATAAAAAGATCGAACGCTACTGCATTGATCGTAAACGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACGGTGAAATATCCCTTcgcg                              <  6:228871/151‑1 (MQ=255)
                                                                                                                    aaGCGCCAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACCGTGAAATATCCCTTCgcgc                             >  4:1114667/1‑151 (MQ=255)
                                                                                                                          cAGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACGGTGAAATATCCCTTCGCGCCACTgg                       >  1:1504348/1‑151 (MQ=255)
                                                                                                                          aaGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACGGTGAAATATCCCTTCGCGCCACTgg                       <  7:823481/150‑1 (MQ=255)
                                                                                                                           aGTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACTGCATTCATAGTAACCGTTCGCGGCTGAAGTAAGGACAGAAGATGCTGAACTTACAGAAAACTATTATGCCCGACCGAACCGCGACCGTCAAATATCCCATCGCGCCACTgg                       >  2:526966/1‑150 (MQ=255)
                                                                                                                            gTCAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACGGTGAAATATCCCTTCGCGCCACTGGAg                     <  7:363032/151‑1 (MQ=255)
                                                                                                                              cAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACGGTGAAATATCCCTTCGCGCCACTGGAGGt                   <  8:240514/151‑1 (MQ=255)
                                                                                                                              cAAAAACCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACGGTGAAATATCCCTTCGCGCCACTGGAGGt                   <  3:1574726/151‑1 (MQ=255)
                                                                                                                               aaaaaCCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACGGTGAAATATCCCTTCGCGCCACTGGAGGTc                  >  8:1401500/1‑151 (MQ=255)
                                                                                                                                aaaaCCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATAATGCGCGCCGGAACCGCGACGGTGAAATATCCCTTCGCGCCACAGGAGGTc                  >  3:902339/1‑150 (MQ=255)
                                                                                                                                 aaaCCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACGGTGAAATATCCCTTCGCGCCACTGGAGGTc                  >  1:185243/1‑149 (MQ=255)
                                                                                                                                 aaaCCGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACGGTGAAATATCCCTTCGCGCCACTGGAGGTCAg                <  8:604635/151‑1 (MQ=255)
                                                                                                                                   aCCGTGCCGTATAAAAAGATCGACCGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACGGTGAAATATCCCTTCGCGCCACTGGAGGTCAGcc              <  5:336011/151‑1 (MQ=255)
                                                                                                                                     cGTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACGGTGAAATATCCCTTCGCGCCACTGGAGGTCAGCCCt            <  4:1471683/151‑1 (MQ=255)
                                                                                                                                      gTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACGGTGAAATATCCCTTCGCGCCACTGGAGGTCAGCCCTg           <  3:407/151‑1 (MQ=255)
                                                                                                                                      gTGCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACACAAGATGCTGAAGTTACTGAAAACTATTATTCGAGCCGGAACCGCGACGGTGAAATATCCCTTCGCGCCACTGGAGGTCAGACCTg           >  1:702538/1‑151 (MQ=255)
                                                                                                                                        gCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACGGTGAAATATCCCTTCGCGCCACTGGAGGTCAGCCCTGGc         >  5:163921/1‑151 (MQ=255)
                                                                                                                                        gCCGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACGGTGAAATATCCCTTCGCGCCACTGGAGGTCAGACCTGGc         >  3:799547/1‑151 (MQ=255)
                                                                                                                                         ccGTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACGGTGAAATATCCCTTCGCGCCACTGGAGGTCAGCCCTGGCt        <  6:466749/151‑1 (MQ=255)
                                                                                                                                           gTATAAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACGGTGAAATATCCCTTCGCGCCACTGGAGGTCAGCCCTGGCttt      >  7:1492220/1‑151 (MQ=255)
                                                                                                                                              taAAAAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACGGTGAAATATCCCTTCGCGCCACTGGAGGTCAGCCCTGGCTTTcgc   >  2:823663/1‑151 (MQ=255)
                                                                                                                                               aaaaAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACGGTGAAATATCCCTTCGCGCCACTGGAGGTCAGCCCTGGCTTTcgcg  >  1:629279/1‑151 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GGCCGGTGTTGCGTTACATGCTGCGCGGCAATACCGTTTCTGACGCACCGCTGATTATCGGTAGCCTTGATCCCTGCTACTCCTGTACCGACCGTGTGACGCTGGTAGATGTGCGCAAGCGCCAGTCAAAAACCGTGCCGTATAAAGAGATCGAACGCTACGGCATTGATCGTAACCGTTCGCCGCTGAAGTAAGGACAGAAGATGCTGAAGTTACTGAAAACTATTATGCGCGCCGGAACCGCGACGGTGAAATATCCCTTCGCGCCACTGGAGGTCAGCCCTGGCTTTCGCG  >  NC_000913/2609961‑2610254

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: