Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 518,863 C→T 6.6% 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,8630CT 6.6% 231.5 / 4.6 75G96E (GGA→GAA) ybbONADP(+)‑dependent aldehyde reductase
Reads supporting (aligned to +/- strand):  ref base C (33/37);  new base T (2/3);  total (35/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.34e-01

TGTCATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAGATCGATCAACACGCCGGTAAATCCCATGCTGTTCATGCGCTCAACATCATCCG  >  NC_000913/518718‑519004
                                                                                                                                                 |                                                                                                                                             
tgTCATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGc                                                                                                                                           >  1:1902411/1‑150 (MQ=255)
  tcatcaCAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCAt                                                                                                                                         >  2:782541/1‑150 (MQ=255)
   catcaCAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGGATTCGGCAtg                                                                                                                                        <  2:1923325/150‑2 (MQ=255)
   catcaAAATACGCTCTTCAGCTTTGGGTAACATCGGGGTTAACAGGCGTATTGTGAGCTGGTGTTCTCGGAAAAATTTGTCGTAAAACTGCTTTTCGATCTGCTAACGGCTGATGGTTGAAAGGGGGCCATACATGCCGAATCGGGCAtt                                                                                                                                        <  2:1623143/150‑1 (MQ=255)
    atcaCAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTg                                                                                                                                       <  1:1262012/150‑1 (MQ=255)
     tcaCAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGt                                                                                                                                      >  2:3729/1‑150 (MQ=255)
      cacaATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGtt                                                                                                                                     <  2:1858890/150‑1 (MQ=255)
      cacaATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGt                                                                                                                                      <  1:1259511/149‑1 (MQ=255)
      cacaATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGt                                                                                                                                      >  2:1259511/1‑149 (MQ=255)
       acaATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGGATTCGGCAtgtatg                                                                                                                                    <  1:700240/150‑6 (MQ=255)
       acaATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGGATTCGGCAtgtatg                                                                                                                                    <  1:674862/150‑6 (MQ=255)
        caATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTaa                                                                                                                                   <  1:633038/150‑1 (MQ=255)
          aTACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAg                                                                                                                                 <  2:449615/150‑1 (MQ=255)
          aTACGCCCTTCACCGAGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTACATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAg                                                                                                                                 >  2:957367/1‑150 (MQ=255)
            aCGCCCTTCACCGTGCGGTAACATCGCGGGAAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGAt                                                                                                                               <  1:1744702/150‑1 (MQ=255)
             cGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATc                                                                                                                              <  1:1566399/150‑1 (MQ=255)
               cccTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATccc                                                                                                                            <  2:903102/150‑1 (MQ=255)
                ccTTCACCGTGCGGTAACATCGGGGGTAACAGGCGCATGGTGAGGTGGTGTGCGCCGAAAAAGTTGGCGGAAAAGTTCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCa                                                                                                                           <  2:841075/150‑1 (MQ=255)
                ccTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCa                                                                                                                           >  2:224470/1‑150 (MQ=255)
                ccTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCa                                                                                                                           <  1:1752277/150‑1 (MQ=255)
                 cTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCAt                                                                                                                          <  1:1041009/150‑1 (MQ=255)
                  ttCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATa                                                                                                                         <  2:1464929/150‑1 (MQ=255)
                       cgagaggTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTTGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGAc                                                                                                                    >  2:258515/6‑150 (MQ=255)
                        gTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACa                                                                                                                   >  2:1472665/1‑150 (MQ=255)
                         tGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACaa                                                                                                                  >  2:822843/1‑150 (MQ=255)
                            ggTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGGCATACATGCCGAATCCGGCATTGTTAAAGCTCCCATACAGACAATTa                                                                                                               >  2:842285/1‑150 (MQ=255)
                              tAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATc                                                                                                             <  1:1647835/150‑1 (MQ=255)
                                   tCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCa                                                                                                        >  2:1352577/1‑150 (MQ=255)
                                    cgcgGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAATAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAg                                                                                                       <  1:1209310/150‑1 (MQ=255)
                                           aaCAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGAta                                                                                                >  2:1241074/1‑149 (MQ=255)
                                           aaCAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATc                                                                                                <  1:741744/150‑1 (MQ=255)
                                           aaCAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATc                                                                                                <  1:123384/150‑1 (MQ=255)
                                               ggCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCt                                                                                            >  2:370636/1‑150 (MQ=255)
                                                    aTGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCtcgtcg                                                                                       <  2:531736/150‑1 (MQ=255)
                                                      ggTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGc                                                                                     <  2:41923/150‑1 (MQ=255)
                                                       gTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCt                                                                                    >  1:886254/1‑150 (MQ=255)
                                                        tGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTg                                                                                   <  1:771609/150‑1 (MQ=255)
                                                          aGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTgcg                                                                                 >  1:498300/1‑150 (MQ=255)
                                                            cTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGcgcg                                                                               <  2:515687/150‑1 (MQ=255)
                                                                    cgcCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGCCGATCACATCGTCGGCTGCGCGATCAacaa                                                                       >  2:993357/1‑149 (MQ=255)
                                                                      ccGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACtt                                                                     >  2:838253/1‑150 (MQ=255)
                                                                        gAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTc                                                                   <  1:224470/150‑1 (MQ=255)
                                                                           aaaGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGTTTATTTCTATGATATGGTCTGCTGCGAGATCAACACTTTCTgg                                                                <  2:694782/150‑1 (MQ=255)
                                                                             aGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGTCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCTCGACTTTCTGGTg                                                              >  1:439654/1‑150 (MQ=255)
                                                                                tGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAAt                                                           <  2:1902411/150‑1 (MQ=255)
                                                                                 ggCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATc                                                          >  2:1953359/1‑150 (MQ=255)
                                                                                   cGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCa                                                        >  2:1308212/1‑150 (MQ=255)
                                                                                       aaaCTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAgatc                                                    >  2:1926591/1‑150 (MQ=255)
                                                                                        aaCTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAgatcg                                                   <  1:370636/150‑1 (MQ=255)
                                                                                          ctgctgTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAgatcgat                                                 >  1:1906384/1‑150 (MQ=255)
                                                                                                ttCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAGATCGATCAacac                                           <  1:993357/150‑1 (MQ=255)
                                                                                                 tCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAGATCGATCAACACg                                          >  1:1627507/1‑150 (MQ=255)
                                                                                                  ccATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAGATCGATCAACACGc                                         >  2:242896/1‑150 (MQ=255)
                                                                                                   cATCTGCGCACGGCTGATGGTGGAATGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAGATCGATCAACACGcc                                        <  2:574938/150‑1 (MQ=255)
                                                                                                             cGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAGATCGATCAACACGCCGGTAAATccc                              >  2:393191/1‑150 (MQ=255)
                                                                                                              ggCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAGATCGATCAACACGCCGGTAAATCCCa                             <  1:724229/150‑1 (MQ=255)
                                                                                                               gCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAGATCGATCAACACGCCGGTAAATCCCAt                            >  1:952274/1‑150 (MQ=255)
                                                                                                               gCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAGATCGATCAACACGCCGGTAAATCCCAt                            >  1:511474/1‑150 (MQ=255)
                                                                                                                    tggtggAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAGATCGATCAACACGCCGGTAAATCCCATGCTGt                       >  2:617655/1‑150 (MQ=255)
                                                                                                                     ggtggAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAGATCGATCAACACGCCGGTAAATCCCATGCTGtt                      <  1:1352577/150‑1 (MQ=255)
                                                                                                                     ggtggAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAGATCGATCAACACGCCGGTAAATCCCATGCTGtt                      <  1:1926591/150‑1 (MQ=255)
                                                                                                                           aaGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAGATCGATCAACACGCCGGTAAATCCCATGCTGTTCATgcg                <  1:258515/150‑1 (MQ=255)
                                                                                                                               gggCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAGATCGATCAACACGCCGGTAAATCCCATGCTGTTCATGCGCTCa            <  2:1906384/150‑1 (MQ=255)
                                                                                                                               gggCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAGATCGATCAACACGCCGGTAAATCCCATGCTGTTCATGCGCTCa            <  1:1472665/150‑1 (MQ=255)
                                                                                                                                 gCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAGATCGATCAACACGCCGGTAAATCCCATGCTGTTCATGCGCTCAAc          <  2:511474/150‑1 (MQ=255)
                                                                                                                                 gCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAGATCGATCAACACGCCGGTAAATCCCATGCTGTTCATGCGCTCAAc          <  1:3729/150‑1 (MQ=255)
                                                                                                                                  ttaacaaTGCCGGATTCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAGATCGATCAACACGCCGGTAAATCCCATGCTGTTCATGCGCTCAAca         >  2:1989095/7‑150 (MQ=255)
                                                                                                                                  ccATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAGATCGATCAACACGCCGGTAAATCCCATGCTGTTCATGCGCTCAAca         >  1:1795741/1‑150 (MQ=255)
                                                                                                                                    aacaaTGCCGGATTCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAGATCGATCAACACGCCGGTAAATCCCATGCTGTTCATGCGCTCAAcata       >  1:562757/5‑149 (MQ=255)
                                                                                                                                     tACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAGATCGATCAACACGCCGGTAAATCCCATGCTGTTCATGCGCTCAAcatca      >  2:800204/1‑150 (MQ=255)
                                                                                                                                     tACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAGATCGATCAACACGCCGGTAAATCCCATGCTGTTCATGCGCTCAAcatca      <  2:1512896/150‑1 (MQ=255)
                                                                                                                                     tACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAGATCGATCAACACGCCGGTAAATCCCATGCTGTTCATGCGCTCAAcatca      >  1:1436941/1‑150 (MQ=255)
                                                                                                                                     tACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAGATCGATCAACACGCCGGTAAATCCCATGCTGTTCATGCGCTCAAcatca      >  1:1394568/1‑150 (MQ=255)
                                                                                                                                        aTGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAGATCGATCAACACGCCGGTAAATCCCATGCTGTTCATGCGCTCAACATCATcc   <  2:1394568/150‑1 (MQ=255)
                                                                                                                                         tGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAGATCGATCAACACGCCGGTAAATCCCATGCTGTTCATGCGCTCAACATCATccc  >  2:1778516/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                             
TGTCATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAGATCGATCAACACGCCGGTAAATCCCATGCTGTTCATGCGCTCAACATCATCCG  >  NC_000913/518718‑519004

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