Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,296,292 C→T R162C (CGC→TGC)  diaA → DnaA initiator‑associating factor for replication initiation

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,296,2920CT100.0% 245.2 / NA 65R162C (CGC→TGC) diaADnaA initiator‑associating factor for replication initiation
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (35/30);  total (35/30)

AGCCATTTCCACCCGTGGCAACAGCCGCGATATTGTTAAAGCAGTTGAAGCCGCCGTTACGCGTGATATGACCATTGTGGCATTGACCGGCTATGACGGCGGCGAACTTGCAGGTTTGTTAGGGCCACAGGATGTGGAGATCCGCATTCCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCTGTGCGATCTGATCGATAACACGCTTTTCCCTCACCAGGATGATTAAGGAGAATACATGAAGGCATTATCGCCAATCGCAGTCCTT  >  NC_000913/3296150‑3296438
                                                                                                                                              |                                                                                                                                                  
aGCCATTTCCACCCGTGGCAACAGCCGCGATATTGTTAAAGCAGTTGAAGCCGCCGTTACGCGTGATATGACCATTGTGGCATTGACCGGCTATGACGGCGGCGAACTTGCAGGTTTGTTAGGGCCACAGGATGTGGAGATCTGCATTcc                                                                                                                                             <  2:184213/150‑1 (MQ=255)
aGCCATTTCCACCCGTGGCAACAGCCGCGATATTGTTAAAGCAGTTGAAGCCGCCGTTACGCGTGATATGACCATTGTGGCATTGACCGGCTATGACGGCGGCGAACTTGCAGGTTTGTTAGGGCCACAGGATGTGGAGATCTGCATTcc                                                                                                                                             <  1:626470/150‑1 (MQ=255)
   cATTTCCCCCCGTGGCAACAGCCGCGATATTGTTAAAGCAGTTGAAGCCGCCGTTACGCGTGATATGACCATTGTGGCATTGACCGGCTATGACGGCGGCGAACTTGCAGGTTTGTTAGGGCCACAGGATGTGGAGATCTGCATTCCTTCg                                                                                                                                         <  1:484034/151‑1 (MQ=255)
   cATTTCCACCCGTGGCAACAGCCGCGATATTGTTAAAGCAGTTGAAGCCGCCGTTACGCGTGATATGACCATTGTGGCATTGACCGGCTATGACGGCGGCGAACTTGCAGGTTTGTTAGGGCCACAGGATGTGGAGATCTGCATTCCTTCg                                                                                                                                         >  1:826901/1‑151 (MQ=255)
      tcccACCCGGGGCAACAGCCGCGATATTGTTAAAGCAGTTGAAGCCGCCGTTACGCGTGATATGACCATTGTGGCATTGACCGGCTATGACGGCGGCGAACTTGCAGGTTTGTTAGGGCCACAGGATGTGGAGATCTGCATTCCTTCGCa                                                                                                                                       <  2:872102/148‑1 (MQ=255)
        ccccccGTGGCAACAGCCGCGATATTGTTAAAGCGGTGGAAGCCGCCGTTACGCGTGATATGACCATTGTGGCATTGACCGGCTATGACGGCGGCGAACTTGCAGGTTTGTTAGGGCCACAGGATGTGGAGATCTGCATTCCTTCGCATCg                                                                                                                                    <  2:688734/148‑1 (MQ=255)
          aCCCGTGGCAACAGCCGCGATATTGTTAAAGCAGTTGAAGCCGCCGTTACGCGTGATATGACCATTGTGGCATTGACCGGCTATGACGGCGGCGAACTTGCAGGTTTGTTAGGGCCACAGGATGTGGAGATCTGCATTCCTTCGCATCGt                                                                                                                                   >  2:527635/1‑150 (MQ=255)
            ccGTGGCAACAGCCGCGATATTGTTAAAGCAGTTGAAGCCGCCGTTACGCGTGATATGACCATTGTGGCATTGACCGGCTATGACGGCGGCGAACTTGCAGGTTTGTTAGGGCCACAGGATGTGGAGATCTGCATTCCTTCGCATCGTAGt                                                                                                                                >  2:421239/1‑151 (MQ=255)
               tGGCAACAGCCGCGATATTGTTAAAGCAGTTGAAGCCGCCGTTACGCGTGATATGACCATTGTGGCATTGACCGGCTATGACGGCGGCGAACTTGCAGGTTTGTTAGGGCCACAGGACGTGGAGATCTGCATTCCTTCGCATCGTAGTGCt                                                                                                                             >  2:191884/1‑151 (MQ=255)
                   aaCAGCCGCGATATTGTTAAAGCAGTTGAAGCCGCCGTTACGCGTGATATGACCATTGTGGCATTGACCGGCTATGACGGCGGCGAACTTGCAGGTTTGTTAGGGCCACAGGATGTGGAGATCTGCATTCCTTCGCATCGTAGTGCTCGCa                                                                                                                         <  1:117289/151‑1 (MQ=255)
                      aGCCGCGATATTGTTAAAGCAGTTGAAGCCGCCGTTACGCGTGATATGACCATTGTGGCATTGACCGGCTATGACGGCGGCGAACTTGCAGGTTTGTTAGGGCCACAGGATGTGGAGATCTGCATTCCTTCGCATCGTAGTGCTCGCATTc                                                                                                                      >  2:212291/1‑151 (MQ=255)
                         cgcgATATTGTTAAAGCAGTTGAAGCCGCCGTTACGCGTGATATGACCATTGTGGCATTGACCGGCTATGACGGCGGCGAACTTGCAGGTTTGTTAGGGCCACAGGATGTGGAGATCTGCATTCCTTCGCATCGTAGTGCTCGCATTCAgg                                                                                                                   <  1:102106/151‑1 (MQ=255)
                               atTGTTAAAGCAGTTGAAGCCGCCGTTACGCGTGATATGACCATTGTGGCATTGACCGGCTATGACGGCGGCGAACTTGCAGGTTTGGTAGGGCCACATGATGGGGAGATCTTCATTCCTTCGCATCGCAGTGCGCGCATTCAGGAAATGc                                                                                                             >  2:96896/1‑151 (MQ=255)
                                    tAAAGCAGTTGAAGCCGCCGTTACGCGTGATATGACCATTGTGGCATTGACCGGCTATGACGGCGGCGAACTTGCAGGTTTGTTAGGGCCACAGGATGTGGAGATCTGCATTCCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATg                                                                                                        >  2:201350/1‑151 (MQ=255)
                                     aaaGCAGTTGAAGCCGCCGTTACGCGTGATATGACCATTGTGGCATTGACCGGCTATGACGGCGGCGAACTTGCAGGTTTGTTAGGGCCACAGGATGTGGAGATCTGCATTCCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGc                                                                                                       >  1:666824/1‑151 (MQ=255)
                                      aaGCAGTTGAAGCCGCCGTTACGCGTGATATGACCATTGTGGCATTGACCGGCTATGACGGCGGCGAACTTGCAGGTTTGTTAGGGCCACAGGATGTGGAGATCTGCATTCCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCt                                                                                                      >  2:635477/1‑151 (MQ=255)
                                        gCAGTTGAAGCAGCCGTTACGCGTGATATGACCATTGTGGCATTTACCGGCTTTGACGGCGGCGGACTTGCAGGTTTGTTAGGGCCCCAGGATGTGGATATCTGCATTGCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTg                                                                                                     >  2:601227/1‑150 (MQ=255)
                                          aGTTGAAGCCGCCGTTACGCGTGATATGACCATTGTGGCATTGACCGGCTATGACGGCGGCGAACTTGCAGGTTTGTTAGGGCCACAGGATGTGGAGATCTGCATTCCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGAc                                                                                                   <  1:786936/150‑1 (MQ=255)
                                          aGTTGAAGCCGCCGTTACGCGTGATATGACCATTGTGGCATTGACCGGCTATGACGGCGGCGAACTTGCAGGTTTGTTAGGGCCACAGGATGTGGAGATCTGCATTCCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACg                                                                                                  >  1:202426/1‑151 (MQ=255)
                                              gAAGCCGCCGTTACGCGTGATATGACCATTGTGGCATTGACCGGCTATGACGGCGGCGAACTTGCAGGTTTGTTAGGGCCACAGGATGTGGAGATCTGCATTCCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATACTGACGGTaa                                                                                              <  2:759190/151‑1 (MQ=255)
                                                 gccgccGTTACGCGTGATATGACCATTGTGGCATTGACCGGCTATGACGGCGGCGAACTTGCAGGTTTGTTAGGGCCACAGGATGTGGAGATCTGCATTCCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAAtt                                                                                           >  1:415876/1‑151 (MQ=255)
                                                 gccgccGTTACGCGTGATATGACCACTGTGGCATTGACCGGCTATGACGGCGGCGAACTTGCAGGTTTGTTAGGGCCACAGGATGTGGAGATCTGCATTCCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAAtt                                                                                           >  2:63225/1‑151 (MQ=255)
                                                  ccgccgTTACGCGTGATATGACCATTGTGGCATTGACCGGCTATGACGGCGGCGAACTTGCAGGTTTGTTAGGGCCACAGGATGTGGAGATCTGCATTCCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTg                                                                                          <  2:650653/151‑1 (MQ=255)
                                                      cgTTACGCGTGATATGACCATTGTGGCATTGACCGGCTATGACGGCGGCGACCTTGCAGGTTTGTTAGGGCCACAGGATGTGGAGATCTCCATTCCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGGAAAtt                                                                                           >  1:507023/1‑146 (MQ=255)
                                                      cgTTACGCGTGATATGACCATTGTGGCATTGACCGGCTATGACGGCGGCGAACTTGCAGGTTTGTTAGGGCCACAGGATGTGGAGATCTGCATTCCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCtg                                                                                      <  2:764253/151‑1 (MQ=255)
                                                       gTTACGCGTGATATGACCATTGTGGCATTGACCGGCTATGACGGCGGCGAACTTGCAGGTTTGTTAGGGCCACAGGATGTGGAGATCTGCATTCCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCtgt                                                                                     <  2:239075/151‑1 (MQ=255)
                                                        ttACGCGTGATATGACCATTGTGGCATTGACCGGCTATGACGGCGGCGAACTTGCAGGTTTGTTAGGGCCACAGGATGTGGAGATCTGCATTCCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCtgtg                                                                                    >  1:61516/1‑151 (MQ=255)
                                                             cgTGATATGACCATTGTGGCATTGACCGGCTATGACGGCGGCGAACTTGCAGGTTTGTTAGGGCCACAGGATGTGGAGATCTGCATTCCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCTGTGCGATc                                                                               >  2:794414/1‑151 (MQ=255)
                                                              gTGATATGACCATTGTGGCATTGACCGGCTATGACGGCGGCGAACTTGCAGGTTTGTTAGGGCCACAGGATGTGGAGATCTGCATTCCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCTGTGCGATCt                                                                              <  1:908678/151‑1 (MQ=255)
                                                               tGATATGACCATTGTGGCATTGACCGGCTATGACGGCGGCGAACTTGCAGGTTTGTAGGGCCACAGGATGTGGAGATCTGCATTCCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCTGTGCGATCTg                                                                             >  2:794044/1‑150 (MQ=255)
                                                                   atGACCATTGTGGCATTGACCGGCTATGACGGCGGCGAACTTGCAGGTTTGTTAGGGCCACAGGATGTGGAGATCTGCATTCCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCTGTGCGATCTGATCg                                                                         >  2:811630/1‑151 (MQ=255)
                                                                      aCCATTGTGGCATTGACCGGCTATGACGGCGGCGAACTTGCAGGTTTGTTAGGGCCACAGGATGTGGAGATCTGCATTCCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCTGTGCGATCTGATCGATa                                                                      <  2:229473/151‑1 (MQ=255)
                                                                          ttGTGGCATTGACCGGCTATGACGGCGGCGAACTTGCAGGTTTGTTAGGGCCACAGGATGTGGAGATCTGCATTCCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCTGTGCGATCTGATCGATAacac                                                                  >  1:630755/1‑151 (MQ=255)
                                                                              ggCATTGACCGGCTATGACGGCGGCGAACTTGCAGGTTTGTTAGGGCCACAGGATGTGGAGATCTGCATTCCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCTGTGCGATCTGATCGATAACACGCtt                                                              >  1:397699/1‑151 (MQ=255)
                                                                                   tGACCGGCTATGACGGCGGCGAACTTGCAGGTTTGTTAGGGCCACAGGATGTGGAGATCTGCATTCCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCTGTGCGATCTGATCGATAACACGCTTTTccc                                                         >  2:453218/1‑151 (MQ=255)
                                                                                    gACCGGCTATGACGGCGGCGAACTTGCAGGTTTGTTAGGGCCACAGGATGTGGAGATCTGCATTCCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCTGTGCGATCTGATCGATAACACGCTTTTCCCt                                                        >  2:885547/1‑151 (MQ=255)
                                                                                    gACCGGCTATGACGGCGGCGAACTTGCAGGTTTGTTAGGGCCACAGGATGTGGAGATCTGCATTCCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCTGTGCGATCTGATCGATAACACGCTTTTCCCt                                                        <  2:477576/151‑1 (MQ=255)
                                                                                       cGGCTATGACGGCGGCGAACTTGCAGGTTTGTTAGGGCCACAGGATGTGGAGATCTGCATTCCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCTGTGCGATCTGATCGATAACACGCTTTTCCCTCAc                                                     >  2:391545/1‑151 (MQ=255)
                                                                                          cTATGACGGCGGCGAACTTGCAGGTTTGTTAGGGCCACAGGATGTGGAGATCTGCATTCCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCTGTGCGATCTGATCGATAACACGCTTTTCCCTCAcc                                                    >  2:710296/1‑149 (MQ=255)
                                                                                           tATGACGGCGGCGAACTTGCAGGTTTGTTAGGGCCACAGGATGTGGAGATCTGCATTCCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCTGTGCGATCTGATCGATAACACGCTTTTCCCTCACCAgg                                                 >  2:516410/1‑151 (MQ=255)
                                                                                            aTGACGGCGGCGAACTTGCAGGTTTGTTAGGGCCACAGGATGTGGAGATCTGCATTCCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCTGTGCGATCTGATCGATAACACGCTTTTCCCTCACCAgg                                                 >  1:906792/1‑150 (MQ=255)
                                                                                            aTGACGGCGGCGAACTTGCAGGTTTGTTAGGGCCACAGGATGTGGAGATCTGCATTCCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCTGTGCGATCTGATCGATAACACGCTTTTCCCTCACCAg                                                  <  2:722728/149‑1 (MQ=255)
                                                                                                cggcggCGAACTTGCAGGTTTGTTAGGGCCACAGGATGTGGAGATCTGCATTCCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCTGTGCGATCTGATCGATAACACGCTTTTCCCTCACCAGgatgat                                            <  2:404500/151‑1 (MQ=255)
                                                                                                  gcggcgAACTTGCAGGTTTGTTAGGGCCACAGGATGTGGAGATCTGCATTCCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCTGTGCGATCTGATCGATAACACGCTTTTCCCTCACCAGGATGATTa                                          <  1:687198/151‑1 (MQ=255)
                                                                                                  gcggcgAACTTGCAGGTTTGTTAGGGCCACAGGATGTGGAGATCTGCATTCCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCTGTGCGATCTGATCGATAACACGCTTTTCCCTCACCAGGATGATTa                                          <  2:263288/151‑1 (MQ=255)
                                                                                                     gcgAACTTGCAGGTTTGTTAGGGCCACAGGATGTGGAGATCTGCATTCCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCTGTGCGATCTGATCGATAACACGCTTTTCCCTCACCAGGATGATTAAgg                                       <  2:397673/151‑1 (MQ=255)
                                                                                                       gAACTTGCAGGTTTGTTAGGGCCACAGGATGTGGAGATCTGCATTCCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCTGTGCGATCTGATCGATAACACGCTTTTCCCTCACCAGGATGATTAAgg                                       >  2:406668/1‑149 (MQ=255)
                                                                                                        aaCTTGCAGGTTTGTTAGGGCCACAGGATGTGGAGATCTGCATTCCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCTGTGCGATCTGATCGATAACACGCTTTTCCCTCACCAGGATGATTAAGgaga                                    <  2:574999/151‑1 (MQ=255)
                                                                                                           ttGCAGGTTTGTTAGGGCCACAGGATGTGGAGATCTGCATTCCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCTGTGCGATCTGATCGATAACACGCTTTTCCCTCACCAGGATGATTAAGGAGAAt                                  >  2:517436/1‑150 (MQ=255)
                                                                                                           ttGCAGGTTTGTTAGGGCCACAGGATGTGGAGATCTGCATTCCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCTGTGCGATCTGATCGATAACACGCTTTTCCCTCACCAGGATGATTAAGGAGAAt                                  >  2:856632/1‑150 (MQ=255)
                                                                                                            tGCAGGTTTGTTAGGGCCACAGGATGTGGAGATCTGCATTCCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCTGTGCGATCTGATCGATAACACGCTTTTCCCTCACCAGGATGATTAAGGAGAATa                                 <  1:33320/150‑1 (MQ=255)
                                                                                                              cAGGTTTGTTAGGGCCACAGGATGTGGAGATCTGCATTCCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCTGTGCGATCTGATCGATAACACGCTTTTCCCTCACCAGGATGATTAAGGAGAATACAt                              >  2:106806/1‑151 (MQ=255)
                                                                                                               aGGTTTGTTAGGGCCACAGGATGTGGAGATCTGCATTCCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCTGTGCGATCTGATCGATAACACGCTTTTCCCTCACCAGGATGATTAAGGAGAATACATg                             <  2:462144/151‑1 (MQ=255)
                                                                                                                ggTTTGTTAGGGCCACAGGATGTGGAGATCTGCATTCCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCTGTGCGATCTGATCGATAACACGCTTTTCCCTCACCAGGATGATTAAGGAGAATACATg                             >  2:107512/1‑150 (MQ=255)
                                                                                                                    tGTTAGGGCCACAGGATGTGGAGATCTGCATTCCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCTGTGCGATCTGATCGATAACACGCTTTTCCCTCACCAGGATGATTAAGGAGAATACATGAAGGc                        >  2:824449/1‑151 (MQ=255)
                                                                                                                     gTTAGGGCCACAGGATGTGGAGATCTGCATTCCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCTGTGCGATCTGATCGATAACACGCTTTTCCCTCACCAGGATGATTAAGGAGAATACATGAAGGCa                       <  2:783503/151‑1 (MQ=255)
                                                                                                                        aGGGCCACAGGATGTGGAGATCTGCATTCCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCTGTGCGATCTGATCGATAACACGCTTTTCCCTCACCAGGATGATTAAGGAGAATACATGAAGGCAtt                     <  1:404971/150‑1 (MQ=255)
                                                                                                                                aGGATGTGGAGATCTGCATTCCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACTGTAAATTGCCTGTGCGATCTGATCGATAACCCGCTTTTCCCTCACCAGGATGATTAAGGAGAATACATGAAGGCATTATCGCCAAt            <  2:476030/151‑1 (MQ=255)
                                                                                                                                aGGATGTGGAGATCTGCATTCCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCTGTGCGATCTGATCGATAACACGCTTTTCCCTCACCAGGATGATTAAGGAGAATACATGAAGGCATTATCGCCAAt            >  1:897015/1‑151 (MQ=255)
                                                                                                                                  gATGTGGAGATCTGCATTCCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCTGTGCGATCTGATCGATAACACGCTTTTCCCTCACCAGGATGATTAAGGAGAATACATGAAGGCATTATCGCCAATCg          >  1:516318/1‑151 (MQ=255)
                                                                                                                                     gtgGAGATCTGCATTCCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCTGTGCGATCTGATCGATAACACGCTTTTCCCTCACCAGGATGATTAAGGAGAATACATGAAGGCATTATCGCCAATCGCa        <  2:43449/150‑1 (MQ=255)
                                                                                                                                        gagaTCTGCATTCCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCTGTGCGATCTGATCGATAACACGCTTTTCCCTCACCAGGATGATTAAGGAGAATACATGAAGGCATTATCGCCAATCGCAGTcc    <  1:66197/151‑1 (MQ=255)
                                                                                                                                        gagaTCTGCATTCCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCTGTGCGATCTGATCGATAACACGCTTTTCCCTCACCAGGATGATTAAGGAGAATACATGAAGGCATTATCGCCAATCGCAGTcc    <  2:641041/151‑1 (MQ=255)
                                                                                                                                         agaTCTGCATTCCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCTGTGCGATCTGATCGATAACACGCTTTTCCCTCACCAGGATGATTAAGGAGAATACATGAAGGCATTATCGCCAATCGCAGTCCt   <  2:490811/151‑1 (MQ=255)
                                                                                                                                          gaTCTGCATTCCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCTGTGCGATCTGATCGATAACACGCTTTTCCCTCACCAGGATGATTAAGGAGAATACATGAAGGCATTATCGCCAATCGCAGTCCtt  >  1:108431/1‑151 (MQ=255)
                                                                                                                                              |                                                                                                                                                  
AGCCATTTCCACCCGTGGCAACAGCCGCGATATTGTTAAAGCAGTTGAAGCCGCCGTTACGCGTGATATGACCATTGTGGCATTGACCGGCTATGACGGCGGCGAACTTGCAGGTTTGTTAGGGCCACAGGATGTGGAGATCCGCATTCCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCTGTGCGATCTGATCGATAACACGCTTTTCCCTCACCAGGATGATTAAGGAGAATACATGAAGGCATTATCGCCAATCGCAGTCCTT  >  NC_000913/3296150‑3296438

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