Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 518,863 C→T 10.4% G96E (GGA→GAA)  ybbO ← NADP(+)‑dependent aldehyde reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913518,8630CT10.4% 217.6 / 11.0 76G96E (GGA→GAA) ybbONADP(+)‑dependent aldehyde reductase
Reads supporting (aligned to +/- strand):  ref base C (31/37);  new base T (2/6);  total (33/43)
Fisher's exact test for biased strand distribution p-value = 4.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATGTCATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAGATCGATCAACACGCCGGTAAATCCCATGCTGTTCATGCGCTCAACATCATCCGGTTTCCGG  >  NC_000913/518716‑519012
                                                                                                                                                   |                                                                                                                                                     
gatgTCATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGGATTCg                                                                                                                                                     <  1:756760/150‑1 (MQ=255)
  tgTCATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGc                                                                                                                                                   >  2:122810/1‑150 (MQ=255)
  tgTCATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGc                                                                                                                                                   <  2:547571/150‑1 (MQ=255)
    tcatcaCAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCAt                                                                                                                                                 <  1:319531/150‑1 (MQ=255)
      atcaCAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGGATTCGGCAtgt                                                                                                                                               <  2:1520365/150‑3 (MQ=255)
      atcaCAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTg                                                                                                                                               <  1:1025401/150‑1 (MQ=255)
       tcaCAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGt                                                                                                                                              >  2:348777/1‑150 (MQ=255)
       tcaCAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGt                                                                                                                                              >  2:557921/1‑150 (MQ=255)
        cacaATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGtt                                                                                                                                             >  2:1246008/1‑150 (MQ=255)
        cacaATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGtt                                                                                                                                             <  2:500144/150‑1 (MQ=255)
           aaTACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGGATTCGGCAtgtatggc                                                                                                                                          <  1:1519932/150‑8 (MQ=255)
            aTACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAg                                                                                                                                         <  1:939208/150‑1 (MQ=255)
            aTACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAg                                                                                                                                         >  2:55334/1‑150 (MQ=255)
               cGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATc                                                                                                                                      <  1:1224662/150‑1 (MQ=255)
               cGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATc                                                                                                                                      <  1:571576/150‑1 (MQ=255)
                 cccTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGTTGGTGTAAAGGTGGCCATACATGCCGATTTTGGCATTGTTTAAGATccc                                                                                                                                    <  2:1356148/150‑1 (MQ=255)
                 cccTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATccc                                                                                                                                    >  1:1460736/1‑150 (MQ=255)
                   cTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCAt                                                                                                                                  <  1:415812/150‑1 (MQ=255)
                          gTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACa                                                                                                                           <  1:1314210/150‑1 (MQ=255)
                            gCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAAt                                                                                                                         >  2:194391/1‑150 (MQ=255)
                                tAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATc                                                                                                                     <  1:348777/150‑1 (MQ=255)
                                    aTCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTc                                                                                                                 >  1:414251/1‑150 (MQ=255)
                                      cgcgGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAg                                                                                                               >  2:1481475/1‑150 (MQ=255)
                                         gggTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGc                                                                                                            <  2:261465/150‑1 (MQ=255)
                                            tAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGAt                                                                                                         <  2:889598/150‑1 (MQ=255)
                                               cAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCAc                                                                                                      >  1:317297/1‑150 (MQ=255)
                                               cAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCAc                                                                                                      <  1:781263/150‑1 (MQ=255)
                                                 ggCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCt                                                                                                    <  1:55334/150‑1 (MQ=255)
                                                 ggCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCt                                                                                                    >  1:1299713/1‑150 (MQ=255)
                                                        ggTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGc                                                                                             >  2:561442/1‑150 (MQ=255)
                                                           gAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTgc                                                                                          >  1:186091/1‑150 (MQ=255)
                                                              cTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGcgcg                                                                                       <  1:1301016/150‑1 (MQ=255)
                                                              cTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGcgcg                                                                                       <  1:1176324/150‑1 (MQ=255)
                                                                ggTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGAt                                                                                     >  1:737274/1‑150 (MQ=255)
                                                                 gtgtGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATc                                                                                    <  2:533006/150‑1 (MQ=255)
                                                                 gtgtGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATc                                                                                    <  2:1055083/150‑1 (MQ=255)
                                                                  tgtgCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCa                                                                                   <  2:317297/150‑1 (MQ=255)
                                                                   gtgCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCaa                                                                                  <  2:407412/150‑1 (MQ=255)
                                                                      cgcCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATAGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAacac                                                                               >  2:849552/1‑150 (MQ=255)
                                                                           aaaaaGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCt                                                                          <  2:834975/150‑1 (MQ=255)
                                                                            aaaaGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTg                                                                         <  2:414251/150‑1 (MQ=255)
                                                                              aaGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGt                                                                       <  2:1299713/150‑1 (MQ=255)
                                                                              aaGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGt                                                                       <  2:88146/150‑1 (MQ=255)
                                                                                gTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCt                                                                                                    >  2:114765/1‑119 (MQ=255)
                                                                                gTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCt                                                                                                    <  1:114765/119‑1 (MQ=255)
                                                                                    gCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAac                                                                                 <  1:506441/134‑1 (MQ=255)
                                                                                    gCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAac                                                                                 >  2:506441/1‑134 (MQ=255)
                                                                                     cGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCa                                                                >  2:1029442/1‑150 (MQ=255)
                                                                                       gAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAga                                                              >  2:724958/1‑150 (MQ=255)
                                                                                         aaaCTGCTGTTCCATCTGCGCACGGCTGATTGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAgatc                                                            <  1:849552/150‑1 (MQ=255)
                                                                                              gctgTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTAATGGTTATCACCTCGTCTGCTGCGCGATCAACACTTTCTGGTGAATCCAGATCTATCa                                                       <  2:1267273/150‑1 (MQ=255)
                                                                                                tgTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAGATCGATCAac                                                     <  1:662203/150‑1 (MQ=255)
                                                                                                  ttCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAGATCGATCAacac                                                   >  1:1144753/1‑150 (MQ=255)
                                                                                                         tGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGAGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAGATCGATCAACACGCCGGTa                                            >  1:1515085/1‑150 (MQ=255)
                                                                                                           cgcACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAGATCGATCAACACGCCGGTaaa                                          <  1:1246008/150‑1 (MQ=255)
                                                                                                              aCGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAGATCGATCAACACGCCGGTAAATcc                                       <  1:989702/150‑1 (MQ=255)
                                                                                                              aCGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAGATCGATCAACACGCCGGTAAATcc                                       <  1:194391/150‑1 (MQ=255)
                                                                                                              aCGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAGATCGATCAACACGCCGGTAAATcc                                       <  1:557921/150‑1 (MQ=255)
                                                                                                                  cTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAGATCGATCAACACGCCGGTAAATCCCATg                                   >  1:214107/1‑150 (MQ=255)
                                                                                                                  cTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAGATCGATCAACACGCCGGTAAATCCCATg                                   >  1:178045/1‑150 (MQ=255)
                                                                                                                        gtggAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAgatcgatc                                                        <  1:1310737/123‑1 (MQ=255)
                                                                                                                        gtggAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAgatcgatc                                                        >  2:1310737/1‑123 (MQ=255)
                                                                                                                        gtggAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAGATCGATCAACACGCCGGTAAATCCCATGCTGTTc                             >  2:587544/1‑150 (MQ=255)
                                                                                                                         tggAAAGGGGGCCATACATGGCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGTTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAGATCGATCAACACGCCGGTAAATCCCATGCTGTTCa                            <  1:1175788/150‑1 (MQ=255)
                                                                                                                          ggAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTACTGGTGAATCAAGATCGATCAACACGCCGGTAAATCCCATGCTGTTCAt                           >  1:527563/1‑150 (MQ=255)
                                                                                                                           gAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAGATCGATCAACACGCCGGTAAATCCCATGCTGTTCATg                          <  1:122810/150‑1 (MQ=255)
                                                                                                                            aaaGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAGATCGATCAACACGCCGGTAAATCCCATGCTGTTCATgc                         <  1:76921/150‑1 (MQ=255)
                                                                                                                               gggggCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAGATCGATCAACACGCCGGTAAATCCCATGCTGTTCATGCGCt                      >  2:453410/1‑150 (MQ=255)
                                                                                                                                ggggCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAGATCGATCAACACGCCGGTAAATCCCATGCTGTTCATGCGCTc                     >  2:1239526/1‑150 (MQ=255)
                                                                                                                                     taacaaTGCCGGATTCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAGATCGATCAACACGCCGGTAAATCCCATGCTGtt                              <  2:1260592/131‑1 (MQ=255)
                                                                                                                                     taacaaTGCCGGATTCGGCAGTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAGATCGATCAACACGCCGGTAAATCCCATGCTGtt                              >  1:1260592/6‑136 (MQ=255)
                                                                                                                                     cATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAGATCGATCAACACGCCGGTAAATCCCATGCTGTTCATGCGCTCAAcat                >  1:84012/1‑150 (MQ=255)
                                                                                                                                      aacaaTGCCGGATTCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAGATCGATCAACACGCCGGTAAATCCCATGCTGTTCATGCGCTCAAcatc               >  1:1240662/5‑150 (MQ=255)
                                                                                                                                        aCATGCCGAATCCGGCAGTGTTAAAGCTGCCATACAGACAATTATCGGTCAGGGCGATCGCCTCGTCGGCTGCGCGATCACCAGTTTCTGGTGAATCCAGATCGCTGAACGGGCCGGTATATCCCATGCTGTTCGTGCGCTGAAcatcat             >  2:587164/1‑150 (MQ=255)
                                                                                                                                          aTGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCACCACTTTCTGGTGAATCCAGATCGATCAACACGCCGGTAAATCCCATGCTGTTCATGCGCTCAACATCATcc           >  2:1029339/1‑150 (MQ=255)
                                                                                                                                            gCCGGATTCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAGATCGATCAACACGCCGGTAAATCCCATGCTGTTCATGCGCTCAACATCATCCgg         <  1:1469502/150‑1 (MQ=255)
                                                                                                                                                   ccGGCATTGTTAAAGATCCCATACAGACAATTATCGTTCAGGGCTATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAGATCGATCAACACGCCGGTAAATCCCATGCTGTTCATGCGCTCAACATCATCCGGTTTCCgg  <  2:349768/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GATGTCATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAGATCGATCAACACGCCGGTAAATCCCATGCTGTTCATGCGCTCAACATCATCCGGTTTCCGG  >  NC_000913/518716‑519012

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