Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,603,573 C→T P586S (CCG→TCG)  hyfB → hydrogenase 4, membrane subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,603,5730CT100.0% 306.7 / NA 80P586S (CCG→TCG) hyfBhydrogenase 4, membrane subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (44/36);  total (44/36)

ACTACTGGCGATGCCTTTGCTGCCTGGCCTTTACTGGCTGTGGTGTCGTTCGCGCCGCGCAGCGTTTCGTCGCACAGGAGATGCCTGGGCATGCGGCTACGGCTGGGAAAATGCGATGGCCCCGTCAGGCAATGGCGTGATGCAGCCGCTGCGTGTGGTCTTTTCTGCGCTATTTCGTCTACGACAACAGCTCGACCCTACGCTGAGGCTAAATAAAGGTCTTGCGCACGTCACCGCCAGGGCTCAGAGCACAGAACCCTTCTGGGATGAGCGGGTGATCCGCCCCATCGTG  >  NC_000913/2603428‑2603719
                                                                                                                                                 |                                                                                                                                                  
actactGGCGATGCCTTTGCTGCCTGGCCTTTACTGGCTGTGGTGTCGTTCGCGCCGCGCAGCGTTTCGTCGCACAGGAGATGCCTGGGCATGCGGCTACGGCTGGGAAAATGCGATGGCCCCGTCAGGCAATGGCGTGATGCAGTCGCt                                                                                                                                                <  1:186365/150‑1 (MQ=255)
 ctactGGCGATGCCTTTGCTGCCTGGCCTTTACTGGCTGTGGTGTCGTTCGCGCCGCGCAGCGTTTCGTCGCACAGGAGATGCCTGGGCATGCGGCTACGGCTGGGAAAATGCGATGGCCCCGTCAGGCAATGGCGTGATGCAGTCGCt                                                                                                                                                <  1:763314/149‑1 (MQ=255)
   actGGCGATGCCTTTGCTGCCTGGCCTTTACTGGCTGTGGTGTCGTTCGCGCCGCGCAGCGTTTCGTCGCACAGGAGATGCCTGGGCATGCGGCTACGGCTGGGAAAATGCGATGGCCCCGTCAGGCAATGGCGTGATGCAGTCGCTGCgt                                                                                                                                            <  2:499749/151‑1 (MQ=255)
   actGGCGATGCCTTTGCTGCCTGGCCTTTACTGGCTGTGGTGTCGTTCGCGCCGCGCAGCGTTTCGTCGCACAGGAGATGCCTGGGCATGCGGCTACGGCTGGGAAAATGCGATGGCCCCGTCAGGCAATGGCGTGATGCAGTCGCTGCg                                                                                                                                             >  1:328432/1‑150 (MQ=255)
     tGGCGATGCCTTTGCTGCCTGGCCTTTACTGGCTGTGGTGTCGTTCGCGCCGCGCAGCGTTTCGTCGCACAGGAGATGCCTGGGCATGCGGCTACGGCTGGGAAAATGCGATGGCCCCGTCAGGCAATGGCGTGATGCAGTCGCTGCgtgt                                                                                                                                          >  2:1256916/1‑151 (MQ=255)
      ggCGATGCCTTTGCTGCCTGGCCTTTACTGGCTGTGGTGTCGTTCGCGCCGCGCAGCGTTTCGTCGCACAGGAGATGCCTGGGCATGCGGCTACGGCTGGGAAAATGCGATGGCCCCGTCAGGCAATGGCGTGATGCAGTCGCTGCgtgt                                                                                                                                          <  2:390775/150‑1 (MQ=255)
      ggCGATGCCTTTGCTGCCTGGCCTTTACTGGCTGTGGTGTCGTTCGCGCCGCGCAGCGTTTCGTCGCACAGGAGATGCCTGGGCATGCGGCTACGGCTGGGAAAATGCGATGGCCCCGTCAGGCAATGGCGTGATGCAGTCGCTGCGtgtg                                                                                                                                         <  2:944446/151‑1 (MQ=255)
                   ctgcCTGGCCTTTACTGGCTGTGGTGTCGTTCGCGCCGCGCAGCGTTTCGTCGCACAGGAGATGCCTGGGCATGCGGCTACGGCTGGGAAAATGCGATGGCCCCGTCAGGCAATGGCGTGATGCAGTCGCTGCGTGTGGTCTTTTc                                                                                                                                 >  2:1169158/1‑146 (MQ=255)
                   ctgcCTGGCCTTTACTGGCTGTGGTGTCGTTCGCGCCGCGCAGCGTTTCGTCGCACAGGAGATGCCTGGGCATGCGGCTACGGCTGGGAAAATGCGATGGCCCCGTCAGGCAATGGCGTGATGCAGTCGCTGCGTGTGGTCTTTTCTgcgc                                                                                                                            >  2:891985/1‑151 (MQ=255)
                    tgcCTGGCCTTTACTGGCTGTGGTGTCGTTCGCGCCGCGCAGCGTTTCGTCGCACAGGAGATGCCTGGGCATGCGGCTACGGCTGGGAAAATGCGATGGCCCCGTCAGGCAATGGCGTGATGCAGTCGCTGCGTGTGGTCTTTTCTGCGCt                                                                                                                           >  2:1194760/1‑151 (MQ=255)
                    tgcCTGGCCTTTACTGGCTGTGGTGTCGTTCGCGCCGCGCAGCGTTTCGTCGCACAGGAGATGCCTGGGCATGCGGCTACGGCTGGGAAAATGCGATGGCCCCGTCAGGCAATGGCGTGATGCAGTCGCTGCGTGTGGTCTTTTCTGCGCt                                                                                                                           >  2:737074/1‑151 (MQ=255)
                     gcCTGGCCTTTACTGGCTGTGGTGTCGTTCGCGCCGCGCAGCGTTTCGTCGCACAGGAGATGCCTGGGCATGCGGCTACGGCTGGGAAAATGCGATGGCCCCGTCAGGCAATGGCGTGATGCAGTCGCTGCGTGTGGTCTTTTCTGCGCTa                                                                                                                          <  2:710477/151‑1 (MQ=255)
                      ccTGGCCTTTACTGGCTGTGGTGTCGTTCGCGCCGCGCAGCGTTTCGTCGCACAGGAGATGCCTGGGCATGCGGCTACGGCTGGGAAAATGCGATGGCCCCGTCAGGCAATGGCGTGATGCAGTCGCTGCGTGTGGTCTTTTCTGCGCTAt                                                                                                                         <  2:247644/151‑1 (MQ=255)
                         ggCCTTTACTGGCTGTGGTGTCGTTCGCGCCGCGCAGCGTTTCGTCGCGCAGGAGATGCCTGGGCATGCGGCTACGGCTGGGAAAATGCGATGGCCCCGTCAGGCAATGGCGTGATGCAGTCGCTGCGTGTGGTCTTTTCTGCGCTAtt                                                                                                                        <  2:423905/149‑1 (MQ=255)
                         ggCCTTTACTGGCTGTGGTGTCGTTCGCGCCGCGCAGCGTTTCGTCGCACAGGAGATGCCTGGGCATGCGGCTACGGCTGGGAAAATGCGATGGCCCCGTCAGGCAATGGCGTGATGCAGTCGCTGCGTGTGGTCTTTTCTGCGCTAttt                                                                                                                       <  1:418948/150‑1 (MQ=255)
                            cTTTACTGGCTGTGGTGTCGTTCGCGCCGCGCAGCGTTTCGTCGCACAGGAGATGCCTGGGCATGCGGCTACGGCTGGGAAAATGCGATGGCCCCGTCAGGCAATGGCGTGATGCAGTCGCTGCGGGGGGTCTTTTCTGCGCTATTTCg                                                                                                                     >  1:70771/1‑149 (MQ=255)
                              ttACTGGCTGTGGTGTCGTTCGCGCCGCGCAGCGTTTCGTCGCACAGGAGATGCCTGGGCATGCGGCTACGGCTGGGAAAATGCGATGGCCCCGTCAGGCAATGGCGTGATGCAGTCGCTGCGTGTGGTCTTTTCTGCGCTATTTCGTc                                                                                                                   <  2:923589/149‑1 (MQ=255)
                              ttACTGGCTGTGGTGTCGTTCGCGCCGCGCAGCGTTTCGTCGCACAGGAGATGCCTGGGCATGCGGCTACGGCTGGGAAAATGCGATGGCCCCGTCAGGCAATGGCGTGATGCAGTCGCTGCGTGTGGTCTTTTCTGCGCTATTTCGTCt                                                                                                                  >  2:598797/1‑150 (MQ=255)
                                   ggCTGTGGTGTCGTTCGCGCCGCGCAGCGTTTCGTCGCACAGGAGATGCCTGGGCATGCGGCTACGGCTGGGAAAATGCGATGGCCCCGTCAGGCAATGGCGTGATGCAGTCGCTGCGTGTGGTCTTTTCTGCGCTATTTCGTCTACGaca                                                                                                            <  2:684129/151‑1 (MQ=255)
                                   ggCTGTGGTGTCGTTCGCGCCGCGCAGCGTTTCGTCGCACAGGAGATGCCTGGGCATGCGGCTACGGCTGGGAAAATGCGATGGCCCCGTCAGGCAATGGCGTGATGCAGTCGCTGCGTGTGGTCTTTTCTGCGCTATTTCGTCTACGa                                                                                                              <  1:1070030/149‑1 (MQ=255)
                                       gtggtgTCGTTCGCGCCGCGCAGCGTTTCGTCGCACAGGAGATGCCTGGGCATGCGGCTACGGCTGGGAAAATGCGATGGCCCCGTCAGGCAATGGCGTGATGCAGTCGCTGCGTGTGGTCTTTTCTGCGCTATTTCGTCTACGacaac                                                                                                          <  2:1244507/149‑1 (MQ=255)
                                          gtgtCGTTCGCGCCGCGCAGCGTTTCGTCGCACAGGAGATGCCTGGGCATGCGGCTACGGCTGGGAAAATGCGATGGCCCCGTCAGGCAATGGCGTGATGCAGTCGCTGCGTGTGGTCTTTTCTGCGCTATTTCGTCTACGACAACAGCTc                                                                                                     <  1:725200/151‑1 (MQ=255)
                                            gtCGTTCGCGCCGCGCAGCGTTTCGTCGCACAGGAGATGCCTGGGCATGCGGCTACGGCTGGGAAAATGCGATGGCCCCGTCAGGCAATGGCGTGATGCAGTCGCTGCGTGTGGTCTTTTCTGCGCTATTTCGTCTACGACAACAGCTCg                                                                                                    >  2:87155/1‑150 (MQ=255)
                                              cgatCGCGCCGCGCAGCGTTTCGTCGCACAGGAGATGCCTGGGCATGCGGCTACGGCTGGGAAAATGCGATGGCCCCGTCAGGCAATGGCGTGATGCAGTCGCTGCGTGTGGTCTTTTCTGCGCTATTTCGTCTACGACAACAGCTCGAcc                                                                                                 <  2:691081/148‑1 (MQ=255)
                                                    cgccgcgccgCGTTTCGTCGCACAGGAGATGCCTGGGCATGCGGCTACGGCTGGGAAAATGCGATGGCCCCGTCAGGCAATGGCGAGATGCAGTCGCTGCGTGTGGTCTTTTCTGCGCTATTTCGTCTACGACAACAGCTCGACCCTACGc                                                                                           >  1:322845/1‑151 (MQ=255)
                                                    cgccgcgcAGCGTTTCGTCGCACAGGAGATGCCTGGGCATGCGGCTACGGCTGGGAAAATGCGATGGCCCCGTCAGGCAATGGCGTGATGCAGTCGCTGCGTGTGGTCTTTTCTGCGCTATTTCGTCTACGACAACAGCTCGACCCTACGc                                                                                           >  2:118727/1‑151 (MQ=255)
                                                       cgcgcAGCGTTTCGTCGCACAGGAGATGCCTGGGCATGCGGCTACGGCTGGGAAAATGCGATGGCCCCGTCAGGCAATGGCGTGATGCAGTCGCTGCGTGTGGTCTTTTCTGCGCTATTTCGTCTACGACAACAGCTCGACCCTACGCTGa                                                                                        <  2:265282/151‑1 (MQ=255)
                                                        gcgcAGCGTTTCGTCGCACAGGAGATGCCTGGGCATGCGGCTACGGCTGGGAAAATGCGATGGCCCCGTCAGGCAATGGCGTGATGCAGTCGCTGCGTGTGGTCTTTTCTGCGCTATTTCGTCTACGACAACAGCTCGACCCTACGCTg                                                                                         >  2:305715/1‑149 (MQ=255)
                                                           cAGCGTTTCGTCGCACAGGAGATGCCTGGGCATGCGGCTACGGCTGGGAAAATGCGATGGCCCCGTCAGGCAATGGCGTGATGCAGTCGCTGCGTGTGGTCTTTTCTGCGCTATTTCGTCTACGACAACAGCTCGACCCTACGCTGAGGCt                                                                                    <  2:1035707/151‑1 (MQ=255)
                                                            cgCGTTTCGTCGCAAAGGAGATGCCTGGGCATGCGGCTACGGCTGGGAAAATGCGATCGCCCCGTCAGGCAATGGCGCGATGCAGTCGCCGAGTGTGGTCTTTTCTGCGCTATTTCGTCTACGACAACCGCTCGACCCTACGCTGAg                                                                                       >  1:331247/2‑147 (MQ=255)
                                                             gCGTTTCGTCGCACAGGAGATGCCTGGGCATGCGGCTACGGCTGGGAAAATGCGATGGCCCCGTCAGGCAATGGCGTGATGCAGTCGCTGCGTGTGGTCTTTTCTGCGCTATTTCGTCTACGACAACAGCTCGACCCTACGCTGAGGCtaa                                                                                  <  1:50203/151‑1 (MQ=255)
                                                              cGTTTCGTCGCACAGGAGATGCCTGGGCATGCGGCTACGGCTGGGAAAATGCGATGGCCCCGTCAGGCAATGGCGTGATGCAGTCGCTGCGTGTGGTCTTTTCTGCGCTATTTCGTCTACGACAACAGCTCGACCCTACGCTGAGGCtaaa                                                                                 <  1:1049978/151‑1 (MQ=255)
                                                               gTTTCGTCGCACAGGAGATGCCTGGGCATGCGGCTACGGCTGGGAAAATGCGATGGCCCCGTCAGGCAATGGCGTGATGCAGTCGCTGCGTGTGGTCTTTTCTGCGCTATTTCGTCTACGACAACAGCTCGACCCTACGCTGAGGCtaaat                                                                                >  1:1283515/1‑151 (MQ=255)
                                                                 ttCGTCGCACAGGAGATGCCTGGGCATGCGGCTACGGCTGGGAAAATGCGATGGCCCCGTCAGGCAATGGCGTGATGCAGTCGCTGCGTGTGGTCTTTTCTGCGCTATTTCGTCTACGACAACAGCTCGACCCTACGCTGAGGCtaaata                                                                               <  2:1153537/150‑1 (MQ=255)
                                                                  tcgtcgCACAGGAGATGCCTGGGCATGCGGCTACGGCTGGGAAAATGCGATGGCCCCGTCAGGCAATGGCGTGATGCAGTCGCTGCGTGTGGTCTTTTCTGCGCTATTTCGTCTACGACAACAGCTCGACCCTACGCTGAGGCtaaataa                                                                              <  1:1006151/150‑1 (MQ=255)
                                                                   cgtcgCACAGGAGATGCCTGGGCATGCGGCTACGGCTGGGAAAATGCGATGGCCCCGTCAGGCAATGGCGTGATGCAGTCGCTGCGTGTGGTCTTTTCTGCGCTATTTCGTCTACGACAACAGCTCGACCCTACGCTGAGGCTAAATAAAg                                                                            <  2:1046241/151‑1 (MQ=255)
                                                                    gtcgCACAGGAGATGCCTGGGCATGCGGCTACGGCTGGGAAAATGCGATGGCCCCGTCAGGCAATGGCGTGATGCAGTCGCTGCGTGTGGTCTTTTCTGCGCTATTTCGTCTACGACAACAGCTCGACCCTACGCTGAGGCTAAATAAAgg                                                                           <  2:636714/151‑1 (MQ=255)
                                                                    gtcgCACAGGAGATGCCTGGGCATGCGGCTACGGCTGGGAAAATGCGATGGCCCCGTCAGGCAATGGCGTGATGCAGTCGCTGCGTGTGGTCTTTTCTGCGCTATTTCGTCTACGACAACAGCTCGACCCTACGCTGAGGCTAAATAAAgg                                                                           >  2:17107/1‑151 (MQ=255)
                                                                     tcgCACAGGAGATGCCTGGGCATGCGGCTACGGCTGGGAAAATGCGATGGCCCCGTCAGGCAATGGCGTGATGCAGTCGCTGCGTGTGGTCTTTTCTGCGCTATTTCGTCTACGACAACAGCTCGACCCTACGCTGAGGCTAAATAAAGGt                                                                          >  2:432530/1‑151 (MQ=255)
                                                                      cgCACAGGAGATGCCTGGGCATGCGGCTACGGCTGGGAAAATGCGATGGCCCCGTCAGGCAATGGCGTGATGCAGTCGCTGCGTGTGGTCTTTTCTGCGCTATTTCGTCTACGACAACAGCTCGACCCTACGCTGAGGCTAAATAAAGGTc                                                                         >  2:228018/1‑151 (MQ=255)
                                                                           aGGAGATGCCTGGGCATGCGGCTACGGCTGGGAAAATGCGATGGCCCCGTCAGGCAATGGCGTGATGCAGTCGCTGCGTGTGGTCTTTTCTGCGCTATTTCGTCTACGACAACAGCTCGACCCTACGCTGAGGCTAAATAAAGGTCTTgcg                                                                    >  1:928302/1‑151 (MQ=255)
                                                                           aGGAGATGCCTGGGCATGCGGCTACGGCTGGGAAAATGCGATGGCCCCGTCAGGCAATGGCGTGATGCAGTCGCTGCGTGTGGTCTTTTCTGCGCTATTTCGTCTACGACAACAGCTCGACCCTACGCTGAGGCTAAATAAAGGTCTTgcg                                                                    <  1:73355/151‑1 (MQ=255)
                                                                            ggAGATGCCTGGGCATGCGGCTTCGGCTGGGAAAATGCGATGGCCCCGTCAGGCAATGGCGTGATGCAGTCGCTGCGTGTGGTCTTTTCTGCGCTATTTCGTCTACGACAACAGCTCGACCCTACGCTGAGGCTAAATAAAGGTCTTgcgc                                                                   >  1:246379/1‑151 (MQ=255)
                                                                             gagaTGCCTGGGCATGCGGCTACGGCTGGGAAAATGCGATGGCCCCGTCAGGCAATGGCGTGATGCAGTCGCTGCGTGTGGTCTTTTCTGCGCTATTTCGTCTACGACAACAGCTCGACCCTACGCTGAGGCTAAATAAAGGTCTTGCGCa                                                                  <  2:743421/151‑1 (MQ=255)
                                                                              agaTGCCTGGGCATGCGGCTACGGCTGGGAAAATGCGATGGCCCCGTCAGGCAATGGCGTGATGCAGTCGCTGCGTGTGGTCTTTTCTGCGCTATTTCGTCTACGACAACAGCTCGACCCTACGCTGAGGCTAAATAAAGGTCTTGCGCAc                                                                 >  1:1078868/1‑151 (MQ=255)
                                                                                 tGCCTGGGCATGCGGCTACGGCTGGGAAAATGCGATGGCCCCGTCAGGCAATGGCGTGATGCAGTCGCTGCGTGTGGTCTTTTCTGCGCTATTTCGTCTACGACAACAGCTCGACCCTACGCTGAGGCTAAATAAAGGTCTTGCGCACGTc                                                              >  1:1083696/1‑151 (MQ=255)
                                                                                    cTGGGCATGCGGCTACGGCTGGGAAAATGCGATGGCCCCGTCAGGCAATGGCGTGATGCAGTCGCTGCGTGTGGTCTTTTCTGCGCTATTTCGTCTACGACAACAGCTCGACCCTACGCTGAGGCTAAATAAAGGTCTTGCGCACGTCAcc                                                           >  1:399867/1‑151 (MQ=255)
                                                                                     tGGGCATGCGGCTACGGCTGGGAAAATGCGATGGCCCCGTCAGGCAATGGCGTGATGCAGTCGCTGCGTGTGGTCTTTTCTGCGCTATTTCGTCTACGACAACAGCTCGACCCTACGCTGAGGCTAAATAAAGGTCTTGCGCACGTCACCg                                                          >  2:1235529/1‑151 (MQ=255)
                                                                                     tGGGCATGCGGCTACGGCTGGGAAAATGCGATGGCCCCGTCAGGCAATGGCGTGATGCAGTCGCTGCGTGTGGCCTTTTCTGCCCTATTTCGTCTACGACAACAGCTCGACCCTACGCTGAGGCTAAATAAAGGTCTTGCGCACGTCACCg                                                          >  1:96263/1‑151 (MQ=255)
                                                                                        gCATGCGGCTGCGGCTGGGAAAATGCGATGGCCCCGTCAGGCAATGGCGTGATGCAGTCGCTGCGTGTGGTCTTTTCTGCGCTATTTCGTCTACGACAACAGCTCGACCCTACGCTGAGGCTAAATAAAGGTCTTGCGCACGTCACCGCCa                                                       <  1:34043/151‑1 (MQ=255)
                                                                                        gCATGCGGCTACGGCTGGGAAAATGCGATGGCCCCGTCAGGCAATGGCGTGATGCAGTCGCTGCGTGTGGTCTTTTCTGCGCTATTTCGTCTACGACAACAGCTCGACCCTACGCTGAGGCTAAATAAAGGTCTTGCGCACGTCACCGCCa                                                       <  2:66914/151‑1 (MQ=255)
                                                                                         cATGCGGCTACGGCTGGGAAAATGCGATGGCCCCGTCAGGCAATGGCGTGATGCAGTCGCTGCGTGTGGTCTTTTCTGCGCTATTTCGTCTACGACAACAGCTCGACCCTACGCTGAGGCTAAATAAAGGTCTTGCGCACGTCACCGcc                                                        >  1:1016832/1‑149 (MQ=255)
                                                                                           tGCGGCTACGGCTGGGAAAATGCGATGGCCCCGTCAGGCAATGGCGTGATGCAGTCGCTGCGTGTGGTCTTTTCTGCGCTATTTCGTCTACGACAACAGCTCGACCCTACGCTGAGGCTAAATAAAGGTCTTGCGCACGTCACCGCCAg                                                      <  2:331273/149‑1 (MQ=255)
                                                                                             cGGCTACGGCTGGGAAAATGCGATGGCCCCGTCAGGCAATGGCGTGATGCAGTCGCTGCGTGTGGTCTTTTCTGCGCTATTTCGTCTACGACAACAGCTCGACCCTACGCTGAGGCTAAATAAAGGTCTTGCGCACGTCACCGCCAggg                                                    >  1:333778/1‑149 (MQ=255)
                                                                                             cGGCTACGGCTGGGAAAATGCGATGGCCCCGTCAGGCAATGGCGTGATGCAGTCGCTGCGTGTGGTCTTTTCTGCGCTATTTCGTCTACGACAACAGCTCGACCCTACGCTGAGGCTAAATAAAGGTCTTGCGCACGTCACCGCCAGGGCt                                                  >  2:251130/1‑151 (MQ=255)
                                                                                                cTACGGCTGGGAAAATGCGATGGCCCCGTCAGGCAATGGCGTGATGCAGTCGCTGCGTGTGGTCTTTTCTGCGCTATTTCGTCTACGACAACAGCTCGACCCTACGCTGAGGCTAAATAAAGGTCTTGCGCACGTCACCGCCAGGGCTc                                                 >  2:690236/1‑149 (MQ=255)
                                                                                                 tACGGCTGGGAAAATGCGATGGCCCCGTCAGGCAATGGCGTGATGCAGTCGCTGCGTGTGGTCTTTTCTGCGCTATTTCGTCTACGACAACAGCTCGACCCTACGCTGAGGCTAAATAAAGGTCTTGCGCACGTCACCGCCAGGGCTc                                                 >  1:404108/1‑148 (MQ=255)
                                                                                                  aCGGCTGGGAAAATGCGATGGCCCCGTCAGGCAATGGCGTGATGCAGTCGCTGCGTGTGGTCTTTTCTGCGCTATTTCGTCTACGACAACAGCTCGACCCTACGCTGAGGCTAAATAAAGGTCTTGCGCACGTCACCGCCAGGGCTCag                                               <  2:349289/149‑1 (MQ=255)
                                                                                                      cTGGGAAAATGCGATGGCCCCGTCAGGCAATGGCGTGATGCAGTCGCTGCGTGTGGTCTTTTCTGCGCTATTTCGTCTACGACAACAGCTCGACCCTACGCTGAGGCTAAATAAAGGTCTTGCGCACGTCACCGCCAGGGCTCAGAGcac                                          >  2:1200224/1‑150 (MQ=255)
                                                                                                       tGGGAAAATGCGATGGCCCCGTCAGGCAATGGCGTGATGCAGTCGCTGCGTGTGGTCTTTTCTGCGCTATTTCGTCTACGACAACAGCTCGACCCTACGCTGAGGCTAAATAAAGGTCTTGCGCACGTCACCGCCAGGGCTCAGAGcac                                          >  1:44928/1‑149 (MQ=255)
                                                                                                        gggAAAATGCGATGGCCCCGTCAGGCAATGGCGTGATGCAGTCGCTGCGTGTGGTCTTTTCTGCGCTATTTCGTCTACGACAACAGCTCGACCCTACGCTGAGGCTAAATAAAGGTCTTGCGCACGTCACCGCCAGGGCTCAGAGCACAg                                        <  2:821298/150‑1 (MQ=255)
                                                                                                         ggAAAATGCGATGGCCCCGTCAGGCAATGGCGTGATGCAGTCGCTGCGTGTGGTCTTTTCTGCGCTATTTCGTCTACGACAACAGCTCGACCCTACGCTGAGGCTAAATAAAGGTCTTGCGCACGTCACCGCCAGGGCTCAGAGCACAGa                                       >  1:479519/1‑150 (MQ=255)
                                                                                                          gAAAATGCGATGGCCCCGTCAGGCAATGGCGTGATGCAGTCGCTGCGTGTGGTCTTTTCTGCGCTATTTCGTCTACGACAACAGCTCGACCCTACGCTGAGGCTAAATAAAGGTCTTGCGCACGTCACCGCCAGGGCTCAGAGCACAGAAc                                     >  2:413467/1‑151 (MQ=255)
                                                                                                          gAAAATGCGATGGCCCCGTCAGGCAATGGCGTGATGCAGTCGCTGCGTGTGGTCTTTTCTGCGCTATTTCGTCTACGACAACAGCTCGACCCTACGCTGAGGCTAAATAAAGGTCTTGCGCACGTCACCGCCAGGGCTCAGAGCACAGAAc                                     <  1:6883/151‑1 (MQ=255)
                                                                                                            aaaTGCGATGGCCCCGTCAGGCAATGGCGTGATGCAGTCGCTGCGTGTGGTCTTTTCTGCGCTATTTCGTCTACGACAACAGCTCGACCCTACGCTGAGGCTAAATAAAGGTCTTGCGCACGTCACCGCCAGGGCTCAGAGCACAGAAccc                                   >  1:217906/1‑151 (MQ=255)
                                                                                                                   aTGGCCCCGTCAGGCAATGGCGTGATGCAGTCGCTGCGTGTGGTCTTTTCTGCGCTATTTCGTCTACGACAACAGCTCGACCCTACGCTGAGGCTAAATAAAGGTCTTGCGCACGTCACCGCCAGGGCTCAGAGCACAGAACCCTTCTggg                            >  2:744277/1‑151 (MQ=255)
                                                                                                                       ccccGTCAGGCAATGGCGTGATGCAGTCGCTGCGTGTGGTCTTTTCTGCGCTATTTCGTCTACGACAACAGCTCGACCCTACGCTGAGGCTAAATAAAGGTCTTGCGCACGTCACCGCCAGGGCTCAGAGCACAGAACCCTTCTGGGATGa                        <  2:893198/151‑1 (MQ=255)
                                                                                                                        cccGTCAGGCAATGGCGTGATGCAGTCGCTGCGTGTGGTCTTTTCTGCGCTATTTCGTCTACGACAACAGCTCGACCCTACGCTGAGGCTAAATAAAGGTCTTGCGCACGTCACCGCCAGGGCTCAGAGCACAGAACCCTTCTGGGATGAg                       <  1:525959/151‑1 (MQ=255)
                                                                                                                            tCAGGCAATGGCGTGATGCAGTCGCTGCGTGTGGTCTTTTCTGCGCTATTTCGTCTACGACAACAGCTCGACCCTACGCTGAGGCTAAATAAAGGTCTTGCGCACGTCACCGCCAGGGCTCAGAGCACAGAACCCTTCTGGGATGAGCggg                   >  2:464666/1‑151 (MQ=255)
                                                                                                                             cAGGCAATGGCGTGATGCAGTCGCTGCGTGTGGTCTTTTCTGCGCTATTTCGTCTACTACAACAGCTCGACCCTACGCTGAGGCTAAATAAAGGTCTTGCGCACGTCACCGCCAGGGCTCAGAGCACAGAACCCTTCTGGGATGATCGGGt                  >  1:1248575/1‑151 (MQ=255)
                                                                                                                               ggCAATGGCGTGATGCAGTCGCTGCGTGTGGTCTTTTCTGCGCTATTTCGTCTACGACAACAGCTCGACCCTACGCTGAGGCTAAATAAAGGTCTTGCGCACGTCACCGCCAGGGCTCAGAGCACAGAACCCTTCTGGGATGAGCGGGTGa                <  2:719330/151‑1 (MQ=255)
                                                                                                                                 cAATGGCGTGATGCAGTCGCTGCGTGTGGTCTTTTCTGCGCTATTTCGTCTACGACAACAGCTCGACCCTACGCTGAGGCTAAATAAAGGTCTTGCGCACGTCACCGCCAGGGCTCAGAGCACAGAACCCTTCTGGGATGAGCGGGTGATc              >  1:905336/1‑151 (MQ=255)
                                                                                                                                  aaTGGCGTGATGCAGTCGCTGCGTGTGGTCTTTTCTGCGCTATTTCGTCTACGACAACAGCTCGACCCTACGCTGAGGCTAAATAAAGGTCTTGCGCACGTCACCGCCAGGGCTCAGAGCACAGAACCCTTCTGGGATGAGCGGGTGATcc             >  1:523608/1‑151 (MQ=255)
                                                                                                                                  aaTGGCGTGATGCAGTCGCTGCGTGTGGTCTTTTCTGCGCTATTTCGTCTACGACAACAGCTCGACCCTACGCTGAGGCTAAATAAAGGTCTTGCGCACGTCACCGCCAGGGCTCAGAGCACAGAACCCTTCTGGGATGAGCGGGTGATcc             >  1:212133/1‑151 (MQ=255)
                                                                                                                                       cGTGATGCAGTCGCTGCGTGTGGTCTTTTCTGCGCTATTTCGTCTACGACAACAGCTCGACCCTACGCTGAGGCTAAATAAAGGTCTTGCGCACGTCACCGCCAGGGCTCAGAGCACAGAACCCTTCTGGGATGAGCGGGTGATCCGcccc        >  2:256780/1‑151 (MQ=255)
                                                                                                                                         tGATGCAGTCGCTGCGTGTGGTCTTTTCTGCGCTATTTCGTCTACGACAACAGCTCGACCCTACGCTGAGGCTAAATAAAGGTCTTGCGCACGTCACCGCCAGGGCTCAGAGCACAGAACCCTTCTGGGATGAGCGGGTGATCCGCCCCAt      >  2:734943/1‑151 (MQ=255)
                                                                                                                                             gCAGTCGCTGCGTGTGGTCTTTTCTGCGCTATTTCGTCTACGACAACAGCTCGACCCTACGCTGAGGCTAAATAAAGGTCTTGCGCACGTCACCGCCAGGGCTCAGAGCACAGAACCCTTCTGGGATGAGCGGGTGATCCGCCCCATCGTg  <  2:351888/151‑1 (MQ=255)
                                                                                                                                              cAGTCGCTGCGTGTGGTCTTTTCTGCGCTATTTCGTCTACGACAACAGCTCGACCCTACGCTGAGGCTAAATAAAGGTCTTGCGTACGTCACCGCCAGGGCTCAGAGCACAGAGCCCTTCTGGGATCCGCGGGTGATCCCCACCATCGt   <  1:82516/149‑1 (MQ=255)
                                                                                                                                              cAGTCGCTGCGTGTGGTCTTTTCTGCGCTATTTCGTCTACGACAACAGCTCGACCCTACGCTGAGGCTAAATAAAGGTCTTGCGCACGTCACCGCCAGGGCTCAGAGCACAGAACCCTTCTGGGATGAGCGGGTGATCCGCCCCATCGt   >  1:640468/1‑149 (MQ=255)
                                                                                                                                              cAGTCGCTGCGTGTGGTCTTTTCTGCGCTATTTCGTCTACGACAACAGCTCGACCCTACGCTGAGGCTAAATAAAGGTCTTGCGCACGTCACCGCCAGGGCTCAGAGCACAGAACCCTTCTGGGATGAGCGGGTGATCCGCCCCATCGTg  >  2:336346/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
ACTACTGGCGATGCCTTTGCTGCCTGGCCTTTACTGGCTGTGGTGTCGTTCGCGCCGCGCAGCGTTTCGTCGCACAGGAGATGCCTGGGCATGCGGCTACGGCTGGGAAAATGCGATGGCCCCGTCAGGCAATGGCGTGATGCAGCCGCTGCGTGTGGTCTTTTCTGCGCTATTTCGTCTACGACAACAGCTCGACCCTACGCTGAGGCTAAATAAAGGTCTTGCGCACGTCACCGCCAGGGCTCAGAGCACAGAACCCTTCTGGGATGAGCGGGTGATCCGCCCCATCGTG  >  NC_000913/2603428‑2603719

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: