Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,821,424	G→C	G196G (GGC→GGG)	alaS ←	alanyl‑tRNA synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,821,424	0	G	C	98.0%	342.0 / ‑3.3	98	G196G (GGC→GGG)	alaS	alanyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base G (0/2);  new base C (43/53);  total (43/55)
Fisher's exact test for biased strand distribution p-value = 5.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTCCAGACCCATACCGGTATCTACAGACGGCTTCGGCAGCGGTTCCATCGTGCCATCGGCCTGGCGGTTGAACTGCATGAAGACGATGTTCCAGATCTCAATGTAGCGGTCGCCGTCTTCTTCCGGGCTTCCCGGAGGGCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTGCACGGGCCGCACGGACCAGTGTCACCCATCTGCCAGAAGTTGTCAGATGCGTATGGCGCACCTTTGTTATCGCCGATGCGAATAATACGTTCGCGCGGGATCCC  >  NC_000913/2821285‑2821570                                                                                                                                            |                                                                                                                                                   gCTCCAGACCCATACCGGTATCTACAGACGGCTTCGGCAGCGGTTCCATCGTGCCATCGGCCTGGCGGTTGAACTGCATGAAGACGATGTTCCAGATCTCAATGTAGCGGTCGCCGTCTTCTTCCGGGCTTCCCGGAGGCCCCCCCCAAAt                                                                                                                                         <  4:893722/151‑1 (MQ=255)  cTCCAGACCCATACCGGTATCTACAGACGGCTTCGGCAGCGGTTCCATCGTGCCATCGGCCTGGCGGTTGAACTGCATGAAGACGATGTTCCAGATCTCAATGTAGCGGTCGCCGTCTTCTTCCGGGCTTCCCGGAGGCCCCCCCCAAAt                                                                                                                                         <  2:1463720/150‑1 (MQ=255)   tCCAGACCCATACCGGTATCTACAGACGGCTTCGGCAGCGGTTCCATCGTGCCATCGGCCTGGCGGTTGAACTGCATGAAAACGGTGTTCCAGATCTCAATGTAGCGGGCGCCGTCCTCTTCCCGGCTTCCCGGAGGCCCCCCCCCAAtc                                                                                                                                        >  2:747145/1‑149 (MQ=255)         cccATACCGGTATCTACAGACGGCTTCGGCAGCGGTTCCATCGTGCCATCGGCCTGGCGGTTGAACTGCATGAAGACGATGTTCCAGATCTCAATGTAGCGGTCGCCGTCTTCTTCCGGGCTTCCCGGAGGCCCCCCCCAAATGTGGTCGc                                                                                                                                 <  4:1304604/151‑1 (MQ=255)          ccATACCGGTATCTACAGACGGCTTCGGCAGCGGTTCCATCGTGCCATCGGCCTGGCGGTTGAACTGCATGAAGACGATGTTCCAGATCTCAATGTAGCGGTCGCCGTCTTCTTCCGGGCTTCCCGGAGGCCCCCCCCAAATGTGGTCGcc                                                                                                                                >  1:214176/1‑151 (MQ=255)          ccATACCGGTATCTACAGACGGCTTCGGCAGCGGTTCCATCGTGCCATCGGCCTGGCGGTTGAACTGCATGAAGACGATGTTCCAGATCTCAATGTAGCGGTCGCCGTCTTCTTCCGGGCTTCCCGGAGGCCCCCCCCAAATGTGGTCGcc                                                                                                                                <  2:1095802/151‑1 (MQ=255)            aTACCGGTATCTACAGACGGCTTCGGCAGCGGTTCCATCGTGCCATCGGCCTGGCGGTTGAACTGCATGAAGACGATGTTCCAGATCTCAATGTAGCGGTCGCCGTCTTCTTCCGGGCTTCCCGGAGGCCCCCCCCAAATGTGGTCGCCg                                                                                                                               >  2:557419/1‑150 (MQ=255)              aCCGGTATCTACAGACGGCTTCGGCAGCGGTTCCATCGTGCCATCGGCCTGGCGGTTGAACTGCATGAAGACGATGTTCCAGATCTCAATGTAGCGGTCGCCGTCTTCTTCCGGGCTTCCCGGAGGCCCCCCCCAAATGTGGTCGCCGTg                                                                                                                             >  2:1004920/1‑150 (MQ=255)                cGGTATCTACAGACGGCTTCGGCAGCGGTTCCATCGTGCCATCGGCCTGGCGGTTGAACTGCATGAAGACGATGTTCCAGATCTCAATGTAGCGGTCGCCGTCTTCTTCCGGACTTCCCGGAGGCCCCCCCCAAATGTGGTCGCCGTGATc                                                                                                                          <  2:328607/151‑1 (MQ=255)                  gTATCTACAGACGGCTTCGGCAGCGGTTCCATCGTGCCATCGGCCTGGCGGTTGAACTGCATGAAGACGATGTTCCAGATCTCAATGTAGCGGTCGCCGTCTTCTTCCGGGCTTCCCGGAGGCCCCCCCCAAATGTGGTCGCCGTGATCGt                                                                                                                        >  4:134656/1‑151 (MQ=255)                    aTCTACAGACGGCTTCGGCAGCGGTTCCATCGTGCCATCGGCCTGGCGGTTGAACTGCATGAAGACGATGTTCCAGATCTCAATGTAGCGGTCGCCGTCTTCTTCCGGGCTTCCCGGAGGCCCCCCCCAAATGTGGTCGCCGTGATCGTa                                                                                                                       <  2:1425186/150‑1 (MQ=255)                     tCTACAGACGGCTTCGGCAGCGGTTCCATCGTGCCATCGGCCTGGCGGTTGAACTGCATGAAGACGATGTTCCAGATCTCAATGTAGCGGTCGCCGTCTTCTTCCGGGCTTCCCGGAGGCCCCCCCCAAATGTGGTCGCCGTGATCGt                                                                                                                        >  2:563588/1‑148 (MQ=255)                        aCAGACGGCTTCGGCAGCGGTTCCATCGTGCCATCGGCCTGGCGGTTGAACTGCATGAAGACGATGTTCCAGATCTCAATGTAGCGGTCGCCGTCTTCTTCCGGGCTTCCCGGAGGCCCCCCCCAAATGTGGTCGCCGTGATCGTagaa                                                                                                                    <  1:69917/149‑1 (MQ=255)                         cAGACGGCTTCGGCAGCGGTTCCATCGTGCCATCGGCCTGGCGGTTGAACTGCATGAAGACGATGTTCCAGATCTCAATGTAGCGGTCGCCGTCTTCTTCCGGGCTTCCCGGAGGCCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGAt                                                                                                                 <  2:961544/151‑1 (MQ=255)                          aGACGGCTTCGGCAGCGGTTCCATCGTGCCATCGGCCTGGCGGTTGAACTGCATGAAGACGATGTTCCAGATCTCAATGTAGCGGTCGCCGTCTTCTTCCGGGCTTCCCGGAGGCCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGAtt                                                                                                                <  4:239946/151‑1 (MQ=255)                          aGACGGCTTCGGCAGCGGTTCCATCGTGCCATCGGCCTGGCGGTTGAACTGCATGAAGACGATGTTCCAGATCTCAATGTAGCGGTCGCCGTCTTCTTCCGGGCTTCCCGGAGGCCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGAtt                                                                                                                >  2:515454/1‑151 (MQ=255)                            aCGGCTTCGGCAGCGGTTCCATCGTGCCATCGGCCTGGCGGTTGAGCTGCATGAAGACGATGTTCCAGATCTCAATGTAGCGGTCGCCGTCTTCTTCCGGGCTTCCCGGAGGCCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGAttt                                                                                                               >  4:345239/1‑150 (MQ=255)                            aCGGCTTCGGCAGCGGTTCCATCGTGCCATCGGCCTGGCGGTTGAACTGCATGAAGACGATGTTCCAGATCTCAATGTAGCGGTCGCCGTCTTCTTCCGGGCTTCCCGGAGGCCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTc                                                                                                              <  1:125803/151‑1 (MQ=255)                              ggCTTCGGCAGCGGTTCCATCGTGCCATCGGCCTGGCGGTTGAACTGCATGAAGACGATGTTCCAGATCTCAATGTAGCGGTCGCCGTCTTCTTCCGGGCTTCCCGGAGGCCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCgg                                                                                                            <  4:413497/151‑1 (MQ=255)                                cTTCGGCAGCGGTTCCATCGTGCCATCGGCCTGGCGGTTGAACTGCATGAAGACGATGTTCCAGATCTCAATGTAGCGGTCGCCGTCTTCTTCCGGGCTTCCCGGGGGCCCCCCCCAAATGTGGTCGCCGGGATCGTAGAAGATTTCGGTg                                                                                                          >  2:882784/1‑151 (MQ=255)                                cTTCGGCAGCGGTTCCATCGTGCCATCGGCCTGGCGGTTGAACTGCATGAAGACGATGTTCCAGATCTCAATGTAGCGGTCGCCGTCTTCTTCCGGGCTTCCCGGAGGCCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTg                                                                                                          >  2:1019215/1‑151 (MQ=255)                                  tCGGCAGCGGTTCCATCGTGCCATCGGCCTGGCGGTTGAACTGCATGAAGACGATGTTCCAGATCTCAATGTAGCGGTCGCCGTCTTCTTCCGGGCTTCCCGGAGGCCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTGc                                                                                                         <  2:552986/150‑1 (MQ=255)                                     gCAGCGGTTCCATCGTGCCATCGGCCTGGCGGTTGAACTGCATGAAGACGATGTTCCAGATCTCAATGTAGCGGTCGCCGTCTTCTTCCGGGCTTCCCGGAGGCCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTGCACgg                                                                                                     >  2:1325284/1‑151 (MQ=255)                                      cAGCGGTTCCATCGTGCCATCGGCCTGGCGGTTGAACTGCATGAAGACGATGTTCCAGATCTCAATGTAGCGGTCGCCGTCTTCTTCCGGGCTTCCCGGAGGCCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTGCACggg                                                                                                    <  4:917564/151‑1 (MQ=255)                                          ggTTCCATCGTGCCATCGGCCTGGCGGTTGAACTGCATGAAGACGATGTTCCAGATCTCAATGTAGCGGTCGCCGTCTTCTTCCGGGCTTCCCGGAGGCCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTGCACGGGCCGc                                                                                                <  4:259596/151‑1 (MQ=255)                                              ccATCGTGCCATCGGCCTGGCGGTTGAACTGCATGAAGACGATGTTCCAGATCTCAATGTAGCGGTCGCCGTCTTCTTCCGGGCTTCCCGGAGGCCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTGCACGGGCCGCACgg                                                                                            >  2:1432129/1‑151 (MQ=255)                                                 tCGTGCCATCGGCCTGGCGGTTGAACTGCATGAAGACGATGTTCCAGATCTCAATGTAGCGGTCGCCGTCTTCTTCCGGGCTTCCCGGAGGCCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTGCACGGGCCGCACGGAcc                                                                                         >  3:472338/1‑151 (MQ=255)                                                 tCGTGCCATCGGCCTGGCGGTTGAACTGCATGAAGACGATGTTCCAGATCTCAATGTAGCGGTCGCCGTCTTCTTCCGGGCTTCCCGGAGGCCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTGCACGGGCCGCACGGAcc                                                                                         >  4:182761/1‑151 (MQ=255)                                                  cGTGCCATCGGCCTGGCGGTTGAACTGCATGAAGACGATGTTCCAGATCTCAATGTAGCGGTCGCCGTCTTCTTCCGGGCTTCCCGGAGGCCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTGCACGGGCCGCACGGACCa                                                                                        <  3:57559/151‑1 (MQ=255)                                                  cGTGCCATCGGCCTGGCGGTTGAACTGCATGAAGACGATGTTCCAGATCTCAATGTAGCGGTCGCCGTCTTCTTCCGGGCTTCCCGGAGGCCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTGCACGGGCCGCACGGACCa                                                                                        <  4:115015/151‑1 (MQ=255)                                                  cGTGCCATCGGCCTGGCGGTTGAACTGCATGAAGACGATGTTCCAGATCTCAATGTAGCGGTCGCCGTCTTCTTCCGGGCTTCCCGGAGGCCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTGCACGGGCCGCACGGACCa                                                                                        <  2:111953/151‑1 (MQ=255)                                                   gTGCCATCGGCCTGGCGGTTGGACTGCATGAAGACGATGTTCCAAATCTCAATGTAGCGGTCGCCGTCTTCTTCCGGGCTTTCCGGAGGCCCCCCCCAAATGTGGTCGCCCTGAACGTAGAAGACTTCGGTGCACGGGCCGCACgg                                                                                            >  2:1066726/1‑146 (MQ=255)                                                   gTGCCATCGGCCTGGCGGTTGAACTGCATGAAGACGATGTTCCAGATCTCAATGTAGCGGTCGCCGTCTTCTTCCGGGCTTCCCGGAGGCCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTGCACGGGCCGCACGGACCa                                                                                        <  4:468129/150‑1 (MQ=255)                                                     gCCATCGGCCTGGCGGTTGAACTGCATGAAGACGATGTTCCAGATCTCAATGTAGCGGTCGCCGTCTTCTTCCGGGCTTCCCGGAGGCCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTGCACGGGCCGCACGGACCAgtg                                                                                     <  2:1462483/151‑1 (MQ=255)                                                     gCCATCGGCCTGGCGGTTGAACTGCATGAAGACGATGTTCCAGATCTCAATGTAGCGGTCGCCGTCTTCTTCCGGGCTTCCCGGAGGCCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTGCACGGGCCGCACGGACCAgt                                                                                      <  4:898073/150‑1 (MQ=255)                                                     gCCATCGGCCTGGCGGTTGAACTGCATGAAGACGATGTTCCAGATCTCAATGTAGCGGTCGCCGTCTTCTTCCGGGCTTCCCGGAGGCCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTGCACGGGCCGCACGGACCAgt                                                                                      <  2:1547786/150‑1 (MQ=255)                                                     gCCATCGGCCTGGCGGTTGAACTGCATGAAGACGATGTTCCAGATCTCAATGTAGCGGTCGCCGTCTTCTTCCGGGCTTCCCGGAGGCCCCCCCCAAATGTGGTCGCCGCGATCGTAGAAGATTTCGGTGCACGGGCCGCACCGGCCAgtg                                                                                     >  4:1034886/1‑151 (MQ=255)                                                          cGGCCTGGCGGTTGAACTGCATGAAGACGATGTTCCAGATCTCAATGTAGCGGTCGCCGTCTTCTTCCGGGCTTCCCGGAGGCCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTGCACGGGCCGCACGGACCAGTGTCAcc                                                                                >  4:1201895/1‑151 (MQ=255)                                                           ggCCTGGCGGTTGAACTGCATGAAGACGATGTTCCAGATCTCAATGTAGCGGTCGCCGTCTTCTTCCGGGCTTCCCGGAGGCCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTGCACGGGCCGCACGGACCAGTGTCAccc                                                                               <  2:134430/151‑1 (MQ=255)                                                               tGGCGGTTGAACTGCATGAAGACGATGTTCCAGATCTCAATGTAGCGGTCGCCGTCTTCTTCCGGGCTTCCCGGAGGCCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTGCACGGGCCGCACGGACCAGTGTCACCCATCt                                                                           >  4:296753/1‑151 (MQ=255)                                                                ggCGGTTGAACTGCATGAAGACGATGTTCCAGATCTCAATGTAGCGGTCGCCGTCTTCTTCCGGGCTTCCCGGAGGCCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTGCACGGGCCGCACGGACCAGTGTCACCCATCTg                                                                          <  3:387087/151‑1 (MQ=255)                                                                 gCGGTTGAACTGCATGAAGACGATGTTCCAGATCTCAATGTAGCGGTCGCCGTCTTCTTCCGGGCTTCCCGGAGGCCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTGCACGGGCCGCACGGACCAGTGTCACCCATCTGc                                                                         <  2:784400/151‑1 (MQ=255)                                                                 gCGGTTGAACTGCATGAAGACGATGTTCCAGATCTCAATGTAGCGGTCGCCGTCTTCTTCCGGGCTTCCCGGAGGCCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTGCACGGGCCGCACGGACCAGTGTCACCCATCTGc                                                                         <  4:330275/151‑1 (MQ=255)                                                                 gCGGTTGAACTGCATGAAGACGATGTTCCAGATCTCAATGTAGCGGTCGCCGTCTTCTTCCGGGCTTCCCGGAGGCCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTGCACGGGCCGCACGGACCAGTGTCACCCATCTGc                                                                         >  4:690953/1‑151 (MQ=255)                                                                   ggTTGAACTGCATGAAGACGATGTTCCAGATCTCAATGTAGCGGTCGCCGTCTTCTTCCGGGCTTCCCGGAGGCCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTGCACGGGCCGCACGGACCAGTGTCACCCATCTGcc                                                                        >  4:284187/1‑150 (MQ=255)                                                                   ggTTGAACTGCATGAAGACGATGTTCCAGATCTCAATGTAGCGGTCGCCGTCTTCTTCCGGGCTTCCCGGAGGCCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTGCACGGGCCGCACGGACCAGTGTCACCCATCTGCCa                                                                       <  2:247842/151‑1 (MQ=255)                                                                      tGAACTGCATGAAGACGATGTTCCAGATCTCAATGTAGCGGTCGCCGTCTTCTTCCGGGCTTCCCGGAGGCCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTGCACGGGCCGCACGGACCAGTGTCACCCATCTGCCAGa                                                                     >  3:384896/1‑150 (MQ=255)                                                                      tGAACTGCATGAAGACGATGTTCCAGATCTCAATGTAGCGGTCGCCGTCTTCTTCCGGGCTTCCCGGAGGCCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTGCACGGGCCGCACGGACCAGTGTCACCCATCTGCCAGa                                                                     >  2:200013/1‑150 (MQ=255)                                                                       gAACTGCATGAAGACGATGTTCCAGATCTCAATGTAGCGGTCGCCGTCTTCTTCCGGGCTTCCCGGAGGCCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTGCACGGGCCGCACGGACCAGTGTCACCCATCTGCCAGaa                                                                    <  3:407050/150‑1 (MQ=255)                                                                       gAACTGCATGAAGACGATGTTCCAGATCTCAATGTAGCGGTCGCCGTCTTCTTCCGGGCTTCCCGGAGGCCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTGCACGGGCCGCACGGACCAGTGTCACCCATCTGCCAGAAg                                                                   <  4:1255838/151‑1 (MQ=255)                                                                        aaCTGCATGAAGACGATGTTCCAGATCTCAATGTAGCGGTCGCCGTCTTCTTCCGGGCTTCCCGGAGGCCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTGCACGGGCCGCACGGACCAGTGTCACCCATCTGCCAGAAGt                                                                  <  4:1475949/151‑1 (MQ=255)                                                                            gCATGAAGACGATGTTCCAGATCTCAATGTAGCGGTCGCCGTCTTCTTCCGGGCTTCCCGGAGGCCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTGCACGGGCCGCACGGACCAGTGTCACCCATCTGCCAGAAGTTGTc                                                              >  1:60892/1‑151 (MQ=255)                                                                             caagAAGACGATGTTCCAAATCTCAATGTAGCGGTCGCCGTCCTCTTCCGGGCTTCCCGGAGGGCCCCCCCAAATTTGGTCGCCGTGAACGTAGAATATTTCGGTGCACGGGCCGCACCGACCAGTGTCACCCATCTGCCAGAAGTTGTCa                                                             <  4:636706/148‑1 (MQ=255)                                                                             cATGAAGACGATGTTCCAGATCTCAATGTAGCGGTCGCCGTCTTCTTCCGGGCTTCCCGGAGGCCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTGCACGGGCCGCACGGACCAGTGTCACCCATCTGCCAGAAGTTGTc                                                              <  2:1053483/150‑1 (MQ=255)                                                                                gAAGACGATGTTCCAGATCTCAATGTAGCGGTCGCCGTCTTCTTCCGGGCTTCCCGGAGGCCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTGCACGGGCCGCACGGACCAGTGTCACCCATCTGCCAGAAGTTGTCAGAt                                                          <  4:81449/151‑1 (MQ=255)                                                                                       aTGTTCCAGATCTCAATGTAGCGGTCGCCGTCTTCTTCCGGGCTTCCCGGAGGCCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTGCACGGGCCGCACGGACCAGTGTCACCCATCTGCCAGAAGTTGTCAGATGCGTATg                                                   >  2:1047757/1‑151 (MQ=255)                                                                                         gTTCCAGATCTCAATGTAGCGGTCGCCGTCTTCTTCCGGGCTTCCCGGAGGCCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTGCACGGGCCGCACGGACCAGTGTCACCCATCTGCCAGAAGTTGTCAGATGCGTATgg                                                  <  4:689548/150‑1 (MQ=255)                                                                                         gTTCCAGATCTCAATGTAGCGGTCGCCGTCTTCTTCCGGGCTTCCCGGAGGCCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTGCACGGGCCGCACGGACCAGTGTCACCCATCTGCCAGAAGTTGTCAGATGCGTATGgc                                                 >  4:382671/1‑151 (MQ=255)                                                                                            ccAGATCTCAATGTAGCGGTCGCCGTCTTCTTCCGGGCTTCCCGGAGGCCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTGCACGGGCCGCACGGACCAGTGTCACCCATCTGCCAGAAGTTGTCAGATGCGTATGgcgc                                               <  4:1677021/150‑1 (MQ=255)                                                                                             cAGATCTCAATGTAGCGGTCGCCGTCTTCTTCCGGGCTTCCCGGAGGCCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTGCACGGGCCGCACGGACCAGTGTCACCCATCTGCCAGAAGTTGTCAGATGCGTATGGCGCAc                                             >  2:1550692/1‑151 (MQ=255)                                                                                             cAGATCCCAATGGAACCGGTGCCGTCTTCTTCCGGGTTTCCCGGAGGCCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTGCACGGGCCGCACGGACCAGTGTCACCCATCTGCCAGAAGTTGTCAGATGCGTATGGCGCAc                                             <  4:1639538/151‑1 (MQ=255)                                                                                              aGATCTCAATGTAGCGGTCGCCGTCTTCTTCCGGGCTTCCCGGAGGCCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTGCACGGGCCGCACGGACCAGTGTCACCCATCTGCCAGAAGTTGTCAGATGCGTATGGCGCAcc                                            >  3:474593/1‑151 (MQ=255)                                                                                               gATCTCAATGTAGCGGTCGCCGTCTTCTTCCGGGCTTCCCGGAGGCCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTGCACGGGCCGCACGGACCAGTGTCACCCATCTGCCAGAAGTTGTCAGATGCGTATGGCGCACCt                                           >  4:1534565/1‑151 (MQ=255)                                                                                                    cAATGTAGCGGTCGCCGTCTTCTTCCGGGCTTCCCGGAGGCCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTGCACGGGCCGCACGGACCAGTGTCACCCATCTGCCAGAAGTTGTCAGATGCGTATGGCGCACCttt                                         <  3:247252/148‑1 (MQ=255)                                                                                                    cAATGTAGCGGTCGCCGTCTTCTTCCGGGCTTCCCGGAGGCCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTGCACGGGCCGCACGGACCAGTGTCACCCATCTGCCAGAAGTTGTCAGATGCGTATGGCGCACCt                                           <  4:144993/146‑1 (MQ=255)                                                                                                    cAATGTAGCGGTCGCCGTCTTCTTCCGGGCTTCCCGGAGGCCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTGCACGGGCCGCACGGACCAGTGTCACCCATCTGCCAGAAGTTGTCAGATGCGTATGGCGCACCTTTGt                                       <  4:417821/150‑1 (MQ=255)                                                                                                    cAATGTAGCGGTCGCCGTCTTCTTCCGGGCTTCCCGGAGGCCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTGCACGGGCCGCACGGACCAGTGTCACCCATCTGCCAGAAGTTGTCAGATGCGTATGGCGCACCTTTGt                                       <  3:41749/150‑1 (MQ=255)                                                                                                      aTGTAGCGGTCGCCGTCTTCTTCCGGGCTTCCCGGAGGCCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTGCACGGGCCGCACGGACCAGTGTCACCCATCTGCCAGAAGTTGTCAGATGCGTATGGCGCACCTTTGTTAt                                    >  3:270220/1‑151 (MQ=255)                                                                                                            cGGTCGCCGTCTTCTTCCGGGCTTCCCGGAGGCCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTGCACGGGCCGCACGGACCAGTGTCACCCATCTGCCAGAAGTTGTCAGATGCGTATGGCGCACCTTTGTTATCGCCGa                              <  3:225775/151‑1 (MQ=255)                                                                                                            cGGTCGCCGTCTTCTTCCGGGCTTCCCGGAGGCCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTGCACGGGCCGCACGGACCAGTGTCACCCATCTGCCAGAAGTTGTCAGATGCGTATGGCGCACCTTTGTTATCGCCGa                              <  3:279039/151‑1 (MQ=255)                                                                                                            cGGTCGCCGTCTTCTTCCGGGCTTCCCGGAGGCCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTGCACGGGCCGCACGGACCAGTGTCACCCATCTGCCAGAAGTTGTCAGATGCGTATGGCGCACCTTTGTTATCGCCGa                              <  4:748143/151‑1 (MQ=255)                                                                                                             ggTCGCCGTCTTCTTCCGGGCTTCCCGGAGGCCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTGCACGGGCCGCACGGACCAGTGTCACCCATCTGCCAGAAGTTGTCAGATGCGTATGGCGCACCTTTGTTATCGCCGAt                             >  4:1383163/1‑151 (MQ=255)                                                                                                                cGCCGTCTTCTTCCGGGCTTCCCGGAGGCCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTGCACGGGCCGCACGGACCAGTGTCACCCATCTGCCAGAAGTTGTCAGATGCGTATGGCGCACCTTTGTTATCGCCGATGc                           <  4:2214/150‑1 (MQ=255)                                                                                                                cGCCGTCTTCTTCCGGGCTTCCCGGAGGCCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTGCACGGGCCGCACGGACCAGTGTCACCCATCTGCCAGAAGTTGTCAGATGCGTATGGCGCACCTTTGTTATCGCCGATGCg                          >  1:47521/1‑151 (MQ=255)                                                                                                                 gCCGTCTTCTTCCGGGCTTCCCGGAGGCCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTGCACGGGCCGCACGGACCAGTGTCACCCATCTGCCAGAAGTTGTCAGATGCGTATGGCGCACCTTTGTTATCGCCGATGCg                          >  1:372720/1‑150 (MQ=255)                                                                                                                 gCCGTCTTCTTCCGGGCTTCCCGGAGGCCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTGCACGGGCCGCACGGACCAGTGTCACCCATCTGCCAGAAGTTGTCAGATGCGTATGGCGCACCTTTGTTATCGCCGATGCg                          >  4:1626566/1‑150 (MQ=255)                                                                                                                 gCCGTCTTCTTCCGGGCTTCCCGGAGGCCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTGCACGGGCCGCACGGACCAGTGTCACCCATCTGCCAGAAGTTGTCAGATGCGTATGGCGCACCTTTGTTATCGCCGATGCGa                         <  4:1570038/151‑1 (MQ=255)                                                                                                                   cGTCTTCTTCCGGGCTTCCCGGAGGCCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTGCACGGGCCGCACGGACCAGTGTCACCCATCTGCCAGAAGTTGTCAGATGCGTATGGCGCACCTTTGTTATCGCCGATGCGaat                       >  4:370434/1‑151 (MQ=255)                                                                                                                    gTCTTCTTCCGGGCTTCCCGGAGGCCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTGCACGGGCCGCACGGACCAGTGTCACCCATCTGCCAGAAGTTGTCAGATGCGTATGGCGCACCTTTGTTATCGCCGATGCGaat                       >  2:43109/1‑150 (MQ=255)                                                                                                                     tcttctTCCGGGCTTCCCGGAGGCCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTGCACGGGCCGCACGGACCAGTGTCACCCATCTGCCAGAAGTTGTCAGATGCGTATGGCGCACCTTTGTTATCGCCGATGCGaata                      >  2:284962/1‑150 (MQ=255)                                                                                                                       ttcttcCGGGCTTCCCGGAGGCCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTGCACGGGCCGCACGGACCAGTGTCACCCATCTGCCAGAAGTTGTCAGATGCGTATGGCGCACCTTTGTTATCGCCGATGCGaataat                    >  4:191506/1‑150 (MQ=255)                                                                                                                        tcttcCGGGCTTCCCGGAGGCCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTGCACGGGCCGCACGGACCAGTGTCACCCATCTGCCAGAAGTTGTCAGATGCGTATGGCGCACCTTTGTTATCGCCGATGCGaataat                    <  4:748216/149‑1 (MQ=255)                                                                                                                        tcttcCGGGCTTCCCGGAGGCCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTGCACGGGCCGCACGGACCAGTGTCACCCATCTGCCAGAAGTTGTCAGATGCGTATGGCGCACCTTTGTTATCGCCGATGCGAataata                   <  4:1669244/150‑1 (MQ=255)                                                                                                                           tcCGGGCTTCCCGGAGGCCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTGCACGGGCCGCACGGACCAGTGTCACCCATCTGCCAGAAGTTGTCAGATGCGTATGGCGCACCTTTGTTATCGCCGATGCGAATAATAc                  <  2:676596/148‑1 (MQ=255)                                                                                                                           tcCGGGCTTCCCGGAGGCCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTGCACGGGCCGCACGGACCAGTGTCACCCATCTGCCAGAAGTTGTCAGATGCGTATGGCGCACCTTTGTTATCGCCGATGCGAATAATAc                  <  3:288283/148‑1 (MQ=255)                                                                                                                              gggCTTCCCGGA‑GGCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTGCACGGGCCGCACGGACCAGTGTCACCCATCTGCCAGAAGTTGTCAGATGCGTATGGCGCACCTTTGTTATCGCCGATGCGAATAATACGTTcgcg           <  2:1582004/151‑1 (MQ=255)                                                                                                                                  ttCCCGGAGGCCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTGCACGGGCCGCACGGACCAGTGTCACCCATCTGCCAGAAGTTGTCAGATGCGTATGGCGCACCTTTGTTATCGCCGATGCGAATAATACGTTCGCGCgg        >  2:364105/1‑151 (MQ=255)                                                                                                                                    cccGGAGGCCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTGCACGGGCCGCATGGACCAGTGTCACCCATCTGCCAGAAGTTGTCAGATGCGTATGGCGCACCTTTGTTATCGCCGATGCGAATAATACGTTCGCGCggg       <  4:352919/150‑1 (MQ=255)                                                                                                                                    cccGGAGGCCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTGCACGGGCCGCACGGACCAGTGTCACCCATCTGCCAGAAGTTGTCAGATGCGTATGGCGCACCTTTGTTATCGCCGATGCGAATAATACGTTCGCGCggg       <  4:738997/150‑1 (MQ=255)                                                                                                                                    cccGGAGGCCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTGCACGGGCCGCACGGACCAGTGTCACCCATCTGCCAGAAGTTGTCAGATGCGTATGGCGCACCTTTGTTATCGCCGATGCGAATAATACGTTCGCGCggg       >  2:1239578/1‑150 (MQ=255)                                                                                                                                    cccGGAGGCCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTGCACGGGCCGCACGGACCAGTGTCACCCATCTGCCAGAAGTTGTCAGATGCGTATGGCGCACCTTTGTTATCGCCGATGCGAATAATACGTTCGCGCGGGa      <  4:165053/151‑1 (MQ=255)                                                                                                                                    cccGGAGGCCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTGCACGGGCCGCACGGACCAGTGTCACCCATCTGCCAGAAGTTGTCAGATGCGTATGGCGCACCTTTGTTATCGCCGATGCGAATAATACGTTCGCGCGGGa      <  4:1326941/151‑1 (MQ=255)                                                                                                                                    cccGGAGGCCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTGCACAGGCCGCACGGACCAGTGTCACCCATCTGCCAGAAGTTGTCAGATGCGTATGGCGCACCTTTGTTATCGCCGATGCGAATAATACGTTCGCGCggg       >  3:258077/1‑150 (MQ=255)                                                                                                                                     ccGGAGGCCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTGCACGGGCCGCACGGACCAGTGTCACCCATCTGCCAGAAGTTGTCAGATGCGTATGGCGCACCTTTGTTATCGCCGATGCGAATAATACGTTCGCGCGGGAt     <  2:1304771/151‑1 (MQ=255)                                                                                                                                     ccGGAGGCCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTGCACGGGCCGCACGGACCAGTGTCACCCATCTGCCAGAAGTTGTCAGATGCGTATGGCGCACCTTTGTTATCGCCGATGCGAATAATACGTTCGCGCGGGAt     >  2:256154/1‑151 (MQ=255)                                                                                                                                     ccGGAGGCCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTGCACGGGCCGCACGGACCAGTGTCACCCATCTGCCAGAAGTTGTCAGATGCGTATGGCGCACCTTTGTTATCGCCGATGCGAATAATACGTTCGCGCGGGAt     >  4:103820/1‑151 (MQ=255)                                                                                                                                     ccGGAGGCCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTGCACGGGCCGCACGGACCAGTGTCACCCATCTGCCAGAAGTTGTCAGATGCGTATGGCGCACCTTTGTTATCGCCGATGCGAATAATACGTTCGCGCGGGAt     >  3:384793/1‑151 (MQ=255)                                                                                                                                        gAGGCCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTGCACGGGCCGCACGGACCAGTGTCACCCATCTGCCAGAAGTTGTCAGATGCGTATGGCGCACCTTTGTTATCGCCGATGCGAATAATACGTTCGCGCGGGATccc  <  4:176027/151‑1 (MQ=255)                                                                                                                                            |                                                                                                                                                   GCTCCAGACCCATACCGGTATCTACAGACGGCTTCGGCAGCGGTTCCATCGTGCCATCGGCCTGGCGGTTGAACTGCATGAAGACGATGTTCCAGATCTCAATGTAGCGGTCGCCGTCTTCTTCCGGGCTTCCCGGAGGGCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTGCACGGGCCGCACGGACCAGTGTCACCCATCTGCCAGAAGTTGTCAGATGCGTATGGCGCACCTTTGTTATCGCCGATGCGAATAATACGTTCGCGCGGGATCCC  >  NC_000913/2821285‑2821570

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