Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,842,330	T→C	R25G (AGA→GGA)	yicN ←	DUF1198 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,842,330	0	T	C	100.0%	296.1 / NA	80	R25G (AGA→GGA)	yicN	DUF1198 family protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (37/43);  total (37/43)

GCGCGGCGTTTTGCAGGTGCGACTCATCCGGACGATGCAAATAACGTAAAAATTCGTCACCAGCTGACTTTCCCATTTGATCGACCATCGACTCCACGGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCTGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTACTAATCCACCAGAATAATCTT  >  NC_000913/3842183‑3842475                                                                                                                                                    |                                                                                                                                                  gcgcGGCGTTTTGCAGGTGCGACTCATCCGGACGATGCAAATAACGTAAAAATTCGTCACCAGCTGACTTTCCCATTTGATCGACCATCGACTCCACGGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCCGGc                                                                                                                                                <  1:247264/151‑1 (MQ=255)     ggcgTTTTGCAGGTGCGACTCATCCGGACGATGCAAATAACGTAAAAATTCGTCACCAGCTGACTTTCCCATTTGATCGACCATCGACTCCACGGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCCGGCCCCg                                                                                                                                            <  1:552168/151‑1 (MQ=255)      gcgTTTTGCAGGTGCGACTCATCCGGACGATGCAAATAACGTAAAAATTCGTCACCAGCTGACTTTCCCATTTGATCGACCATCGACTCCACGGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCCGGCCCCgg                                                                                                                                           <  1:1069784/151‑1 (MQ=255)       cgTTTTGCAGGTGCGACTCATCCGGACGATGCAAATAACGTAAAAATTCGTCACCAGCTGACTTTCCCATTTGATCGACCATCGACTCCACGGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCCGGCCCCGGa                                                                                                                                          <  2:912013/151‑1 (MQ=255)        gTTTTGCAGGTGCGACTCATCCGGACGATGCAAATAACGTAAAAATTCGTCACCAGCTGACTTTCCCATTTGATCGACCATCGACTCCACGGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCCGGCCCCGGAt                                                                                                                                         <  1:1140673/151‑1 (MQ=255)         ttttGCAGGTGCGACTCATCCGGACGATGCAAATAACGTAAAAATTCGTCACCAGCTGACTTTCCCATTTGATCGACCATCGACTCCACGGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCCGGCCCCgg                                                                                                                                           >  1:228292/1‑148 (MQ=255)                ggTGCGACTCATCCGGACGATGCAAATAACGTAAAAATTCGTCACCAGCTGACTTTCCCATTTGATCGACCATCGACTCCACGGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCCGGCCCCGGATGTCAGCAc                                                                                                                                 <  2:711141/151‑1 (MQ=255)                 gTGCGACTCATCCGGACGATGCAAATAACGTAAAAATTCGTCACCAGCTGACTTTCCCATTTGATCGACCATCGACTCCACGGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCCGGCCCCGGATGTCAGCAcc                                                                                                                                >  2:810900/1‑151 (MQ=255)                   gccaCTCATCCGGACGATGCAAATAACGTAAAAATTCGTCACCAGCTGACTTTCCCATTTGATCGACCATCGACTCCACGGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCCGGCCCCGGATGTCAGCACCCg                                                                                                                              >  1:702919/4‑151 (MQ=255)                    cGACTCATCCGGACGATGCAAATAACGTAAAAATTCGTCACCAGCTGACTTTCCCATTTGATCGACCATCGACTCCACGGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCCGGCCCCGGATGTCAGCACCCGa                                                                                                                             <  2:706033/151‑1 (MQ=255)                     gtcTCATCCGGACGATGCAAATAACGTAAAAATTCGTCACCAGCTGACTTTCCCATTTGATCGACCATCGACTCCACGGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCCGGCCCCGGATGTCAGCACCCGa                                                                                                                             >  2:613133/3‑150 (MQ=255)                     gACTCATCCGGACGATGCAAATAACGTAAAAATTCGTCACCAGCTGACTTTCCCATTTGATCGACCATCGACTCCACGGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCCGGCCCCGGATGTCAGCACCCGa                                                                                                                             >  2:613138/1‑150 (MQ=255)                     gACTCATCCGGACGATGCAAATAACGTAAAAATTCGTCACCAGCTGACTTTCCCATTTGATCGACCATCGACTCCACGGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCCGGCCCCGGATGTCAGCACCCGa                                                                                                                             >  2:558945/1‑150 (MQ=255)                      aCTCATCCGGACGATGCAAATAACGTAAAAATTCGTCACCAGCTGACTTTCCCATTTGATCGACCATCGACTCCACGGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCCGGCCCCGGATGTCAGCCCCCGaaa                                                                                                                           <  2:974618/151‑1 (MQ=255)                       cTCATCCGGACGATGCAAATAACGTAAAAATTCGTCACCAGCTGACTTTCCCATTTGATCGACCATCGACTCCACGGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCCGGCCCCGGATGTCAGCACCCGaaa                                                                                                                           >  2:358523/1‑150 (MQ=255)                        tCATCCGGACGATGCAAATAACGTAAAAATTCGTCACCAGCTGACTTTCCCATTTGATCGACCATCGACTCCACGGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCCGGCCCCGGATGTCAGCACCCGAAaac                                                                                                                         >  2:887934/1‑151 (MQ=255)                           tCCGGACGATGCAAATAACGTAAAAATTCGTCACCAGCTGACTTTCCCATTTGATCGACCATCGACTCCACGGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCCGGCCCCGGATGTCAGCACCCGAAaacca                                                                                                                       >  2:594807/1‑150 (MQ=255)                             cGGACGATGCAAATAACGTAAAAATTCGTCACCAGCTGACTTTCCCATTTGATCGACCATCGACTCCACGGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCCGGCCCCGGATGTCAGCACCCGAAaaccaacc                                                                                                                    >  2:123571/1‑151 (MQ=255)                                aCGATGCAAATAACGTAAAAATTCGTCACCAGCTGACTTTCCCATTTGATCGACCATCGACTCCACGGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCCGGCCCCGGATGTCAGCACCCGAAAACCAACcaca                                                                                                                 <  2:840117/151‑1 (MQ=255)                                aCGATGCAAATAACGTAAAAATTCGTCACCAGCTGACTTTCCCATTTGATCGACCATCGACTCCACGGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCCGGCCCCGGATGTCAGCACCCGAAAACCAACcac                                                                                                                  >  1:777705/1‑150 (MQ=255)                                   aTGCAAATAACGTAAAAATTCGTCACCAGCTGACTTTCCCATTTGATCGACCATCGACTCCACGGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCCGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAac                                                                                                              <  2:138558/151‑1 (MQ=255)                                      cAAATAACGTAAAAATTCGTCACCAGCTGACTTTCCCATTTGATCGACCATCGACTCCACGGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCCGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACt                                                                                                           >  1:142275/1‑151 (MQ=255)                                      cAAATAACGTAAAAATTCGTCACCAGCTGACTTTCCCATTTGATCGACCATCGACTCCACGGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCCGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACt                                                                                                           >  2:165060/1‑151 (MQ=255)                                         aTAACGTAAAAATTCGTCACCAGCTGACTTTCCCATTTGATCGACCATCGACTCCACGGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCCGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTAcc                                                                                                        >  2:1175091/1‑151 (MQ=255)                                         aTAACGTAAAAATTCGTCACCAGCTGACTTTCCCATTTGATCGACCATCGACTCCACGGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCCGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTAcc                                                                                                        >  1:722024/1‑151 (MQ=255)                                              gTAAAAATTCGTCACCAGCTGACTTTCCCATTTGATCGACCATCGACTCCACGGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCCGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAg                                                                                                   <  1:972372/151‑1 (MQ=255)                                                aaaaaTTCGTCACCAGCTGACTTTCCCATTTGATCGACCATCGACTCCACGGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCCGGCCCCGGATGTCAGCACCCGAAAACCAACCACaaa                                                                                                               >  2:186803/1‑137 (MQ=255)                                                aaaaaTTCGTCACCAGCTGACTTTCCCATTTGATCGACCATCGACTCCACGGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCCGGCCCCGGATGTCAGCACCCGAAAACCAACCACaaa                                                                                                               <  1:186801/137‑1 (MQ=255)                                                aaaaaTTCGTCACCAGCTGACTTTCCCATTTGATCGACCATCGACTCCACGGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCCGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCg                                                                                                 <  2:743854/151‑1 (MQ=255)                                                 aaaaTTCGTCACCAGCTGACTTTCCCATTTGATCGACCATCGACTCCACGGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCCGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGt                                                                                                <  1:591282/151‑1 (MQ=255)                                                  aaaTTCGTCACCAGCTGACTTTCCCATTTGATCGACCATCGACTCCACGGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCCGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGt                                                                                                >  2:1014747/1‑150 (MQ=255)                                                  aaaTTCGTCACCAGCTGACTTTCCCATTTGATCGACCATCGACTCCACGGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCCGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTg                                                                                               <  1:507873/151‑1 (MQ=255)                                                  aaaTTCGTCACCAGCTGACTTTCCCATTTGATCGACCATCGACTCCACGGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCCGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTg                                                                                               <  2:263984/151‑1 (MQ=255)                                                    atacGTCACCAGCTGACTTTCCCATTTGATCGACCATCGACTCCACGGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCCGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGc                                                                                             >  2:1134319/4‑151 (MQ=255)                                                           aCCAGCTGACTTTCCCATTTGATCGACCATCGACTCCACGGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCCGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACCCTACCGCCAGCGTGGCGAGCAt                                                                                       <  1:747165/150‑1 (MQ=255)                                                           aCCAGCTGACTTTCCCATTTGATCGACCATCGACTCCACGGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCCGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCa                                                                                        <  1:1112375/149‑1 (MQ=255)                                                           aCCAGCTGACTTTCCCATTTGATCGACCATCGACTCCACGGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCCGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCAtt                                                                                      >  1:1177498/1‑151 (MQ=255)                                                            ccAGCTGACTTTCCCATTTGATCGACCATCGACTCCACGGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCCGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTa                                                                                     <  2:854738/151‑1 (MQ=255)                                                             cAGCTGACTTTCCCATTTAATCGACCATCGACTCCACGGTTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCCGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTAt                                                                                    <  1:355391/151‑1 (MQ=255)                                                               gCTGACTTTCCCATTTGATCGACCATCGACTCCACGGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCCGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATcc                                                                                  <  2:898706/151‑1 (MQ=255)                                                                cTGACTTTCCCATTTGATCGACCATCGACTCCACGGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCCGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATcc                                                                                  >  2:1132047/1‑150 (MQ=255)                                                                     tttCCCATTTGATCGACCATCGACTCCACGGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCCGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATc                                                                             >  2:803132/1‑150 (MQ=255)                                                                     tttCCCATTTGATCGACCATCGACTCCACGGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCCGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATc                                                                             >  1:734301/1‑150 (MQ=255)                                                                     tttCCCATTTGATCGACCATCGACTCCACGGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCCGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATc                                                                             >  1:730407/1‑150 (MQ=255)                                                                     tttCCCATTTGATCGACCATCGACTCCACGGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCCGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATc                                                                             >  2:90667/1‑150 (MQ=255)                                                                           aTTTGATCGACCATCGACTCCACGGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCCGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGtt                                                                      <  2:102304/151‑1 (MQ=255)                                                                           aTTTGATCGACCATCGACTCCACGGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCCGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGtt                                                                      <  1:529661/151‑1 (MQ=255)                                                                            tttGATCGACCATCGACTCCACGGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCCGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTtc                                                                     >  1:1049676/1‑151 (MQ=255)                                                                                  cGACCATCGACTCCACGGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCCGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGc                                                               <  1:910588/151‑1 (MQ=255)                                                                                  cGACCATCGACTCCACGGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCCGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGc                                                               <  2:410935/151‑1 (MQ=255)                                                                                    aCCATCGACTCCACGGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCCGGCCCCGGATGTCAGCACCCGAAACCCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCg                                                              <  1:463367/150‑1 (MQ=255)                                                                                      cATCGACTCCACGGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCCGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTaa                                                           <  1:481582/151‑1 (MQ=255)                                                                                          gACTCCACGGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCCGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATc                                                         >  2:9581/1‑149 (MQ=255)                                                                                              ccACGGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCCGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACt                                                   <  2:1149428/151‑1 (MQ=255)                                                                                              ccACGGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCCGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACt                                                   <  2:862929/151‑1 (MQ=255)                                                                                               cACGGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCCGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAAc                                                    <  1:822559/149‑1 (MQ=255)                                                                                               cACGGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCCGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTc                                                  <  2:265671/151‑1 (MQ=255)                                                                                                  gggTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCCGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATg                                               <  1:1086529/151‑1 (MQ=255)                                                                                                  gggTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCCGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATg                                               <  2:887021/151‑1 (MQ=255)                                                                                                           aTCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCCGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTgg                                       >  1:1189547/1‑150 (MQ=255)                                                                                                            tCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCCGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGa                                      >  2:922952/1‑150 (MQ=255)                                                                                                             cGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCCGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGaa                                     <  2:858226/150‑1 (MQ=255)                                                                                                             cGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCCGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGaa                                     <  1:703634/150‑1 (MQ=255)                                                                                                              gATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCCGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAAc                                   >  1:16770/1‑151 (MQ=255)                                                                                                                  ttCAGCCGATCGCTGAGACGGCGAATCGCTTTTCCGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCt                                <  1:558942/150‑1 (MQ=255)                                                                                                                   tCAGCCGATCGCTGAGACGGCGAATCGCTTTTCCGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTc                               <  1:386990/150‑1 (MQ=255)                                                                                                                    cAGCCGATCGCTGAGACGGCGAATCGCTTTTCCGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGc                             <  1:931440/151‑1 (MQ=255)                                                                                                                    cAGCCGATCGCTGAGACGGCGAATCGCTTTTCCGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGc                             >  2:1118666/1‑151 (MQ=255)                                                                                                                        cGATCGCTGAGACGGCGAATCGCTTTTCCGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGt                         <  2:672664/151‑1 (MQ=255)                                                                                                                           tCGCTGAGACGGCGAATCGCTTTTCCGGCCCCGGATGTCAGCACCCGAAACCCAACCACAAACCCTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTa                       <  2:585121/150‑1 (MQ=255)                                                                                                                            cGCTGAGACGGCGAATCGCTTTTCCGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTAc                      <  2:16770/150‑1 (MQ=255)                                                                                                                            cGCTGAGACGGCGAATCGCTTTTCCGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTACt                     >  1:711283/1‑151 (MQ=255)                                                                                                                             gCTGAGACGGCGAATCGCTTTTCCGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTACt                     >  2:670033/1‑150 (MQ=255)                                                                                                                              cTGAGACGGCGAATCGCTTTTCCGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTACTa                    >  2:413653/1‑150 (MQ=255)                                                                                                                                gagaCGGCGAATCGCTTTTCCGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTACTAAt                  <  2:524781/150‑1 (MQ=255)                                                                                                                                 agaCGGCGAATCGCTTTTCCGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTACTAATcc                >  1:958245/1‑151 (MQ=255)                                                                                                                                     ggCGAATCGCTTTTCCGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTACTAATccacc             >  2:842659/1‑150 (MQ=255)                                                                                                                                              cTTTTCCGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTACTAATCCACCAGAATAATCt   >  2:217115/1‑151 (MQ=255)                                                                                                                                               ttttCCGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTACTAATCCACCAGAATAATc    <  2:777715/149‑1 (MQ=255)                                                                                                                                               ttttCCGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTACTAATCCACCAGAATAATCtt  >  1:360343/1‑151 (MQ=255)                                                                                                                                                    |                                                                                                                                                  GCGCGGCGTTTTGCAGGTGCGACTCATCCGGACGATGCAAATAACGTAAAAATTCGTCACCAGCTGACTTTCCCATTTGATCGACCATCGACTCCACGGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCTGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTACTAATCCACCAGAATAATCTT  >  NC_000913/3842183‑3842475

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