Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,649,566 C→T intergenic (+132/‑43) dicF → / → dicB Qin prophage; small regulatory RNA DicF/Qin prophage; cell division inhibition protein DicB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,649,5660CT100.0% 1162.4 / NA 285intergenic (+132/‑43)dicF/dicBQin prophage; small regulatory RNA DicF/Qin prophage; cell division inhibition protein DicB
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (151/134);  total (151/134)

Only 100 of 285 total aligned reads displayed.

GGCTGAGCGCACGCGGAACAGTTAAAACCAAAAACAGTGTTATGGGTGGATTCTCTGTATCCGGCGTTAATTGTTAACTGGTTAACGTCACCTGGAGGCACCAGGCACTGCATCACAAAATTCATTGTTGAGGACGCGATAATGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTG  >  NC_000913/1649468‑1649663
                                                                                                  |                                                                                                 
ggCTGAGCGCACGCGGAACAGTTAAAACCAAAAACAGTGTTATGGGTGGATTCTCTGTATCCGGCGTTAATTGTTAACTGGTTAACGTCACCTGGAGGTAc                                                                                                 >  2:2385820/1‑101 (MQ=255)
  cTGAGCGCACGCGGAACAGTTAAAACCAAAAACAGTGTTATGGGTGGATTCTCTGTATCCGGCGTTAATTGTTAACTGGTTAACGTCACCTGGAGGTACCa                                                                                               >  1:2385822/1‑101 (MQ=255)
   tGAGCGCACGCGGAACAGTTAAAACCAAAAACAGTGTTATGGGTGGATTCTCTGTATCCGGCGTTAATTGTTAACTGGTTAACGTCACCTGGAGGTACCAg                                                                                              <  2:2385504/101‑1 (MQ=255)
    gAGCGCACGCGGAACAGTTAAAACCAAAAACAGTGTTATGGGTGGATTCTCTGTATCCGGCGTTAATTGTTAACTGGTTAACGTCACCTGGAGGTACCAgg                                                                                             >  2:2385824/1‑101 (MQ=255)
     aGCGCACGCGGAACAGTTAAAACCAAAAACAGTGTTATGGGTGGATTCTCTGTATCCGGCGTTAATTGTTAACTGGTTAACGTCACCTGGAGGTACCAGGc                                                                                            <  1:2385129/101‑1 (MQ=255)
     aGCGCACGCGGAACAGTTAAAACCAAAAACAGTGTTATGGGTGGATTCTCTGTATCCGGCGTTAATTGTTAACTGGTTAACGTCACCTGGAGGTACCAGGc                                                                                            <  2:2385406/101‑1 (MQ=255)
     aGCGCACGCGGAACAGTTAAAACCAAAAACAGTGTTATGGGTGGATTCTCTGTATCCGGCGTTAATTGTTAACTGGTTAACGTAACCTGGAGGTACCAGGc                                                                                            <  1:2385513/101‑1 (MQ=255)
      gcgcACGCGGAACAGTTAAAACCAAAAACAGTGTTATGGGTGGATTCTCTGTATCCGGCGTTAATTGTTAACTGGTTAACGTCACCTGGAGGTACCAGGCa                                                                                           <  1:2385461/101‑1 (MQ=255)
      gcgcACGCGGAACAGTTAAAACCAAAAACAGTGTTATGGGTGGATTCTCTGTATCCGGCGTTAATTGTTAACTGGTTAACGTCACCTGGAGGTACCAGGCa                                                                                           <  2:2385636/101‑1 (MQ=255)
       cgcACGCGGAACAGTTAAAACCAAAAACAGTGTTATGGGTGGATTCTCTGTATCCGGCGTTAATTGTTAACTGGTTAACGTCACCTGGAGGTACCAGGCAc                                                                                          >  2:2385828/1‑101 (MQ=255)
        gcACGCGGAACAGTTAAAACCAAAAACAGTGTTATGGGTGGATTCTCTGTATCCGGCGTTAATTGTTAACTGGTTAACGTCACCTGGAGGTACCAGGCACt                                                                                         >  1:2385829/1‑101 (MQ=255)
         cACGCGGAACAGTTAAAACCAAAAACAGTGTTATGGGTGGATTCTCTGTATCCGGCGTTAATTGTTAACTGGTTAACGTCACCTGGAGGTACCAGGCACTg                                                                                        <  2:2385598/101‑1 (MQ=255)
         cACGCGGAACAGTTAAAACCAAAAACAGTGTTATGGGTGGATTCTCTGTATCCGGCGTTAATTGTTAACTGGTTAACGTCACCTGGAGGTACCAGGCACTg                                                                                        >  1:2385832/1‑101 (MQ=255)
         cACGCGGAACAGTTAAAACCAAAAACAGTGTTATGGGTGGATTCTCTGTATCCGGCGTTAATTGTTAACTGGTTAACGTCACCTGGAGGTACCAGGCACTg                                                                                        <  2:2385431/101‑1 (MQ=255)
           cgcgGAACAGTTAAAACCAAAAACAGTGTTATGGGTGGATTCTCTGTATCCGGCGTTAATTGTTAACTGGTTAACGTCACCTGGAGGTACCAGGCACTGCa                                                                                      >  1:2385838/1‑101 (MQ=255)
           cgcgGAACAGTTAAAACCAAAAACAGTGTTATGGGTGGATTCTCTGTATCCGGCGTTAATTGTTAACTGGTTAACGTCACCTGGAGGTACCAGGCACTGCa                                                                                      >  2:2385837/1‑101 (MQ=255)
           cgcgGAACAGTTAAAACCAAAAACAGTGTTATGGGTGGATTCTCTGTATCCGGCGTTAATTGTTAACTGGTTAACGTCACCTGGAGGTACCAGGCACTGCa                                                                                      <  1:2385578/101‑1 (MQ=255)
              ggAACAGTTAAAACCAAAAACAGTGTTATGGGTGGATTCTCTGTATCCGGCGTTAATTGTTAACTGGTTAACGTCACCTGGAGGTACCAGGCACTGCATca                                                                                   <  2:2385262/101‑1 (MQ=255)
              ggAACAGTTAAAACCAAAAACAGTGTTATGGGTGGATTCTCTGTATCCGGCGTTAATTGTTAACTGGTTAACGTCACCTGGAGGTACCAGGCACTGCATca                                                                                   <  1:2385410/101‑1 (MQ=255)
                aaCAGTTAAAACCAAAAACAGTGTTATGGGTGGATTCTCTGTATCCGGCGTTAATTGTTAACTGGTTAACGTCACCTGGAGGTACCAGGCACTGCATcaca                                                                                 >  2:2385844/1‑101 (MQ=255)
                aaCAGTTAAAACCAAAAACAGTGTTATGGGTGGATTCTCTGTATCCGGCGTTAATTGTTAACTGGTTAACGTCACCTGGAGGTACCAGGCACTGCATcaca                                                                                 >  1:2385843/1‑101 (MQ=255)
                 aCAGTTAAAACCAAAAACAGTGTTATGGGTGGATTCTCTGTATCCGGCGTTAATTGTTAACTGGTTAACGTCACCTGGAGGTACCAGGCACTGCATCACaa                                                                                <  1:2385289/101‑1 (MQ=255)
                 aCAGTTAAAACCAAAAACAGTGTTATGGGTGGATTCTCTGTATCCGGCGTTAATTGTTAACTGGTTAACGTCACCTGGAGGTACCAGGCACTGCATCACaa                                                                                >  1:2385848/1‑101 (MQ=255)
                  cAGTTAAAACCAAAAACAGTGTTATGGGTGGATTCTCTGTATCCGGCGTTAATTGTTAACTGGTTAACGTCACCTGGAGGTACCAGGCACTGCATCACaaa                                                                               >  2:2385850/1‑101 (MQ=255)
                   aGTTAAAACCAAAAACAGTGTTATGGGTGGATTCTCTGTATCCGGCGTTAATTGTTAACTGGTTAACGTCACCTGGAGGTACCAGGCACTGCATCACaaaa                                                                              >  1:2385853/1‑101 (MQ=255)
                   aGTTAAAACCAAAAACAGTGTTATGGGTGGATTCTCTGTATCCGGCGTTAATTGTTAACTGGTTAACGTCACCTGGAGGTACCAGGCACTGCATCACaaaa                                                                              <  1:2385387/101‑1 (MQ=255)
                     ttAAAACCAAAAACAGTGTTATGGGTGGATTCTCTGTATCCGGCGTTAATTGTTAACTGGTTAACGTCACCTGGAGGTACCAGGCACTGCATCACAAAAtt                                                                            >  1:2385857/1‑101 (MQ=255)
                      tAAAACCAAAAACAGTGTTATGGGTGGATTCTCTGTATCCGGCGTTAATTGTTAACTGGTTAACGTCACCTGGAGGTACCAGGCACTGCATCACAAAATTc                                                                           >  2:2385859/1‑101 (MQ=255)
                        aaaCCAAAAACAGTGTTATGGGTGGATTCTCTGTATCCGGCGTTAATTGTTAACTGGTTAACGTCACCTGGAGGTACCAGGCACTGCATCACAAAATTCAt                                                                         >  1:2385864/1‑101 (MQ=255)
                        aaaCCAAAAACAGTGTTATGGGTGGATTCTCTGTATCCGGCGTTAATTGTTAACTGGTTAACGTCACCTGGAGGTACCAGGCACTGCATCACAAAATTCAt                                                                         >  2:2385865/1‑101 (MQ=255)
                         aaCCAAAAACAGTGTTATGGGTGGATTCTCTGTATCCGGCGTTAATTGTTAACTGGTTAACGTCACCTGGAGGTACCAGGCACTGCATCACAAAATTCAtt                                                                        <  1:2385526/101‑1 (MQ=255)
                         aaCCAAAAACAGTGTTATGGGTGGATTCTCTGTATCCGGCGTTAATTGTTAACTGGTTAACGTCACCTGGAGGTACCAGGCACTGCATCACAAAATTCAtt                                                                        <  2:2385465/101‑1 (MQ=255)
                          aCCAAAAACAGTGTTATGGGTGGATTCTCTGTATCCGGCGTTAATTGTTAACTGGTTAACGTCACCTGGAGGTACCAGGCACTGCATCACAAAATTCAttg                                                                       >  2:2385868/1‑101 (MQ=255)
                           ccAAAAACAGTGTTATGGGTGGATTCTCTGTATCCGGCGTTAATTGTTAACTGGTTAACGTCACCTGGAGGTACCAGGCACTGCATCACAAAATTCAttgt                                                                      <  1:2385437/101‑1 (MQ=255)
                            cAAAAACAGTGTTATGGGTGGATTCTCTGTATCCGGCGTTAATTGTTAACTGGTTAACGTCACCTGGAGGTACCAGGCACTGCATCGCAAAATTCAttgtt                                                                     >  1:2385871/1‑101 (MQ=255)
                            cAAAAACAGTGTTATGGGTGGATTCTCTGTATCCGGCGTTAATTGTTAACTGGTTAACGTCACCTGGAGGTACCAGGCACTGCATCACAAAATTCAttgtt                                                                     <  1:2385352/101‑1 (MQ=255)
                            cAAAAACAGTGTTATGGGTGGATTCTCTGTATCCGGCGTTAATTGTTAACTGGTTAACGTCACCTGGAGGTACCAGGCACTGCATCACAAAATTCAttgtt                                                                     <  2:2385557/101‑1 (MQ=255)
                               aaaCAGTGTTATGGGTGGATTCTCTGTATCCGGCGTTAATTGTTAACTGGTTAACGTCACCTGGAGGTACCAGGCACTGCATCACAAAATTCATTGTTGAg                                                                  >  2:2385874/1‑101 (MQ=255)
                                aaCAGTGTTATGGGTGGATTCTCTGTATCCGGCGTTAATTGTTAACTGGTTAACGTCACCTGGAGGTACCAGGCACTGCATCACAAAATTCATTGTTGAgg                                                                 >  2:2385877/1‑101 (MQ=255)
                                    gtgtTATGGGTGGATTCTCTGTATCCGGCGTTAATTGTTAACTGGTTAACGTCACCTGGAGGTACCAGGCACTGCATCACAAAATTCATTGTTGAGGAcgc                                                             >  1:2385880/1‑101 (MQ=255)
                                     tgtTATGGGTGGATTCTCTGTATCCGGCGTTAATTGTTAACTGGTTAACGTCACCTGGAGGTACCAGGCACTGCATCACAAAATTCATTGTTGAGGAcgcg                                                            <  2:2385480/101‑1 (MQ=255)
                                      gtTATGGGTGGATTCTCTGTATCCGGCGTTAATTGTTAACTGGTTAACGTCACCTGGAGGTACCAGGCACTGCATCACAAAATTCATTGTTGAGGACGCGa                                                           >  1:2385883/1‑101 (MQ=255)
                                      gtTATGGGTGGATTCTCTGTATCCGGCGTTAATTGTTAACTGGTTAACGTCACCTGGAGGTACCAGGCACTGCATCACAAAATTCATTGTTGAGGACGCGa                                                           <  1:2385398/101‑1 (MQ=255)
                                        tATGGGTGGATTCTCTGTATCCGGCGTTAATTGTTAACTGGTTAACGTCACCTGGAGGTACCAGGCACTGCATCACAAAATTCATTGTTGAGGACGCGATa                                                         >  1:2385887/1‑101 (MQ=255)
                                        tATGGGTGGAATCTCTGTATCCGGCGTTAATTGTTAACTGGTTAACGTCACCTGGAGGTACCAGGCACTGCATCACAAAATTCATTGTTGAGGACGCGATa                                                         >  1:2385889/1‑101 (MQ=255)
                                         aTGGGTGGATTCTCTGTATCCGGCGTTAATTGTTAACTGGTTAACGTCACCTGGAGGTAc                                                                                                 >  2:2385890/1‑60 (MQ=255)
                                         aTGGGTGGACTCTCTGTATCCGGCGTTAATTGTTAACTGGTTAACGTCACCTGGAGGTACCAGGCACTGCATCACAAAATTCATTGTTGAGGACGCGATaa                                                        <  2:2385336/101‑1 (MQ=255)
                                          tGGGTGGATTCTCTGTATCCGGCGTTAATTGTTAACTGGTTAACGTCACCTGGAGGTACCAGGCACTGCATCACAAAATTCATTGTTGAGGACGCGATAAt                                                       <  2:2385493/101‑1 (MQ=255)
                                           gggTGGATTCTCTGTATCCGGCGTTAATTGTTAACTGGTTAACGTCACCTGGAGGTACCAGGCACTGCATCACAAAATTCATTGTTGAGGACGCGATAATg                                                      <  1:2385363/101‑1 (MQ=255)
                                           gggTGGATTCTCTGTATCCGGCGTTAATTGTTAACTGGTTAACGTCACCTGGAGGTACCAGGCACTGCATCACAAAATTCATTGTTGAGGACGCGATAATg                                                      <  2:2385673/101‑1 (MQ=255)
                                           gggTGGATTCTCTGTATCCGGCGTTAATTGTTAACTGGTTAACGTCACCTGGAGGTACCAGGCACTGCATCACAAAATTCATTGTTGAGGACGCGATAATg                                                      <  2:2385516/101‑1 (MQ=255)
                                           gggTGGATTCTCTGTATCCGGCGTTAATTGTTAACTGGTTAACGTCACCTGGAGGTACCAGGCACTGCATCACAAAATTCATTGTTGAGGACGCGATAATg                                                      <  1:2385498/101‑1 (MQ=255)
                                           gggTGGATTCTCTGTATCCGGCGTTAATTGTTAACTGGTTAACGTCACCTGGAGGTACCAGGCACTGCATCACAAAATTCATTGTTGAGGACGCGATAATg                                                      <  1:2385537/101‑1 (MQ=255)
                                            ggTGGATTCTCTGTATCCGGCGTTAATTGTTAACTGGTTAACGTCACCTGGAGGTACCAGGCACTGCATCACAAAATTCATTGTTGAGGACGCGATAATGa                                                     >  2:2385895/1‑101 (MQ=255)
                                             gTGGATTCTCTGTATCCGGCGTTAATTGTTAACTGGTTAACGTCACCTGGAGGTACCAGGCACTGCATCACAAAATTCATTGTTGAGGACGCGATAATGaa                                                    >  1:2385897/1‑101 (MQ=255)
                                              tGGATTCTCTGTATCCGGCGTTAATTGTTAACTGGTTAACGTCACCTGGAGGTACCAGGCACTGCATCACAAAATTCATTGTTGAGGACGCGATAATGaaa                                                   <  1:2385561/101‑1 (MQ=255)
                                               ggATTCTCTGTATCCGGCGTTAATTGTTAACTGGTTAACGTCACCTGGAGGTACCAGGCACTGCATCACAAAATTCATTGTTGAGGACGCGATAATGaaaa                                                  >  2:2385899/1‑101 (MQ=255)
                                               ggATTCTCTGTATCCGGCGTTAATTGTTAACTGGTTAACGTCACCTGGAGGTACCAGGCACTGCATCACAAAATTCATTGTTGAGGACGCGATAATGaaaa                                                  <  2:2385385/101‑1 (MQ=255)
                                                gATTCTCTGTATCCGGCGTTAATTGTTAACTGGTTAACGTCACCTGGAGGTACCAGGCACTGCATCACAAAATTCATTGTTGAGGACGCGATAATGAAAAc                                                 >  2:2385901/1‑101 (MQ=255)
                                                    ctctGTATCCGGCGTTAATTGTTAACTGGTTAACGTCACCTGGAGGTACCAGGCACTGCATCACAAAATTCATTGTTGAGGACGCGATAATGAAAACGtta                                             >  2:2385906/1‑101 (MQ=255)
                                                    ctctGTATCCGGCGTTAATTGTTAACTGGTTAACGTCACCTGGAGGTACCAGGCACTGCATCACAAAATTCATTGTTGAGGACGCGATAATGAAAACGtta                                             <  2:2385448/101‑1 (MQ=255)
                                                     tctGTATCCGGCGTTAATTGTTAACTGGTTAACGTCACCTGGAGGTACCAGGCACTGCATCACAAAATTCATTGTTGAGGACGCGATAATGAAAACGttat                                            >  1:2385909/1‑101 (MQ=255)
                                                     tctGTATCCGGCGTTAATTGTTAACTGGTTAACGTCACCTGGAGGTACCAGGCACTGCATCACAAAATTCATTGTTGAGGACGCGATAATGAAAACGttat                                            <  1:2385470/101‑1 (MQ=255)
                                                      ctGTATCCGGCGTTAATTGTTAACTGGTTAACGTCACCTGGAGGTACCAGGCACTGCATCACAAAATTCATTGTTGAGGACGCGATAATGAAAACGttatt                                           <  1:2385601/101‑1 (MQ=255)
                                                      ctGTATCCGGCGTTAATTGTTAACTGGTTAACGTCACCTGGAGGTACCAGGCACTGCATCACAAAATTCATTGTTGAGGACGCGATAATGAAAACGttatt                                           <  2:2385653/101‑1 (MQ=255)
                                                        gTATCCGGCGTTAATTGTTAACTGGTTAACGTCACCTGGAGGTACCAGGCACTGCATCACAAAATTCATTGTTGAGGACGCGATAATGAAAACGTTATTAc                                         <  2:2385215/101‑1 (MQ=255)
                                                         tATCCGGCGTTAATTGTTAACTGGTTAACGTCACCTGGAGGTACCAGGCACTGCATCACAAAATTCATTGTTGAGGACGCGATAATGAAAACGTTATTAcc                                        >  2:2385912/1‑101 (MQ=255)
                                                         tATCCGGCGTTAATTGTTAACTGGTTAACGTCACCTGGAGGTACCAGGCACTGCATCACAAAATTCATTGTTGAGGACGCGATAATGAAAACGTTATTAcc                                        >  1:2385913/1‑101 (MQ=255)
                                                           tCCGGCGTTAATTGTTAACTGGTTAACGTCACCTGGAGGTACCAGGCACTGCATCACAAAATTCATTGTTGAGGACGCGATAATGAAAACGTTATTACCaa                                      <  2:2385573/101‑1 (MQ=255)
                                                                cGTTAATTGTTAACTGGTTAACGTCACCTGGAGGTACCAGGCACTGCATCACAAAATTCATTGTTGAGGACGCGATAATGAAAACGTTATTACCAAACGtt                                 <  1:2385695/101‑1 (MQ=255)
                                                                cGTTAATTGTTAACTGGTTAACGTCACCTGGAGGTACCAGGCACTGCATCACAAAATTCATTGTTGAGGACGCGATAATGAAAACGTTATTACCAAACGtt                                 >  1:2385918/1‑101 (MQ=255)
                                                                cGTTAATTGTTAACTGGTTAACGTCACCTGGAGGTACCAGGCACTGCATCACAAAATTCATTGTTGAGGACGCGATAATGAAAACGTTATTACCAAACGtt                                 >  2:2385919/1‑101 (MQ=255)
                                                                 gTTAATTGTTAACTGGTTAACGTCACCTGGAGGTACCAGGCACTGCATCACAAAATTCATTGTTGAGGACGCGATAATGAAAACGTTATTACCAAACGTTa                                >  2:2385922/1‑101 (MQ=255)
                                                                 gTTAATTGTTAACTGGTTAACGTCACCTGGAGGTACCAGGCACTGCATCACAAAATTCATTGTTGAGGACGCGATAATGAAAACGTTATTACCAAACGTTa                                <  2:2385576/101‑1 (MQ=255)
                                                                   tAATTGTTAACTGGTTAACGTCACCTGGAGGTACCAGGCACTGCATCACAAAATTCATTGTTGAGGACGCGATa                                                         >  1:2385931/1‑74 (MQ=255)
                                                                   tAATTGTTAACTGGTTAACGTCACCTGGAGGTACCAGGCACTGCATCACAAAATTCATTGTTGAGGACGCGATAATGAAAACGTTATTACCAAACGTTAAt                              >  2:2385930/1‑101 (MQ=255)
                                                                   tAATTGTTAACTGGTTAACGTCACCTGGAGGTACCAGGCACTGCATCACAAAATTCATTGTTGAGGACGCGATAATGAAAACGTTATTACCAAACGTTAAt                              >  2:2385927/1‑101 (MQ=255)
                                                                    aaTTGTTAACTGGTTAACGTCACCTGGAGGTACCAGGCACTGCATCACAAAATTCATTGTTGAGGACGCGATAATGAAAACGTTATTACCAAACGTTAATa                             >  2:2385934/1‑101 (MQ=255)
                                                                      ttGTTAACTGGTTAACGTCACCTGGAGGTACCAGGCACTGCATCACAAAATTCATTGTTGAGGACGCGATAATGAAAACGTTATTACCAAACGTTAATACg                           >  1:2385936/1‑101 (MQ=255)
                                                                        gTTAACTGGTTAACGTCACCTGGAGGTACCAGGCACTGCATCACAAAATTCATTGTTGAGGACGCGATAATGAAAACGTTATTACCAAACGTTAATACGTc                         >  1:2385939/1‑101 (MQ=255)
                                                                        gTTAACTGGTTAACGTCACCTGGAGGTACCAGGCACTGCATCACAAAATTCATTGTTGAGGACGCGATAATGAAAACGTTATTACCAAACGTTAATACGTc                         <  2:2385369/101‑1 (MQ=255)
                                                                         ttAACTGGTTAACGTCACCTGGAGGTACCAGGCACTGCATCACAAAATTCATTGTTGAGGACGCGATAATGAAAACGTTATTACCAAACGTTAATACGTCt                        <  2:2385692/101‑1 (MQ=255)
                                                                          tAACTGGTTAACGTCACCTGGAGGTACCAGGCACTGCATCACAAAATTCATTGTTGAGGACGCGATAATGAAAACGTTATTACCAAACGTTAATACGTCTg                       >  1:2385942/1‑101 (MQ=255)
                                                                          tAACTGGTTAACGTCACCTGGAGGTACCAGGCACTGCATCACAAAATTCATTGTTGAGGACGCGATAATGAAAACGTTATTACCAAACGTTAATACGTCTg                       <  2:2385530/101‑1 (MQ=255)
                                                                           aaCTGGTTAACGTCACCTGGAGGTACCAGGCACTGCATCACAAAATTCATTGTTGAGGACGCGATAATGAAAACGTTATTACCAAACGTTTATACGTCTGa                      <  1:2385484/101‑1 (MQ=255)
                                                                             cTGGTTAACGTCACCTGGAGGTACCAGGCACTGCATCACAAAATTCATTGTTGAGGACGCGATAATGAAAACGTTATTACCAAACGTTAATACGTCTGAAg                    >  2:2385946/1‑101 (MQ=255)
                                                                             cTGGTTAACGTCACCTGGAGGTACCAGGCACTGCATCACAAAATTCATTGTTGAGGACGCGATAATGAAAACGTTATTACCAAACGTTAATACGTCTGAAg                    <  1:2385245/101‑1 (MQ=255)
                                                                              tGGTTAACGTCACCTGGAGGTACCAGGCACTGCATCACAAAATTCATTGTTGAGGACGCGATAATGAAAACGTTATTACCAAACGTTAATACGTCTGAAgg                   <  1:2385630/101‑1 (MQ=255)
                                                                                gTTAACGTCACCTGGAGGTACCAGGCACTGCATCACAAAATTCATTGTTGAGGACGCGATAATGAAAACGTTATTACCAAACGTTAATACGTCTGAAGgtt                 >  1:2385948/1‑101 (MQ=255)
                                                                                 ttAACGTCACCTGGAGGTACCAGGCACTGCATCACAAAATTCATTGTTGAGGACGCGATAATGAAAACGTTATTACCAAACGTTAATACGTCTGAAGgttg                <  2:2385512/101‑1 (MQ=255)
                                                                                  tAACGTCACCTGGAGGTACCAGGCACTGCATCACAAAATTCATTGTTGAGGACGCGATAATGAAAACGTTATTACCAAACGTTAATACGTCTGAAGgttgt               >  2:2385952/1‑101 (MQ=255)
                                                                                      gTCACCTGGAGGTACCAGGCACTGCATCACAAAATTCATTGTTGAGGACGCGATAATGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTg           <  2:2385543/101‑1 (MQ=255)
                                                                                      gTCACCTGGAGGTACCAGGCACTGCATCACAAAATTCATTGTTGAGGACGCGATAATGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTg           >  1:2385955/1‑101 (MQ=255)
                                                                                      gTCACCTGGAGGTACCAGGCACTGCATCACAAAATTCATTGTTGAGGACGCGATAATGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTg           >  2:2385956/1‑101 (MQ=255)
                                                                                       tCACCTGGAGGTACCAGGCACTGCATCACAAAATTCATTGTTGAGGACGCGATAATGAAAACGTTATTACCaaa                                     <  2:2385297/74‑1 (MQ=255)
                                                                                       tCACCTGGAGGTACCAGGCACTGCATCACAAAATTCATTGTTGAGGACGCGATAATGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGa          >  1:2385960/1‑101 (MQ=255)
                                                                                           cTGGAGGTACCAGGCACTGCATCACAAAATTCATTGTTGAGGACGCGATAATGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAAtt      >  2:2385964/1‑101 (MQ=255)
                                                                                            tGGAGGTACCAGGCACTGCATCACAAAATTCATTGTTGAGGACGCGATAATGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTg     <  2:2385550/101‑1 (MQ=255)
                                                                                            tGGAGGTACCAGGCACTGCATCACAAAATTCATTGTTGAGGACGCGATAATGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTg     >  1:2385967/1‑101 (MQ=255)
                                                                                            tGGAGGTACCAGGCACTGCATCACAAAATTCATTGTTGAGGACGCGATAATGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTg     >  2:2385968/1‑101 (MQ=255)
                                                                                                  |                                                                                                 
GGCTGAGCGCACGCGGAACAGTTAAAACCAAAAACAGTGTTATGGGTGGATTCTCTGTATCCGGCGTTAATTGTTAACTGGTTAACGTCACCTGGAGGCACCAGGCACTGCATCACAAAATTCATTGTTGAGGACGCGATAATGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTG  >  NC_000913/1649468‑1649663

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