Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,483,658	A→T	5.8%	D60E (GAT→GAA) ‡	valS ←	valyl‑tRNA synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,483,658	0	A	T	5.8%	320.9 / 5.8	97	D60E (GAT→GAA) ‡	valS	valyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base A (43/48);  new base T (3/3);  total (46/51)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.55e-01

AAGCTTCGCGGCCGTAGTCGTGACGGGTTTTACCTTCTTCTGCGGCAATCTTGCGCTCAACGACCATCTGGGTAGCGATCCCGGCGTGGTCAGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGTTCGGCGGCGGGATCATGATGCAGAAACTTTCCTGGCTTTCATCGCCATTAGGCTTAAAGTAGCCCTGCTTTTCCCAGTGCTCGTAAAGC  >  NC_000913/4483510‑4483799                                                                                                                                                     |                                                                                                                                              aaGCTTCGCGGCCGTAGTCGTGACGGGTTTTACCTTCTTCTGCGGCAATCTTGCGCTCAACGACCATCTGGGTAGCGATCCCGGCGTGGTCAGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTa                                                                                                                                               >  4:28079/1‑149 (MQ=255)   gCTTCGCGGCCGTAGTCGTGACGGGTTTTACCTTCTTCTGCGGCAATCTTGCGCTCAACGACCATCTGGGTAGCGATCCCGGCGTGGTCAGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATc                                                                                                                                             <  5:369648/149‑1 (MQ=255)   gCTTCGCGGCCGTAGTCGTGACGGGTTTTACCTTCTTCTGCGGCAATCTTGCGCTCAACGACCATCTGGGTAGCGATCCCGGCGTGGTCAGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATc                                                                                                                                             <  4:110591/149‑1 (MQ=255)     ttCGCGGCCGTAGTCGTGACGGGTTTTACCTTCTTCTGCGGCAATCTTGCGCTCAACGACCATCTGGGTAGCGATCCCGGCGTGGTCAGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATcaa                                                                                                                                           >  4:286256/1‑147 (MQ=255)       cgcgGCCGTAGTCGTGACGGGTTTTACCTTCTTCTGCGGCAATCTTGCGCTCAACGACCATCTGGGTAGCGATCCCGGCGTGGTCAGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCa                                                                                                                                           <  4:354003/147‑1 (MQ=255)        gcgGCCGTAGTCGTGACGGGTTTTACCTTCTTCTGCGGCAATCTTGCGCTCAACGACCATCTGGGTAGCGATCCCGGCGTGGTCAGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCatga                                                                                                                                        >  3:117055/1‑149 (MQ=255)         cgGCCGTAGTCGTGACGGGTTTTACCTTCTTCTGCGGCAATCTTGCGCTCAACGACCATCTGGGTAGCGATCCCGGCGTGGTCAGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCatgat                                                                                                                                       <  3:28079/149‑1 (MQ=255)          ggCCGTAGTCGTGACGGGTTTTACCTTCTTCTGCGGCAATCTTGCGCTCAACGACCATCTGGGTAGCGATCCCGGCGTGGTCAGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCatgatg                                                                                                                                      <  7:220840/149‑1 (MQ=255)          ggCCGTAGTCGTGACGGGTTTTACCTTCTTCTGCGGCAATCTTGCGCTCAACGACCATCTGGGTAGCGATCCCGGCGTGGTCAGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCatgatg                                                                                                                                      <  5:376579/149‑1 (MQ=255)            ccGTAGTCGTGACGGGTTTTACCTTCTTCTGCGGCAATCTTGCGCTCAACGACCATCTGGGTAGCGATCCCGGCGTGGTCAGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGt                                                                                                                                    <  7:219587/149‑1 (MQ=255)             cGTAGTCGTGACGGGTTGTTCCTTCTTCTGCGGCTATCTTGAGCTCAACGACCATGTGGTTAGCGATCCCGGCGTGGTCAGGACCGACCTGCCACAGGGTGTTTGTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGtt                                                                                                                                   <  7:369720/149‑1 (MQ=255)               tAGTCGTGACGGGTTTTACCTTCTTCTGCGGCAATCTTGCGCTCAACGACCATCTGGGTAGCGATCCCGGCGTGGTCAGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTg                                                                                                                                 >  6:341555/1‑149 (MQ=255)                 gTCGTGACGGGTTTTACCTTCTTCTGCGGCAATCTTGCGCTCAACGACCATCTGGGTAGCGATCCCGGCGTGGTCAGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCt                                                                                                                               >  8:249627/1‑149 (MQ=255)                      gACGGGTTTTACCTTCTTCTGCGGCAATCTTGCGCTCAACGACCATCTGAGTAGCGATCCCGGCGTGGTCAGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCa                                                                                                                                           >  8:62508/1‑132 (MQ=255)                      gACGGGTTTTACCTTCTTCTGCGGCAATCTTGCGCTCAACGACCATCTGAGTAGCGATCCCGGCGTGGTCAGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCa                                                                                                                                           <  7:62508/132‑1 (MQ=255)                        cGGGTTTTACCTTCTTCTGCGGCAATCTTGCGCTCAACGACCATCTGGGTAGCGATCCCGGCGTGGTCAGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGc                                                                                                                        >  3:75454/1‑149 (MQ=255)                        cGGGTTTTACCTTCTTCTGCGGCAATCTTGCGCTCAACGACCATCTGGGTAGCGATCCCGGCGTGGTCAGTACCGACCTGCCACAGGGTGGTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGc                                                                                                                        <  5:224729/149‑1 (MQ=255)                        cGGGTTTTACCTTCTTCTGCGGCAATCTTGCGCTCAACGACCATCTGGGTAGCGATCCCGGCGTGGGCAGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGc                                                                                                                        <  1:324420/149‑1 (MQ=255)                          ggTTTTACCTTCTTCTGCGGCAATCTTGCGCTCAACGACCATCTGGGTAGCGATCCCGGCGTGGTCAGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGt                                                                                                                      <  1:237281/149‑1 (MQ=255)                            ttttACCTTCTTCTGCGGCAATCTTGCGCTCAACGACCATCTGGGTAGCGATCCCGGCGTGGTCAGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGa                                                                                                                    <  6:90814/149‑1 (MQ=255)                                 ccTTCTTCTGCGGCAATCTTGCGCTCAACGACCATCTGGGTAGCGATCCCGGCGTGGTCAGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCat                                                                                                               <  1:116448/149‑1 (MQ=255)                                  cttcttCTGCGGCAATCTTGCGCTCAACGACCATCTGGGTAGCGATCCCGGCGTGGTCAGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCat                                                                                                                                          <  5:13724/121‑1 (MQ=255)                                  cttcttCTGCGGCAATCTTGCGCTCAACGACCATCTGGGTAGCGATCCCGGCGTGGTCAGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCat                                                                                                                                          >  6:13724/1‑121 (MQ=255)                                     cttctGCGGCAATCTTGCGCTCAACGACCATCTGGGTAGCGATCCCGGCGTGGTCAGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCatat                                                                                                             <  5:320088/147‑1 (MQ=255)                                     cttctGCGGCAATCTTGCGCTCAACGACCATCTGGGTAGCGATCCCGGCGTGGTCAGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCATATg                                                                                                            >  2:98527/1‑148 (MQ=255)                                      ttctGCGGCAATCTTGCGCTCAACGACCATCTGGGTAGCGATCCCGGCGTGGTCAGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCATATGc                                                                                                           >  8:121708/1‑148 (MQ=255)                                         tGCGGCAATCTTGCGCTCAACGACCATCTGGGTAGCGATCCCGGCGTGGTCAGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGa                                                                                                                            >  4:205652/1‑128 (MQ=255)                                            ggCAATCTTGCGCTCAACGACCATCTGGGTAGCGATCCCGGCGTGGTCAGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGc                                                                                                    >  1:185231/1‑149 (MQ=255)                                             gCAATCTTGCGCTCAACGACCATCTGGGTAGCGATCCCGGCGTGGTCAGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGc                                                                                                                        >  1:181679/1‑128 (MQ=255)                                             gCAATCTTGCGCTCAACGACCATCTGGGTAGCGATCCCGGCGTGGTCAGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGc                                                                                                                        <  2:181679/128‑1 (MQ=255)                                             gCAATCTTGCGCTCAACGACCATCTGGGTAGCGATCCCGGCGTGGTCAGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGcc                                                                                                   >  4:149945/1‑149 (MQ=255)                                                      cgcTCAACGACCATCTGGGTAGCGATCCCGGCGTGGTCAGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGt                                                                                           <  3:149945/148‑1 (MQ=255)                                                       gcTCAACGACCATCTGGGTAGCGATCCCGGCGTGGTCAGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGTTc                                                                                         >  6:218155/1‑149 (MQ=255)                                                        cTCAACGACCATCTGGGTAGCGATCCCGGCGTGGTCAGTACCGACCTGCCACAGGGTGTTTTTGACCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGTTcg                                                                                        >  8:81528/1‑149 (MQ=255)                                                         tCAACGACCATCTGGGTAGCGATCCCGGCGTGGTCAGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGTTcgg                                                                                       >  2:233363/1‑149 (MQ=255)                                                         tCAACGACCATCTGGGTAGCGATACCGGCGTGGTCAGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGTTcgg                                                                                       >  7:32313/1‑149 (MQ=255)                                                               aCCATCTGGGTAGCGATCCCGGCGTGGTCAGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGTTcggcggcgg                                                                                 <  4:117055/149‑1 (MQ=255)                                                                  aTCTGGGTAGCGATCCCGGCGTGGTCAGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGAccc                                                                                                                 <  8:244700/114‑1 (MQ=255)                                                                  aTCTGGGTAGCGATCCCGGCGTGGTCAGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGAccc                                                                                                                 >  7:244700/1‑114 (MQ=255)                                                                       ggTAGCGATCCCGGCGTGGTCAGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGTTCGGCGGCGGGATCatga                                                                         <  4:276611/149‑1 (MQ=255)                                                                       ggTAGCGATCCCGGCGTGGTCAGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCACATGCAAACTGCCGGTGACGTTCGGCGGCGGGATc                                                                             <  4:257295/145‑1 (MQ=255)                                                                        gTAGCGATCCCGGCGTGGTCAGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGTTCGGCGGCGGGATCatgat                                                                        <  7:249627/149‑1 (MQ=255)                                                                        gTAGCGATCCCGGCGTGGTCAGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGGTCGGCGGCGGGATCatgat                                                                        <  1:233363/149‑1 (MQ=255)                                                                        gTAGCGATCCCGGCGTGGTCAGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTATGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGTTCGGCGGCGGGATCatga                                                                         <  4:275975/148‑1 (MQ=255)                                                                        gTAGCGATCCCGGCGTGGTCAGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTATGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGTTCGGCGGCGGGATCatga                                                                         >  3:275975/1‑148 (MQ=255)                                                                         tAGCGATCCCGGCGTGGTCAGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGTTCGGCGGCGGGATCatgatg                                                                       >  6:295689/1‑149 (MQ=255)                                                                           gCGATCCCGGCGTGGTCAGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGTTCGGCGGCGGGATCATGATGc                                                                      <  4:75454/148‑1 (MQ=255)                                                                             gATCCCGGCGTGGTCAGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGTTCGGCGGCGGGATCATGATGCAg                                                                    >  4:36461/1‑148 (MQ=255)                                                                                    gCGTGGTCAGTACCGACCTGCGACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGTTCGGCGGCGGTATCATGATGCAGAAACTTTc                                                            <  8:32313/149‑1 (MQ=255)                                                                                    gCGTGGTCAGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGTTCGGCGGCGGGATCATGATGCAGAAACTTTc                                                            <  4:325663/149‑1 (MQ=255)                                                                                    gCGTGGTCAGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGTTCGGCGGCGGGATCATGATGCAGAAACTTTc                                                            <  3:36461/149‑1 (MQ=255)                                                                                      gTGGTCAGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGTTCGGCGGCGGGATCATGATGCAGAAACTTTCCt                                                          <  2:185231/149‑1 (MQ=255)                                                                                        ggTCAGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTa                                                                                                                                               <  7:91359/62‑1 (MQ=255)                                                                                        ggTCAGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTa                                                                                                                                               >  8:91359/1‑62 (MQ=255)                                                                                         gTCAGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCatgatg                                                                                                                                      <  5:262114/70‑1 (MQ=255)                                                                                         gTCAGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCatgatg                                                                                                                                      >  6:262114/1‑70 (MQ=255)                                                                                         gTCAGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGTTCGGCGGCGGGATCATGATGCAGAAACTTTCCTGGc                                                       >  2:95617/1‑149 (MQ=255)                                                                                         gTCAGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGTTCGGCGGCGGGATCATGATGCAGAAACTTTCCTGGc                                                       >  3:226260/1‑149 (MQ=255)                                                                                          tCAGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGa                                                                                              >  4:295710/1‑109 (MQ=255)                                                                                          tCAGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGa                                                                                              <  3:295710/109‑1 (MQ=255)                                                                                            aGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGTTCGGCGGCGGGa                                                                               <  5:298076/122‑1 (MQ=255)                                                                                            aGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGTTCGGCGGCGGGATCATGATGCAGAAACTTTCCTGGCttt                                                    <  5:295689/149‑1 (MQ=255)                                                                                            aGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACACATATGCAAACTGCCGGTGACGTTCGGCGGCGGGa                                                                               >  6:298076/1‑122 (MQ=255)                                                                                             gTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGTTCGGCGGCGGGATCATGATGCAGAAACTTTCCTGGCTTTc                                                   <  5:341555/149‑1 (MQ=255)                                                                                                ccGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGTTCGGCGGCGGGATCATGATGCAGAAACTTTCCTGGCTTTCATc                                                <  1:147920/149‑1 (MQ=255)                                                                                                ccGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGTTCGGCGGCGGGATCATGATGCAGAAACTTTCCTGGCTTTCATc                                                <  1:95617/149‑1 (MQ=255)                                                                                                ccGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGTTCGGCGGCGGGATCATGATGCAGAAACTTTCCTGGCTTTCATc                                                >  2:373524/1‑149 (MQ=255)                                                                                                     cTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGTTCGGCGGCGGGATCATGATGCAGAAACTTTCCTGGCTTTCATCGCCAt                                           <  3:286256/149‑1 (MQ=255)                                                                                                     cTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGTTCGGCGGCGGGATCATGATGCAGAAACTTTCCTGGCTTTCATCGCCAt                                           >  1:238740/1‑149 (MQ=255)                                                                                                     cTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCAATATGCAAACTGCCGGTGACGTTCGGCGGCGGGATCATGATGTAGAAACTTTCCTGGCTTTCATCGCCAt                                           >  6:200814/1‑149 (MQ=255)                                                                                                       gCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGGGACGCATATGCAAACTGCCGGGGCCGTTCGGCGGCGGGATCATGATGCAGAAACTTTCCTGTCTTTCATCGCCAt                                           >  3:38665/1‑147 (MQ=255)                                                                                                        ccACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGTTCTTGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGTTCGGCGGCGGTATCATGATCCAGAAACTTTCCTGGCTTTCATCGCCATTa                                         <  2:7967/148‑1 (MQ=255)                                                                                                        ccACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGTTCGGCGGCGGGATCATGATGCAGAAACTTTCCTGGCTTTCATCGCCATTa                                         <  5:142927/148‑1 (MQ=255)                                                                                                        ccACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGTTCGGCGGCGGGATCATGATGCAGAAACTTTCCTGGCTTTCATCGCCATTa                                         <  7:121708/148‑1 (MQ=255)                                                                                                         cacaGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGTTCGGCGGCGGGATCATGATGCAGAAACTTTCCTGGCTTTCATCGCCATTAg                                        <  1:274264/148‑1 (MQ=255)                                                                                                         cacaGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGTTCGGCGGCGGGATCATGATGCAGAAACTTTCCTGGCTTTCATCGCCATTAg                                        >  2:274264/1‑148 (MQ=255)                                                                                                          acaGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGTTCGGCGGCGGGATCATGATGCAGAAACTTTCCTGGCTTTCATCGCCATTAGGc                                      <  5:218155/149‑1 (MQ=255)                                                                                                          acaGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGTTCGGCGGCGGGATCATGATGCAGAAACTTTCCTGGCTTTCATCGCCATTAGGc                                      >  1:373046/1‑149 (MQ=255)                                                                                                          acaGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGTTCGGCGGCGGGATCATGATGCAGAAACTTTCCTGGCTTTCATCGCCATTAGGc                                      >  4:30767/1‑149 (MQ=255)                                                                                                          acaGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGTTCGGCGGCGGGATCATGATGCAGAAACTTTCCTGGCTTTCATCGCCATTAGGc                                      >  5:86587/1‑149 (MQ=255)                                                                                                              ggTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGTTCGGCGGCGGGATCATGATGCAGAAACTTTCCTGGCTTTCATCGCCATTAGGCTTa                                   >  8:87081/1‑148 (MQ=255)                                                                                                               gtgtTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTCTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGTTCGGCGGCGGGATCATGATGCAGAAACTTTCCTGGCTTTCATCACCATTAGGCTTa                                   >  6:290348/1‑147 (MQ=255)                                                                                                                tgtTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGTTCGGCGGCGGGATCATGATGCAGAAACTTTCCTGGCTTTCATCGCCATTAGGCTTAAAg                                <  7:81528/149‑1 (MQ=255)                                                                                                                         ccTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGTTCGGCGGCGGGATCATGATGCAGAAACTTTCCTGGCTTTCATCGCCATTAGGCTTAAAGTAGCCCTGc                       >  6:65410/1‑149 (MQ=255)                                                                                                                         ccTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGTTCGGCGGCGGGATCATGATGCAGAAACTTTCCTGGCTTTCATCGCCATTAGGCTTAAAGTAGCCCTGc                       >  4:85706/1‑149 (MQ=255)                                                                                                                           tGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGTTCGGCGGCGGGATCATGATGCAGAAACTTTc                                                            >  2:435/1‑110 (MQ=255)                                                                                                                           tGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGTTCGGCGGCGGGATCATGATGCAGAAACTTTc                                                            <  1:435/110‑1 (MQ=255)                                                                                                                            gCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGTTCGGCGGCGGGATCATGATGCAGAAACTTTCCTGGCTTTCATCGCCATTAGGCTTAAAGTAGCCCTGCttt                    >  5:29942/1‑149 (MQ=255)                                                                                                                            gCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGTTCGGCGGCGGGATCATGATGCAGAAACTTTCCTGGCTTTCATCGCCATTAGGCTTAAAGTAGCCCTGCttt                    >  7:213303/1‑149 (MQ=255)                                                                                                                                 cgcTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGTTCGGCGGCGGGATCATGATGCAGAAACTTTCCTGGCTTTCATCGCCATTAGGCTTAAAGTAGCCCTGCTTTTCCCa               <  7:87081/149‑1 (MQ=255)                                                                                                                                 cgcTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGTTCGGCGGCGGGATCATGATGCAGAAACTTTCCTGGCTTTCATCGCCATTAGGCTTAAAGTAGCCCTGCTTTTCCCa               <  4:226260/149‑1 (MQ=255)                                                                                                                                    tGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGTTCGGCGGCGGGATCATGATGCAGAAACTTTCCTGGCTTTCATCGCCATTAGGCTTAAAGTAGCCCTGCTTTTCCCAGTg            >  5:59951/1‑149 (MQ=255)                                                                                                                                    tGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGTTCGGCGGCGGGATCATGATGCAGAAACTTTCCTGGCTTTCATCGCCATTACGCTTAAAGTAGCCCTGCTTTTCCCAGTg            >  1:209045/1‑149 (MQ=255)                                                                                                                                         gCGGATCATGGTATCCATGATGGTCTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGTTCGGCGGCGGGATCATGATGCAGAAACTTTCCTGGCTTTCATCACCATTAGGCTTAAAGTAGCCCTGCTTTTCCCAGTGCTCGt       <  5:290348/149‑1 (MQ=255)                                                                                                                                         aaccaTCATGGATACCATGTTGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGTTCGGCGGCGGGATCATGATGCAGAAACTTTCCTTGc                                                       <  5:328475/97‑1 (MQ=255)                                                                                                                                         aaccaTCATGGATACCATGATGGTTTGCTGGAAGGCGTGACACATATGCAAACAGCCGGTCACGAACGGCGGCGGGATCATGATACACAAACTAACCTGGc                                                       >  6:328475/5‑101 (MQ=255)                                                                                                                                          accaTCATGGATACCATGATGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGTTCGGCGGCGGGATCATGATGCAGAAACTTTCCt                                                          <  2:189459/94‑1 (MQ=255)                                                                                                                                          accaTCATGGATACCATGATGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGTTCGGCGGCGGGATCATGATGCAGAAACTTTCCt                                                          >  1:189459/4‑97 (MQ=255)                                                                                                                                          accaTCATGGATACCATGATGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGTTCGGCGGCGGGATCATGATGCAGAAACTTTCCTGGCTTTCATCGCCATTAGGCTTAAAGTAGCCCTGCTTTTCCCAGTGCTCGt       >  4:113470/4‑148 (MQ=255)                                                                                                                                              tCATGGTTTCCATGATGGTTTGCTGGATGGCGTGACCCATATGCTAACTGCCGGTGACGTTCGGCGGCGGGATCATGATGCAGAAACTTTCCTGGCTTTCATCGCCATTAGGCTTAAAGTAGCCCTGCTTTTCCCAGTGCTCGTAAAGc  <  2:209045/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                              AAGCTTCGCGGCCGTAGTCGTGACGGGTTTTACCTTCTTCTGCGGCAATCTTGCGCTCAACGACCATCTGGGTAGCGATCCCGGCGTGGTCAGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGTTCGGCGGCGGGATCATGATGCAGAAACTTTCCTGGCTTTCATCGCCATTAGGCTTAAAGTAGCCCTGCTTTTCCCAGTGCTCGTAAAGC  >  NC_000913/4483510‑4483799

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑