Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,608,367 T→C I484T (ATC→ACC)  hyfF → hydrogenase 4, membrane subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,608,3670TC100.0% 273.8 / NA 75I484T (ATC→ACC) hyfFhydrogenase 4, membrane subunit
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (35/40);  total (35/40)

CTGCGCGGGTGTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCATCTGTTTCGGAGAAACACAGTGAACGTTAATTCA  >  NC_000913/2608235‑2608501
                                                                                                                                    |                                                                                                                                      
cTGCGCGGGTGTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCACCAGGAt                                                                                                                                  >  1:696820/1‑139 (MQ=255)
cTGCGCGGGTGTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCACCAGGAt                                                                                                                                  >  1:185993/1‑139 (MQ=255)
     cgGGTGTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCACCAGGATCCTgg                                                                                                                             >  1:1534835/1‑139 (MQ=255)
      gggTGTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTG‑‑‑ACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCACCAGGATCCTggcggg                                                                                                                         <  2:741379/139‑1 (MQ=255)
         tgtTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCACCAGGATCCTggcggg                                                                                                                         <  1:1016559/139‑1 (MQ=255)
         tgtTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCACCAGGATCCTggcggg                                                                                                                         <  1:861302/139‑1 (MQ=255)
         tgtTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCACCAGGATCCTggcggg                                                                                                                         <  2:696820/139‑1 (MQ=255)
         tgtTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGACGATGGGAACGCATATTCCACAACCTGTCACCAGGATCCTggcggg                                                                                                                         <  2:371382/139‑1 (MQ=255)
          gtTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCACCAGGATCCTggcgggc                                                                                                                        <  2:691733/139‑1 (MQ=255)
           ttAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCACCAGGATCCTggcgggcg                                                                                                                       >  2:825848/1‑139 (MQ=255)
            tAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCACCAGGATCCTGGCGGgcgc                                                                                                                      <  2:1534835/139‑1 (MQ=255)
             aaTGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCACCAGGATCCTGGCGGGCGCt                                                                                                                     <  1:325211/139‑1 (MQ=255)
              aTGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCACCAGGATCCTGGCGGGCGCtt                                                                                                                    <  2:204443/139‑1 (MQ=255)
               tGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCACCAGGATCCTGGCGGGCGCTTc                                                                                                                   <  2:1214471/139‑1 (MQ=255)
               tGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCACCAGGATCCTGGCGGGCGCTTc                                                                                                                   <  1:825848/139‑1 (MQ=255)
               tGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCCGTCACCAGGATCCTGGCGGGCGCTTc                                                                                                                   <  1:745608/139‑1 (MQ=255)
                 gCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCACCAGGATCCTGGCGGGCGCTTCCa                                                                                                                 <  1:1633143/139‑1 (MQ=255)
                  cGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCACCAGGATCCTggcggg                                                                                                                         <  2:1232556/130‑1 (MQ=255)
                  cGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCACCAGGATCCTggcggg                                                                                                                         >  1:1232556/1‑130 (MQ=255)
                   gAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCACCAGGATCCTggcg                                                                                                                           <  2:1595997/127‑1 (MQ=255)
                   gAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCACCAGGATCCTggcg                                                                                                                           >  1:1595997/1‑127 (MQ=255)
                        cgccgcAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCACCAGGATCCTGGCGGGCGCTTCCACTATAGt                                                                                                          <  1:1242355/139‑1 (MQ=255)
                           cgcAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCACCAGGATCCTGGCGGGCGCTTc                                                                                                                   >  2:528439/1‑127 (MQ=255)
                           cgcAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCACCAGGATCCTGGCGGGCGCTTc                                                                                                                   <  1:528439/127‑1 (MQ=255)
                                            ggTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCACCAGGATCCTGGCGGgcgc                                                                                                                      >  1:869353/1‑107 (MQ=255)
                                            ggTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCACCAGGATCCTGGCGGgcgc                                                                                                                      <  2:869353/107‑1 (MQ=255)
                                                 tctcGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCACCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATctgcct                                                                                 >  1:1345777/1‑139 (MQ=255)
                                                     ggCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCACCAGGATCCTGGCGGGCGCTTCCACTATAGt                                                                                                          <  2:1598065/110‑1 (MQ=255)
                                                     ggCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCACCAGGATCCTGGCGGGCGCTTCCACTATAGt                                                                                                          >  1:1598065/1‑110 (MQ=255)
                                                          gTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCACCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGt                                                                        >  2:248134/1‑139 (MQ=255)
                                                          gTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCACCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGt                                                                        >  1:751134/1‑139 (MQ=255)
                                                             gACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCACCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGc                                                                     >  2:250334/1‑139 (MQ=255)
                                                                caccTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCACCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCAcc                                                                  <  2:751134/139‑1 (MQ=255)
                                                                   cTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCACCAGGATCCTGGCGGGCGCTTCCACTATAGTcc                                                                                                        >  2:1439521/1‑98 (MQ=255)
                                                                   cTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCACCAGGATCCTGGCGGGCGCTTCCACTATAGTcc                                                                                                        <  1:1439521/98‑1 (MQ=255)
                                                                     cGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCACCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCa                                                             >  2:212273/1‑139 (MQ=255)
                                                                     cGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCACCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCa                                                             >  1:286378/1‑139 (MQ=255)
                                                                             gTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCACCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGAttttt                                                     <  1:212273/139‑1 (MQ=255)
                                                                                aTTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCACCAGGATCCTGGCGGGCGCTTCCACTATAGTCctct                                                                                                     <  2:24005/88‑1 (MQ=255)
                                                                                aTTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCACCAGGATCCTGGCGGGCGCTTCCACTATAGTCctct                                                                                                     >  1:24005/1‑88 (MQ=255)
                                                                                aTTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCACCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGAt                                                                                       >  1:133846/1‑102 (MQ=255)
                                                                                aTTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCACCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGAt                                                                                       <  2:133846/102‑1 (MQ=255)
                                                                                 ttCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCACCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAg                                                                                                  >  2:17562/1‑90 (MQ=255)
                                                                                 ttCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCACCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAg                                                                                                  <  1:17562/90‑1 (MQ=255)
                                                                                   ctgctgGGCATGATGCTGGCGATGGGGACCCATATTCCACAGCCTGTCACCAGGATCCTGGCGGGCGCTTCCACTAAAGACCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCtt                                               <  2:520218/139‑1 (MQ=255)
                                                                                     gctgGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCACCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCa                                             >  2:1436211/1‑139 (MQ=255)
                                                                                     gctgGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCACCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCa                                             <  2:214887/139‑1 (MQ=255)
                                                                                     gctgGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCACCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCa                                             >  1:214887/1‑139 (MQ=255)
                                                                                      ctgGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCACCAGGATCCTGGCGGGCGCTTCCACtata                                                                                                            <  2:1336352/75‑1 (MQ=255)
                                                                                      ctgGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCACCAGGATCCTGGCGGGCGCTTCCACtata                                                                                                            >  1:1336352/1‑75 (MQ=255)
                                                                                          tCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCACCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCTGGCAc                                        <  2:519832/139‑1 (MQ=255)
                                                                                                   tGGCGATGGGAACGCATATTCCACAACCTGTCACCAGGATCCTGGCGGGCGCTTCCACTATAGTCctct                                                                                                     <  2:149362/69‑1 (MQ=255)
                                                                                                   tGGCGATGGGAACGCATATTCCACAACCTGTCACCAGGATCCTGGCGGGCGCTTCCACTATAGTCctct                                                                                                     >  1:149362/1‑69 (MQ=255)
                                                                                                     gCGATGGGAACGCATATTCCACAACCTGTCACCAGGATCCTGGCGGgcgc                                                                                                                      <  1:803314/50‑1 (MQ=255)
                                                                                                     gCGATGGGAACGCATATTCCACAACCTGTCACCAGGATCCTGGCGGgcgc                                                                                                                      >  2:803314/1‑50 (MQ=255)
                                                                                                          gggAACGCATATTCCACAACCTGTCACCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCATCTGTTTCGgag                        >  2:1661290/1‑139 (MQ=255)
                                                                                                          gggAACGCATATTCCACAACCTGTCACCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCATCTGTTTCGgag                        >  1:1127923/1‑139 (MQ=255)
                                                                                                            gAACGCATATTCCACAACCTGTCACCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAg                                                                                                  <  1:436192/63‑1 (MQ=255)
                                                                                                            gAACGCATATTCCACAACCTGTCACCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAg                                                                                                  >  2:436192/1‑63 (MQ=255)
                                                                                                                    atTCCACAACCTGTCACCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCATCTGTTTCGGAGAAACACAGTg              <  1:1661290/139‑1 (MQ=255)
                                                                                                                    atTCCACAACCTGTCACCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCATCTGTTTCGGAGAAACACAGTg              <  2:463415/139‑1 (MQ=255)
                                                                                                                             ccTGTCGCCAGGATCCTGGCGGGCGCTTCCACTATAGt                                                                                                          <  1:1074193/38‑1 (MQ=255)
                                                                                                                             ccTGTCGCCAGGATCCTGGCGGGCGCTTCCACTATAGt                                                                                                          >  2:1074193/1‑38 (MQ=255)
                                                                                                                              cTGTCACCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCATCTGTTTCGGAGAAAcaca                 <  2:65422/126‑1 (MQ=255)
                                                                                                                              cTGTCACCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCATCTGTTTCGGAGAAAcaca                 >  1:65422/1‑126 (MQ=255)
                                                                                                                                gTCACCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGAttttt                                                     >  2:59104/1‑88 (MQ=255)
                                                                                                                                gTCACCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGAttttt                                                     <  1:59104/88‑1 (MQ=255)
                                                                                                                                gTCACCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCa                                    <  2:1472706/105‑1 (MQ=255)
                                                                                                                                gTCACCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCa                                    >  1:1472706/1‑105 (MQ=255)
                                                                                                                                gTCACCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCATCTGttt                             <  2:1243729/112‑1 (MQ=255)
                                                                                                                                gTCACCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCATCTGttt                             >  1:1243729/1‑112 (MQ=255)
                                                                                                                                gTCACCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCATCTGTTTCGGAGAAACACAGTg              <  2:654326/127‑1 (MQ=255)
                                                                                                                                gTCACCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCATCTGTTTCGGAGAAACACAGTg              >  1:654326/1‑127 (MQ=255)
                                                                                                                                gTCACCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCATCTGTTTCGGAGAAACACAGTGAACGTTAATTCa  >  2:907976/1‑139 (MQ=255)
                                                                                                                                gTCACCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCATCTGTTTCGGAGAAACACAGTGAACGTTAATTCa  >  2:907977/1‑139 (MQ=255)
                                                                                                                                 tcaccaGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTc                                              >  2:330218/1‑94 (MQ=255)
                                                                                                                                 tcaccaGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTc                                              <  1:330218/94‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                      
CTGCGCGGGTGTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCATCTGTTTCGGAGAAACACAGTGAACGTTAATTCA  >  NC_000913/2608235‑2608501

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