Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,104,750 T→C H112R (CAC→CGC)  wbbJ ← putative lipopolysaccharide biosynthesis O‑acetyl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,104,7500TC100.0% 234.1 / NA 66H112R (CAC→CGC) wbbJputative lipopolysaccharide biosynthesis O‑acetyl transferase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (35/31);  total (35/31)

ACGACTCCATTACCAATAATTGTTCCAGGCAAAACCGTCACATTCTCACCCAACCAAACCCTCTGGCCAATTACAACAGCTGAAGATTCCAACGTGCGCATGTCTGGAGGTATATTTGGCGAACTCATTGGATCAGAGTGCTTAAAGGAACCGTGATTATGATCGGTAATAAATACTTTACTTGCAATAAGCGTATCCCGACCTATCGTAACGCTCTCAATTGAGGCGATATGAACATAGTCGTTAACTTGCACATTATCGGAAAAAAAAATCACGCCACGT  >  NC_000913/2104612‑2104893
                                                                                                                                          |                                                                                                                                               
aCGACTCCATTACCAATAATTGTTCCAGGCAAAACCGTCACATTCTCACCCAACCAAACCCTCTGGCCAATTACAACAGCTGAAGATTCCAACGTGCGCATGTCTGGAGGTATATTTGGCGAACTCATTGGATCAGAGCGCTTAAAGGAAc                                                                                                                                     <  1:618324/151‑1 (MQ=255)
   aCTCCATTACCAATAATTGTTCCAGGCAAAACCGTCACATTCTCACCCAACCAAACCCTCTGGCCAATTACAACAGCTGAAGATTCCAACGTGCGCATGTCTGGAGGTATATTTGGCGAACTCATTGGATCAGAGCGCTTAAAGGAACCGt                                                                                                                                  >  1:265277/1‑151 (MQ=255)
      ccATTACCAATAATTGTTCCAGGCAAAACCGTCACATTCTCACCCAACCAAACCCTCTGGCCAATTACAACAGCTGAAGATTCCAACGTGCGCATGTCTGGAGGTATATTTGGCGAACTCATTGGATCAGAGCGCTTAAAGGAACCGTGa                                                                                                                                <  2:265251/150‑1 (MQ=255)
      ccATTACCAATAATTGTTCCAGGCAAAACCGTCACATTCTCACCCAACCAAACCCTCTGGCCAATTACAACAGCTGAAGATTCCAACGTGCGCATGTCTGGAGGTATATTTGGCGAACTCATTGGATCAGAGCGCTTAAAGGAACCGTGAt                                                                                                                               <  1:653946/151‑1 (MQ=255)
        aTTACCAATAATTGTTCCAGGCAAAACCGTCACATTCTCACCCAACCAAACCCTCTGGCCAATTACAACAGCTGAAGATTCCAACGTGCGCATGTCTGGAGGTATATTTGGCGAACTCATTGGATCAGAGCGCTTAAAGGAACCGTGAtt                                                                                                                              >  1:376334/1‑150 (MQ=255)
             cAATAATTGTTCCAGGCAAAACCGTCACATTCTCACCCAACCAAACCCTCTGGCCAATTACAACAGCTGAAGATTCCAACGTGCGCATGTCTGGAGGTATATTTGGCGAACTCATTGGATCAGAGCGCTTAAAGGAACCGTGATTATGATc                                                                                                                        >  2:379728/1‑151 (MQ=255)
             cAATAATTGTTCCAGGCAAAACCGTCACATTCTCACCCAACCAAACCCTCTGGCCAATTACAACAGCTGAAGATTCCAACGTGCGCATGTCTGGAGGTATATTTGGCGAACTCATTGGATCAGAGCGCTTAAAGGAACCGTGATTATGATc                                                                                                                        >  1:519683/1‑151 (MQ=255)
              aataatTGTTCCAGGCAAAACCGTCACATTCTCACCCAACCAAACCCTCTGGCCAATTACAACAGCTGAAGATTCCAACGTGCGCATGTCTGGAGGTATATTTGGCGAACTCATTGGATCAGAGCGCTTAAAGGAACCGTGATTATGATCg                                                                                                                       >  2:900523/1‑151 (MQ=255)
              aataatTGTTCCAGGCAAAACCGTCACATTCTCACCCAACCAAACCCTCTGGCCAATTACAACAGCTGAAGATTCCAACGTGCGCATGTCTGGAGGTATATTTGGCGAACTCATTGGATCAGAGCGCTTAAAGGAACCGTGATTATGATCg                                                                                                                       >  2:17829/1‑151 (MQ=255)
              aataatTGTTCCAGGCAAAACCGTCACATTCTCACCCAACCAAACCCTCTGGCCAATTACAACAGCTGAAGATTCCAACGTGCGCATGTCTGGAGGTATATTTGGCGAACTCATTGGATCAGAGCGCTTAAAGGAACCGTGATTATGATCg                                                                                                                       >  1:786070/1‑151 (MQ=255)
              aataatTGTTCCAGGCAAAACCGTCACATTCTCACCCAACCAAACCCTCTGGCCAATTACAACAGCTGAAGATTCCAACGTGCGCATGTCTGGAGGTATATTTGGCGAACTCATTGGATCAGAGCGCTTAAAGGAACCGTGATTATGATCg                                                                                                                       >  1:750973/1‑151 (MQ=255)
                 aatTGTTCCAGGCAAAACCGTCACATTCTCACCCAACCAAACCCTCTGGCCAATTACAACAGCTGAAGATTCCAACGTGCGCATGTCTGGAGGTATATTTGGCGAACTCATTGGATCAGAGCGCTTAAAGGAACCGTGATTATGATCGGt                                                                                                                     >  2:615477/1‑150 (MQ=255)
                  atTGTTCCAGGCAAAACCGTCACATTCTCACCCAACCAAACCCTCTGGCCAATTACAACAGCTGAAGATTCCAACGTGCGCATGTCTGGAGGTATATTTGGCGAACTCATTGGATCAGAGCGCTTAAAGGAACCGTGATTATGATCGGta                                                                                                                    >  2:546309/1‑150 (MQ=255)
                     gTTCCAGGCAAAACCGTCACATTCTCACCCAACCAAACCCTCTGGCCAATTACAACAGCTGAAGATTCCAACGTGCGCATGTCTGGAGGTATATTTGGCGAACTCATTGGATCAGAGCGCTTAAAGGAACCGTGATTATGATCGGtaata                                                                                                                 >  1:627563/1‑150 (MQ=255)
                     gTTCCAGGCAAAACCGTCACATTCTCACCCAACCAAACCCTCTGGCCAATTACAACAGCTGAAGATTCCAACGTGCGCATGTCTGGAGGTATATTTGGCGAACTCATTGGATCAGAGCGCTTAAAGGAACCGTGATTATGATCGGtaata                                                                                                                 >  1:544685/1‑150 (MQ=255)
                           ggCAAAACCGTCACATTCTCACCCAACCAAACCCTCTGGCCAATTACAACAGCTGAAGATTCCAACGTGCGCATGTCTGGAGGTATATTTGGCGAACTCATTGGATCAGAGCGCTTAAAGGAACCGTGATTATGATCGGTAATAAAtactt                                                                                                          <  1:39167/151‑1 (MQ=255)
                             cAAAACCGTCACATTCTCACCCAACCAAACCCTCTGGCCAATTACAACAGCTGAAGATTCCAACGTGCGCATGTCTGGAGGTATATTTGGCGAACTCATTGGATCAGAGCGCTTAAAGGAACCGTGATTATGATCGGTAATAAAtacttta                                                                                                        <  2:396821/151‑1 (MQ=255)
                              aaaaCCGTCACATTCTCACCCAACCAAACCCTCTGGCCAATTACAACAGCTGAAGATTCCAACGTGCGCATGTCTGGAGGTATATTTGGCGAACTCATTGGATCAGAGCGCTTAAAGGAACCGTGATTATGATCGGTAATAAAtactttac                                                                                                       <  1:257092/151‑1 (MQ=255)
                               aaaCCGTCACATTCTCACCCAACCAAACCCTCTGGCCAATTACAACAGCTGAAGATTCCAACGTGCGCATGTCTGGAGGTATATTTGGCGAACTCATTGGATCAGAGCGCTTAAAGGAACCGTGATTATGATCGGTAATAAAtactttact                                                                                                      <  2:662347/151‑1 (MQ=255)
                                aaCCGTCACATTCTCACCCAACCAAACCCTCTGGCCAATTACAACAGCTGAAGATTCCAACGTGCGCATGTCTGGAGGTATATTTGGCGAACTCATTGGATCAGAGCGCTTAAAGGAACCGTGATTATGATCGGTAATAAAtactttactt                                                                                                     <  1:562565/151‑1 (MQ=255)
                                  ccGTCACATTCTCACCCAACCAAACCCTCTGGCCAATTACAACAGCTGAAGATTCCAACGTGCGCATGTCTGGAGGTATATTTGGCGAACTCATTGGATCAGAGCGCTTAAAGGAACCGTGATTATGATCGGTAATAAATACTTTACTTGc                                                                                                   <  1:227427/151‑1 (MQ=255)
                                   cGTCACATTCTCACCCAACCAAACCCTCTGGCCAATTACAACAGCTGAAGATTCCAACGTGCGCATGTCTGGAGGTATATTTGGCGAACTCATTGGATCAGAGCGCTTAAAGGAACCGTGATTATGATCGGTAATAAATACTTTACTTGCa                                                                                                  <  2:740680/151‑1 (MQ=255)
                                      cacaTTCTCACCCACCCAAACCCTCTGGCCAATTACAACAGCTGAAGATTCCAACGTGCGCATGTCTGGAGGTATATTTGGCGAACTCATTGGATCAGAGCGCTTAAAGGAACCGTGATTATGATCGGTAATAAATACTTTACTTGCAATa                                                                                               <  2:89614/151‑1 (MQ=255)
                                          ttCTCACCCAACCAAACCCTCTGGCCAATTACAACAGCTGAAGATTCCAACGTGCGCATGTCTGGAGGTATATTTGGCGAACTCATTGGATCAGAGCGCTTAAAGGAACCGTGATTATGATCGGTAATAAATACTTTACTTGCAATAAGCg                                                                                           >  1:244654/1‑151 (MQ=255)
                                               aCCCAACCAAACCCTCTGGCCAATTACAACAGCTGAAGATTCCAACGTGCGCATGTCTGGAGGTATATTTGGCGAACTCATTGGATCAGAGCGCTTAAAGGAACCGTGATTATGATCGGTAATAAATACTTTACTTGCAATAAGCGTATcc                                                                                      >  1:5427/1‑151 (MQ=255)
                                                cccAACCAAACCCTCTGGCCAATTACAACAGCTGAAGATTCCAACGTGCGCATGTCTGGAGGTATATTTGGCGAACTCATTGGATCAGAGCGCTTAAAGGAACCGTGATTATGATCGGTAATAAATACTTTACTTGCAATAAGCGTATccc                                                                                     >  1:382731/1‑151 (MQ=255)
                                                    accaaACCCTCTGGCCAATTACAACAGCTGAAGATTCCAACGTGCGCATGTCTGGAGGTATATTTGGCGAACTCATTGGATCAGAGCGCTTAAAGGAACCGTGATTATGATCGGTAATAAATACTTTACTTGCAATAAGCGTATCCCGAcc                                                                                 <  2:827742/151‑1 (MQ=255)
                                                     ccaaACCCTCTGGCCAATTACAACAGCTGAAGATTCCAACGTGCGCATGTCTGGAGGTATATTTGGCGAACTCATTGGATCAGAGCGCTTAAAGGAACCGTGATTATGATCGGTAATAAATACTTTACTTGCAATAAGCGTATCCCGACCt                                                                                >  1:20309/1‑151 (MQ=255)
                                                      caaACCCTCTGGCCAATTACAACAGCTGAAGATTCCAACGTGCGCATGTCTGGAGGTATATTTGGCGAACTCATTGGATCAGAGCGCTTAAAGGAACCGTGATTATGATCGGTAATAAATACTTTACTTGCAATAAGCGTATCCCGACCTa                                                                               <  2:552270/151‑1 (MQ=255)
                                                        aaCCCTCTGGCCAATTACAACAGCTGAAGATTCCAACGTGCGCATGTCTGGAGGTATATTTGGCGAACTCATTGGATCAGAGCGCTTAAAGGAACCGTGATTATGATCGGTAATAAATACTTTACTTGCAATAAGCGTATCCCGACCTATc                                                                             <  1:403560/151‑1 (MQ=255)
                                                        aaCCCTCTGGCCAATTACAACAGCTGAAGATTCCAACGTGCGCATGTCTGGAGGTATATTTGGCGAACTCATTGGATCAGAGCGCTTAAAGGAACCGTGATTATGATCGGTAATAAATACTTTACTTGCAATAAGCGTATCCCGACCTATc                                                                             >  2:111556/1‑151 (MQ=255)
                                                         aCCCTCTGGCCAATTACAACAGCTGAAGATTCCAACGTCCGCATGTCTGGAGGTATATTTGGCCAACTCATTGGATCAGAGCGCTTAAAGTAACCGTGATTATGATCGGTAATAAATACTTTACTTGCAAGAAGCGTATCCCGACCTACCg                                                                            <  2:167062/151‑1 (MQ=255)
                                                          cccTCTGGCCAATTACAACAGCTGAAGATTCCAACGTGCGCATGTCTGGAGGTATATTTGGCGAACTCATTGGATCAGAGCGCTTAAAGGAACCGTGATTATGATCGGTAATAAATACTTTACTTGCAATAAGCGTATCCCGACCTATCg                                                                            <  1:232992/150‑1 (MQ=255)
                                                          cccTCTGGCCAATTACAACAGCTGAAGATTCCAACGTGCGCATGTCTGGAGGTATATTTGGCGAACTCATTGGATCAGAGCGCTTAAAGGAACCGTGATTATGATCGGTAATAAATACTTTACTTGCAATAAGCGTATCCCGACCTATCGt                                                                           <  1:244009/151‑1 (MQ=255)
                                                           ccTCTGGCCAATTACAACAGCTGAAGATTCCAACGTGCGCATGTCTGGAGGTATATTTGGCGAACTCATTGGATCAGAGCGCTTAAAGGAACCGTGATTATGATCGGTAATAAATACTTTACTTGCAATAAGCGTATCCCGACCTATCGt                                                                           >  1:324983/1‑150 (MQ=255)
                                                            ctctGGCCAATTACAACAGCTGAAGATTCCAACGTGCGCATGTCTGGAGGTATATTTGGCGAACTCATTGGATCAGAGCGCTTAAAGGAACCGTGATTATGATCGGTAATAAATACTTTACTTGCAATAAGCGTATCCCGACCTATCGTa                                                                          >  1:499205/1‑150 (MQ=255)
                                                                ggCCAATTACAACAGCTGAAGATTCCAACGTGCGCATGTCTGGAGGTATATTTGGCGAACTCATTGGATCAGAGCGCTTAAAGGAACCGTGATTATGATCGGTAATAAATACTTTACTTGCAATAAGCGTATCCCGACCTATCGTAAc                                                                        <  1:442892/148‑1 (MQ=255)
                                                                  ccAATTACAACAGCTGAAGATTCCAACGTGCGCATGTCTGGAGGTATATTTGGCGAACTCATTGGATCAGAGCGCTTAAAGGAACCGTGATTATGATCGGTAATAAATACTTTACTTGCAATAAGCGTATCCCGACCTATCGTAACGctct                                                                   >  1:279372/1‑151 (MQ=255)
                                                                   cAATTACAACAGCTGAAGATTCCAACGTGCGCATGTCTGGAGGTATATTTGGCGAACTCATTGGATCAGAGCGCTTAAAGGAACCGTGATTATGATCGGTAATAAATACTTTACTTGCAATAAGCGTATCCCGACCTATCGTAACGCtctc                                                                  >  1:106912/1‑151 (MQ=255)
                                                                    aaTTACAACAGCTGAAGATTCCAACGTGCGCATGTCTGGAGGTATATTTGGCGAACTCATTGGATCAGAGCGCTTAAAGGAACCGTGATTATGATCGGTAATAAATACTTTACTTGCAATAAGCGTATCCCGACCTATCGTAACGCtctc                                                                  >  2:763756/1‑150 (MQ=255)
                                                                            caGCTGAAGATTCCAACGTGCGCATGTCTGGAGGTATATTTGGCGAACTCATTGGATCAGAGCGCTTAAAGGAACCGTGATTATGATCGGTAATAAATACTTTACTTGCAATAAGCGTATCCCGACCTATCGTAACGCTCTCAATTGAGGc                                                         >  2:166686/1‑151 (MQ=255)
                                                                                  aaGATTCCAACGTGCGCATGTCTGGAGGTATATTTGGCGAACTCATTGGATCAGAGCGCTTAAAGGAACCGTGATTATGATCGGTAATAAATACTTTACTTGCAATAAGCGTATCCCGACCTATCGTAACGCTCTCAATTGAGGCGATATg                                                   <  1:162549/151‑1 (MQ=255)
                                                                                  aaGATTCCAACGTGCGCATGTCTGAGGTATATTTGGCGAACTCATTGGATCAGAGCGCTTAAAGGAACCGTGATTATGATCGGTAATAAATACTTTACTTGCAATAAGCGTATCCCGACCTATCGTAACGCTCTCAATTGAGGCGATATGa                                                  <  1:321354/151‑1 (MQ=255)
                                                                                   aGATTCCAACGTGCGCATGTCTGGAGGTATATTTGGCGAACTCATTGGATCAGAGCGCTTAAAGGAACCGTGATTATGATCGGTAATAAATACTTTACTTGCAATAAGCGTATCCCGACCTATCGTAACGCTCTCAATTGAGGCGATATg                                                   >  2:797913/1‑150 (MQ=255)
                                                                                   aGATTCCAACGTGCGCATGTCTGGAGGTATATTTGGCGAACTCATTGGATCAGAGCGCTTAAAGGAACCGTGATTATGATCGGTAATAAATACTTTACTTGCAATAAGCGTATCCCGACCTATCGTAACGCTCTCAATTGAGGCGATATg                                                   >  1:339961/1‑150 (MQ=255)
                                                                                     aTTCCAACGTGCGCATGTCTGGAGGTATATTTGGCGAACTCATTGGATCAGAGCGCTTAAAGGAACCGTGATTATGATCGGTAATAAATACTTTACTTGCAATAAGCGTATCCCGACCTATCGTAACGCTCTCAATTGAGGCGATATGaa                                                 <  2:627514/150‑1 (MQ=255)
                                                                                       tCCAACGTGCGCATGTCTGGAGGTATATTTGGCGAACTCATTGGATCAGAGCGCTTAAAGGAACCGTGATTATGATCGGTAATAAATACTTTACTTGCAATAAGCGTATCCCGACCTATCGTAACGCTCTCAATTGAGGCGATATGAACa                                               <  2:668341/150‑1 (MQ=255)
                                                                                       tCCAACGTGCGCATGTCTGGAGGTATATTTGGCGAACTCATTGGATCAGAGCGCTTAAAGGAACCGTGATTATGATCGGTAATAAATACTTTACTTGCAATAAGCGTATCCCGACCTATCGTAACGCTCTCAATTGAGGCGATATGAACa                                               <  1:652254/150‑1 (MQ=255)
                                                                                        ccAACGTGCGCATGTCTGGAGGTATATTTGGCGAACTCATTGGATCAGAGCGCTTAAAGGAACCGTGATTATGATCGGTAATAAATACTTTACTTGCAATAAGCGTATCCCGACCTATCGTAACGCTCTCAATTGAGGCGATATGAACAt                                              <  2:422695/150‑1 (MQ=255)
                                                                                            cGTGCGCATGTCTGGAGGTATATTTGGCGAACTCATTGGATCAGAGCGCTTAAAGGAACCGTGATTATGATCGGTAATAAATACTTTACTTGCAATAAGCGTATCCCGACCTATCGTAACGCTCTCAATTGAGGCGATATGAACATAGTCg                                         >  2:883066/1‑151 (MQ=255)
                                                                                             gTGCGCATGTCTGGAGGTATATTTGGCGAACTCATTGGATCAGAGCGCTTAAAGGAACCGTGATTATGATCGGTAATAAATACTTTACTTGCAATAAGCGTATCCCGACCTATCGTAACGCTCTCAATTGAGGCGATATGAACATAGTCGt                                        <  1:546359/151‑1 (MQ=255)
                                                                                                cgcATGTCTGGAGGTATATTTGGCGAACTCATTGGATCAGAGCGCTTAAAGGAACCGTGATTATGATCGGTAATAAATACTTTACTTGCAATAAGCGTATCCCGACCTATCGTAACGCTCTCAATTGAGGCGATATGAACATAGTCGTTa                                      >  2:74686/1‑150 (MQ=255)
                                                                                                  cATGTCTGGAGGTATATTTGGCGAACTCATTGGATCAGAGCGCTTAAAGGAACCGTGATTATGATCGGTAATAAATACTTTACTTGCAATAAGCGTATCCCGACCTATCGTAACGCTCTCAATTGAGGCGATATGAACATAGTCGTTAACt                                   >  1:500543/1‑151 (MQ=255)
                                                                                                   aTGTCTGGAGGTATATTTGACGAACTCATTGGATCAGAGCGCTTAAAGGAACCGTGATTATGATCGGTAATAAATACTTTACTTGCAAGAAGCGTATCCCGACCTGTCGTAACGCTCTCAATTGAGGCGATATGAACATAGTCGTTAAc                                    <  1:841919/149‑1 (MQ=255)
                                                                                                     gTCTGGAGGTATATTTGGCGAACTCATTGGATCAGAGCGCTTAAAGGAACCGTGATTATGATCGGTAATAAATACTTTACTTGCAATAAGCGTATCCCGACCTATCGTAACGCTCTCAATTGAGGCGATATGAACATAGTCGTTAACTTGc                                >  1:666152/1‑151 (MQ=255)
                                                                                                             gTATATTTGGCGAACTCATTGGATCAGAGCGCTTAAAGGAACCGTGGTTATGATCGGTAATAAATACTTTACTTGCAATAAGCGTATCCCGACCTATCGTAACGCTCTCAATTGAGGCGATATGAACATAGTCGTTAACTTGCACATTATc                        <  1:344169/151‑1 (MQ=255)
                                                                                                             gTATATTTGGCGAACTCATTGGATCAGAGCGCTTAAAGGAACCGTGATTATGATCGGTAATAAATACTTTACTTGCAATAAGCGTATCCCGACCTATCGTAACGCTCTCAATTGAGGCGATATGAACATAGTCGTTAACTTGCACATTATc                        >  1:355271/1‑151 (MQ=255)
                                                                                                               atatTTGGCGAACTCATTGGATCAGAGCGCTTAAAGGAACCGTGATTATGATCGGTAATAAATACTTTACTTGCAATAAGCGTATCCCGACCTATCGTAACGCTCTCAATTGAGGCGATATGAACATAGTCGTTAACTTGCACATTATCgg                      >  1:403228/1‑151 (MQ=255)
                                                                                                                  tttGGCGAACTCATTGGATCAGAGCGCTTAAAGGAACCGTGATTATGATCGGTAATAAATACTTTACTTGCAATAAGCGTATCCCGACCTATCGTAACGCTCTCAATTGAGGCGATATGAACATAGTCGTTAACTTGCACATTATCGGaa                    >  2:319760/1‑150 (MQ=255)
                                                                                                                   ttGGCGAACTCATTGGATCAGAGCGCTTAAAGGAACCGTGATTATGATCGGTAATAAATACTTTACTTGCAATAAGCGTATCCCGACCTATCGTAACGCTCTCAATTGAGGCGATATGAACATAGTCGTTAACTTGCACATTATCGGaaa                   <  2:223112/150‑1 (MQ=255)
                                                                                                                      gCGAACTCATTGGATCAGAGCGCTTAAAGGAACCGTGATTATGATCGGTAATAAATACTTTACTTGCAATAAGCGTATCCCGACCTATCGTAACGCTCTCAATTGAGGCGATATGAACATAGTCGTTAACTTGCACATTATCGGaaaaaa                >  2:446513/1‑150 (MQ=255)
                                                                                                                         aaCTCATTGGATCAGAGCGCTTAAAGGAACCGTGATTATGATCGGTAATAAATACTTTACTTGCAATAAGCGTATCCCGACCTATCGTAACGCTCTCAATTGAGGCGATATGAACATAGACGTTAACTTGCACATTATCGGAAAAAAAAAt            <  1:845257/151‑1 (MQ=255)
                                                                                                                             cATTGGATCAGAGCGCTTAAAGGAACCGTGATTATGATCGGTAATAAATACTTTACTTGCAATAAGCGTATCCCGACCTATCGTAACGCTCTCAATTGAGGCGATATGAACATAGTCGTTAACTTGCACATTATCGGAAAAAAAAATCACg        <  1:158771/151‑1 (MQ=255)
                                                                                                                                tGGATCAGAGCGCTTAAAGGAACCGTGATTATGATCGGTAATAAATACTTTACTTGCAATAAGCGTATCCCGACCTATCGTAACGCTCTCAATTGAGGCGATATGAACATAGTCGTTAACTTGCACATTATCGGAAAAAAAAATCACGcc      >  1:650544/1‑150 (MQ=255)
                                                                                                                                  gATCAGAGCGCTTAAAGGAACCGTGATTATGATCGGTAATAAATACTTTACTTGCAATAAGCGTATCCCGACCTATCGTAACGCTCTCAATTGAGGCGATATGAACATAGTCGTTAACTTGCACATTATCGGAAAAAAAAATCACGCCACg   <  1:338050/151‑1 (MQ=255)
                                                                                                                                   aTCAGAGCGCTTAAAGGAACCGTGATTATGATCGGTAATAAATACTTTACTTGCAATAAGCGTATCCCGACCTATCGTAACGCTCTCAATTGAGGCGATATGAACATAGTCGTTAACTTGCACATTATCGGAAAAAAAAATCACGCCACGt  >  2:35509/1‑151 (MQ=255)
                                                                                                                                          |                                                                                                                                               
ACGACTCCATTACCAATAATTGTTCCAGGCAAAACCGTCACATTCTCACCCAACCAAACCCTCTGGCCAATTACAACAGCTGAAGATTCCAACGTGCGCATGTCTGGAGGTATATTTGGCGAACTCATTGGATCAGAGTGCTTAAAGGAACCGTGATTATGATCGGTAATAAATACTTTACTTGCAATAAGCGTATCCCGACCTATCGTAACGCTCTCAATTGAGGCGATATGAACATAGTCGTTAACTTGCACATTATCGGAAAAAAAAATCACGCCACGT  >  NC_000913/2104612‑2104893

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