Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,451,123	C→A	V22V (GTC→GTA)	ytfR →	putative sugar ABC transporter ATPase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,451,123	0	C	A	100.0%	264.2 / NA	67	V22V (GTC→GTA)	ytfR	putative sugar ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (32/35);  total (32/35)

GCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTCAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGCGGGAAAATCAACGCTAATCAAAGCATTAACTGGTGTATACCACGCCGATCGCGGCACCATCTGGCTGGAAGGCC  >  NC_000913/4450976‑4451268                                                                                                                                                    |                                                                                                                                                  gCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTaaaa                                                                                                                                                <  1:527707/151‑1 (MQ=255) gCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTaaa                                                                                                                                                 >  2:762744/1‑150 (MQ=255)  cGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTaaa                                                                                                                                                 >  1:756105/1‑149 (MQ=255)  cGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTaaa                                                                                                                                                 >  1:512010/1‑149 (MQ=255)    aTGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCg                                                                                                                                             >  1:905668/1‑151 (MQ=255)      ggTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCGtt                                                                                                                                           <  1:146257/151‑1 (MQ=255)       gTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCGtt                                                                                                                                           >  2:83889/1‑150 (MQ=255)         gCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCGtt                                                                                                                                           >  2:85300/1‑148 (MQ=255)         gCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAg                                                                                                                                         <  2:331241/150‑1 (MQ=255)         gCACCGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAg                                                                                                                                         <  1:543580/150‑1 (MQ=255)           aCTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAGACa                                                                                                                                      <  1:568024/151‑1 (MQ=255)                aGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAGACAACGtt                                                                                                                                 <  2:706618/151‑1 (MQ=255)                aGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAGACAACGtt                                                                                                                                 <  2:706181/151‑1 (MQ=255)                   cggATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAACCCATGACGACCTACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTGCCCCGGCGTAAAAGCGTTAGACAACGTTGa                                                                                                                               >  1:581201/3‑150 (MQ=255)                                tACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAGACAACGTTGATTTCAGCCTgcgc                                                                                                                  <  2:209403/150‑1 (MQ=255)                                tACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAGACAACGTTGATTTCAGCCTGCGcc                                                                                                                 >  1:690929/1‑151 (MQ=255)                                 aCCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCg                                                                                                                <  1:838459/151‑1 (MQ=255)                                  ccGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGt                                                                                                               <  1:626643/151‑1 (MQ=255)                                    gcgTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAGACAACGTTGATTTCAGCCTgcgc                                                                                                                  <  1:240775/146‑1 (MQ=255)                                    gcgTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTgg                                                                                                             >  2:902279/1‑151 (MQ=255)                                          tATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGGCCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGaaa                                                                                                        <  2:346108/150‑1 (MQ=255)                                          tATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAAt                                                                                                       >  1:521767/1‑151 (MQ=255)                                           aTCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTCTCCCCGGCGTAAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATc                                                                                                      <  1:476490/151‑1 (MQ=255)                                                gCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAGACAACGTGGATTTCAGCCGGCGCCGTGGCGAAATCATGgc                                                                                                 <  1:577719/151‑1 (MQ=255)                                                 ccTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCAttg                                                                                                  >  2:527270/1‑147 (MQ=255)                                                     cTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCgctgc                                                                                            >  2:960751/1‑151 (MQ=255)                                                        ggAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCg                                                                                         <  1:483289/151‑1 (MQ=255)                                                          aGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGt                                                                                       >  2:607325/1‑151 (MQ=255)                                                           gCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTg                                                                                      <  1:713517/151‑1 (MQ=255)                                                           gCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTg                                                                                      >  2:221840/1‑151 (MQ=255)                                                             aCGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGa                                                                                     >  2:690915/1‑150 (MQ=255)                                                                  ttACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACgg                                                                               >  1:557028/1‑151 (MQ=255)                                                                    ccacCGGGGGAAACAATGACGACCGACCAACACCAGGAGATCCTCCGCCCCGAAGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACgggg                                                                             >  1:529995/2‑151 (MQ=255)                                                                    acacCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACggg                                                                              <  2:867243/150‑1 (MQ=255)                                                                     cacCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGc                                                                            <  1:113622/151‑1 (MQ=255)                                                                      acCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGCg                                                                           <  2:450698/151‑1 (MQ=255)                                                                      acCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGCg                                                                           >  2:402718/1‑151 (MQ=255)                                                                        cGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGCggg                                                                         >  1:172291/1‑151 (MQ=255)                                                                        cGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGCggg                                                                         >  1:347298/1‑151 (MQ=255)                                                                         gggggAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGCGGGa                                                                        <  2:168701/151‑1 (MQ=255)                                                                          ggggAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGCggg                                                                         >  1:5083/1‑149 (MQ=255)                                                                          ggggAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGCGGGaa                                                                       <  2:539182/151‑1 (MQ=255)                                                                           gggAAACCATGACGACCGACCAACGCCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGCGGGa                                                                        <  1:42159/149‑1 (MQ=255)                                                                                  cATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGCGGGAAAATCAACg                                                               >  2:470456/1‑151 (MQ=255)                                                                                    tGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGCGGGAAAATCAACGc                                                              >  1:874689/1‑150 (MQ=255)                                                                                               cAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGCGGGAAAATCAACGCTAATCAAAGCAt                                                  <  1:987100/151‑1 (MQ=255)                                                                                               cAACACCAGGAGATCCTCCACACCGAAGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGCGGGAAAATCAACGCTAATCAAAGCAt                                                  >  2:141014/1‑151 (MQ=255)                                                                                                   acCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGCGGGAAAATCAACGCTAATCAAAGCATTaa                                               >  1:869526/1‑150 (MQ=255)                                                                                                      aGGAGATCCCCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGCGGGAAAATCAACGCTAATCAAAGCATTAACTgg                                           <  2:809119/151‑1 (MQ=255)                                                                                                       ggAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGCGGGAAAATCAACGCTAATCAAAGCATTAACTGgt                                          <  2:541875/151‑1 (MQ=255)                                                                                                        gagaTCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGCGGGAAAATCAACGCTAATCAAAGCATTAACTGgtg                                         >  2:572426/1‑151 (MQ=255)                                                                                                        gagaTCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGCGGGAAAATCAACGCTAATCAAAGCATTAACTGgtg                                         <  1:314395/151‑1 (MQ=255)                                                                                                        gagaTCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGCGGGAAAATCAACGCTAATCAAAGCATTAACTGgtg                                         <  2:683137/151‑1 (MQ=255)                                                                                                          gaTCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGCGGGAAAATCAACGCTAATCAAAGCATTAACTGgtgt                                        >  2:31558/1‑150 (MQ=255)                                                                                                            tcctccGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGCGGGAAAATCAACGCTAATCAAAGCATTAACTGGTGtat                                      >  1:751554/1‑150 (MQ=255)                                                                                                               tccGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGCGGGAAAATCAACGCTAATCAAAGCATTAACTGGTGTATAcc                                   >  1:346080/1‑150 (MQ=255)                                                                                                                ccGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGCGGGAAAATCAACGCTAATCAAAGCATTAACTGGTGGATACCa                                  <  1:765041/150‑1 (MQ=255)                                                                                                                 cGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGCGGGAAAATCAACGCTAATCAAAGCATTAACTGGTGTATACCACg                                <  1:85263/151‑1 (MQ=255)                                                                                                                      cGAAGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGCGGGAAAATCAACGCTAATCAAAGCATTAACTGGTGTATACCACGCCGAt                           <  2:48132/151‑1 (MQ=255)                                                                                                                        aaGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGCGGGAAAATCAACGCTAATCAAAGCATTAACTGGTGTATACCACGCCGATcg                         >  2:470964/1‑151 (MQ=255)                                                                                                                        aaGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGCGGGAAAATCAACGCTAATCAAAGCATTAACTGGTGTATACCACGCCGATcg                         <  1:855245/151‑1 (MQ=255)                                                                                                                        aaGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGCGGGAAAATCAACGCTAATCAAAGCATTAACTGGTGTATACCACGCCGATcg                         <  1:604053/151‑1 (MQ=255)                                                                                                                         aGGATTAAGTAAATTTTTCCCCGGAGGAAAAGCGTTAGACAACGTTGAATACAGCCTGCGCCGTGCCGAAATCATGGCGCTGCTCAGTGAAAACCGGGCGTGAAAATCAACGCGATACAAAGCATTAACTTGTGGACTCCACGc                               <  1:344585/144‑1 (MQ=255)                                                                                                                          ggATTAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGCGGGAAAATCAACGCTAATCAAAGCATTAACTGGTGTATACCACGCCGATcgcg                       >  1:453190/1‑151 (MQ=255)                                                                                                                          ggATTAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGCGGGAAAATCAACGCTAATCAAAGCATTAACTGGTGTATACCACGCCGATcgcg                       <  1:919895/151‑1 (MQ=255)                                                                                                                                aGTAAATTTTTCCCCGGCGTAAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGCGGGAAAATCAACGCTAATCAAAGCATTAACTGGTGTATACCACGCCGATCGCGGCAcc                  >  1:566070/1‑150 (MQ=255)                                                                                                                                  tAAATTTTTCCCCGGCGTAAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGCGGGAAAATCAACGCTAATCAAAGCATTAACTGGTGTATACCACGCCGATCGCGGCACCATc               >  2:31510/1‑151 (MQ=255)                                                                                                                                     aTTTTTCCCCGGCGTAAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGCGGGAAAATCAACGCTAATCAAAGCATTAACTGGTGTATACCACGCCGATCGCGGCACCATctgg            <  2:4780/151‑1 (MQ=255)                                                                                                                                         ttCCCCGGCGTAAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGCGGGAAAATCAACGCTAATCAAAGCATTAACTGGTGTATACCACGCCGATCGCGGCACCATctggctg         <  1:598854/150‑1 (MQ=255)                                                                                                                                          tCCCCGGCGTAAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGCGGGAAAATCAACGCTAATCAAAGCATTAACTGGTGTATACCACGCCGATCGCGGCACCCTctggctgg        >  1:92625/1‑150 (MQ=255)                                                                                                                                               ggCGTAAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGCGGGAAAATCAACGCTAATCAAAGCATTAACTGGTGTATACCACGCCGATCGCGGCACCATCTGGCTGGAAGGcc  >  2:77911/1‑151 (MQ=255)                                                                                                                                                    |                                                                                                                                                  GCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTCAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGCGGGAAAATCAACGCTAATCAAAGCATTAACTGGTGTATACCACGCCGATCGCGGCACCATCTGGCTGGAAGGCC  >  NC_000913/4450976‑4451268

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑