Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,756,553 A→T K286* (AAA→TAA)  pykF → pyruvate kinase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,756,5530AT100.0% 256.8 / NA 63K286* (AAA→TAA) pykFpyruvate kinase I
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (35/28);  total (35/28)

TCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGA  >  NC_000913/1756406‑1756695
                                                                                                                                                   |                                                                                                                                              
tctGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCtaaa                                                                                                                                             >  2:615013/1‑151 (MQ=255)
  tGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAAc                                                                                                                                           <  2:253623/151‑1 (MQ=255)
  tGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAAc                                                                                                                                           <  2:38269/151‑1 (MQ=255)
  tGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAAc                                                                                                                                           <  2:419431/151‑1 (MQ=255)
    aCGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAAcc                                                                                                                                          <  1:429380/150‑1 (MQ=255)
    aCGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAAcc                                                                                                                                          <  2:94444/150‑1 (MQ=255)
         atcatGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGcc                                                                                                                                    <  2:598813/151‑1 (MQ=255)
          tcatGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGcc                                                                                                                                    >  2:85183/1‑150 (MQ=255)
             tGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGa                                                                                                                                 >  1:205067/1‑150 (MQ=255)
             tGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGa                                                                                                                                 >  1:123054/1‑150 (MQ=255)
             tGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGa                                                                                                                                 >  2:182858/1‑150 (MQ=255)
                                gggTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGa                                                                                                              >  1:614207/1‑150 (MQ=255)
                                gggTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGa                                                                                                              >  1:522410/1‑150 (MQ=255)
                                 ggTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACCCGCGCAGAAACCGGGGACg                                                                                                            >  2:276808/1‑151 (MQ=255)
                                         aaTCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAAc                                                                                                    <  1:85171/151‑1 (MQ=255)
                                                    aagaagTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGa                                                                                         >  1:637274/1‑151 (MQ=255)
                                                      gaagTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACCCCATCCTCGACg                                                                                       >  2:483099/1‑151 (MQ=255)
                                                       aagTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGAc                                                                                        <  1:212825/149‑1 (MQ=255)
                                                       aagTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACgg                                                                                      <  1:473019/151‑1 (MQ=255)
                                                          ttATCTTCGCCCAGAAGATGATGGACGAAAAGTGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCTTCGACGGTAc                                                                                   <  1:103807/151‑1 (MQ=255)
                                                              cTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGAc                                                                               <  1:598821/151‑1 (MQ=255)
                                                               ttCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACg                                                                              <  1:300225/151‑1 (MQ=255)
                                                               ttCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACg                                                                              <  2:165512/151‑1 (MQ=255)
                                                                tCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGc                                                                             <  2:348588/151‑1 (MQ=255)
                                                                tCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGc                                                                             <  1:325702/151‑1 (MQ=255)
                                                                tCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGc                                                                             <  1:69759/151‑1 (MQ=255)
                                                                   cccAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAg                                                                           <  1:314355/150‑1 (MQ=255)
                                                                   cccAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAg                                                                           <  2:29100/150‑1 (MQ=255)
                                                                        aagaTGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGa                                                                        <  1:4271/148‑1 (MQ=255)
                                                                        aagaTGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGc                                                                     <  1:515963/151‑1 (MQ=255)
                                                                          gatgatgatCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATg                                                                      >  1:328692/1‑148 (MQ=255)
                                                                                 atCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTg                                                            <  1:338817/151‑1 (MQ=255)
                                                                                   cGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGaa                                                          >  2:387616/1‑151 (MQ=255)
                                                                                       aaaTGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCg                                                      >  2:246775/1‑151 (MQ=255)
                                                                                         aTGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCa                                                    >  1:148568/1‑151 (MQ=255)
                                                                                         aTGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCa                                                    >  1:309502/1‑151 (MQ=255)
                                                                                         aTGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCa                                                    >  1:404737/1‑151 (MQ=255)
                                                                                         aTGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCa                                                    >  1:571840/1‑151 (MQ=255)
                                                                                                     aCGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGATGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATAc                                        <  1:260931/151‑1 (MQ=255)
                                                                                                     aCGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATAc                                        <  2:553640/151‑1 (MQ=255)
                                                                                                     aCGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATAc                                        <  2:440332/151‑1 (MQ=255)
                                                                                                     aCGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATAc                                        <  2:299213/151‑1 (MQ=255)
                                                                                                      cGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATAcc                                       <  2:204882/151‑1 (MQ=255)
                                                                                                      cGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATAc                                        >  1:170666/1‑150 (MQ=255)
                                                                                                      cGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATAc                                        >  1:72509/1‑150 (MQ=255)
                                                                                                      cGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATAc                                        >  1:305286/1‑150 (MQ=255)
                                                                                                      cGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGGGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATAc                                        >  2:17618/1‑150 (MQ=255)
                                                                                                       gTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATAcc                                       >  2:280230/1‑150 (MQ=255)
                                                                                                            gTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTgg                                 >  1:560376/1‑151 (MQ=255)
                                                                                                            gTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTgg                                 >  1:182291/1‑151 (MQ=255)
                                                                                                              cGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGaa                               >  1:621317/1‑151 (MQ=255)
                                                                                                              cGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGaa                               >  1:293560/1‑151 (MQ=255)
                                                                                                                ttATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAg                              >  1:313901/1‑150 (MQ=255)
                                                                                                                ttATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAg                              >  1:169322/1‑150 (MQ=255)
                                                                                                                ttATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAg                              >  1:24896/1‑150 (MQ=255)
                                                                                                                           aCCCAGATGCTGGGTTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATc                   <  1:572273/150‑1 (MQ=255)
                                                                                                                           aCCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATc                   <  2:382919/150‑1 (MQ=255)
                                                                                                                            cccAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCAt                 >  2:182504/1‑151 (MQ=255)
                                                                                                                               aGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATgg               >  2:388293/1‑150 (MQ=255)
                                                                                                                                 aTGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGa            >  2:649939/1‑151 (MQ=255)
                                                                                                                                      ggATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCAt        >  1:240253/1‑150 (MQ=255)
                                                                                                                                      ggATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCAt        >  2:338639/1‑150 (MQ=255)
                                                                                                                                           ccATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGa  >  1:59306/1‑151 (MQ=255)
                                                                                                                                           ccATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGa  >  2:322672/1‑151 (MQ=255)
                                                                                                                                                   |                                                                                                                                              
TCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGA  >  NC_000913/1756406‑1756695

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