Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,152,439 C→T G12R (GGA→AGA)  yegK ← ser/thr phosphatase‑related protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,152,4390CT100.0% 244.7 / NA 61G12R (GGA→AGA) yegKser/thr phosphatase‑related protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (36/25);  total (36/25)

CATCGCGCCTTCGCCACCCTGCGAGACACTGCCTGCGCCATCAGCGACAAATACCGACAATAATGGCTGCTGGTCATTTAACCAGGCAATCTGCATCTGGCAGGCATCCTGACAAGGCAAGTCCGCGCTGATGTGCGACGTTCCAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCT  >  NC_000913/2152296‑2152586
                                                                                                                                               |                                                                                                                                                   
cATCGCGCCTTCGCCACCCTGCGAGACACTGCCTGCGCCATCAGCGACAAATACCGACAATAATGGCTGCTGGTCATTTAACCAGGCAATCTGCATCTGGCAGGCATCCTGACAAGGCAAGTCCGCGCTGATGTGCGACGTTCTAACGGtt                                                                                                                                              <  2:15785/151‑1 (MQ=255)
      gcCTTCGCCACCCTGCGAGACACTGCCTGCGCCATCAGCGACAAATACCGACAATAATGGCTGCTGGTCATTTAACCAGGCAATCTGCATCTGGCAGGCATCCTGACAAGGCAAGTCCGCGCTGATGTGCGACGTTCTAACGGTTGATGCa                                                                                                                                        <  2:137135/151‑1 (MQ=255)
                cccTGCGAGACACTGCCTGCGCCATCAGCGACAAATACCGACAATAATGGCTGCTGGTCATTTAACCAGGCAATCTGCATCTGGCAGGCATCCTGACAAGGCAAGTCCGCGCTGATGTGCGACGTTCTAACGGTTGATGCATAGACCAGa                                                                                                                               <  2:776026/150‑1 (MQ=255)
                   tGCGAGACACTGCCTGCGCCATCAGCGACAAATACCGACAATAATGGCTGCTGGTCATTTAACCAGGCAATCTGCATCTGGCAGGCATCCTGACAAGGCAAGTCCGCGCTGATGTGCGACGTTCTAACGGTTGATGCATAGACCAGACGcc                                                                                                                           >  2:435037/1‑151 (MQ=255)
                         acacTGCCTGCGCCATTAGCGACAAATACCGACAATAATGGCTGCTCGTCATTTAACCCGGCAATCTGCATCTGGCAGGCATCCTGACAAGGCAAGTCCGCGCTAATGTGCGGCGTTCTAACGGTTGATGCATAGAACAGACGCCAGc                                                                                                                        >  2:254579/1‑148 (MQ=255)
                                  gcgcCATCAGCGACAAATACCGACAATAATGGCTGCTGGTCATTTAACCAGGCAATCTGCATCTGGCAGGCATCCTGACAAGGCAAGTCCGCGCTGATGTGCGACGTTCTAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGAc                                                                                                              <  1:506675/149‑1 (MQ=255)
                                           gCGACAAATACCGACAATAATGGCTGCTGGTCATTTAACCAGGCAATCTGCATCTGGCAGGCATCCTGACAAGGCAAGTCCGCGCTGATGTGCGACGTTCTAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCttt                                                                                                   <  2:41270/151‑1 (MQ=255)
                                           gCGACAAATACCGACAATAATGGCTGCTGGTCATTTAACCAGGCAATCTGCATCTGGCAGGCATCCTGACAAGGCAAGTCCGCGCTGATGTGCGACGTTCTAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCttt                                                                                                   >  2:409792/1‑151 (MQ=255)
                                             gACAAATACCGACAATAATGGCTGCTGGTCATTTAACCAGGCAATCTGCATCTGGCAGGCATCCTGACAAGGCAAGTCCGCGCTGATGTGCGACGTTCTAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTg                                                                                                 >  1:798510/1‑151 (MQ=255)
                                              aCAAATACCGACAATAATGGCTGCTGGTCATTTAACCAGGCAATCTGCATCTGGCAGGCATCCTGACAAGGCAAGTCCGCGCTGATGTGCGACGTTCTAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTgg                                                                                                >  1:575994/1‑151 (MQ=255)
                                               cAAATACCGACAATAATGGCTGCTGGTCATTTAACCAGGCAATCTGCATCTGGCAGGCATCCTGACAAGGCAAGTCCGCGCTGATGTGCGACGTTCTAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGgc                                                                                               >  2:133771/1‑151 (MQ=255)
                                               cAAATACCGACAATAATGGCTGCTGGTCATTTAACCAGGCAATCTGCATCTGGCAGGCATCCTGACAAGGCAAGTCCGCGCTGATGTGCGACGTTCTAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGgc                                                                                               <  1:576798/151‑1 (MQ=255)
                                                aaaTACCGACAATAATGGCTGCTGGTCATTTAACCAGGCAATCTGCATCTGGCAGGCATCCTGACAAGGCAAGTCCGCGCTGATGTGCGACGTTCTAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGgcg                                                                                              >  2:180516/1‑151 (MQ=255)
                                                aaaTACCGACAATAATGGCTGCTGGTCATTTAACCAGGCAATCTGCATCTGGCAGGCATCCTGACAAGGCAAGTCCGCGCTGATGTGCGACGTTCTAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGgcg                                                                                              >  1:776724/1‑151 (MQ=255)
                                                 aaTACCGACAATAATGGCTGCTGGTCATTTAACCAGGCAATCTGCATCTGGCAGGCATCCTGACAAGGCAAGTCCGCGCTGATGTGCGACGTTCTAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGgcgc                                                                                             <  2:544017/151‑1 (MQ=255)
                                                     ccGACAATAATGGCTGCTGGTCATTTAACCAGGCAATCTGCATCTGGCAGGCATCCTGACAAGGCAAGTCCGCGCTGATGTGCGACGTTCTAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTcc                                                                                         >  1:87385/1‑151 (MQ=255)
                                                      cGACAATAATGGCTGCTGGTCATTTAACCAGGCAATCTGCATCTGGCAGGCATCCTGACAAGGCAAGTCCGCGCTGATGTGCGACGTTCTAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTcc                                                                                         >  1:380594/1‑150 (MQ=255)
                                                      cGACAATAATGGCTGCTGGTCATTTAACCAGGCAATCTGCATCTGGCAGGCATCCTGACAAGGCAAGTCCGCGCTGATGTGCGACGTTCTAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTcc                                                                                         >  2:818512/1‑150 (MQ=255)
                                                      cGACAATAATGGCTGCTGGTCATTTAACCAGGCAATCTGCATCTGGCAGGCATCCTGACAAGGCAAGTCCGCGCTGATGTGCGACGTTCTAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCa                                                                                        <  2:397227/151‑1 (MQ=255)
                                                        aCAATAATGGCTGCTGGTCATTTAACCAGGCAATCTGCATCTGGCAGGCATCCTGACAAGGCAAGTCCGCGCTGATGTGCGACGTTCTAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGc                                                                                      >  1:63509/1‑151 (MQ=255)
                                                           ataatGGCTGCTGGTCATTTAACCAGGCAATCTGCATCTGGCAGGCATCCTGACAAGGCAAGTCCGCGCTGATGTGCGACGTTCTAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGcaca                                                                                   <  1:895243/151‑1 (MQ=255)
                                                           ataatGGCTGCTGGTCATTTAACCAGGCAATCTGCATCTGGCAGGCATCCTGACAAGGCAAGTCCGCGCTGATGTGCGACGTTCTAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGcac                                                                                    >  1:462756/1‑150 (MQ=255)
                                                            taatGGCTGCTGGTCATTTAACCAGGCAATCTGCATCTGGCAGGCATCCTGACAAGGCAAGTCCGCGCTGATGTGCGACGTTCTAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGcaca                                                                                   >  2:837437/1‑150 (MQ=255)
                                                             aatGGCTGCTGGTCATTTAACCAGGCAATCTGCATCTGGCAGGCATCCTGACAAGGCAAGTCCGCGCTGATGTGCGACGTTCTAACGGTTGATGCATAGACCAGACGCCAGCTCACCCTGACGTCCAGCCTTTTGGCGCTTCCAGCACaa                                                                                  >  2:723447/1‑150 (MQ=255)
                                                             aatGGCTGCTGGTCATTTAACCAGGCAATCTGCATCTGGCAGGCATCCTGACAAGGCAAGTCCGCGCTGATGTGCGACGTTCTAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCGAGCACAAc                                                                                 <  2:78330/151‑1 (MQ=255)
                                                                ggCTGCTGGTCATTTAACCAGGCAATCTGCATCTGGCAGGCATCCTGACAAGGCAAGTCCGCGCTGATGTGCGACGTTCTAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTc                                                                              >  2:735466/1‑151 (MQ=255)
                                                                   tgctgGTCATTTAACCAGGCAATCTGCATCTGGCAGGCATCCTGACAAGGCAAGTCCGCGCTGATGTGCGACGTTCTAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGt                                                                           >  1:576395/1‑151 (MQ=255)
                                                                      tgGTCATTTAACCAGGCAATCTGCATCTGGCAGGCATCCTGACAAGGCAAGTCCGCGCTGATGTGCGACGTTCTAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTcc                                                                        >  2:728792/1‑151 (MQ=255)
                                                                                aCCAGGCAATCTGCATCTGGCAGGCATCCTGACAAGGCAAGTCCGCGCTGATGTGCGACGTTCTAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAg                                                              <  2:759068/151‑1 (MQ=255)
                                                                                   aGGCAATCTGCATCTGGCAGGCATCCTGACAAGGCAAGTCCGCGCTGATGTGCGACGTTCTAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCgg                                                           >  1:843550/1‑151 (MQ=255)
                                                                                   aGGCAATCTGCATCTGGCAGGCATCCTGACAAGGCAAGTCCGCGCTGATGTGCGACGTTCTAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCgg                                                           <  2:146412/151‑1 (MQ=255)
                                                                                     gCAATCTGCATCTGGCAGGCATCCTGACAAGGCAAGTCCGCGCTGATGTGCGACGTTCTAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCggg                                                          >  1:164534/1‑150 (MQ=255)
                                                                                      cAATCTGCATCTGGCAGGCATCCTGACAAGGCAAGTCCGCGCTGATGTGCGACGTTCTAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGa                                                         >  2:741192/1‑150 (MQ=255)
                                                                                       aaTCTGCATCTGGCAGGCATCCTGACAAGGCAAGTCCGCGCTGATGTGCGACGTTCTAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGaa                                                        >  2:153716/1‑150 (MQ=255)
                                                                                       aaTCTGCATCTGGCAGGCATCCTGACAAGGCAAGTCCGCGCTGATGTGCGACGTTCTAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGaa                                                        <  2:771912/150‑1 (MQ=255)
                                                                                       aaTCTGCATCTGGCAGGCATCCTGACAAGGCAAGTCCGCGCTGATGTGCGACGTTCTAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGaa                                                        >  2:840599/1‑150 (MQ=255)
                                                                                         tCTGCATCTGGCAGGCATCCTGACAAGGCAAGTCCGCGCTGATGTGCGACGTTCTAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAAcc                                                     >  2:656322/1‑151 (MQ=255)
                                                                                            gCATCTGGCAGGCATCCTGACAAGGCAAGTCCGCGCTGATGTGCGACGTTCTAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAAcc                                                     <  1:163082/148‑1 (MQ=255)
                                                                                                  ggcaggcaTCCTGACAAGGCAAGTCCGCGCTGATGTGCGACGTTCTAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGc                                            >  1:716899/1‑151 (MQ=255)
                                                                                                   gcaggcaTCCTGACAAGGCAAGTCCGCGCTGATGTGCGACGTTCTAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCg                                           <  1:212499/151‑1 (MQ=255)
                                                                                                   gcaggcaTCCTGACAAGGCAAGTCCGCGCTGATGTGCGACGTTCTAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCg                                           <  2:247326/151‑1 (MQ=255)
                                                                                                   gcaggcaTCCTGACAAGGCAAGTCCGCGCTGATGTGCGACGTTCTAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCg                                           >  1:348916/1‑151 (MQ=255)
                                                                                                    caggcaTCCTGACAAGGCAAGTCCGCGCTGATGTGCGACGTTCTAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGa                                          <  2:355467/151‑1 (MQ=255)
                                                                                                     aggcaTCCTGACAAGGCAAGTCCGCGCTGATGTGCGACGTTCTAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGa                                          >  2:628071/1‑150 (MQ=255)
                                                                                                      ggcaTCCTGACAAGGCAAGTCCGCGCTGATGTGCGACGTTCTAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAAc                                        >  2:773145/1‑151 (MQ=255)
                                                                                                        caTCCTGACAAGGCAAGTCCGCGCTGATGTGCGACGTTCTAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTg                                      >  2:75869/1‑151 (MQ=255)
                                                                                                           ccTGACAAGGCAAGTCCGCGCTGATGTGCGACGTTCTAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGa                                    <  1:839392/150‑1 (MQ=255)
                                                                                                                 aaGGCAAGTCCGCGCTGATGTGCGACGTTCTAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAg                             <  2:414707/151‑1 (MQ=255)
                                                                                                                 aaGGCAAGTCCGCGCTGATGTGCGACGTTCTAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAg                             <  1:229078/151‑1 (MQ=255)
                                                                                                                      aaGTCCGCGCTGATGTGCGACGTTCTAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTaa                         >  1:365468/1‑150 (MQ=255)
                                                                                                                          ccGCGCTGATGTGCGACGTTCTAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGt                    >  1:738536/1‑151 (MQ=255)
                                                                                                                            gcgcTGATGTGCGACGTTCTAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTc                  >  1:407845/1‑151 (MQ=255)
                                                                                                                              gcTGATGTGCGACGTTCTAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCAc                >  2:17684/1‑151 (MQ=255)
                                                                                                                                tGATGTGCGACGTTCTAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACgg              >  1:64620/1‑151 (MQ=255)
                                                                                                                                 gATGTGCGACGTTCTAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGa             <  2:284331/151‑1 (MQ=255)
                                                                                                                                  aTGTGCGACGTTCTAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGa             >  1:756851/1‑150 (MQ=255)
                                                                                                                                    gtgCGACGTTCTAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTACCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAAtt          <  1:466090/151‑1 (MQ=255)
                                                                                                                                       cGACGTTCTAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTa       <  1:95362/151‑1 (MQ=255)
                                                                                                                                        gACGTTCTAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTaa      <  2:731165/151‑1 (MQ=255)
                                                                                                                                          cGTTCTAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAAt     <  1:17706/150‑1 (MQ=255)
                                                                                                                                            ttCTAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCt  >  2:901733/1‑151 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
CATCGCGCCTTCGCCACCCTGCGAGACACTGCCTGCGCCATCAGCGACAAATACCGACAATAATGGCTGCTGGTCATTTAACCAGGCAATCTGCATCTGGCAGGCATCCTGACAAGGCAAGTCCGCGCTGATGTGCGACGTTCCAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCT  >  NC_000913/2152296‑2152586

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