Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,421,553 G→A E77K (GAA→AAA)  folX → D‑erythro‑7,8‑dihydroneopterin triphosphate 2'‑epimerase and dihydroneopterin aldolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,421,5530GA100.0% 347.6 / NA 98E77K (GAA→AAA) folXD‑erythro‑7,8‑dihydroneopterin triphosphate 2'‑epimerase and dihydroneopterin aldolase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (50/48);  total (50/48)

GTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGCACGTGAACATCACTGGGTGACGTATGCTGAAGTGGAGATCGATAAACTGCACGCGCTGCGCTACGCCGATTCGGTATCCATGACCTTAAGCTGGCAGCGTT  >  NC_000913/2421421‑2421685
                                                                                                                                    |                                                                                                                                    
gTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGAAAAAAt                                                                                                                                >  2:1035226/1‑139 (MQ=255)
 ttATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGAAAAAAtt                                                                                                                               >  2:1086141/1‑139 (MQ=255)
    tCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGAAAAAATTAAc                                                                                                                            >  1:829683/1‑139 (MQ=255)
    tCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGAAAAAATTAAc                                                                                                                            >  1:51004/1‑139 (MQ=255)
     cAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGAAAAAATTAACt                                                                                                                           <  2:53442/139‑1 (MQ=255)
     cAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGAAAAAATTAACt                                                                                                                           <  2:390958/139‑1 (MQ=255)
        tgtgACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGAAAAAATTAACTCAg                                                                                                                        >  2:842835/1‑139 (MQ=255)
         gtgACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGAAAAAATTAACTCAgg                                                                                                                       <  2:936018/139‑1 (MQ=255)
         ataaCGACCCACTACCCCGCAGATAAAGCGCTCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTTAGCATGTAGAGAATAACCGTTTCTCTTTGCTGAAAAAATTAACTCAgg                                                                                                                       <  1:210790/136‑1 (MQ=255)
           gACGATCCATTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGAAAAAATTAACTCAGGAt                                                                                                                     <  1:943932/139‑1 (MQ=255)
           gACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGAAAAAATTAACTCAGGAt                                                                                                                     <  2:167318/139‑1 (MQ=255)
            aCGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGAAAAAATTAACTCAGGAtg                                                                                                                    >  2:250138/1‑139 (MQ=255)
                tCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGAAAAAATTAACTCAGGATGTGCt                                                                                                                <  1:842835/139‑1 (MQ=255)
                 ccACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGAAAAAATTAACTCAGGATGTGCTc                                                                                                               <  1:241394/139‑1 (MQ=255)
                  cACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCCTGTAGAGAATAACCGTTTCTCTTTGCTGAAAAAATTAACTCAGGAt                                                                                                                     >  1:480939/1‑132 (MQ=255)
                  cACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGaaaaaa                                                                                                                                 >  1:1016349/1‑120 (MQ=255)
                  cACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGaaaaaa                                                                                                                                 <  2:1016349/120‑1 (MQ=255)
                  cACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGAAAAAATTAACTCAGGAt                                                                                                                     <  1:33524/132‑1 (MQ=255)
                  cACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGAAAAAATTAACTCAGGAt                                                                                                                     >  2:33524/1‑132 (MQ=255)
                  cACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGAAAAAATTAACTCAGGAt                                                                                                                     <  2:480939/132‑1 (MQ=255)
                   aCTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGTACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGAAAAAATTAACTCAGGATGTGCTCGa                                                                                                             >  2:833852/1‑139 (MQ=255)
                       ccccGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGAAAAAATTAACTCAGGATGTGCTCGATATc                                                                                                         >  2:560921/1‑139 (MQ=255)
                          cgccgATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGAAAAAATTAACTCAGGAtgt                                                                                                                   >  2:281442/1‑126 (MQ=255)
                          cgccgATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGAAAAAATTAACTCAGGAtgt                                                                                                                   <  1:281442/126‑1 (MQ=255)
                          cgccgATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGAAAAAATTAACTCAGGATGTGCTCGATATCGCa                                                                                                      >  2:357623/1‑139 (MQ=255)
                          cgccgATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGAAAAAATTAACTCAGGATGTGCTCGATATCGCa                                                                                                      >  1:887814/1‑139 (MQ=255)
                          cgccgATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGAAAAAATTAACTCAGGATGTGCTCGATATCGCa                                                                                                      >  2:791361/1‑139 (MQ=255)
                          cgccgATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGAAAAAATTAACTCAGGATGTGCTCGATATCGCa                                                                                                      >  2:1050472/1‑139 (MQ=255)
                          cgccgATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGAAAAAATTAACTCAGGATGTGCTCGATATCGCa                                                                                                      >  2:275986/1‑139 (MQ=255)
                          cgccgATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGAAAAAATTAACTCAGGATGTGCTCGATATCGCa                                                                                                      >  1:299772/1‑139 (MQ=255)
                           gccgATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGAAAAAATTGACTCAGGATGTGCTCGATATCGCAc                                                                                                     <  1:102821/139‑1 (MQ=255)
                           gccgATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGAAAAAATTAACTCAGGATGTGCTCGATATCGCAc                                                                                                     <  2:985577/139‑1 (MQ=255)
                                tAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGTACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGAAAAAATTAACTCAGGATGTGCTCGATATCGCACGTGaa                                                                                                <  1:833852/139‑1 (MQ=255)
                                tAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGAAAAAATTAACTCAGGATGTGCTCGATATCGCACGTGaa                                                                                                <  2:436567/139‑1 (MQ=255)
                                     cgcgcACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGAAAAAATTAACTCAGGATGTGCTCGATATCGCACGTGAACATCa                                                                                           <  2:299772/139‑1 (MQ=255)
                                     cgcgcACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGAAAAAATTAACTCAGGATGTGCTCGATATCGCACGTGAACATCa                                                                                           <  1:1086141/139‑1 (MQ=255)
                                     cgcgcACTAGCGAAGATATCAACGATGCGCTGAATTATCGCAACGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGAAAAAATTAACTCAGGATGTGCTCGATATCGCACGTGAACATCa                                                                                           <  1:787008/139‑1 (MQ=255)
                                       cgcACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGAAAAAATTAACTCAGGATGTGCTCGATATCGCACGTGAACATCACt                                                                                         >  2:999195/1‑139 (MQ=255)
                                         cACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGAAAAAATTAACTCAGGATGTGCTCGATATCGCACGTGAACATCACTgg                                                                                       >  2:422041/1‑139 (MQ=255)
                                         cACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGAAAAAATTAACTCAGGATGTGCTCGATATCGCACGTGAACATCACTgg                                                                                       <  2:559228/139‑1 (MQ=255)
                                               cGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGAAAAAATTAACTCAGGAt                                                                                                                     >  1:1019225/1‑103 (MQ=255)
                                               cGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGAAAAAATTAACTCAGGAt                                                                                                                     <  2:1019225/103‑1 (MQ=255)
                                                 aaGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGAAAAAATTAACTCAGGATGTGCTCGATATCGCACGTGAACATCACTGGGTGACGTa                                                                               <  1:560921/139‑1 (MQ=255)
                                                  aGATATCAACGATGCGCGGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGAAAAAATTAACTCAGGATGTGCTCGATATCGCACGTGAACATCACTGGGTGACGTAt                                                                              >  2:512725/1‑139 (MQ=255)
                                                      atCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGAAAAAATTAACTCAGGATGTGCTCGATATCGCACGTGaa                                                                                                >  1:445568/1‑117 (MQ=255)
                                                      atCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGAAAAAATTAACTCAGGATGTGCTCGATATCGCACGTGaa                                                                                                <  2:445568/117‑1 (MQ=255)
                                                      atCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGAAAAAATTAACTCAGGATGTGCTCGATATCGCACGTGAACATCACTGGGTGACGTATGCTg                                                                          >  2:424945/1‑139 (MQ=255)
                                                      atCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGAAAAAATTAACTCAGGATGTGCTCGATATCGCACGTGAACATCACTGGGTGACGTATGCTg                                                                          >  2:81028/1‑139 (MQ=255)
                                                        cAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGAAAAAATTAACTCAGGATGTGCTCGATATCGCACGTGAACATCACTGGGTGACGTATGCTg                                                                          >  1:968278/1‑137 (MQ=255)
                                                        cAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGAAAAAATTAACTCAGGATGTGCTCGATATCGCACGTGAACATCACTGGGTGACGTATGCTg                                                                          <  2:968278/137‑1 (MQ=255)
                                                        cAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGAAAAAATTAACTCAGGATGTGCTCGATATCGCACGTGAACATCACTGGGTGACGTATGCTGaa                                                                        >  1:74341/1‑139 (MQ=255)
                                                         aaCGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGAAAAAATTAACTCAGGATGTGCTCGATATCGCACGTGAACATc                                                                                            <  2:113896/118‑1 (MQ=255)
                                                         aaCGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGAAAAAATTAACTCAGGATGTGCTCGATATCGCACGTGAACATc                                                                                            >  1:113896/1‑118 (MQ=255)
                                                           cGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGAAAAAATTAACTCAGGATGTGCTCGATATCGCACGTGAACATCACTGGGTGACGTATGCTGAAGTg                                                                     <  1:999195/139‑1 (MQ=255)
                                                             aTGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGAAAAAATTAACTCAGGATGTGCTCGATATCGCACGTGAACATCACTGGGTGACGTATGCTGAAGTGga                                                                   <  2:597864/139‑1 (MQ=255)
                                                              tGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGAAAAAATTAACTCAGGATGTGCTCGATATCGCACGTGAACATCACTGGGTGACGTATGCTGAAGTGgag                                                                  >  1:681261/1‑139 (MQ=255)
                                                               gcgcTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGAAAAAATTAACTCAGGATGTGCTCGATATCGCACGTGAACATCACTGGGTGACGTATGCTGAAGTGgaga                                                                 <  2:983352/139‑1 (MQ=255)
                                                                cgcTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGAAAAAATTAACTCAGGATGTGCTCGATATCGCACGTGAACATCACTgg                                                                                       >  2:888328/1‑116 (MQ=255)
                                                                cgcTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGAAAAAATTAACTCAGGATGTGCTCGATATCGCACGTGAACATCACTgg                                                                                       <  1:888328/116‑1 (MQ=255)
                                                                     aaTTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGAAAAAATTAACTCAGGATGTGCTCGATATCGCACGTGaa                                                                                                <  1:238002/102‑1 (MQ=255)
                                                                     aaTTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGAAAAAATTAACTCAGGATGTGCTCGATATCGCACGTGaa                                                                                                >  2:238002/1‑102 (MQ=255)
                                                                      aTTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGAAAAAATTAACTCAGGATGTGCTCGATATCGCACGTGAACATCACTGGGTGACGTATGCTGAAGTGGAGATCGATaa                                                          >  2:31313/1‑139 (MQ=255)
                                                                           cGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGAAAAAATTAACTCAgg                                                                                                                       <  2:961613/73‑1 (MQ=255)
                                                                           cGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGAAAAAATTAACTCAgg                                                                                                                       >  1:961613/1‑73 (MQ=255)
                                                                               ccGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGAAAAAATTAACTCAGGATGTGCTCGATATCGCACGTGAACATCACTGGGTGACGTATGCTGAAGTGGAGATCGATAAACTGCAcgc                                                 >  2:41024/1‑139 (MQ=255)
                                                                               ccGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGAAAAAATTAACTCAGGATGTGCTCGATATCGCACGTGAACATCACTGGGTGACGTATGCTGAAGTGGAGATCGATAAACTGCAcgc                                                 >  2:290190/1‑139 (MQ=255)
                                                                                cgtagccgataaCATAATTCAGCATGTAGAGAATAACCGTGTCTCTTTGCTGAAAAGATTAACTCATGGTGTTCTCGATAACAGACGTGAAGGTCACTGGGGGACGTGTGCTGAAGTGGAGATCGATAAACTGGAcgcg                                                <  1:77580/129‑1 (MQ=255)
                                                                                cGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGAAAAAATTAACTCAGGATGTGCTCGATATCGCACGTGAACATCACTGGGTGACGTATGCTGAAGTGGAGATCGATAAACTGCAcgcg                                                <  1:424945/139‑1 (MQ=255)
                                                                                        aaaaCATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGAAAAAATTAACTCAGGATGTGCTCGATATCGCACGTGAACATCACTGGGTGACGTATGCTGAAGTGGAGATCGATAAACTGCACGCGCTGCGCTa                                        <  1:791361/139‑1 (MQ=255)
                                                                                          aaCATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGAAAAAATTAACTCAGGATGTGCTCGATATCGCACGTGAACATCACTGGGTGACGTATGCTGAAGTGGAGATCGATAAACTGCACGCGCTGCGCTACg                                      <  2:966936/139‑1 (MQ=255)
                                                                                          aaCATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGAAAAAATTAACTCAGGATGTGCTCGATATCGCACGTGAACATCACTGGGTGACGTATGCTGAAGTGGAGATCGATAAACTGCACGCGCTGCGCTACg                                      <  1:512725/139‑1 (MQ=255)
                                                                                               tattCAGCATGTAGAGAATAACCGTTTCTCTATGCTGAAAAAATTAACTCAGGATGTGCTCGATATCGCACGTGAACATCACTGGGTGACGTATGCTGAAGTGGAGATCGATAAACTGCACGCGCTGCGCTACGCCGAt                                 <  2:829683/139‑1 (MQ=255)
                                                                                                 ttCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGAAAAAATTAACTCAGGATGTGCTCGATATCGCACGTGAACATCACTGGGTGACGTATGCTGaa                                                                        <  1:823335/98‑1 (MQ=255)
                                                                                                 ttCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGAAAAAATTAACTCAGGATGTGCTCGATATCGCACGTGAACATCACTGGGTGACGTATGCTGaa                                                                        >  2:823335/1‑98 (MQ=255)
                                                                                                  tCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGAAAAAATTAACTCAGGATGTGCTCGATATCGCACGTGAACATCACTGGGTGACGTATGCTGAAGTGGAGATCGATAAACTGCACGCGCTGCGCTACGCCGATTCg                              <  1:275986/139‑1 (MQ=255)
                                                                                                      cATGTAGAGAATAACCGTTTCTCTTTGCTGAAAAAATTAACTCAGGATGTGCTCGATATCGCACGTGAAc                                                                                               >  1:261395/1‑70 (MQ=255)
                                                                                                      cATGTAGAGAATAACCGTTTCTCTTTGCTGAAAAAATTAACTCAGGATGTGCTCGATATCGCACGTGAAc                                                                                               <  2:261395/70‑1 (MQ=255)
                                                                                                         gTAGAGAATAACCGTTTCTCTTTGCTGAAAAAATTAACTCAGGATGTGCTCGATATCGCACGTGAACATCACTGGGTGACGTATGCTGAAGTGGAGATCGATAAACTGCACGCGCTGCGCTACGCCGATTCGGTATCCa                       <  1:250138/139‑1 (MQ=255)
                                                                                                           agagAATAACCGTTTCTCTTTGCTGAAAAAATTAACTCAGGATGTGCTCGATATCGCACGTGAACATCACTGGGTGACGTATGCTGAAGTGGAGATCGATAAACTGCACGCGCTGCGCTACGCCGATTCGGTATCCATg                     >  1:395277/1‑139 (MQ=255)
                                                                                                               aaTAACCGTTTCTCTTTGCTGAAAAAATTAACTCAGGATGTGCTCGATATCGCACGTGAACATCACTGGGTGACGTATGCTGAAGTGGAGATCGATAAACTGCACGCGCTGCGCTACGCCGATTCGGTATCCATGACCt                 >  1:126990/1‑139 (MQ=255)
                                                                                                                aTAACCGTTTCCCTTTGCTGAAAAATTTAACTCAGGATGTGCTCGATATCGCCCGTGAACATCACTGGGTGACGTATGCTGAAGTGGAGATCGATAAACTGCACGTGCTGCGCTACGCCGAATCGGTATCCATGACCtt                >  1:146048/1‑139 (MQ=255)
                                                                                                                    ccGTTTCTCTTTGCTGAAAAAATTAACTCAGGATGTGCTCGATATCGCACGTGAACATCAc                                                                                          >  2:175736/1‑61 (MQ=255)
                                                                                                                    ccGTTTCTCTTTGCTGAAAAAATTAACTCAGGATGTGCTCGATATCGCACGTGAACATCAc                                                                                          <  1:175736/61‑1 (MQ=255)
                                                                                                                    ccGTTTCTCTTTGCTGAAAAAATTAACTCAGGATGTGCTCGATATCGCACGTGAACATCATTGGGTGACGTATGCTGAAGTGGAGATCGATAAACTGCACGCGCTGCGCTACGCCGATTCGGTATCCATGAc                   >  2:300985/1‑132 (MQ=255)
                                                                                                                    ccGTTTCTCTTTGCTGAAAAAATTAACTCAGGATGTGCTCGATATCGCACGTGAACATCACTGGGTGACGTa                                                                               >  2:784561/1‑72 (MQ=255)
                                                                                                                    ccGTTTCTCTTTGCTGAAAAAATTAACTCAGGATGTGCTCGATATCGCACGTGAACATCACTGGGTGACGTATGCTGAAGTGGAGATCGATAAACTGCAcgcg                                                >  2:900647/1‑103 (MQ=255)
                                                                                                                    ccGTTTCTCTTTGCTGAAAAAATTAACTCAGGATGTGCTCGATATCGCACGTGAACATCACTGGGTGACGTATGCTGAAGTGGAGATCGATAAACTGCAcgcg                                                <  1:900647/103‑1 (MQ=255)
                                                                                                                    ccGTTTCTCTTTGCTGAAAAAATTAACTCAGGATGTGCTCGATATCGCACGTGAACATCACTGGGTGACGTATGCTGAAGTGGAGATCGATAAACTGCACGCGCTGCGCTACGCCGATTCGGTATCCATGAc                   <  1:300985/132‑1 (MQ=255)
                                                                                                                    ccGTTTCTCTTTGCTGAAAAAATTAACTCAGGATGTGCTCGATATCGCACGTGAACAACACTGGGTGACGTATGCTGAAGTGGAGATCGATAAACTGCACGCGCTGCGCTACGCCGATTCGGTATCCATGACCTTAAGc            >  2:712816/1‑139 (MQ=255)
                                                                                                                    ccGTTTCTCTTTGCTGAAAAAATTAACACAGGATGTGCTCGATATCGCACGTGAACATCACTGGGTGACGTa                                                                               <  1:784561/72‑1 (MQ=255)
                                                                                                                       tttCTCTTTGCTGAAAAAATTAACTCAGGATGTGCTCGATATCGCACGTGAACATCACTGGGTGACGTATGCTGAAGTGGAGATCGATAAACTGCACGCGCTGCGCTACGCCGATTCGGTATCCATGACCTTAAGCTgg         >  1:78370/1‑139 (MQ=255)
                                                                                                                         tctcTTTGCTGAAAAAATTAACTCAGGATGTGCTCGATATCGCACGTGAACATCACTGGGTGACGTATGCTGAAGTGga                                                                   >  1:222280/1‑79 (MQ=255)
                                                                                                                         tctcTTTGCTGAAAAAATTAACTCAGGATGTGCTCGATATCGCACGTGAACATCACTGGGTGACGTATGCTGAAGTGga                                                                   <  2:222280/79‑1 (MQ=255)
                                                                                                                         tctcTTTGCTGAAAAAATTAACTCAGGATGTGCTCGATATCGCACGTGAACATCACTGGGTGACGTATGCTGAAGTGGAGATCGATAAACTGCACGCGCTGCGCTACGCCGATTCGGTATCCa                       >  2:779611/1‑123 (MQ=255)
                                                                                                                         tctcTTTGCTGAAAAAATTAACTCAGGATGTGCTCGATATCGCACGTGAACATCACTGGGTGACGTATGCTGAAGTGGAGATCGATAAACTGCACGCGCTGCGCTACGCCGATTCGGTATCCa                       <  1:779611/123‑1 (MQ=255)
                                                                                                                            ctTTGCTGAAAAAATTAACTCAGGATGTGCTCGATATCGCACGTGAACATCACTGGGTGACGTATGCTGAAGTGGAGATCGATAAACTGCACGCGCTGCGCTACGCCGATTCgg                             >  1:992990/1‑114 (MQ=255)
                                                                                                                            ctTTGCTGAAAAAATTAACTCAGGATGTGCTCGATATCGCACGTGAACATCACTGGGTGACGTATGCTGAAGTGGAGATCGATAAACTGCACGCGCTGCGCTACGCCGATTCgg                             <  2:992990/114‑1 (MQ=255)
                                                                                                                              ttGCTGAAAAAATTAACTCAGGATGTGCTCGATATCGCACGTGAACATCACTGGGTGACGTATGCTGAAGTGGAGATCGATAAACTGCACGCGCTGCGCTACGCCGATTCGGTATCCATGACCTTAAGCTGGCAGCGtt  >  2:75579/1‑139 (MQ=255)
                                                                                                                               tGCTGAAAAAATTAACTCAGGATGTGCTCGATATCGCAc                                                                                                     <  1:117816/39‑1 (MQ=255)
                                                                                                                               tGCTGAAAAAATTAACTCAGGATGTGCTCGATATCGCAc                                                                                                     >  2:117816/1‑39 (MQ=255)
                                                                                                                                    |                                                                                                                                    
GTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGCACGTGAACATCACTGGGTGACGTATGCTGAAGTGGAGATCGATAAACTGCACGCGCTGCGCTACGCCGATTCGGTATCCATGACCTTAAGCTGGCAGCGTT  >  NC_000913/2421421‑2421685

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