Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,508,007	T→A	5.2%	R409R (CGT→CGA)	ydcP →	putative peptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,508,007	0	T	A	5.2%	261.0 / 3.6	88	R409R (CGT→CGA)	ydcP	putative peptidase
Reads supporting (aligned to +/- strand):  ref base T (39/43);  new base A (2/3);  total (42/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.64e-01

TACGGTCGAGAAAACCGGAGAAAATCAGTACCGCGTCTGGCCCAATGAAATGCCAGCAGATTTGCACAAAATTCGTCCACATCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTCGACGAAGCCAATAACGCCGAAAAAGCAATGAACAATCTGAA  >  NC_000913/1507860‑1508153                                                                                                                                                    |                                                                                                                                                   tACGGTCGAGAAAACCGGAGAAAATCAGTACCGCGTCTGGCCCAATGAAATGCCAGCAGATTTGCACAAAATTCGTCCACATCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTc                                                                                                                                                   >  7:373084/1‑149 (MQ=255) tACGGTCGAGAAAACCGGAGAAAATCAGTACCGCGTCTGGCCCAATGAAATGCCAGCAGATTTGCACAAAATTCGTCCACATCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCGCCAGCGAACGTc                                                                                                                                                   >  3:272642/1‑149 (MQ=255)   cGGTCGAGAAAAGCGGAGTAAATCAGTACCGCGTCTGGGCCAATGAAATGCCAGCAGATTTGCACAAAATTCGTCCACATCACCCACTAAAGCGTAATCTTGATCATAACTGGGAGCAGGCACTGACAAAAACCTCCAGCGAACGTCgg                                                                                                                                                 <  7:11601/149‑1 (MQ=255)     gTCGAGAAAACCGGAGAAAATCAGTACCGCGTCTGGCCCAATGAAATGCCAGCAGATTTGCACAAAATTCGTCCACATCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGt                                                                                                                                               >  4:188695/1‑149 (MQ=255)      tCGAGAAAACCGGAGAAAATCAGTACCGCGTCTGGCCCAATGAAATGCCAGCAGATTTGCACAAAATTCGTCCACATCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTg                                                                                                                                              <  2:273801/149‑1 (MQ=255)        gagaAAACCGGAGAAAATCAGTACCGCGTCTGGCCCAATGCAATGCCAGCAGATTTGCACAAAATTCGTCCACATCCCCCCCTAACCCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAACACCCTCCAGAGCACGGCGGGTGGc                                                                                                                                            >  3:276942/1‑149 (MQ=255)        gagaAAACCGGAGAAAATCAGTACCGCGTCTGGCCCAATGAAATGCCAGCAGATTTGCACAAAATTCGTCCACATCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTgg                                                                                                                                             <  7:20096/148‑1 (MQ=255)         agaAAACCGGAGAAAATCAGTACCGCGTCTGGCCCAATGAAATGCCAGCAGATTTGCACAAAATTCGTCCACATCCCCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCg                                                                                                                                           >  8:119175/1‑149 (MQ=255)                cGGAGAAAATCAGTACCGCGTCTGGCCCAATGAAATGCCAGCAGATTTGCACAAAATTCGTCCACATCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACa                                                                                                                                    <  1:5316/149‑1 (MQ=255)                cGGAGAAAATCAGTACCGCGTCTGGCCCAATGAAATGCCAGCAGATTTGCACAAAATTCGTCCACATCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACa                                                                                                                                    <  8:316681/149‑1 (MQ=255)                  gagaAAATCAGTACCGCGTCTGGCCCAATGAAATGCCAGCAGATTTGCACAAAATTCGTCCACATCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACAtt                                                                                                                                  >  4:252955/1‑149 (MQ=255)                     aaaaTCAGTACCGCGTCTGGCCCAATGAAATGCCAGCAGATTTGCACAAAATTCGTCCACATCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGaa                                                                                                                               <  1:241971/149‑1 (MQ=255)                     aaaaTCAGTACCGCGTCTGGCCCAATGAAATGCCAGCAGATTTGCACAAAATTCGTCCACATCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGACGGGTGGCGGTAGACATTGaa                                                                                                                               <  2:358035/149‑1 (MQ=255)                      aaaTCAGTACCGCGTCTGGCCCAATGAAATGCCAGCAGATTTGCACAAAATTCGTCCACATCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAAc                                                                                                                              <  7:85430/149‑1 (MQ=255)                      aaaTCAGTACCGCGTCTGGCCCAATGAAATGCCAGCAGATTTGCAAAAAATTCGTCCACATCACCCACTATACCGTAATCTTGATCATAACTTGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCTGTAGACATTGAAc                                                                                                                              <  5:120816/149‑1 (MQ=255)                               ccGCGTCTGGCCCAATGAAATGCCAGCATATTTGCACAACATTCGTCCACATCACCCACTAAACAGTAATCTTTATCATAACTGGCAGCAGGCACTTACAAAAACCTCCAGAGAAAGTCGGGGGGCGGTATACATGTAACTGGGCGGCt                                                                                                                     >  6:354453/1‑149 (MQ=255)                               ccGCGTCTGGCCCAATGAAATGCCAGCAGATTTGCACAAAATTCGTCCACATCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCt                                                                                                                     <  8:267231/149‑1 (MQ=255)                               ccGCGTCTGGCCCAATGAAATGCCAGCAGATTTGCACAAAATTCGTCCACATCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCt                                                                                                                     >  3:212546/1‑149 (MQ=255)                               ccGCGTCTGGCCCAATGAAATGCCAGCAGATTTGCACAAAATTCGTCCACATCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCt                                                                                                                     >  4:312926/1‑149 (MQ=255)                                cgcgTCTGGCCCAATGAAATGCCAGCAGATTTGCACAAAATTCGTCCACATCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTg                                                                                                                    <  7:313402/149‑1 (MQ=255)                                cgcgTCTGGCCCAATGAAATGCCAGCAGATTTGCACAAAATTCGTCCACATCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGACGGGTGGCGGTAGACATTGAACTGGGCGGCTg                                                                                                                    >  2:58693/1‑149 (MQ=255)                                cgcgTCTGGCCCAATGAAATGCCAGCAGATTTGCACAAAATTCGTCCACATCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGACGGGTGGCGGTAGACATTGAACTGGGCGGCTg                                                                                                                    >  2:64467/1‑149 (MQ=255)                                 gcgTCTGGCCCAATGAAATGCCAGCAGATTTGCACAAAATTCGTCCACATCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTgg                                                                                                                   <  3:188695/149‑1 (MQ=255)                                  cgTCTGGCCCAATGACATGCCAGCAGATTTGCACAAAATGCGTCCACATCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGc                                                                                                                  <  5:354453/149‑1 (MQ=255)                                   gTCTGGCCCAATGAAATGCCAGCAGATTTGCACAAAATTCGTCCACATCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCa                                                                                                                 <  8:45171/149‑1 (MQ=255)                                   gTCTGGCCCAATGAAATGCCAGCAGATTTGCACAAAATTCGTCCACATCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTAGGGTGGCGGTAGACATTGAACTGGGCGGCTgg                                                                                                                   >  1:161305/1‑147 (MQ=255)                                        gCCCAATGAAATGCCAGCAGATTTGCACAAACTTCGTCCACCTCACCCACTACACCATAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGa                                                                                                              <  2:351989/147‑1 (MQ=255)                                        gCCCAATGAAATGCCAGCAGATTTGCACAAAATTCGTCCACATCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGa                                                                                                              <  8:26247/147‑1 (MQ=255)                                          ccAATGAAATGCCAGCAGATTTGCACAAAATTCGTCCACATCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGaacaa                                                                                                          >  3:330652/1‑149 (MQ=255)                                          ccAATGAAATGCCAGCAGATTTGCACAAAATTCGTCCACATCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGaacaa                                                                                                          <  1:77602/149‑1 (MQ=255)                                           cAATGAAATGCCAGCAGATTTGCACAAAATTCGTCCACATCACCCACTAAACCGTAATCTTGATCATACCTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGaacaac                                                                                                         >  8:302903/1‑149 (MQ=255)                                           cAATGAAATGCCAGCAGATTTGCACAAAATTCGTCCACATCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGaacaac                                                                                                         <  1:410422/149‑1 (MQ=255)                                              tGAAATGCCAGCAGATTTGCACAAAATTCGTCCACATCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGa                                                                                                      <  3:252955/149‑1 (MQ=255)                                              tGAAATGCCAGCAGATTTGCACAAAATTCGTCCACATCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGa                                                                                                      <  8:373084/149‑1 (MQ=255)                                                aaaTGCCAGCAGATTTGCACAAAATTCGTCCACATCACCCACTAAACCGTAATCTTGATAATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATAGAACTCGGCGGCTGGCAGGAACAACAGAt                                                                                                     >  4:113959/1‑148 (MQ=255)                                                 aaTGCCAGCAGATTTGCACAAAATTCGTCCACATCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTc                                                                                                   >  4:88807/1‑149 (MQ=255)                                                 aaTGCCAGCAGATTTGCACAAAATTCGTCCACATCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTAGCAGGAACAACTGATTc                                                                                                   <  2:299280/149‑1 (MQ=255)                                                   tGCCAGCAGATTTGCACAAAATTCGTCCACATCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGACCGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGCTTCt                                                                                                  >  6:181517/1‑148 (MQ=255)                                                     ccAGCAGATTTGCACAAAATTCGTCCACATCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGAc                                                                                               <  4:212546/149‑1 (MQ=255)                                                     ccAGCAGATTTGCACAAAATTCGTCCACATCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGAc                                                                                               <  8:53765/149‑1 (MQ=255)                                                      cagcagATTTGCACAAAATTCGTCCACATCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGAcc                                                                                              >  4:162237/1‑149 (MQ=255)                                                       agcagATTTGCACAAAATTCGTCCACATCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGAccc                                                                                             >  7:114004/1‑149 (MQ=255)                                                           gATTTGCACAAAATTCGTCCACATCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCAc                                                                                         <  3:187580/149‑1 (MQ=255)                                                            aTTTGCACAAAATTCGTCCACATCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCCAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCAc                                                                                         >  8:216076/1‑148 (MQ=255)                                                                 cacaAAATTCGTCCACATCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTg                                                                                    >  5:226282/1‑148 (MQ=255)                                                                    aaaaTTCGTCCACATCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAaga                                                                                <  8:405971/149‑1 (MQ=255)                                                                           gTCCACATCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGgtgt                                                                          >  2:13517/1‑148 (MQ=255)                                                                           gTCCACATCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGACGTGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTc                                                                         <  1:64467/149‑1 (MQ=255)                                                                           gTCCACATCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGACGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTc                                                                         <  1:58693/149‑1 (MQ=255)                                                                            tCCACATCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGgtgt                                                                          >  5:400192/1‑147 (MQ=255)                                                                              cacaTCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAg                                                                       >  7:94846/1‑148 (MQ=255)                                                                                     cccACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGc                                                               <  6:400192/149‑1 (MQ=255)                                                                                            aaCCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCt                                                        <  6:226282/149‑1 (MQ=255)                                                                                             aCCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTg                                                       >  1:382883/1‑149 (MQ=255)                                                                                              ccGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGaa                                                                                  <  5:100988/121‑1 (MQ=255)                                                                                              ccGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGaa                                                                                  >  6:100988/1‑121 (MQ=255)                                                                                              ccGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAATAGGGTGTCAGCATCACGCATACGCTg                                                       >  1:263140/1‑148 (MQ=255)                                                                                               cGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGa                                                     >  7:178532/1‑149 (MQ=255)                                                                                                  aaTCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACgg                                                  >  1:277851/1‑149 (MQ=255)                                                                                                      ttGATCATAACTGGCAGCAGGCACTGCCAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGCACTGGGCGGCGGGCAGGAACAACTGATTATGACCCTCACAAGTGAAGAGGGTGTCAGCATCACGCATACTCTGGACGTGCAg                                              >  6:329927/1‑149 (MQ=255)                                                                                                         aTCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTTGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTc                                           <  4:300321/149‑1 (MQ=255)                                                                                                          tCATAACTGGCAGCGGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTc                                           >  3:197033/1‑148 (MQ=255)                                                                                                                   gcagcaGGCACTGACAAAAACCTCCAGCGCACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTCGACGAAGcc                                  >  1:26735/1‑148 (MQ=255)                                                                                                                   gcagcaGGCACTGACAAAAACCTCCAGCGACCGTCGGTGGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTCGACGAAGc                                   <  3:113959/147‑1 (MQ=255)                                                                                                                    cagcagGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTCGACGAAGCCaa                                <  3:29253/149‑1 (MQ=255)                                                                                                                     agcagGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGGGTCAGCATCACGCATACGCTGGACGTGCAGTTCGACGAAGCCaa                                >  3:42652/1‑148 (MQ=255)                                                                                                                      gcagGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTCGACGAAGCCAAt                               >  8:121044/1‑148 (MQ=255)                                                                                                                      gcagGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTCGACGAAGCCAATa                              <  2:382883/149‑1 (MQ=255)                                                                                                                       cagGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTCGACGAAGCCAATaa                             <  3:6064/149‑1 (MQ=255)                                                                                                                        agGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTCGACGAAGCCAATAAc                            <  3:312926/149‑1 (MQ=255)                                                                                                                           cACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTCGACGAAGCCAATAACGcc                         >  1:250341/1‑149 (MQ=255)                                                                                                                           cACTGACAAAAACCTCCAGCGAACGTCGGGGGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTCGACGAAGCCAATAACGcc                         <  4:272642/149‑1 (MQ=255)                                                                                                                             atGACAAAAACCTCCAGCGAACGTAGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGAGCAGTACGACGAAGCCAAGAACGCCg                        >  8:304492/2‑148 (MQ=255)                                                                                                                              tGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTCGCCGAAGCCAATAACGCCGa                       >  4:356110/1‑148 (MQ=255)                                                                                                                                aCAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTCGACGAAGCCAATAACGCCGaaa                     >  8:259786/1‑148 (MQ=255)                                                                                                                                  aaaaaCCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTCGACGAAGCCAATAACGCCGAAAAAg                  >  5:310246/1‑149 (MQ=255)                                                                                                                                  aaaaaCCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTCGACGAAGCCAATAACGCCGAAAAAg                  <  3:88807/149‑1 (MQ=255)                                                                                                                                  aaaaaCCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCAGACGCTGGACGGGCAGTTCGACGAAGCCAAATACGCCGAAAAAg                  >  1:61309/1‑149 (MQ=255)                                                                                                                                     aaCCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTCGACGAAGCCAATAACGCCGAAAAAGCaa               <  4:330652/149‑1 (MQ=255)                                                                                                                                      aCCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTCGACGAAGCCAATAACGCCGAAAAAGCAAt              <  7:302903/149‑1 (MQ=255)                                                                                                                                      aCCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTCGACGAAGCCAATAACGCCGAAAAAGCAAt              <  8:175437/149‑1 (MQ=255)                                                                                                                                         tCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCa                                                              <  4:100465/98‑1 (MQ=255)                                                                                                                                         tCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCa                                                              >  3:100465/1‑98 (MQ=255)                                                                                                                                         tCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTCGACGAAGCCAATAACGCCGAAAAAGCAATGaa           >  2:129222/1‑149 (MQ=255)                                                                                                                                              cGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTCGACGAAGCCAATAACGCCGAAAAAGCAATGAACAATc      <  2:250341/149‑1 (MQ=255)                                                                                                                                                  cGTCGGGTGGCGGTAGACATTTAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTCGACGAAGCCAATAACGCCGAAAAAGCAATGAACAATCTGaa  <  8:282607/149‑1 (MQ=255)                                                                                                                                                  cGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTCGACGAAGCCAATAACGCCGAAAAAGCAATGAACAATCTGaa  <  7:216076/149‑1 (MQ=255)                                                                                                                                                  cGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTCGACGAAGCCAATAACGCCGAAAAAGCAATGAACAATCTGaa  <  8:94846/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                   TACGGTCGAGAAAACCGGAGAAAATCAGTACCGCGTCTGGCCCAATGAAATGCCAGCAGATTTGCACAAAATTCGTCCACATCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTCGACGAAGCCAATAACGCCGAAAAAGCAATGAACAATCTGAA  >  NC_000913/1507860‑1508153

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