Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,346,086 A→G I414T (ATC→ACC)  fumB ← anaerobic class I fumarate hydratase (fumarase B)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,346,0860AG100.0% 334.8 / NA 93I414T (ATC→ACC) fumBanaerobic class I fumarate hydratase (fumarase B)
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (49/44);  total (49/44)

ATGCTGCCGCCGTGGGATTGCAGCAGATCCACGTAGGAGTCCATACGGCCTGCGGTGGTTGGG‑CCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGATGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACGAGTGGATACCGGGTATTGCGAAAGCTGGGCGAGG  >  NC_000913/4345954‑4346217
                                                                                                                                     |                                                                                                                                   
aTGCTGCCGCCGTGGGATTGCAGCAGATCCACGTAGGAGTCCATACGGCCTGCGGTGGTTGGGCCCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTAc                                                                                                                                <  1:1188433/139‑1 (MQ=255)
   ctgcCGCCGTGGGATTGCAGCAGATCCACGTAGGAGTCCATACGGCCTGCGGTGGTTGGG‑CCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCg                                                                                                                            >  2:876001/1‑139 (MQ=255)
    tgcCGCCGTGGGATTGCAGCAGATCCACGTAGGAGTCCATACGGCCTGCGGTGGTTGGG‑CCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCgg                                                                                                                           <  2:1619976/139‑1 (MQ=255)
      ccgccgTGGGATTGCAGCAGATCCACGTAGGAGTCCATACGGCCTGCGGTGGTTGGG‑CCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGaa                                                                                                                         >  1:1616606/1‑139 (MQ=255)
      ccgccgTGGGATTGCAGCAGATCCACGTAGGAGTCCATACGGCCTGCGGTGGTTGGG‑CCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGaa                                                                                                                         >  1:354089/1‑139 (MQ=255)
       cgccgTGGGATTGCAGCAGATCCACGTAGGAGTCCATACGGCCTGCGGTGGTTGGG‑CCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAg                                                                                                                        <  1:806127/139‑1 (MQ=255)
        gccgTGGGATTGCAGCAGATCCACGTAGGAGTCCATACGGCCTGCGGTGGTTGGG‑CCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGt                                                                                                                       >  1:1273333/1‑139 (MQ=255)
        gccgTGGGATTGCAGCAGATCCACGTAGGAGTCCATACGGCCTGCGGTGGTTGGG‑CCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCATTGGTGTACTGCGGAAGt                                                                                                                       >  2:500961/1‑139 (MQ=255)
         ccgTGGGATTGCAGCAGATCCACGTAGGAGTCCATACGGCCTGCGGTGGTTGGG‑CCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGtt                                                                                                                      >  2:1226611/1‑139 (MQ=255)
           gTGGGATTGCAGCAGATCCACGTAGGAGTCCATACGGCCTGCGGTGGTTGGG‑CCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCt                                                                                                                    >  2:679144/1‑139 (MQ=255)
               gATTGCAGCAGATCCACGTAGGAGTCCATACGGCCTGCGGTGGTTGGG‑CCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTAc                                                                                                                <  2:592499/139‑1 (MQ=255)
                      gcagATCCACGTAGGAGTCCATACGGCCTGCGGTGGTTGGG‑CCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTc                                                                                                         >  2:1563921/1‑139 (MQ=255)
                           tCCACGTAGGAGTCCATACGGCCTGCGGTGGTTGGG‑CCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCa                                                                                                    <  2:1645201/139‑1 (MQ=255)
                            ccACGTAGGAGTCCATACGGCCTGCGGTGGTTGGG‑CCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAg                                                                                                   <  1:1437459/139‑1 (MQ=255)
                                   ggAGTCCATACGGCCTGCGGTGGTTGGG‑CCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTc                                                                                            >  2:1459147/1‑139 (MQ=255)
                                    gAGTCCATACGGCCTGCGGTGGTTGGG‑CCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCa                                                                                           >  1:395268/1‑139 (MQ=255)
                                    gAGTCCATACGGCCTGCGGTGGTTGGG‑CCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCa                                                                                           >  1:635071/1‑139 (MQ=255)
                                     aGTCCATACGGCCTGCGGTGGTTGGG‑CCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAg                                                                                                                        >  1:225718/1‑109 (MQ=255)
                                     aGTCCATACGGCCTGCGGTGGTTGGG‑CCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAg                                                                                                                        <  2:225718/109‑1 (MQ=255)
                                         cATACGGCCTGCGGTGGTTGGG‑CCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTg                                                                                      >  1:1187351/1‑139 (MQ=255)
                                           tACGGCCTGCGGTGGTTGGG‑CCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGc                                                                                    <  1:1164030/139‑1 (MQ=255)
                                           tACGGCCTGCGGTGGTTGGG‑CCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGc                                                                                    <  2:635071/139‑1 (MQ=255)
                                           tACGGCCTGCGGTGGTTGGG‑CCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGc                                                                                    >  2:1322759/1‑139 (MQ=255)
                                             cGGCCTGCGGTGGTTGGG‑CCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCgt                                                                                  <  2:410378/139‑1 (MQ=255)
                                             cGGCCTGCGGTGGTTGGG‑CCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCgt                                                                                  <  2:509127/139‑1 (MQ=255)
                                                 cTGCGGTGGTTGGG‑CCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGc                                                                              >  2:317861/1‑139 (MQ=255)
                                                  tGCGGTGGTTGGG‑CCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCa                                                                             <  2:354089/139‑1 (MQ=255)
                                                   gCGGTGGTTGGG‑CCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTAc                                                                                                                >  2:1208800/1‑103 (MQ=255)
                                                   gCGGTGGTTGGG‑CCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTAc                                                                                                                <  1:1208800/103‑1 (MQ=255)
                                                   gCGGTGGTTGGG‑CCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCgtg                                                                                 <  1:1322118/134‑1 (MQ=255)
                                                   gCGGTGGTTGGG‑CCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCgtg                                                                                 >  2:1322118/1‑134 (MQ=255)
                                                        ggtTGGG‑CCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATAtct                                                                       <  2:1616606/139‑1 (MQ=255)
                                                        ggtTGGG‑CCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATAtct                                                                       <  2:693707/139‑1 (MQ=255)
                                                          ttGGG‑CCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCg                                                                     <  2:878598/139‑1 (MQ=255)
                                                          ttGGG‑CCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCg                                                                     >  1:858745/1‑139 (MQ=255)
                                                          ttGGG‑CCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCg                                                                     <  1:317861/139‑1 (MQ=255)
                                                             gg‑CCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGcc                                                                  <  1:19370/139‑1 (MQ=255)
                                                             gg‑CCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGcc                                                                  <  1:1563921/139‑1 (MQ=255)
                                                              g‑CCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCgg                                                                                                                           <  1:993996/81‑1 (MQ=255)
                                                               gacAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCgg                                                                                                                           >  2:993996/3‑81 (MQ=255)
                                                                    gTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTAc                                                                                                                <  1:81419/87‑1 (MQ=255)
                                                                    gTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTAc                                                                                                                >  2:81419/1‑87 (MQ=255)
                                                                    gTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGcc                                                                  >  2:1200638/1‑133 (MQ=255)
                                                                    gTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGcc                                                                  <  1:1200638/133‑1 (MQ=255)
                                                                    gTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGAt                                                            >  1:1597576/1‑139 (MQ=255)
                                                                         ccTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGcc                                                                  <  2:1544534/128‑1 (MQ=255)
                                                                         ccTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGcc                                                                  >  1:1544534/1‑128 (MQ=255)
                                                                             aTGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCgg                                                                                                                           <  2:1302284/67‑1 (MQ=255)
                                                                             aTGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCgg                                                                                                                           >  1:1302284/1‑67 (MQ=255)
                                                                               ggATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCgg                                                 >  2:1374893/1‑139 (MQ=255)
                                                                                 aTAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCgt                                                                                  <  2:427463/104‑1 (MQ=255)
                                                                                 aTAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCgt                                                                                  >  1:427463/1‑104 (MQ=255)
                                                                                  tAACCGGCAGGGGTTTTCTCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCtt                                                                                              >  2:919321/1‑91 (MQ=255)
                                                                                  tAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCtt                                                                                              <  1:919321/91‑1 (MQ=255)
                                                                                  tAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGa                                              >  2:1276593/1‑139 (MQ=255)
                                                                                    aCCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCgtg                                                                                 >  2:686809/1‑102 (MQ=255)
                                                                                    aCCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCgtg                                                                                 <  1:686809/102‑1 (MQ=255)
                                                                                    aCCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGc                                            >  1:192929/1‑139 (MQ=255)
                                                                                        gCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACa                                        <  1:679144/139‑1 (MQ=255)
                                                                                        gCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACa                                        >  2:1658113/1‑139 (MQ=255)
                                                                                         cAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTg                                                                                                                                   >  1:130412/1‑47 (MQ=255)
                                                                                            gggTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCa                                                                                                    >  1:1190536/1‑75 (MQ=255)
                                                                                            gggTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCa                                                                                                    <  2:1190536/75‑1 (MQ=255)
                                                                                            gggTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACg                                    >  1:714024/1‑139 (MQ=255)
                                                                                             ggTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTAc                                                                                                                <  2:1373770/62‑1 (MQ=255)
                                                                                             ggTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTAc                                                                                                                >  1:1373770/1‑62 (MQ=255)
                                                                                               ttttCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACGAGt                                 <  1:1276593/139‑1 (MQ=255)
                                                                                               ttttCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACGAGt                                 <  2:395268/139‑1 (MQ=255)
                                                                                               ttttCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACGAGt                                 <  1:1374893/139‑1 (MQ=255)
                                                                                               ttttCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGACGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACGAGt                                 <  2:1187351/139‑1 (MQ=255)
                                                                                                tttCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAg                                                                                                                        >  2:184801/1‑51 (MQ=255)
                                                                                                tttCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAg                                                                                                                        <  1:184801/51‑1 (MQ=255)
                                                                                                      cGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACGAGTGGATAcc                          <  1:560215/139‑1 (MQ=255)
                                                                                                         aCCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACGAGTGGATACCggg                       >  1:36606/1‑139 (MQ=255)
                                                                                                         aCCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACGAGTGGATACCggg                       >  1:1154345/1‑139 (MQ=255)
                                                                                                           ccGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACGAGTGGATACCGGGTa                     <  2:1273333/139‑1 (MQ=255)
                                                                                                              cgTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGcc                                                                  <  2:869386/91‑1 (MQ=255)
                                                                                                              cgTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGcc                                                                  >  1:869386/1‑91 (MQ=255)
                                                                                                                        cGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACGAGTGGATACCGGGTATTGCGAAAGCTgg        >  1:951803/1‑139 (MQ=255)
                                                                                                                          ggTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACGAGTGGATACCGGGTATTGCGAAAGCTGGGc      >  2:140834/1‑139 (MQ=255)
                                                                                                                           gTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACGAGTGGATACCGGGTATTGCGAAAGCTGGGCg     >  1:1646572/1‑139 (MQ=255)
                                                                                                                             gATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGc                                                                                    >  1:407548/1‑58 (MQ=255)
                                                                                                                             gATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGc                                                                                    <  2:407548/58‑1 (MQ=255)
                                                                                                                             gATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACGAGTGGATACCGGGTATTGCGAAAGCTGGGCg     <  2:109433/137‑1 (MQ=255)
                                                                                                                             gATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACGAGTGGATACCGGGTATTGCGAAAGCTGGGCg     >  1:109433/1‑137 (MQ=255)
                                                                                                                             gATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACGAGTGGATACCGGGTATTGCGAAAGCTGGGCGAg   >  1:1213338/1‑139 (MQ=255)
                                                                                                                             gATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACGAGTGGATACCGGGTATTGCGAAAGCTGGGCGAg   >  2:1200635/1‑139 (MQ=255)
                                                                                                                             gATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACGAGTGGATACCGGGTATTGCGAAAGCTGGGCGAg   >  1:369412/1‑139 (MQ=255)
                                                                                                                              aTCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACGAGTGGATACCGGGTATTGCGAAAGCTGGGCGAgg  <  2:858745/139‑1 (MQ=255)
                                                                                                                                cTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGc                                                                                    >  2:100787/1‑55 (MQ=255)
                                                                                                                                cTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGc                                                                                    <  1:100787/55‑1 (MQ=255)
                                                                                                                                cTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCa                                                                             >  2:667589/1‑62 (MQ=255)
                                                                                                                                cTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCa                                                                             <  1:667589/62‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                   
ATGCTGCCGCCGTGGGATTGCAGCAGATCCACGTAGGAGTCCATACGGCCTGCGGTGGTTGGG‑CCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGATGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACGAGTGGATACCGGGTATTGCGAAAGCTGGGCGAGG  >  NC_000913/4345954‑4346217

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: ‑