Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,343,880 A→G V296A (GTG→GCG)  yciR ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,343,8800AG100.0% 294.0 / NA 94V296A (GTG→GCG) yciRhypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (42/52);  total (42/52)

GAGGTGTTGGATGCCAGTACCAGAAACTCATCCCCACCTGGACGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACACGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGCGTCGTTGACCTTTTTGAAATTATCCAAATCAAGATACACAACCCCAACT  >  NC_000913/1343790‑1343969
                                                                                          |                                                                                         
gAGGTGTTGGATGCCAGTACCAGAAACTCATCCTCACCTGGACGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGAcgcg                                                                                         <  2:1483999/93‑1 (MQ=255)
 aGGTGTTGGATGCCAGTACCAGAAACTCATCCTCACCTGGACGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGt                                                                                        >  1:4283511/1‑93 (MQ=255)
  ggTGTTGGATGCCAGTACCAGAAACTCATCCTCACCTGGACGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTc                                                                                       <  2:194597/93‑1 (MQ=255)
    tgtTGGATGCCAGTACCAGAAACTCATCCTCACCTGGACGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTcgc                                                                                     <  2:3259813/93‑1 (MQ=255)
     gtTGGATGCCAGTACCAGAAACTCATCCTCACCTGGACGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTcgcg                                                                                    >  2:4652116/1‑93 (MQ=255)
     gtTGGATGCCAGTACCAGAAACTCATCCTCACCTGGACGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTcgcg                                                                                    <  1:4077417/93‑1 (MQ=255)
      ttGGATGCCAGTACCAGAAACTCATCCTCACCTGGACGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGAcgcgtcgcgt                                                                                   <  2:1640585/93‑1 (MQ=255)
      gtagATCCCAGTACCAGATACTCATCCTCACCTGGACGCGCCAACGCCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGAcgcgtcgcgt                                                                                   <  2:1169160/90‑1 (MQ=255)
         gATGCCAGTACCAGAAACTCATCCTCACCTGGACGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGtaat                                                                                >  1:3643468/1‑93 (MQ=255)
          aTGCCAGTACCAGAAACTCATCCTCACCTGGACGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGtaata                                                                               >  1:3824169/1‑93 (MQ=255)
           tGCCAGTACCAGAAACTCATCCTCACCTGGACGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGtaataa                                                                              <  1:1557696/93‑1 (MQ=255)
            gCCAGTACCAGAAACTCATCCTCACCTGGACGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAAc                                                                             <  1:2038003/93‑1 (MQ=255)
            gCCAGTACCAGAAACTCATCCTCACCTGGACGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAAc                                                                             <  2:710744/93‑1 (MQ=255)
            gCCAGTACCAGAAACTCATCCTCACCTGGACGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAAc                                                                             <  1:36772/93‑1 (MQ=255)
             ccAGTACCAGAAACTCATCCTCACCTGGACGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACt                                                                            <  1:2183982/93‑1 (MQ=255)
             ccAGTACCAGAAACTCATCCTCACCTGGACGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACt                                                                            <  1:3614782/93‑1 (MQ=255)
             ccAGTACCAGAAACTCATCCTCACCTGGACGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACt                                                                            <  2:697344/93‑1 (MQ=255)
               aGTACCAGAAACTCATCCTCACCTGGACGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTgg                                                                          >  1:2711184/1‑93 (MQ=255)
                gTACCAGAACCTCATCCTCACCTGGACGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGt                                                                         <  1:2906447/93‑1 (MQ=255)
                gTACCAGAAACTCATCCTCACCTGGACGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGt                                                                         >  2:3857466/1‑93 (MQ=255)
                gTACCAGAAACTCATCCTCACCTGGACGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGt                                                                         >  1:876897/1‑93 (MQ=255)
                gTACCAGAAACTCATCCTCACCTGGACGCGCCAACACCTGGGCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGCCGCGTAATAACTGGt                                                                         >  2:3631014/1‑93 (MQ=255)
                 tACCAGAAACTCATCCTCACCTGGACGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTc                                                                        >  1:334758/1‑93 (MQ=255)
                      gAAACTCATCCTCACCTGGACGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACcaa                                                                   <  2:3814811/93‑1 (MQ=255)
                        aaCTCATCCTCACCTGGACGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACcaaac                                                                 >  2:1832632/1‑93 (MQ=255)
                         aCTCATCCTCACCTGGACGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACcaaaca                                                                >  1:3791716/1‑93 (MQ=255)
                           tCATCCTCACCTGGACGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACcaaacaaa                                                              <  2:1610304/93‑1 (MQ=255)
                           tCATCCTCACCTGGACGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACcaaacaaa                                                              <  2:47347/93‑1 (MQ=255)
                            cATCCTCACCTGGACGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAAt                                                             <  2:3057090/93‑1 (MQ=255)
                             aTCCTCACCTGGACGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATg                                                            <  2:536740/93‑1 (MQ=255)
                               cctcACCTGGACGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGcc                                                          <  1:3304381/90‑1 (MQ=255)
                                ctcACCTGGACGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGccc                                                         >  1:2310516/3‑93 (MQ=255)
                                   aCCTGGACGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATa                                                      <  1:3788011/93‑1 (MQ=255)
                                     cTGGACGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAgg                                                    <  2:4107521/93‑1 (MQ=255)
                                      tGGACGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGc                                                   >  2:148001/1‑93 (MQ=255)
                                      tGGACGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGc                                                   >  2:1886431/1‑93 (MQ=255)
                                       ggACGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGcg                                                  <  1:2217436/93‑1 (MQ=255)
                                       ggACGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGcg                                                  >  2:684962/1‑93 (MQ=255)
                                        gACGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGcgt                                                 >  2:3668091/1‑93 (MQ=255)
                                        gACGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGcgt                                                 >  2:3823685/1‑93 (MQ=255)
                                         aCGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGcgtc                                                <  2:2190006/93‑1 (MQ=255)
                                          cgcgCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGcgtcg                                               <  1:3071000/93‑1 (MQ=255)
                                          cgcgCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGcgtcg                                               >  2:3196977/1‑93 (MQ=255)
                                           gcgcCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGcgtcgt                                              <  2:1925168/93‑1 (MQ=255)
                                            cgcCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGCGTCGtt                                             <  2:3920362/93‑1 (MQ=255)
                                            cgcCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGCGTCGtt                                             <  1:86678/93‑1 (MQ=255)
                                            cgcCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGCGTCGtt                                             <  2:4181907/93‑1 (MQ=255)
                                            cgcCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGCGTCGtt                                             <  1:3427286/93‑1 (MQ=255)
                                              ccAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGCGTCGTTGa                                           <  2:2006102/93‑1 (MQ=255)
                                              ccAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGCGTCGTTGa                                           <  1:27778/93‑1 (MQ=255)
                                               cAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGCGTCGTTGAc                                          <  1:1283368/93‑1 (MQ=255)
                                                 acacCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGCGTCGTTGACCt                                        >  1:3571366/1‑93 (MQ=255)
                                                 acacCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGCGTCGTTGACCt                                        >  2:4080150/1‑93 (MQ=255)
                                                 acacCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGCGTCGTTGACCt                                        <  2:880561/93‑1 (MQ=255)
                                                   acCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGCGTCGTTGACCttt                                      >  2:3641896/1‑93 (MQ=255)
                                                    ccTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGCGTCGTTGACCtttt                                     <  1:2143919/93‑1 (MQ=255)
                                                    ccTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGCGTCGTTGACCtttt                                     <  1:2599749/93‑1 (MQ=255)
                                                      tGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGCGTCGTTGACCttttta                                   >  2:171502/1‑92 (MQ=255)
                                                      tGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGCGTCGTTGACCTTTTTg                                   >  1:4513169/1‑93 (MQ=255)
                                                      tGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGCGTCGTTGACCTTTTTg                                   >  1:3923421/1‑93 (MQ=255)
                                                       ggTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGCGTCGTTGACCTTTTTGa                                  <  2:1777122/93‑1 (MQ=255)
                                                       ggTCATGTTCGAGACAGCTTAAAATAGCCAAAGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGCGTCGTTGACCTTTTTGa                                  >  1:4053390/1‑93 (MQ=255)
                                                         tCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGCGTCGTTGACCTTTTTGaaa                                >  1:222123/1‑93 (MQ=255)
                                                         tCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGCGTCGTTGACCTTTTTGaaa                                >  1:75674/1‑93 (MQ=255)
                                                         tCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGCGTCGTTGACCTTTTTGaaa                                >  1:3720057/1‑93 (MQ=255)
                                                          cATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGCGTCGTTGACCTTTTTGAAAt                               >  1:2945574/1‑93 (MQ=255)
                                                           aTGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGCGTCGTTGACCTTTTTGAAAtt                              >  2:909467/1‑93 (MQ=255)
                                                            tGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGCGTCGTTGAGCTTTTTGAAATTa                             >  1:2679425/1‑93 (MQ=255)
                                                            tGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGCGTCGTTGACCTTTTTGAAATTa                             >  2:2782668/1‑93 (MQ=255)
                                                             gTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGCGTCGTTGACCTTTTTGAAATTAt                            <  1:456979/93‑1 (MQ=255)
                                                              ttCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGCGTCGTTGACCTTTTTGAAATTATc                           <  1:1933947/93‑1 (MQ=255)
                                                              ttCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGCGTCGTTGACCTTTTTGAAATTATc                           <  1:2110612/93‑1 (MQ=255)
                                                              ttCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGCGTCGTTGACCTTTTTGAAATTATc                           >  2:1285229/1‑93 (MQ=255)
                                                               tCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGCGTCGTTGACCTTTTTGAAATTATcc                          <  2:1410568/93‑1 (MQ=255)
                                                                cGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGCGTCGTTGACCTTTTTGAAATTATCCa                         <  2:2715675/93‑1 (MQ=255)
                                                                cGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGCGTCGTTGACCTTTTTGAAATTATCCa                         <  2:2917741/93‑1 (MQ=255)
                                                                cGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGCGTCGTTGACCTTTTTGAAATTATCCa                         <  2:3948821/93‑1 (MQ=255)
                                                                 gagaCAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGCGTCGTTGACCTTTTTGAAATTATCCaa                        >  2:2682679/1‑93 (MQ=255)
                                                                  agaCAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGCGTCGTTGACCTTTTTGAAATTATCCaaa                       <  1:1262611/93‑1 (MQ=255)
                                                                  agaCAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGCGTCGTTGACCTTTTTGAAATTATCCaaa                       <  1:4348255/93‑1 (MQ=255)
                                                                   gaCAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGCGTCGTTGACCTTTTTGAAATTATCCAAAt                      >  2:29266/1‑93 (MQ=255)
                                                                    aCAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGCGTCGTTGACCTTTTTGAAATTATCCAAATc                     <  1:2719308/93‑1 (MQ=255)
                                                                      aGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGCGTCGTTGACCTTTTTGAAATTATCCAAATCaa                   <  2:2170203/93‑1 (MQ=255)
                                                                        cTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGCGTCGTTGACCTTTTTGAAATTATCCAAATCAAGa                 <  1:3555258/93‑1 (MQ=255)
                                                                        cTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGCGTCGTTGACCTTTTTGAAATTATCCAAATCAAGa                 <  2:2017982/93‑1 (MQ=255)
                                                                         ttAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGCGTCGTTGACCTTTTTGAAATTATCCAAATCAAGAt                >  2:127719/1‑93 (MQ=255)
                                                                           aaaaTAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGCGTCGTTGACCTTTTTGAAATTATCCAAATCAAGATac              <  1:2343940/93‑1 (MQ=255)
                                                                            aaaTAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGCGTCGTTGACCTTTTTGAAATTATCCAAATCAAGATaca             >  1:3696262/1‑93 (MQ=255)
                                                                             aaTAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGCGTCGTTGACCTTTTTGAAATTATCCAAATCAAGATacac            <  2:3698493/93‑1 (MQ=255)
                                                                             aaTAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGCGTCGTTGACCTTTTTGAAATTATCCAAATCAAGATacac            >  1:177199/1‑93 (MQ=255)
                                                                               tAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGCGTCGTTGACCTTTTTGAAATTATCCAAATCAAGATACACaa          <  1:48702/93‑1 (MQ=255)
                                                                               tAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGCGTCGTTGACCTTTTTGAAATTATCCAAATCAAGATACACaa          >  2:23355/1‑93 (MQ=255)
                                                                                   cAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGCGTCGTTGACCTTTTTGAAATTATCCAAATCAAGATACACAAcccc      >  1:3808623/1‑93 (MQ=255)
                                                                                   cAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGCGTCGTTGACCTTTTTGAAATTATCCAAATCAAGATACACAAcccc      <  1:87233/93‑1 (MQ=255)
                                                                                   cAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGCGTCGTTGACCTTTTTGAAATTATCCAAATCAAGATACACAAcccc      >  1:619497/1‑93 (MQ=255)
                                                                                    aaTGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGCGTCGTTGACCTTTTTGAAATTATCCAAATCAAGATACACAACCCCa     >  1:3109563/1‑93 (MQ=255)
                                                                                       gACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGCGTCGTTGACCTTTTTGAAATTATCCAAATCAAGATACACAACCCCAACt  <  1:167211/93‑1 (MQ=255)
                                                                                          |                                                                                         
GAGGTGTTGGATGCCAGTACCAGAAACTCATCCCCACCTGGACGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACACGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGCGTCGTTGACCTTTTTGAAATTATCCAAATCAAGATACACAACCCCAACT  >  NC_000913/1343790‑1343969

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