Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,672,258	T→C	5.1%	L100P (CTG→CCG)	ydgD →	putative peptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,672,258	0	T	C	5.1%	253.2 / 3.2	79	L100P (CTG→CCG)	ydgD	putative peptidase
Reads supporting (aligned to +/- strand):  ref base T (37/38);  new base C (2/2);  total (39/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.30e-01

ACCCAATCACCGTGGGATGCGGTTGGGCAACTGGAAACGGCCAGCGGCAATTTATGTACGGCGACGCTGATTGCACCCAATCTGGCATTAACGGCAGGACACTGTTTATTGACACCTCCAAAGGGTAAAGCCGATAAAGCAGTGGCGCTGCGTTTTGTGTCAAATAAAGGTCTTTGGCGCTATGAGATCCACGACATAGAAGGCCGCGTTGATCCGACACTGGGAAAGCGGTTAAAAGCAGATGGGGATGGTTGGATTGTACCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATT  >  NC_000913/1672110‑1672406                                                                                                                                                     |                                                                                                                                                     aCCCAATCACCGTGGGATGCGGTTGGGCAACTGGAAACGGCCAGCGGCAATTTATGTACGGCGACGCTGATTGCACCCAATCTGGCATTAACGGCAGGACACTGTTTATTGACACCTCCAAAGGGTAAAGCCGATAAAGCAGTGGCGCt                                                                                                                                                      <  2:77660/149‑1 (MQ=255)   ccAATCACCGTGGGATGCGGTTGGGCAACTGGAAACGGCCAGCGGCAATTTATGTACGGCGACGCTGATTGCACCCAATCTGGCATTAACGGCAGGACACTGTTTATTGACACCTCCAAAGGGTAAAGCCGATAAAGCAGTGGCGCTGc                                                                                                                                                    >  1:153938/1‑149 (MQ=255)   ccAATCACCGTGGGATGCGGTTGGGCAACTGGAAACGGCCAGCGGCAATTTATGTACGGCGACGCTGATTGCACCCAATCTGGCATTAACGGCAGGACACTGTTTATTGACACCTCCAAAGGGTAAAGCCGATAAAGCAGTGGCGCTGc                                                                                                                                                    >  1:124759/1‑149 (MQ=255)     aaTCACCGTGGGATGCGGTTGGGCAACTGGAAACGGCCAGCGGCAATTTATGTACGGCGACGCTGATTGCACCCAATCTGGCATTAACGGCAGGACACTGTTTATTGACACCTCCAAAGGGTAAAGCCGATAAAGCAGTGGCGCTGCGt                                                                                                                                                  >  4:239696/1‑149 (MQ=255)      aTCACCGTGGGATGCGGTTGGGCAACTGGAAACGGCCAGCGGCAATTTATGTACGGCGACGCTGATTGCACCCAATCTGGCATTAACGGCAGGACACTGTTTATTGACACCTCCAAAGGGTAAAGCCGATAAAGCAGTGGCGCTGCGtt                                                                                                                                                 <  4:312328/149‑1 (MQ=255)       tCACCGTGGGATGCGGTTGGGCAACTGGAAACGGCCAGCGGCAATTTATGTACGGCGACGCTGATTGCACCCAATCTGGCATTAACGGCAGGACACTGTTTATTGACACCTCCAAAGGGTAAAGCCGATAAAGCAGTGGCGCTGCGttt                                                                                                                                                >  2:591110/1‑149 (MQ=255)       tCACCGTGGGATGCGGTTGGGCAACTGGAAACGGCCAGCGGCAATTTATGTACGGCGACGCTGATTGCACCCAATCTGGCATTAACGGCAGGACACTGTTTATTGACACCTCCAAAGGGTAAAGCCGATAAAGCAGTGGCGCTGCGttt                                                                                                                                                >  4:326947/1‑149 (MQ=255)        cACCGTGGGATGCGGTTGGGCAACTGGAAACGGCCAGCGGCAATTTATGTACGGCGACGCTGATTGCACCCAATCTGGCATTAACGGCAGGACACTGTTTATTGACACCTCCAAAGGGTAAAGCCGATAAAGCAGTGGCGCTGCGtttt                                                                                                                                               <  2:662777/149‑1 (MQ=255)         aCCGTGGGATGCGGTTGGGCAACTGGAAACGGCCAGCGGCAATTTATGTACGGCGACGCTGATTGCACCCAATCTGGCATTAACGGCAGGACACTGTTTATTGACACCTCCAAAGGGTAAAGCCGATAAAGCAGTGGCGCTGCGTTTtg                                                                                                                                              >  4:747405/1‑149 (MQ=255)         aCCGTGGGATGCGGTTGGGCAACTGGAAACGGCCAGCGGCAATTTATGTACGGCGACGCTGATTGCACCCAATCTGGCATTAACGGCAGGACACTGTTTATTGACACCTCCAAAGGGTAAAGCCGATAAAGCAGTGGCGCTGCGTTTtg                                                                                                                                              <  4:597565/149‑1 (MQ=255)          ccGTGGGATGCGGTTGGGCAACTGGAAACGGCCAGCGGCAATTTATGTACGGCGACGCTGATTGCACCCAATCTGGCATTAACGGCAGGACACTGTTTATTGACACCTCCAAAGGGTAAAGCCGATAAAGCAGTGGCGCTGCGTTTtgt                                                                                                                                             <  4:591523/149‑1 (MQ=255)           cGTGGGATGCGGTTGGGCAACTGGAAACGGCCAGCGGCAATTTATGTACGGCGACGCTGATTGCACCCAATCTGGCATTAACGGCAGGACACTGTTTATTGACACCTCCAAAGGGTAAAGCCGATAAAGCAGTGGCGCTGCGTTTtgtg                                                                                                                                            >  2:675708/1‑149 (MQ=255)            gTGGGATGCGGTTGGGCAACTGGAAACGGCCAGCGGCAATTTATGTACGGCGACGCTGATTGCACCCAATCTGGCATTAACGGCAGGACACTGTTTATTGACACCTCCAAAGGGTAAAGCCGATAAAGCAGTGGCGCTGCGTTTTgtgt                                                                                                                                           <  3:223576/149‑1 (MQ=255)            gTGGGATGCGGTTGGGCAACTGGAAACGGCCAGCGGCAATTTATGTACGGCGACGCTGATTGCACCCAATCTGGCATTAACGGCAGGACACTGTTTATTGACACCTCCAAAGGGTAAAGCCGATAAAGCAGTGGCGCTGCGTTTTgtgt                                                                                                                                           <  2:813707/149‑1 (MQ=255)             tGGGATGCGGTTGGGCAACTGGAAACGGCCAGCGGCAATTTATGTACGGCGACGCTGATTGCACCCAATCTGGCATTAACGGCAGGACACTGTTTATTGACACCTCCAAAGGGTAAAGCCGATAAAGCAGTGGCGCTGCGTTTTGTGTc                                                                                                                                          <  2:519941/149‑1 (MQ=255)                    cGGTTGGGCAACTGGAAACGGCCAGCGGCAATTTATGTACGGCGACGCTGATTGCACCCAATCTGGCATTAACGGCAGGACACTGTTTATTGACACCTCCAAAGGGTAAAGCCGATAAAGCAGTGGCGCTGCGTTTTGTGTCAAATaaa                                                                                                                                   <  3:514322/149‑1 (MQ=255)                    cGGTTGGGCAACTGGAAACGGCCAGCGGCAATTTATGTACGGCGACGCTGATTGCACCCAATCTGGCATTAACGGCAGGACACTGTTTATTGACACCTCCAAAGGGTAAAGCCGATAAAGCAGTGGCGCTGCGTTTTGTGTCAAATaaa                                                                                                                                   <  1:855791/149‑1 (MQ=255)                    cGGTTGGGCAACTGGAAACGGCCAGCGGCAATTTATGTACGGCGACGCTGATTGCACCCAATCTGGCATTAACGGCAGGACACTGTTTATTGACACCTCCAAAGGGTAAAGCCGATAAAGCAGTGGCGCTGCGTTTTGTGTCAAATaaa                                                                                                                                   <  4:349585/149‑1 (MQ=255)                      gTTGGGCAACTGGAAACGGCCAGCGGCAATTTATGTACGGCGACGCTGATTGCACCCAATCTGGCATTAACGGCAGGACACTGTTTATTGACACCTCCAAAGGGTAAAGCCGATAAAGCAGTGGCGCTGCGTTTTGTGTCAAATAAAg                                                                                                                                  <  2:401681/148‑1 (MQ=255)                      gTTGGGCAACTGGAAACGGCCAGCGGCAATTTATGTACGGCGACGCTGATTGCACCCAATCTGGCATTAACGGCAGGACACTGTTTATTGACACCTCCAAAGGGTAAAGCCGATAAAGCAGTGGCGCTGCGTTTTGTGTCAAATAAAg                                                                                                                                  >  1:401681/1‑148 (MQ=255)                          ggCAACTGGAAACGGCCAGCGGCAATTTATGTACGGCGACGCTGATTGCACCCAATCTGGCATTAACGGCAGGACACTGTTTATTGACACCTCCAAAGGGTAAAGCCGATAAAGCAGTGGCGCTGCGTTTTGTGTCAAATAAAGGTCtt                                                                                                                             <  4:137763/149‑1 (MQ=255)                          ggCAACTGGAAACGGCCAGCGGCAATTTATGTACGGCGACGCTGATTGCACCCAATCTGGCATTAACGGCAGGACACTGTTTATTGACACCTCCAAAGGGTAAAGCCGATAAAGCAGTGGCGCTGCGTTTTGTGTCAAATAAAGGTCtt                                                                                                                             <  1:591110/149‑1 (MQ=255)                           gCAACTGGAAACGGCCAGCGGCAATTTATGTACGGCGACGCTGATTGCACCCAATCTGGCATTAACGGCAGGACACTGTTTATTGACACCTCCAAAGGGTAAAGCCGATAAAGCAGTGGCGCTGCGTTTTGTGTCAAATAAAGGTCttt                                                                                                                            <  1:520385/149‑1 (MQ=255)                             aaCTGGAAACGGCCAGCGGCAATTTATGTACGGCGACGCTGATTGCACCCAATGTGGCATTAACGGCAGGACACTGTTTATTGACACCTCCAAAGGGTAAAGCCGATAAAGCAGTGGCGCTGCGTTTTGTGTCAAATAAAGGTCTTTgg                                                                                                                          >  2:184029/1‑149 (MQ=255)                             aaCTGGAAACGGCCAGCGGCAATTTATGTACGGCGACGCTGATTGCACCCAATCTGGCATTAACGGCAGGACACTGTTTATTGACACCTCCAAAGGGTAAAGCCGATAAAGCAGTGGCGCTGCGTTTTGTGTCAAATAAAGGTCTTTgg                                                                                                                          >  2:578264/1‑149 (MQ=255)                               cTGGAAACGGCCAGCGGCAATTTATGTACGGCGACGCTGATTGCACCCAATCTGGCATTAACGGCAGGACACTGTTTATTGACACCTCCAAAGGGTAAAGCCGATAAAGCAGTGGCGCTGCGTTTTGTGTCAAATAAAGGTCTTTGgcg                                                                                                                        >  3:270134/1‑149 (MQ=255)                               cTGGAAACGGCCAGCGGCAATTTATGTACGGCGACGCTGATTGCACCCAATCTGGCATTAACGGCAGGACACTGTTTATTGACACCTCCAAAGGGTAAAGCCGATAAAGCAGTGGCGCTGCGTTTTGTGTCAAATAAAGGTCTTTGgcg                                                                                                                        >  1:868234/1‑149 (MQ=255)                               cTGGAAACGGCCAGCGGCAATTTATGTACGGCGACGCTGATTGCACCCAATCTGGCATTAACGGCAGGACACTGTTTATTGACACCTCCAAAGGGTAAAGCCGATAAAGCAGTGGCGCTGCGTTTTGTGTCAAATAAAGGGCTTTGgcg                                                                                                                        >  1:599346/1‑149 (MQ=255)                                tGGAAACGGCCAGCGGCAATTTATGTACGGCGACGCTGATTGCACCCAATCTGGCATTAACGGCAGGACACTGTTTATTGACACCTCCAAAGGGTAAAGCCGATAAAGCAGTGGCGCTGCGTTTTGTGTCAAATAAAGGTCTTTGgcgc                                                                                                                       >  1:740071/1‑149 (MQ=255)                                         ccAGCGGCAATTTATGTACGGCGACGCTGATTGCACCCAATCTGGCATTAACGGCAGGACACTGTTTATTGACACCTCCAAAGGGTAAAGCCGATAAAGCAGTGGCGCTGCGTTTTGTGTCAAATAAAGGTCTTTGGCGCTATGAGATc                                                                                                              <  2:288385/149‑1 (MQ=255)                                              ggCAATTTATGTACGGCGACGCTGATTGCACCCAATCTGGCATTAACGGCAGGACACTGTTTATTGACACCTCCAAAGGGTAAAGCCGATAAAGCAGTGGCGCTGCGTTTTGTGTCAAATAAAGGTCTTTGGCGCTATGAGATCCACGa                                                                                                         >  1:560568/1‑149 (MQ=255)                                                  aTTTATGTACGGCGACGCTGATTGCACCCAATCTGGCATTAACGGCAGGACACTGTTTATTGACACCTCCAAAGGGTAAAGCCGATAAAGCAGTGGCGCTGCGTTTTGTGTCAAATAAAGGTCTTTGGCGCTATGAGATCCACGACATa                                                                                                     >  2:248821/1‑149 (MQ=255)                                                           cGGCGACGCTGATTGCACCCAATCTGGCATTAACGGCAGGACACTGTTTATTGACACCTCCAAAGGGTAAAGCCGATAAAGCAGTGGCGCTGCGTTTTGTGTCAAATAAAGGTCTTTGGCGCTATGAGATCCACGACATAGAAGGCcgc                                                                                            <  3:747405/149‑1 (MQ=255)                                                              cGACGCTGATTGCACCCAATCTGGCATTAACGGCAGGACACTGTTTATTGACACCTCCAAAGGGTAAAGCCGATAAAGCAGTGGCGCTGCGTTTTGTGTCAAATAAAGGTCTTTGGCGCTATGAGATCCACGACATAGAAGGCCGCGtt                                                                                         <  2:868234/149‑1 (MQ=255)                                                                aCGCTGATTGCACCCAATCTGGCATTAACGGCAGGACACTGTTTATTGACACCTCCAAAGGGTAAAGCCGATAAAGCAGTGGCGCTGCGTTTTGTGTCAAATAAAGGTCTTTGGCGCTATGAGATCCACGACATAGAAGGCCGCGTTGa                                                                                       <  2:586880/149‑1 (MQ=255)                                                                  gCTGATTGCACCCAATCTGGCATTAACGGCAGGACACTGTTTATTGACACCTCCAAAGGGTAAAGCCGATAAAGCAGTGGCGCTGCGTTTTGTGTCAAATAAAGGTCTTTGGCGCTATGAGATCCACGACATAGAAGGCCGCGTTGATc                                                                                     <  2:686879/149‑1 (MQ=255)                                                                  gCTGATTGCACCCAATCTGCCATTAACGGCAGGACACTGTTTATTGACACCTCCAAAGGGTAAAGCCGATAAAGCAGTGGCGCTGCGTTTTGTGTCAAATAAAGGCCTTTGGCGCTATGAGCTCCACGACATAGAAGGCCGAGTTGATc                                                                                     >  2:809415/1‑149 (MQ=255)                                                                    tGATTGCACCCAATCTGGCATTAACGGCAGGACACTGTTTATTGACACCTCCAAAGGGTAAAGACGATAAAGCAGTGGCGCTGCGTTTTGTGTCAAATAAAGGTCTTTGGCGCTATGAGATCCACGACATAGAAGGCCGCGTTGATCCg                                                                                   <  1:809415/149‑1 (MQ=255)                                                                           aCCCAATCTGGCATTAACGGCAGGACACTGTTTATTGACACCTCCAAAGGGTAAAGCCGATAAAGCAGTGGCGCTGCGTTTTGTGTCAAATAAAGGTCTTTGGCGCTATGAGATCCACGACATAGAAGGCCGCGTTGATCCGTCACTgg                                                                            <  2:623290/149‑1 (MQ=255)                                                                           aCCCAATCTGGCATTAACGGCAGGACACTGTTTATTGACACCTCCAAAGGGTAAAGCCGATAAAGCAGTGGCGCTGCGTTTTGTGTCAAATAAAGGTCTTTGGCGCTATGAGATCCACGACATAGAAGGCCGCGTTGATCCGACACTgg                                                                            <  1:157134/149‑1 (MQ=255)                                                                           aCCCAATCTGGCATTAACGGCAGGACACTGTTTATTGACACCTCCAAAGGGTAAAGCCGATAAAGCAGTGGCGCTGCGTTTTGTGTCAAATAAAGGTCTTTGGCGCTATGAGATCCACGACATAGAAGGCCGCGTTGATCCGACACTgg                                                                            <  4:54570/149‑1 (MQ=255)                                                                           aCCCAATCTGGCATTAACGGCAGGACACTGTTTATTGACACCTCCAAAGGGTAAAGCCGATAAAGCAGTGGCGCTGCGTTTTGTGTCAAATAAAGGTCTTTGGCGCTATGAGATCCACGACATAGAAGGCCGCGTTGATCCGACACTgg                                                                            >  4:24384/1‑149 (MQ=255)                                                                           aCCCAATCTGGCATTAACGGCAGGACACTGTTTATTGACACCTCCAAAGGGTAAAGCCGATAAAGCAGTGGCGCTGCGTTTTGTGTCAAATAAAGGTCTTTGGCGCTATGAGATCCACGACATAGAAGGCCGCGTTGATCCGACACTgg                                                                            >  4:343588/1‑149 (MQ=255)                                                                             ccAATCTGGCATTAACGGCAGGGCACTGTTTATTGACACCTCCAAAGGGTAAAGCCGATAAAGCAGTGGCGCTGCGTTTTGTGTCAAATAAAGGTCTTTGGCGCTAt                                                                                                                    <  3:179999/107‑1 (MQ=255)                                                                             ccAATCTGGCATTAACGGCAGGACACTGTTTATTGACACCTCCAAAGGGTAAAGCCGATAAAGCAGTGGCGCTGCGTTTTGTGTCAAATAAAGGTCTTTGGCGCTAt                                                                                                                    >  4:179999/1‑107 (MQ=255)                                                                               aaTCTGGCATTAACGGCAGGACACTGTTTATTGACACCTCCAAAGGGTAAAGCCGATAAAGCAGTGGCGCTGCGTTTTGTGTCAAATAAAGGTCTTTGGCGCTATGAGATCCACGACATAGAAGGCCGCGTTGATCCGACACTGGGaaa                                                                        <  1:625951/149‑1 (MQ=255)                                                                                 tCTGGCATTAACGGCAGGACACTGTTTATTGACACCTCCAGAGGGTAAAGCCGATAAAGCAGTGGCGCTGCGTTTTGTGTCAAATAAAGGTCTTTGGCGCTATGAGATCCACGACATAGAAGGCCGCGTTGATCCGACACTGGGAAAGc                                                                      <  2:124759/149‑1 (MQ=255)                                                                                  cTGGCATTAACGGCAGGACACTGTTTATTGACACCTCCAAAGGGTAAAGCCGATAAAGCAGTGGCGCTGCGTTTTGTGTCAAATAAAGGTCTTTGGCGCTATGAGATCCACGACATAGAAGGCCGCGTTGATCCGACACTGGGAAAGCg                                                                     <  3:528551/149‑1 (MQ=255)                                                                                  cTGGCATTAACGGCAGGACACTGTTTATTGACACCTCCAAAGGGTAAAGCCGATAAAGCAGTGGCGCTGCGTTTTGTGTCAAATAAAGGTCTTTGGCGCTATGAGATCCACGACATAGAAGGCCGCGTTGATCCGACACTGGGAAAGCg                                                                     <  1:675708/149‑1 (MQ=255)                                                                                       aTTAACGGCAGGACACTGTTTATTGACACCTCCAAAGGGTAAAGCCGATAAAGCAGTGGCGCTGCGTTTTGTGTCAAATAAAGGTCTTTGGCGCTATGAGATCCACGACATAGAAGGCCGCGTTGATCCGACACTGGGAAAGCGGTTaa                                                                >  4:499853/1‑149 (MQ=255)                                                                                         tAACGGCATGACACTGTTTATTGACACCTCCAAAGGGTAAAGCCGATAAAGCAGTGGCGCTGCGTTTTGTGTCAAATAAAGGTCTTTGGCGCTATGAGATCCACGACATAGAAGGCCGCGTTGATCCGACACTGGGAAAGCGGTTaaaa                                                              <  2:740071/149‑1 (MQ=255)                                                                                         tAACGGCAGGACACTGTTTATTGACACCTCCAAAGGGTAAAGCCGATAAAGCAGTGGCGCTGCGTTTTGTGTCAAATAAAGGTCTTTGGCGCTATGAGATCCACGACATAGAAGGCCGCGTTGATCCGACACTGGGAAAGCGGTTaaaa                                                              <  1:759238/149‑1 (MQ=255)                                                                                                     accgTTTATTGACACCTCCAAAGGGTAAAGCCGATAAAACAGTGGCGCTGCGTTTTATGTCAAATAAGGGTCTTTGGCTCTATGATATCCACGACATAGAAGGCCGCGGTGATCCGACACTGGGACAGCGGTTAAAAGCAGATGGGGac                                                  >  2:156861/4‑148 (MQ=255)                                                                                                     acTGTTTATTGACACCTCCAAAGGGTAAAGCCGATAAAGCAGTGGCGCTGCGTTTTGTGTCAAATAAAGGTCTTTGGCGCTATGAGATCCACGACATAGAAGGCCGCGTTGATCCGACACTGGGAAAGCGGTTAAAAGCAGATGGGGAt                                                  <  2:612031/149‑1 (MQ=255)                                                                                                       tGTTTATTGACACCTCCAAAGGGTAAAGCCGATAAAGCAGTGGCGCTGCGTTTTGTGTCAAATAAAGGTCTTTGGCGCTATGAGATCCACGACATAGAAGGCCGCGTTGATCCGACACTGGGAAAGCGGTTAAAAGCAGATGGGGATgg                                                <  2:153938/149‑1 (MQ=255)                                                                                                        gTTTATTGACACCTCCAAAGGGTAAAGCCGATAAAGCAGTGGCGCTGCGTTTTGTGTCAAATAAAGGGCTTTGGCGCTATGAGATCCACGACATAGAAGGCCGCGTTGATCCGACACTGGGAAAGCGGTTAAAAGCAGATGGGGATGGt                                               >  2:862504/1‑149 (MQ=255)                                                                                                             ttGACACCTCCAAAGGGTAAAGCCGATAAAGCAGTGGCGCTGCGTTTTGTGTCAAATAAAGGTCTTTGGCGCTATGAGATCCACGACATAGAAGGCCGCGTTGATCCGACACTGGGAAAGCGGTTAAAAGCAGATGGGGATGGTTGGAt                                          >  4:278956/1‑149 (MQ=255)                                                                                                                acacCTCCAAAGGGTAAAGCCGATAAAGCAGTGGCGCTGCGTTTTGTGTCAAATAAAGGTCTTTGGCGCTATGAGATCCACGACATAGAAGGCCGCGTTGATCCGACACTGGGAAAGCGGTTAAAAGCAGATGGGGATGGTTGGATTGt                                       >  4:782842/1‑149 (MQ=255)                                                                                                                 cacCTCCAAAGGGTAAAGCCGATAAAGCAGTGGCGCTGCGTTTTGTGTCAAATAAAGGTCTTTGGCGCTATGAGATCCACGACATAGAAGGCCGCGTTGATCCGACACTGGGAAAGCGGTTAAAAGCAGATGGGGATGGTTGGATTGTa                                      <  3:343588/149‑1 (MQ=255)                                                                                                                  acCTCCAAAGGGTAAAGCCGATAAAGCAGTGGCGCTGCGTTTTGTGTCAAATAAAGGTCTTTGGCGCTATGAGATCCACGACATAGAAGGCCGCGTTGATCCGACACTGGGAAAGCGGTTAAAAGCAGATGGGGATGGTTGGATTGTaa                                     >  3:270756/1‑148 (MQ=255)                                                                                                                       cAAAGGGTAAAGCCGATAAAGCAGTGGCGCTGCGTTTTGTGTCAAATAAAGGTCTTTGGCGCTATGAGATCCACGACATAGAAGGCCGCGTTGATCCGACACTGGGAAAGCGGTTAAAAGCAGATGGGGATGGTTGGATTGTACCTccc                                <  1:216760/149‑1 (MQ=255)                                                                                                                        aaaGGGTAAAGCCGATAAAGCAGTGGCGCTGCGTTTTGTGTCAAATAAAGGTCTTTGGCGCTATGAGATCCACGACATAGAAGGCCGCGTTGATCCGACACTGGGAAAGCGGTTAAAAGCAGATGGGGATGGTTGGATTGTACCTCCCg                               >  3:424305/1‑149 (MQ=255)                                                                                                                         aaGGGTAAAGCCGATAAAGCAGTGGCGCTGCGTTTTGTGTCAAATAAAGGTCTTTGGCGCTATGAGATCCACGACATAGAAGGCCGCGTTGATCCGACACTGGGAAAGCGGTTAAAAGCAGATGGGGATGGTTGGATTGTACCTCCCGc                              >  1:325722/1‑149 (MQ=255)                                                                                                                           gggTAAAGCCGATAAAGCAGTGGCGCTGCGTTTTGTGTCAAATAAAGGTCTTTGGCGCTATGAGATCCACGACATAGAAGGCCGCGTTGATCCGACACTGGGAAAGCGGTTAAAAGCAGATGGGGATGGTTGGATTGTACCTCCCGCAg                            <  2:262117/149‑1 (MQ=255)                                                                                                                           gggTAAAGCCGATAAAGCAGTGGCGCTGCGTTTTGTGTCAAATAAAGGTCTTTGGCGCTATGAGATCCACGACATAGAAGGCCGCGTTGATCCGACACTGGGAAAGCGGTTAAAAGCAGATGGGGATGGTTGGATTGTACCTCCCGCAg                            >  3:560914/1‑149 (MQ=255)                                                                                                                              tAAAGCCGATAAAGCAGTGGCGCTGCGTTTTGTGTCAAATAAAGGTCTTTGGCGCTATGAGATCCACGACATAGAAGGCCGCGTTGAt                                                                                      >  3:134336/1‑88 (MQ=255)                                                                                                                              tAAAGCCGATAAAGCAGTGGCGCTGCGTTTTGTGTCAAATAAAGGTCTTTGGCGCTATGAGATCCACGACATAGAAGGCCGCGTTGAt                                                                                      <  4:134336/88‑1 (MQ=255)                                                                                                                               aaaGCCGATAAAGCAGTGGCGCTGCGTTTTGTGTCAAATAAAGGTCTTTGGCGCTATGAGATCCACGACATAGAAGGCCGCGTTGATCCGACACTGGGAAAGCGGTTAAAAGCAGATGGGGATGGTTGGATTGTACCTCCCGCAGCcgc                        >  4:762920/1‑149 (MQ=255)                                                                                                                                aaGCCGATAAAGCAGTGGCGCTGCGTTTTGTGTCAAATAAAGGTCTTTGGCGCTATGAGATCCACGACATAGAAGGCCGCGTTGATCCGACACTGGGAAAGCGGTTAAAAGCAGATGGGGATGGTTGGATTGTACCTCCCGCAGCcgcg                       >  2:229837/1‑149 (MQ=255)                                                                                                                                aaGCCGATAAAGCAGTGGCGCTGCGTTTTGTGTCAAATAAAGGTCTTTGGCGCTATGAGATCCACGACATAGAAGGCCGCGTTGATCCGACACTGGGAAAGCGGTTAAAAGCAGATGGGGATGGTTGGATTGTACCTCCCGCAGCcgcg                       >  1:226421/1‑149 (MQ=255)                                                                                                                                  gCCGATAAAGCAGTGGCGCTGCGTTTTGTGTCAAATAAAGGTCTTTGGCGCTATGAGATCCACGACATAGAAGGCCGCGTTGATCCGACACTGGGAAAGCGGTTAAAAGCAGATGGGGATGGTTGGATTGTACCTCCCGCAGCCGCGcc                     >  3:296402/1‑149 (MQ=255)                                                                                                                                    cGATAAAGCAGTGGCGCTGCGTTTTGTGTCAAATAAAGGTCTTTGGCGCTATGAGATCCACGACATAGAAGGCCGCGTTGATCCGACACTGGGAAAGCGGTTAAAAGCAGATGGGGATGGTTGGATTGTACCTCCCGCAGCCGCGCCGt                   <  3:50670/149‑1 (MQ=255)                                                                                                                                     gATAAAGCAGTGGCGCTGCGTTTTGTGTCAAATAAAGGTCTTTGGCGCTATGAGATCCACGACATAGAAGGCCGCGTTGATCCGACACTGGGAAAGCGGTTAAAAGCAGATGGGGATGGTTGGATTGTACCTCCCGCAGCCGCGCCGTg                  >  3:265254/1‑149 (MQ=255)                                                                                                                                     gATAAAGCAGTGGCGCTGCGTTTTGTGTCAAATAAAGGTCTTTGGCGCTATGAGATCCACGACATAGAAGGCCGCGTTGATCCGACACTGGGAAAGCGGTTAAAAGCAGATGGGGATGGTTGGATTGTACCTCCCGCAGCCGCGCCGTg                  >  2:869978/1‑149 (MQ=255)                                                                                                                                     gATAAAGCAGTGGCGCTGCGTTTTGTGTCAAATAAAGGTCTTTGGCGCTATGAGATCCACGACATAGAAGGCCGCGTTGATCCGACACTGGGAAAGCGGTTAAAAGCAGATGGGGATGGTTGGATTGTACCTCCCGCAGCCGCGCCGTg                  >  4:73474/1‑149 (MQ=255)                                                                                                                                           gCAGTGGCGCCGCGTTTTGTGTCAAATAAAGGTCTTTGGCGCTATGAGATCCACGACATAGAAGGCCGCGTTGATCCGACACTGGGAAAGc                                                                      <  1:244968/91‑1 (MQ=255)                                                                                                                                           gCAGTGGCGCCGCGTTTTGTGTCAAATAAAGGTCTTTGGCGCTATGAGATCCACGACATAGAAGGCCGCGTTGATCCGACACTGGGAAAGc                                                                      >  4:582748/1‑91 (MQ=255)                                                                                                                                           gCAGTGGCGCCGCGTTTTGTGTCAAATAAAGGTCTTTGGCGCTATGAGATCCACGACATAGAAGGCCGCGTTGATCCGACACTGGGAAAGc                                                                      >  2:244968/1‑91 (MQ=255)                                                                                                                                           gCAGTGGCGCCGCGTTTTGTGTCAAATAAAGGTCTTTGGCGCTATGAGATCCACGACATAGAAGGCCGCGTTGATCCGACACTGGGAAAGc                                                                      <  3:582748/91‑1 (MQ=255)                                                                                                                                                   gcTGCGTTTTGTGTCAAATAAAGGTCTTTGGCGCTATGAGATCCACGACATAGAAGGCCGCGTTGATCCGACACTGGGAAAGCGGTTAAAAGCAGATGGGGATGGTTGGATTGTACCTCCCGCAGCCGCGCCGTGGGACTTCGgattga    <  1:360078/149‑1 (MQ=255)                                                                                                                                                     tGCGTTTTGTGTCAAATAAAGGTCTTTGGCGCTATGAGATCCACGACATAGAAGGCCGCGTTGATCCGACACTGGGAAAGCGGTTAAAAGCAGATGGGGATGGTTGGATTGTACCTCCCGCAGCCGCGCCGTGGGACTTCGgattgatt  <  3:90501/149‑1 (MQ=255)                                                                                                                                                     tGCGTTTTGTGTCAAATAAAGGTCTTTGGCGCTATGAGATCCACGACATAGAAGGCCGCGTTGATCCGACACTGGGAAAGCGGTTAAAAGCAGATGGGGATGGTTGGATTGTACCTCCCGCAGCCGCGCCGTGGGACTTCGgattgatt  >  4:738949/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                     ACCCAATCACCGTGGGATGCGGTTGGGCAACTGGAAACGGCCAGCGGCAATTTATGTACGGCGACGCTGATTGCACCCAATCTGGCATTAACGGCAGGACACTGTTTATTGACACCTCCAAAGGGTAAAGCCGATAAAGCAGTGGCGCTGCGTTTTGTGTCAAATAAAGGTCTTTGGCGCTATGAGATCCACGACATAGAAGGCCGCGTTGATCCGACACTGGGAAAGCGGTTAAAAGCAGATGGGGATGGTTGGATTGTACCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATT  >  NC_000913/1672110‑1672406

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