Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,156,603 T→C intergenic (+85/‑20) yqhD → / → dkgA aldehyde reductase, NADPH‑dependent/2,5‑diketo‑D‑gluconate reductase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,156,6030TC100.0% 277.5 / NA 71intergenic (+85/‑20)yqhD/dkgAaldehyde reductase, NADPH‑dependent/2,5‑diketo‑D‑gluconate reductase A
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (32/39);  total (32/39)

TGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCAAAAAGCGTTAGAAGTGGGTTATCG  >  NC_000913/3156457‑3156750
                                                                                                                                                  |                                                                                                                                                   
tGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGCTCat                                                                                                                                                 >  2:470152/1‑151 (MQ=255)
 gggCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGCTCata                                                                                                                                                <  2:302951/151‑1 (MQ=255)
 gggCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGCTCat                                                                                                                                                 >  1:334952/1‑150 (MQ=255)
  ggCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGCTCatat                                                                                                                                               <  1:984466/151‑1 (MQ=255)
    cGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGCTCAtata                                                                                                                                              >  1:652900/1‑150 (MQ=255)
      aaaaTCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGCTCATATaaa                                                                                                                                            <  1:872801/150‑1 (MQ=255)
      aaaaTCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGCTCATATaa                                                                                                                                             >  2:305354/1‑149 (MQ=255)
        aaTCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGCTCATATAAagg                                                                                                                                          >  1:86183/1‑150 (MQ=255)
        aaTCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGCTCATATAAagg                                                                                                                                          >  2:321762/1‑150 (MQ=255)
            aTGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGCTCATATAAAGGAGGaa                                                                                                                                     <  1:452574/151‑1 (MQ=255)
               aCATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGCTCATATAAAGGAGGAACGt                                                                                                                                  <  2:133599/151‑1 (MQ=255)
                cATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGCTCATATAAAGGAGGAACGTa                                                                                                                                 <  2:279696/151‑1 (MQ=255)
                 aTTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGCTCATATAAAGGAGGAACGTAt                                                                                                                                <  1:583682/151‑1 (MQ=255)
                   tACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGCTCATATAAAGGAGGAACGTATg                                                                                                                               <  1:308572/150‑1 (MQ=255)
                    aCGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGCTCATATAAAGGAGGAACGTATGGc                                                                                                                             <  2:43757/151‑1 (MQ=255)
                     cattgtgtggcAGCCACCGTAGATACGAAGCCACCCGCTAAGCTATTTGCGCCTCAAACATGCGTGTACGGGCATTTCGTCCAGACTTAAGTACACATCACCACACCGGAGCCTGCTCCGGTGATCTCATATAAAGGAGGAACGTATgg                                                                                                                              <  1:198199/139‑1 (MQ=255)
                     cGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGCTCATATAAAGGAGGAACGTATGGCt                                                                                                                            >  2:665529/1‑151 (MQ=255)
                     cGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGCTCATATAAAGGAGGAACGTATGGCt                                                                                                                            <  2:771282/151‑1 (MQ=255)
                           aTGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGCTCATATAAAGGAGGAACGTATGGCTAATCCa                                                                                                                      <  2:187975/151‑1 (MQ=255)
                           aTGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGCTCATATAAAGGAGGAACGTATGGCTAATCCa                                                                                                                      <  2:267967/151‑1 (MQ=255)
                             gTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGCTCATATAAAGGAGGAACGTATGGCTAATCCaa                                                                                                                     <  2:331727/150‑1 (MQ=255)
                              tCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGCTCATATAAAGGAGGAACGTATGGCTAATCCAAcc                                                                                                                   >  1:838395/1‑151 (MQ=255)
                                 gccgccGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGCTCATATAAAGGAGGAACGTATGGCTAATCCAACCGtt                                                                                                                <  2:555831/151‑1 (MQ=255)
                                    gccgTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGCTCATATAAAGGAGGAACGTATGGCTAATCCAACCGttatt                                                                                                             >  2:859653/1‑151 (MQ=255)
                                      cgTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGCTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTaa                                                                                                           <  2:775608/151‑1 (MQ=255)
                                       gTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGCTCATATAAAGGAGGAACGTATGGCTAATCCAACCGttatta                                                                                                            >  2:694538/1‑149 (MQ=255)
                                       gTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGCTCATATAAAGGAGGAACGTATGGCTAATCCAACCGttatta                                                                                                            >  1:661946/1‑149 (MQ=255)
                                          tataCGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGCTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCt                                                                                                        >  2:992291/1‑150 (MQ=255)
                                               gAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGCTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAgg                                                                                                   >  2:505149/1‑150 (MQ=255)
                                                aaGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGCTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGAt                                                                                                 >  1:811180/1‑151 (MQ=255)
                                                    cgcccgcTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGCTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGc                                                                                              >  1:217786/1‑150 (MQ=255)
                                                    cgcccgcTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGCTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCa                                                                                             <  2:94943/151‑1 (MQ=255)
                                                    cgcccgcTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGCTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCa                                                                                             <  1:378416/151‑1 (MQ=255)
                                                            aaGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGCTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATg                                                                                     >  1:934534/1‑151 (MQ=255)
                                                            aaGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGCTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATg                                                                                     <  1:567661/151‑1 (MQ=255)
                                                                      cGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGCTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTgg                                                                           <  1:291872/151‑1 (MQ=255)
                                                                       gCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGCTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTggg                                                                          >  2:259582/1‑151 (MQ=255)
                                                                              aCTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGCTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGc                                                                   >  1:991081/1‑151 (MQ=255)
                                                                                 ttCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGCTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCg                                                                  >  2:146931/1‑149 (MQ=255)
                                                                                   cGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGCTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTg                                                              <  1:211461/151‑1 (MQ=255)
                                                                                   cGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGCTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTg                                                              <  2:204295/151‑1 (MQ=255)
                                                                                     ttttCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGCTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGgc                                                            >  1:231423/1‑151 (MQ=255)
                                                                                      tttCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGCTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGgc                                                            >  2:281297/1‑150 (MQ=255)
                                                                                        tCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGCTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGgcaa                                                          <  2:612580/150‑1 (MQ=255)
                                                                                         cGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGCTCATATAAAGGAGGAACGTATGGCTAGTCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGgcaag                                                         <  1:956362/150‑1 (MQ=255)
                                                                                             cATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGCTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGt                                                    <  2:312761/151‑1 (MQ=255)
                                                                                             cATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGCTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGt                                                    >  1:597973/1‑151 (MQ=255)
                                                                                              aTTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGCTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTa                                                   <  1:936424/151‑1 (MQ=255)
                                                                                                ttCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGCTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTaa                                                  <  2:1024358/150‑1 (MQ=255)
                                                                                                 tCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGCTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAAt                                                 <  1:828292/150‑1 (MQ=255)
                                                                                                  cGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGCTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATg                                                <  1:638838/150‑1 (MQ=255)
                                                                                                   gTCCAGACGTAAGTACAAAACACCTCACCGGAGCCTGCTCCTGTGAGCTCATATAAAGGACGAACGTATGGATAATCCAACCGTTCTTAAGCTACAGGATGGGAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAAtta                                               <  1:604939/150‑3 (MQ=255)
                                                                                                     ccAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGCTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAgg                                             >  1:53649/1‑150 (MQ=255)
                                                                                                            tAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGCTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCTTCTGGCAAGCAAGCAATGAGGAAGTAAt                                      <  1:375392/150‑1 (MQ=255)
                                                                                                             aaGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGCTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATc                                     >  1:241305/1‑150 (MQ=255)
                                                                                                             aaGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGCTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCa                                    <  2:657566/151‑1 (MQ=255)
                                                                                                              aGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGCTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCAc                                   >  2:398929/1‑151 (MQ=255)
                                                                                                                ttCACAACACCTCACCGGAGCCTGCTCCGGTGAGCTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCg                                 >  1:40330/1‑151 (MQ=255)
                                                                                                                ttCACAACACCTCACCGGAGCCTGCTCCGGTGAGCTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCg                                 >  2:302622/1‑151 (MQ=255)
                                                                                                                 tCACAACACCTCACCGGAGCCTGCTCCGGTGAGCTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCg                                 <  1:323901/150‑1 (MQ=255)
                                                                                                                              ccGGAGCCTGCTCCGGTGAGCTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCAAAAAGc                    <  2:505567/150‑1 (MQ=255)
                                                                                                                              ccGGAGCCTGCTCCGGTGAGCTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCAAAAAGc                    <  1:270177/150‑1 (MQ=255)
                                                                                                                                 gAGCCTGCTCCGGTGAGCTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCAAAAAGCGTTa                <  2:335027/151‑1 (MQ=255)
                                                                                                                                       gCTCCGGTGAGCTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCAAAAAGCGTTAGAAGTg          >  1:485817/1‑151 (MQ=255)
                                                                                                                                        cTCCGGTGAGCTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCAAAAAGCGTTAGAAGTgg         >  2:201255/1‑151 (MQ=255)
                                                                                                                                         tCCGGTGAGCTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCAAAAAGCGTTAGAAGTggg        >  1:244557/1‑151 (MQ=255)
                                                                                                                                          ccGTTGAGCTCATATAATGGAGGACCGTATGGCTAATCCACCCGTTATTAAGCTACAGGATGGCAATTTCTTGCCCAAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAATTAATCACCGCCATTCAAAAAGCGTTAGAAGTGGGt       <  2:939937/151‑1 (MQ=255)
                                                                                                                                          ccGGTGAGCTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCAAAAAGCGTTAGAAGTgg         <  1:953298/149‑1 (MQ=255)
                                                                                                                                            ggTGAGCTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCAAAAAGCGTTAGAAGTGGGTTa     <  2:86190/151‑1 (MQ=255)
                                                                                                                                             gTGAGCTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCAAAAAGCGTTAGAAGTGGGTTa     >  1:572470/1‑150 (MQ=255)
                                                                                                                                               gAGCTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCAAAAAGCGTTAGAAGTGGGTTATCg  >  2:182050/1‑151 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
TGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCAAAAAGCGTTAGAAGTGGGTTATCG  >  NC_000913/3156457‑3156750

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