Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,366,953 G→A G10R (GGA→AGA)  yhcF → putative transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,366,9530GA100.0% 238.6 / NA 75G10R (GGA→AGA) yhcFputative transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (42/33);  total (42/33)

GTAAAAGCGTTAGCCTCACACTGACAATGGCTCCAGCTTTTAAGGATGAAGGGGAACTAACCGACATGACAGATATTACAGGCAATCTGACGGTCCTGGTGGAGGCAAAATGAACAATGTAAAATTACTGATTGCCGGAAGTGCCTTTTTTGCCATGTCAGCGCAAGCCGCTGATAGAGTATCAATTGACGTTAAGGTGACTCTGGAAGCTGCAGCTTGTACTCCAATACTAAGCAACGGCGGAGTTGTTAACTTCGGTAGTCATTCAGTCA  >  NC_000913/3366817‑3367088
                                                                                                                                        |                                                                                                                                       
gTAAAAGCGTTAGCCTCACACTGACAATGGCTCCAGCTTTTAAGGATGAAGGGGAACTAACCGACATGACAGATATTACAGGCAATCTGACGGTCCTGGTGGAGGCAAAATGAACAATGTAAAATTACTGATTGCCAGa                                                                                                                                       >  2:756155/1‑139 (MQ=255)
    aaGCGTTAGCCTCACACTGACAATGGCTCCAGCTTTTAAGGATGAAGGGGAACTAACCGACATGACAGATATTACAGGCAATCTGACGGTCCTGGTGGAGGCAAAATGAACAATGTAAAATTACTGATTGCCAGAAGTg                                                                                                                                   <  1:1246439/139‑1 (MQ=255)
         ttAGCCTCACACTGACAATGGCTCCAGCTTTTAAGGATGAAGGGGAACTAACCGACATGACAGATATTACAGGCAATCTGACGGTCCTGGTGGAGGCAAAATGAACAATGTAAAATTACTGATTGCCAGAAGTGCCttt                                                                                                                              <  1:709171/139‑1 (MQ=255)
         ttAGCCTCACACTGACAATGGCTCCAGCTTTTAAGGATGAAGGGGAACTAACCGACATGACAGATATTACAGGCAATCTGACGGTCCTGGTGGAGGCAAAATGAACAATGTAAAATTACTGATTGCCAGAAGTGCCttt                                                                                                                              <  1:756155/139‑1 (MQ=255)
         ttAGCCTCACACTGACAATGGCTCCAGCTTTTAAGGATGAAGGGGAACTAACCGACATGACAGATATTACAGGCAATCTGACGGTCCTGGTGGAGGCAAAATGAACAATGTAAAATTACTGATTGCCAGAAGTGCCttt                                                                                                                              >  2:709171/1‑139 (MQ=255)
         ttAGCCTCACACTGACAATGGCTCCAGCTTTTAAGGATGAAGGGGAACTAACCGACATGACAGATATTACAGGCAATCTGACGGTCCTGGTGGAGGCAAAATGAACAATGTAAAATTACTGATTGCCAGAAGTGCCttt                                                                                                                              <  1:215326/139‑1 (MQ=255)
                           tGGCTCCAGCTTTTAAGGATGAAGGGGAACTAACCGACATGACAGATATTACAGGCAATCTGACGGTCCTGGTGGAGGCAAAATGAACAATGTAAAATTACTGATTGCCAGAAGTGc                                                                                                                                  >  1:621945/1‑117 (MQ=255)
                           tGGCTCCAGCTTTTAAGGATGAAGGGGAACTAACCGACATGACAGATATTACAGGCAATCTGACGGTCCTGGTGGAGGCAAAATGAACAATGTAAAATTACTGATTGCCAGAAGTGc                                                                                                                                  <  2:621945/117‑1 (MQ=255)
                           tGGCTCCAGCTTTTAAGGATGAAGGGGAACTAACCGACATGACAGATATTACAGGCAATCTGACGGTCCTGGTGGAGGCAAAATGAACAATGTAAAATTACTGATTGCCAGAAGTGCCTTTTTTGCCATGTCAGCGCaa                                                                                                            >  2:1301722/1‑139 (MQ=255)
                             gCTCCAGCTTTTAAGGATGAAGGGGAACTAACCGACATGACAGATATTACAGGCAATCTGACGGTCCTGGTGGAGGCAAAATGAACAATGTAAAATTACTGATTGCCAGAAGTGCCTTTTTTGCCATGTCAGCGCAAGc                                                                                                          <  1:768215/139‑1 (MQ=255)
                             gCTCCAGCTTTTAAGGATGAAGGGGAACTAACCGACATGACAGATATTACAGGCAATCTGACGGTCCTGGTGGAGGCAAAATGAACAATGTAAAATTACTGATTGCCAGAAGTGCCTTTTTTGCCATGTCAGCGCAAGc                                                                                                          <  1:370189/139‑1 (MQ=255)
                             gCTCCAGCTTTTAAGGATGAAGGGGAACTAACCGACATGACAGATATTACAGGCAATCTGACGGTCCTGGTGGAGGCAAAATGAACAATGTAAAATTACTGATTGCCAGAAGTGCCTTTTTTGCCATGTCAGCGCAAGc                                                                                                          <  1:610338/139‑1 (MQ=255)
                             gCTCCAGCTTTTAAGGATGAAGGGGAACTAACCGACATGACAGATATTACAGGCAATCTGACGGTCCTGGTGGAGGCAAAATGAACAATGTAAAATTACTGATTGCCAGAAGTGCCTTTTTTGCCATGTCAGCGCAAGc                                                                                                          >  1:375581/1‑139 (MQ=255)
                             ccTCCGGCCTTTAAGGTTGAGGGGGAACTAACCGACGTAACATATATTACAAGAAATCTGACGGTCCTGATGGACGCAAAATGAACAATTTAAAATTACTGATTGCCAGAAGTGCCTTTTTTGCCAGGTCAGCGCAAgt                                                                                                          >  2:480089/2‑138 (MQ=255)
                                                      aaCTAACCGACATGACAGATATTACAGGCAATCTGACGGTCCTGGTGGAGGCAAAATGAACAATGTAAAATTACTGATTGCCAGAAGTGCCttt                                                                                                                              <  1:782398/94‑1 (MQ=255)
                                                      aaCTAACCGACATGACAGATATTACAGGCAATCTGACGGTCCTGGTGGAGGCAAAATGAACAATGTAAAATTACTGATTGCCAGAAGTGCCttt                                                                                                                              >  2:782398/1‑94 (MQ=255)
                                                          aaCCGACATGACAGATATTACAGGCAATCTGACGGTCCTGGTGGAGGCAAAATGAACAATGTAAAATTTCTGATTGCCAGAAGTGCCTTTTTTGCCATGTCAGCGCAAGCCGCTGATAGAGTATCAATTGACGTTAAgg                                                                             >  1:962945/1‑139 (MQ=255)
                                                          aaCCGACATGACAGATATTACAGGCAATCTGACGGTCCTGGTGGAGGCAAAATGAACAATGTAAAATTACTGATTGCCAGAAGTGCCTTTTTTGCCATGTCAgcgc                                                                                                              >  2:484462/1‑106 (MQ=255)
                                                          aaCCGACATGACAGATATTACAGGCAATCTGACGGTCCTGGTGGAGGCAAAATGAACAATGTAAAATTACTGATTGCCAGAAGTGCCTTTTTTGCCATGTCAgcgc                                                                                                              <  1:484462/106‑1 (MQ=255)
                                                          aaCCGACATGACAGATATTACAGGCAATCTGACGGTCCTGGTGGAGGCAAAATGAACAATGTAAAATTACTGATTGCCAGAAGTGCCTTTTTTGCCATGTCAgcgc                                                                                                              <  1:428745/106‑1 (MQ=255)
                                                          aaCCGACATGACAGATATTACAGGCAATCTGACGGTCCTGGTGGAGGCAAAATGAACAATGTAAAATTACTGATTGCCAGAAGTGCCTTTTTTGCCATGTCAgcgc                                                                                                              >  2:428745/1‑106 (MQ=255)
                                                                         tatTACAGGCAATCTGACGGTCCTGGTGGAGGCAAAATGAACAATGTAAAATTACTGATTGCCAGAAGTGCCTTTTTTGCCATGTCAGCGCAAGCCGCTGATAGAGTATCAATTGACGTTAAGGTGACTCTGGAAGCTg                                                              <  1:371624/139‑1 (MQ=255)
                                                                            tACAGGCAATCTGACGGTCCTGGTGGAGGCAAAATGAACAATGTAAAATTACTGATTGCCAGAAGTGCCTTTTTTGc                                                                                                                         >  1:717991/1‑77 (MQ=255)
                                                                            tACAGGCAATCTGACGGTCCTGGTGGAGGCAAAATGAACAATGTAAAATTACTGATTGCCAGAAGTGCCTTTTTTGc                                                                                                                         <  2:717991/77‑1 (MQ=255)
                                                                             aCAGGCAATCTGACGGTCCTGGTGGAGGCAAAATGAACAATGTAAAATTACTGATTGCCAGAAGTGCCTTTTTTGCCATGTCAGCGCAAGCCGCTGATa                                                                                                  <  2:1221642/99‑1 (MQ=255)
                                                                             aCAGGCAATCTGACGGTCCTGGTGGAGGCAAAATGAACAATGTAAAATTACTGATTGCCAGAAGTGCCTTTTTTGCCATGTCAGCGCAAGCCGCTGATa                                                                                                  >  1:1221642/1‑99 (MQ=255)
                                                                             aCAGGCAATCTGACGGTCCTGGTGGAGGCAAAATGAACAATGTAAAATTACTGATTGCCAGAAGTGCCTTTTTTGCCATGTCAGCGCAAGCCGCTGATAGAGTATCAATTGACGTTAAGGTGACTCTGGAAGCTGCAGc                                                          >  1:860711/1‑139 (MQ=255)
                                                                                ggCAATCTGACGGTTCTGGTGGAGGCAAAATGAACAATGTAAAATTACTGATTGCCAGAAGTGCCTTTTTTGCCATTTCAGCGCAAGCCTCTGATAGAGTATCAATTGACGTTAAGGTGACTCTGTAAGCTGCAGCTTg                                                       <  2:794011/139‑1 (MQ=255)
                                                                                ggCAATCTGACGGTCCTGGTGGAGGCAAAATGAACAATGTAAAATTACTGATTGCCAGAAGTGCCTTTTTTGCCATGTCAGCGCAAGCCGCTGATAGAGTATCAATTGACGTTAAGGTGACTCTGGAAGCTGCAGCTTg                                                       >  2:865558/1‑139 (MQ=255)
                                                                                ggCAATCTGACGGTCCTGGTGGAGGCAAAATGAACAATGTAAAATTACTGATTGCCAGAAGTGCCTTTTTTGCCATGTCAGCGCAAGCCGCTGATAGAGTATCAATTGACGTTAAGGTGACTCTGGAAGCTGCAGCTTg                                                       <  1:755728/139‑1 (MQ=255)
                                                                                ggCAATCTGACGGTCCTGGTGGAGGCAAAATGAACAATGTAAAATTACTGATTGCCAGAAGTGCCTTTTTTGCCATGTCAGCGCAAGCCGCTGATAGAGTATCAATTGACGTTAAGGTGACTCTGGAAGCTGCAGCTTg                                                       <  1:402395/139‑1 (MQ=255)
                                                                                 gCAATCTGACGGTCCTGGTGGAGGCAAAATGAACAATGTAAAATTACTGATTGCCAGAAGTGCCTTTTTTGCCATGTCAGCGCAAGCCGCTGATAGAGTATCAATTGACGTTAAGGTGACTCTGGAAGCTGCAGCTTGt                                                      >  1:913709/1‑139 (MQ=255)
                                                                                  cAATCTGACGGTCCTGGTGGAGGCAAAATGAACAATGTAAAATTACTGATTGCCAGAAGTGCCTTTTTTGCCATGTCAGCGCAAGCCGCTGATAGAGTa                                                                                             >  2:1229949/1‑99 (MQ=255)
                                                                                  cAATCTGACGGTCCTGGTGGAGGCAAAATGAACAATGTAAAATTACTGATTGCCAGAAGTGCCTTTTTTGCCATGTCAGCGCAAGCCGCTGATAGAGTa                                                                                             <  1:1229949/99‑1 (MQ=255)
                                                                                  cAATCTGACGGTCCTGGTGGAGGCAAAATGAACAATGTAAAATTACTGATTGCCAGAAGTGCCTTTTTTGCCATGTCAGCGCAAGCCGCTGATAGAGTATCAATTGACGTTAAGGTGACTCTGGAAGCTGCAGCTTGTa                                                     <  2:706670/139‑1 (MQ=255)
                                                                                   aaTCTGACGGTCCTGGTGGAGGCAAAATGAACAATGTAAAATTACTGATTGCCAGAAGTGCCTTTTTTGCCATGTCAGCGCAAGCCGCTGATAGAGTATCAAt                                                                                        <  2:1051014/103‑1 (MQ=255)
                                                                                   aaTCTGACGGTCCTGGTGGAGGCAAAATGAACAATGTAAAATTACTGATTGCCAGAAGTGCCTTTTTTGCCATGTCAGCGCAAGCCGCTGATAGAGTATCAAt                                                                                        >  1:1051014/1‑103 (MQ=255)
                                                                                    aTCTGACGGTCCTGGTGGAGGCAAAATGAACAATGTAAAATTACTGATTGCCAGAAGTGCCTTTTTTGCCATGTCAGCGCAAGc                                                                                                          <  2:1295298/84‑1 (MQ=255)
                                                                                    aTCTGACGGTCCTGGTGGAGGCAAAATGAACAATGTAAAATTACTGATTGCCAGAAGTGCCTTTTTTGCCATGTCAGCGCAAGc                                                                                                          >  1:1295298/1‑84 (MQ=255)
                                                                                    aTCTGACGGTCCTGGTGGAGGCAAAATGAACAATGTAAAATTACTGATTGCCAGAAGTGCCTTTTTTGCCATGTCAGCGCAAGCCGCTGATAGAGTATCAATTGACGTTAAGGTGACTCTGGAAGCTGc                                                             <  1:942081/129‑1 (MQ=255)
                                                                                    aTCTGACGGTCCTGGTGGAGGCAAAATGAACAATGTAAAATTACTGATTGCCAGAAGTGCCTTTTTTGCCATGTCAACGCAAGCCGCTGATAGAGTATCAATTGACGTTAAGGTGACTCTGGAAGCTGc                                                             >  2:942081/1‑129 (MQ=255)
                                                                                     tCTGACGGTCCTGGTGGAGGCAAAATGAACAATGTAAAATTACTGATTGCCAGAAGTGCCTTTTTTGCCATGTCAGCGCAAGCCGc                                                                                                       <  1:713583/86‑1 (MQ=255)
                                                                                     tCTGACGGTCCTGGTGGAGGCAAAATGAACAATGTAAAATTACTGATTGCCAGAAGTGCCTTTTTTGCCATGTCAGCGCAAGCCGc                                                                                                       >  2:713583/1‑86 (MQ=255)
                                                                                        gACGGTCCTGGTGGAGGCAAAATGAACAATGTAAAATTACTGATTGCCAGAAGTGCCTTTTTTGCCATGTCAGCGCAAGCCGCTGATAGAGTATCAATTGACGTTAAGGTGACTCTGGAAGCTGCAGCTTGTACTCCaa                                               <  2:860711/139‑1 (MQ=255)
                                                                                          cGGTCCTGGTGGAGGCAAAATGAACAATGTAAAATTTCTGATTGCCAGAAGTGCCTTTTTTGCCATGTCAGCGCAAGCCGCTGATAGAGTATCAATTGACGTTAAGGTGACTCTGGAAGCTGCAGCTTGTACTCCAATa                                             <  2:962945/139‑1 (MQ=255)
                                                                                           ggTCCTGGTGGAGGCAAAATGAACAATGTAAAATTACTGATTGCCAGAAGTGCCTTTTTTGCCATGTCAGCGCAAGCCGCTGATAGAGTATCAATTGACGTTAAGGTGACTCTGGAAGCTGCAGCTTGTACTCCAATAc                                            >  2:1205489/1‑139 (MQ=255)
                                                                                           ggTCCTGGTGGAGGCAAAATGAACAATGTAAAATTACTGATTGCCAGAAGTGCCTTTTTTGCCATGTCAGCGCAAGCCGCTGATAGAGTATCAATTGACGTTAAGGTGACTCTGGAAGCTGCAGCTTGTACTCCAATAc                                            >  1:1172308/1‑139 (MQ=255)
                                                                                             tCCTGGTGGAGGCAAAATGAACAATGTAAAATTACTGATTGCCAGAAGTGCCTTTTTTGCCATGTCAGCGCAAGCCGCTGATAGAGTATCAATTGACGTTAAGGTGACTCTGGAAGCTGCAGCTTGTACTCCAATACTa                                          >  2:1272238/1‑139 (MQ=255)
                                                                                               cTGGTGGAGGCAAAATGAACAATGTAAAATTACTGATTGCCAGAAGTGCCTTTTTTGCCATGTCAGCGCAAGCCGc                                                                                                       <  2:485381/76‑1 (MQ=255)
                                                                                               cTGGTGGAGGCAAAATGAACAATGTAAAATTACTGATTGCCAGAAGTGCCTTTTTTGCCATGTCAGCGCAAGCCGc                                                                                                       >  1:485381/1‑76 (MQ=255)
                                                                                                       ggCAAAATGAACAATGTAAAATTACTGATTGCCAGAAGTGCCttt                                                                                                                              >  2:392306/1‑45 (MQ=255)
                                                                                                       ggCAAAATGAACAATGTAAAATTACTGATTGCCAGAAGTGCCttt                                                                                                                              <  1:392306/45‑1 (MQ=255)
                                                                                                       ggCAAAATGAACAATGTAAAATTACTGATTGCCAGAAGTGCCTTTTTTGCCATGTCAGCGCAAGCCGCTGATag                                                                                                 >  2:506620/1‑74 (MQ=255)
                                                                                                       ggCAAAATGAACAATGTAAAATTACTGATTGCCAGAAGTGCCTTTTTTGCCATGTCAGCGCAAGCCGCTGATag                                                                                                 <  1:506620/74‑1 (MQ=255)
                                                                                                       ggCAAAATGAACAATGTAAAATTACTGATTGCCAGAAGTGCCTTTTTTGCCATGTCAGCGCAAGCCGCTGATAGAGTATCAATTGACGTTAAGGTGACTCTGGAAGCTGCAGCTTGTACTCCAATACTAAGCAAcggcg                                >  1:217772/1‑139 (MQ=255)
                                                                                                       ggCAAAATGAACAATGTAAAATTACTGATTGCCAGAAGTGCCTTTTTTGCCATGTCAGCGCAAGCCGCTGATAGAGTATCAATTGACGTTAAGGTGACTCTGGAAGCTGCAGCTTGTACTCCAATACGAAGCAAcggcg                                >  1:854007/1‑139 (MQ=255)
                                                                                                            aaTGAACAATGTAAAATTACTGATTGCCAGAAGTGCCTTTTTTGCCATGTCAGCGCAAGCCGCTGATAGAGTATCAATTGACGTTAAGGTGACTCTGGAAGCTGCAGCTTGTACTCCAATACTAAGCAACGGCGGAgtt                           >  2:561138/1‑139 (MQ=255)
                                                                                                                     tGTAAAATTACTGATTGCCAGAAGTGCCTTTTTTGCCATGTCAGCGCAAGCCGCTGATAGAGTATCAATTGACGTTAAGGTGACTCTGGAAGCTGCAGCTTGTACTCCAATACTAAGCAACGGCGGAGTTGTTAACTTc                  <  1:1272238/139‑1 (MQ=255)
                                                                                                                       tAAAATTACTGATTGCCAGAAGTGCCTTTTTTGCCATGTCAGCGCAAGCCGCTGATAGAGTATCAATTGACGTTAAGGTGACTCTGGAAGCTGCAGCTTGTACTCCAATACTAAGCAACGGCGGAGTTGTTAACTTCgg                >  1:395547/1‑139 (MQ=255)
                                                                                                                       tAAAATTACTGATTGCCAGAAGTGCCTTTTTTGCCATGTCAGCGCAAGCCGCTGATAGAGTATCAATTGACGTTAAGGTGACTCTGGAAGCTGCAGCTTGTACTCCAATACTAAGCAACGGCGGAGTTGTTAACTTCgg                >  1:402615/1‑139 (MQ=255)
                                                                                                                          aaTTACTGATTGCCAGAAGTGCCTTTTTTGCCATGTCAGCGCAAGCCGCTGATAGAGTATCAATTGACGTTAAGGTGACTCTGGAAGCTGCAGCTTGTACTCCAATACTAAGCAACGGCGGAGTTGTTAACTTCGGTAg             >  1:1102323/1‑139 (MQ=255)
                                                                                                                          aaTTACTGATTGCCAGAAGTGCCTTTTTTGCCATGTCAGCGCAAGCCGCTGATAGAGTATCAATTGACGTTAAGGTGACTCTGGAAGCTGCAGCTTGTACTCCAATACTAAGCAACGGCGGAGTTGTTAACTTCGGTAg             >  2:407161/1‑139 (MQ=255)
                                                                                                                               cTGATTGCCAGAAGTGCCTTTTTTGCCATGTCAGCGCAAGCCGCTGATagag                                                                                               >  2:415423/1‑52 (MQ=255)
                                                                                                                               cTGATTGCCAGAAGTGCCTTTTTTGCCATGTCAGCGCAAGCCGCTGATagag                                                                                               <  1:415423/52‑1 (MQ=255)
                                                                                                                               cTGATTGCCAGAAGTGCCTTTTTTGCCATGTCAGCGCAAGCCGCTGATAGAGTATCAATTGACGTTAAGGTGACTCTGGAAGCTGCAGCTTGTACTCCAATACTAAGCAACGGCGGAGTTGTTAACTTCGGTAGTCAtt        >  2:350955/1‑139 (MQ=255)
                                                                                                                               cTGATTGCCAGAAGTGCCTTTTTTGCCATGTCAGCGCAAGCCGCTGATAGAGTATCAATTGACGTTAAGGTGACTCTGGAAGCTGCAGCTTGTACTCCAATACTAAGCAACGGCGGAGTTGTTAACTTCGGTAGTCAtt        >  2:1265419/1‑139 (MQ=255)
                                                                                                                               cTGATTGCCAGAAGTGCCTTTTTTGCCATGTCAGCGCAAGCCGCTGATAGAGTATCAATTGACGTTAAGGTGACTCTGGAAGCTGCAGCTTGTACTCCAATACTAAGCAACGGCGGAGTTGTTAACTTCGGTAGTCAtt        >  1:596317/1‑139 (MQ=255)
                                                                                                                               cTGATTGCCAGAAGGGCCTTTTTTGCCATGTCAGCGCAAGCCGCTGATAGAGTATCAATTGACGTTAAGGTGACTCTGGAAGCTGCAGCTTGTACTCCAATACTAAGCAACGGCGGAGTTGTTAACTTCGGTAGTCAtt        >  2:984780/1‑139 (MQ=255)
                                                                                                                                 gATTGCCAGAAGTGCCTTTTTTGCCATGTCAGCGCAAGCCGCTGATAGAGTATCAATTGACGTTAAGGTGACTCTGGAAGCTGCAGCTTGTACTCCAATACTAAGCAACGGCGGAGTTGTTAACTTCGGTAGTCATTCa      <  1:561138/139‑1 (MQ=255)
                                                                                                                                    tGCCAGAAGTGCCTTTTTTGCCATGTCAGCGCAAGCCGCTGATAGAGTATCAATTGACGTTAAGGTGACTCTGGAAGCTGCAGCTTg                                                       <  2:657940/87‑1 (MQ=255)
                                                                                                                                    tGCCAGAAGTGCCTTTTTTGCCATGTCAGCGCAAGCCGCTGATAGAGTATCAATTGACGTTAAGGTGACTCTGGAAGCTGCAGCTTg                                                       >  1:657940/1‑87 (MQ=255)
                                                                                                                                    tGCCAGAAGTGCCTTTTTTGCCATGTCAGCGCAAGCCGCTGATAGAGTATCAATTGACGTTAAGGTGACTCTGGAAGCTGCAGCTTGTa                                                     >  2:420744/1‑89 (MQ=255)
                                                                                                                                    tGCCAGAAGTGCCTTTTTTGCCATGTCAGCGCAAGCCGCTGATAGAGTATCAATTGACGTTAAGGTGACTCTGGAAGCTGCAGCTTGTa                                                     <  1:420744/89‑1 (MQ=255)
                                                                                                                                    tGCCAGAAGTGCCTTTTTTGCCATGTCAGCGCAAGCCGCTGATAGAGTATCAATTGACGTTAAGGTGACTCTGGAAGCTGCAGCTTGTACTCCAATACTAAGCAACGGCGGAGTTGTTAACTTCGGTAGTCATTCAGTc   >  1:613865/1‑139 (MQ=255)
                                                                                                                                     gCCAGAAGTGCCTTTTTTGCCATGTCAGCGCAAGCCGCTGATAGAGTATCAATTGACGTTAAGGTGACTCTGGAAGCTGCAGCTTGTACTCCAATACTAAGCAACGGCGGAGTTGTTAACTTCGGTAGTCATTCAGTCa  >  2:614643/1‑139 (MQ=255)
                                                                                                                                        |                                                                                                                                       
GTAAAAGCGTTAGCCTCACACTGACAATGGCTCCAGCTTTTAAGGATGAAGGGGAACTAACCGACATGACAGATATTACAGGCAATCTGACGGTCCTGGTGGAGGCAAAATGAACAATGTAAAATTACTGATTGCCGGAAGTGCCTTTTTTGCCATGTCAGCGCAAGCCGCTGATAGAGTATCAATTGACGTTAAGGTGACTCTGGAAGCTGCAGCTTGTACTCCAATACTAAGCAACGGCGGAGTTGTTAACTTCGGTAGTCATTCAGTCA  >  NC_000913/3366817‑3367088

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