Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,566,711	C→T	5.7%	L441L (CTG→CTA)	glgA ←	glycogen synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,566,711	0	C	T	5.7%	229.6 / 2.7	72	L441L (CTG→CTA)	glgA	glycogen synthase
Reads supporting (aligned to +/- strand):  ref base C (30/37);  new base T (2/2);  total (33/39)
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.16e-01

GGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACCACAGTACAAAAGCACGTCGAATAGCCCGTAACAGCGACCAGGCATTACTATCTTCAAAGACAAACCCACTGGCGACGCCATCTGCAAGGTTCTCAAGAGAACAGTCAGAAACCGTATCAGCAAGCCCACCGGTGCGCCGCACTAACGGC  >  NC_000913/3566571‑3566858                                                                                                                                             |                                                                                                                                                    ggAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCCCGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACCACAGTACaa                                                                                                                                              <  2:251773/148‑1 (MQ=255)    gCAATCATATAGGCGGTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACCACAGTACAAAAg                                                                                                                                           <  5:290330/148‑1 (MQ=255)       ttCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACCACAGTACAAAAGCAc                                                                                                                                        >  4:130710/1‑148 (MQ=255)        tCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACCACAGTACAAAAGCACg                                                                                                                                       >  5:316575/1‑148 (MQ=255)          atatAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACCACAGTACAAAAGCACGt                                                                                                                                      <  7:72992/147‑1 (MQ=255)           tataGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACCACAGTACAAAAGCACGTc                                                                                                                                     >  7:12983/1‑147 (MQ=255)            ataGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACCACAGTACAAAAg                                                                                                                                           >  4:233485/1‑140 (MQ=255)            ataGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCAATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACCACAGTACAAAAg                                                                                                                                           <  3:233485/140‑1 (MQ=255)             taGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACCACAGTACAAAAGCACGTCGaa                                                                                                                                  <  4:240933/148‑1 (MQ=255)             taGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACCACAGTACAAAAGCACGTCGa                                                                                                                                   <  8:107794/147‑1 (MQ=255)              aGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACCACAGTACAAAAGCACGTCGaa                                                                                                                                  >  1:94006/1‑147 (MQ=255)                gCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACCACAGTACAAAAGCACGTCGAATAg                                                                                                                               <  2:126069/148‑1 (MQ=255)                 cGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACCACAGTACAAAAGCACGTCGAATAGc                                                                                                                              >  2:102573/1‑148 (MQ=255)                    ttCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACcac                                                                                                                                                     <  8:5058/122‑1 (MQ=255)                    ttCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACcac                                                                                                                                                     >  7:5058/1‑122 (MQ=255)                    ttCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACCACAGTACAAAAGCACGTCGAATAGCCCg                                                                                                                           >  1:147943/1‑148 (MQ=255)                    ttCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACCACAGTACAAAAGCACGTCGAATAGCCCg                                                                                                                           >  2:69361/1‑148 (MQ=255)                       cTGAAAACTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACCACAGTACAAAAGCACGTCGAATAGCCCGTa                                                                                                                         <  1:69361/147‑1 (MQ=255)                        tGAAAACTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACCACAGTACAAAAGCACGTCGAATAGCCCGTaa                                                                                                                        >  4:236287/1‑147 (MQ=255)                        tGAAAACTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACCACAGTACAAAAGCACGTCGAATAGCCCGTaa                                                                                                                        >  3:315265/1‑147 (MQ=255)                          aaaaCTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACCACAGTACAAAAGCACGTCGAATAGCCCGTAACa                                                                                                                      <  1:93885/147‑1 (MQ=255)                           aaaCTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACCACAGTACAAAAGCACGTCGAAt                                                                                                                                 <  2:376723/135‑1 (MQ=255)                           aaaCTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACCACAGTACAAAAGCACGTCGAAt                                                                                                                                 >  1:376723/1‑135 (MQ=255)                             aCTATTTAAAGCGATAGTAAAGCTCACGTTACTACTTCGCCGCCACCTGCCAGCTAAAATCCATTGACATAGCCTGACGTTGCACACACCGCCACAGTCAAGGACGGGACCAGAGTACAAAAGCACGTCGAATCTCCCGTAACAGcca                                                                                                                  >  6:72935/1‑146 (MQ=255)                                 tttCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACCACAGTACAAAAGCACGTCGAATAGCCCGTAACAGCGAcc                                                                                                                <  6:156878/146‑1 (MQ=255)                                    cAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACCACAGTACAAAAGCACGTCGAATAGCCCGTAACAGCGACCAGGc                                                                                                            <  6:84429/147‑1 (MQ=255)                                     aaGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACCACAGTACAAAAGCACGTCGAATAGCCCGTAACAGCg                                                                                                                   <  3:167504/139‑1 (MQ=255)                                     aaGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACCACAGTACAAAAGCACGTCGAATAGCCCGTAACAGCg                                                                                                                   >  4:167504/1‑139 (MQ=255)                                      aGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACCACAGTACAAAAGCACGTCGAATAGCCCGTAACAGCGACCAGGCAtt                                                                                                         <  8:197736/148‑1 (MQ=255)                                      aGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACCACAGTACAAAAGCACGTCGAATAGCCCGTAACAGCGACCAGGCAtt                                                                                                         <  1:52108/148‑1 (MQ=255)                                       gCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACCACAGTACAAAAGCACGTCGAATAGCCCGTAACAGCGACCAGGCATTa                                                                                                        <  4:315265/148‑1 (MQ=255)                                          aTAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACCACAGTACAAAAGCACGTCGAATAGCCCGTAACAGCGACCAGGCATTACTa                                                                                                     <  4:100277/148‑1 (MQ=255)                                          aTAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGGGAAGGACGGGACCACAGTACAAAAGCACGTCGAATAGCCCGTAACAGCGACCAGGCATTACTa                                                                                                     <  6:184539/148‑1 (MQ=255)                                               aaaGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACCACAGTACAAACGCACGTCGAATAGCCCGTAACAGCGACCAGGCATTACTATCtt                                                                                                 >  3:95128/1‑147 (MQ=255)                                               aaaGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACCACAGTACAAAAGCACGTCGAATAGCCCGTAACAGCGACCAGGCATTACTATCtt                                                                                                 >  8:291283/1‑147 (MQ=255)                                               aaaGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACCACAGTACAAAAGCACGTCGAATAGCCCGTAACAGCGACCAGGCATTACTATCtt                                                                                                 >  8:231795/1‑147 (MQ=255)                                                      aCGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACCATAg                                                                                                                                                   >  5:318085/1‑90 (MQ=255)                                                      aCGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACCATAg                                                                                                                                                   >  5:133185/1‑90 (MQ=255)                                                      aCGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACCATAg                                                                                                                                                   <  6:318085/90‑1 (MQ=255)                                                      aCGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACCATAg                                                                                                                                                   <  6:133185/90‑1 (MQ=255)                                                         gTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACCACAGTACAAAAGCACGTCGAATAGCCCGTAACAGCGACCAGGCATTACTATCTTCAAAGACAAAc                                                                                      <  5:287630/148‑1 (MQ=255)                                                           aCGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACCACAGTACAAAAGCACGTCGAATAGCCCGTAACAGCGACCATGCATTACTATCTTCAAAGCCAAAcc                                                                                     >  3:27958/1‑147 (MQ=255)                                                           aCGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACCACAGTACAAAAGCACGTCGAATAGCCCGTAACAGCGACCAGGCATTACTATCTTCAAAGACAAAcc                                                                                     >  7:129503/1‑147 (MQ=255)                                                            cGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACCACAGTACAAAAGCACGTCGAATAGCCCGTAACAGCGACCAGGCATTACTATCTTCAAAGACAAACCCa                                                                                   <  8:129503/148‑1 (MQ=255)                                                               cTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACCACAGTACAAAAGCACGTCGAATAGCCCGTAACAGCGACCAGGCATTACTATCTTCAAAGACAAACCCACTg                                                                                >  2:91508/1‑148 (MQ=255)                                                                   gccgcGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACCACAGTACAAAAGCACGTCGAATAGCCCGTAACAGCGACCAGGCATTACTATCTTCAAAGACAAACCCACTGGCg                                                                             <  2:147943/147‑1 (MQ=255)                                                                      gcgACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACCACAGTACAAAAGCACGTCGAATAGCCCGTAACAGCGACCAGGCATTACTATCTTCAAAGACAAACCCACTGGCGACGc                                                                         <  8:89141/148‑1 (MQ=255)                                                                        gACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACCACAGTACAAAAGCACGGCGAATAGCCCGTAACAGCGACCAGGCATTACTATCTTCAAAGACAAACCCACTGGCGACGCCa                                                                       <  6:316575/148‑1 (MQ=255)                                                                         aCCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACCACAGTACAAAAGCACGTCGAATAGCCCGTAACAGCGACCAGGCATTACTATCTTCAAAGACAAACCCACTGGCGACGCCAt                                                                      <  3:236287/148‑1 (MQ=255)                                                                                aGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACCACAGTACAAAAGCACGTCGAATAGCCCGTAACAGCGACCAGGCATTACTATCTTCAAAGACAAACCCACTGGCGACGCCATCTGCAAg                                                               <  1:102573/148‑1 (MQ=255)                                                                                    aaaaTCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACCACAGTACAAAAGCACGTCGAATAGCCCGTAACAGCGACCAGGCATTACTATCTTCAAAGACAAACCCACTGGCGACGCCATCTGCAAGGTtc                                                           >  4:225659/1‑148 (MQ=255)                                                                                     aaaTCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACCACAGTACAAAAGCACGTCGAATAGCCCGTAACAGCGACCAGGCATTACTATCTTCAAAGACAAACCCACTGGCGACGCCATCTGCAAGGTtct                                                          >  3:137766/1‑148 (MQ=255)                                                                                      aaTCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACCACAGTACAAAAGCACGTCGAATAGCCCGTAACAGCGACCAGGCATTACTATCTTCAAAGACAAACCCACTGGCGACGCCATCTGCAAGGTtct                                                          >  1:100434/1‑147 (MQ=255)                                                                                         ccATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACCACAGTACAAAAGCACGTCGAATAGCCCGTAACAGCGACCAGGCATTACTATCTTCAAAGACAAACCCACTGGCGACGCCATc                                                                     >  5:240512/1‑133 (MQ=255)                                                                                         ccATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACCACAGTACAAAAGCACGTCGAATAGCCCGTAACAGCGACCAGGCATTACTATCTTCAAAGACAAACCCACTGGCGACGCCATc                                                                     <  6:240512/133‑1 (MQ=255)                                                                                         ccATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACCACAGTACAAAAGCACGTCGAATAGCCCGTAACAGCGACCAGGCATTACTATCTTCAAAGACAAACCCACTGGCGACGCCATCTGCAAGGTtctc                                                         >  4:253833/1‑145 (MQ=255)                                                                                            ttGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACCACAGTACAAAAGCACGTCGAATAGCCCGTAACAGCGACCAGGCATTACTATCTTCAAAGACAAACCCACTGGCGACGCCATCTGCAAGGTTCTCAAgaga                                                   <  8:12983/148‑1 (MQ=255)                                                                                              gCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACCACAGTACAAAAGCACGTCGAATAGCCCGTAACAGCGACCAGGCATTACTATCTTCAAAGACAAACCCACTGGCGACGCCATCTGCAAGGTTCTCAAGAGaa                                                  <  6:35769/147‑1 (MQ=255)                                                                                                cATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACCACAGTACAAAAGCACGTCGAATAGCCCGTAACAGCGACCAGGCATTACTATCTTCAAAGACAAACCCACTGGCGACGCCATCTGCAAGGTTCTCAAGAGAACAg                                               <  3:130710/148‑1 (MQ=255)                                                                                                  tAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACCACAGTACAAAAGCACGTCGAATATCCCGTAACAGCGACCAGGCATTACTATCTTCAAAGACAAACCCACTGGCGACGCCATCTGCAAGGTTCTCAAGAGAACAGTc                                             <  3:173019/148‑1 (MQ=255)                                                                                                      cTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACCACAGTACAAAAGCACGTCGAATAGCCCGTAACAGCGACCAGGCATTACTATCTTCAAAGACAAACCCACTGGCGACGCCATCTGCAAGGTTCTCAAGAGAACAGTCAGa                                          >  6:319050/1‑147 (MQ=255)                                                                                                         aCGTTGCACAAACCGCCACAGTGAAGGACGGGACCACAGTACAAAAGCACGTCGAATAGCCCg                                                                                                                           <  4:284259/63‑1 (MQ=255)                                                                                                         aCGTTGCACAAACCGCCACAGTGAAGGACGGGACCACAGTACAAAAGCACGTCGAATAGCCCg                                                                                                                           >  3:284259/1‑63 (MQ=255)                                                                                                             tGCACAAACCGCCACAGTGAAGGACGGGACCACAGTACAAAAGCACGTCGAATAGCCCGTAACAGCGACCAGGCATTACTATCTTCAAAGACAAACCCACTGGCGACGCCATCTGCAAGGTTCTCAAGAGAACAGTCAGAAACCGTAt                                  >  5:307604/1‑148 (MQ=255)                                                                                                                 caAACCGCCACAGTGAAGGACGGGACCACAGTACAAAAGCACGTCGAATAGCCCGTAACAGCGACCAGGCATTACTATCTTCAAAGACAAACCCACTGGCGACGCCATCTGCAAGGTTCTCAAGAGAACAGTCAGAAACCGTATCAGc                              <  7:291283/148‑1 (MQ=255)                                                                                                                  aaaCCGCCACAGTGAAGGACGGGACCACAGTACAAAAGCACGTCGAATAGCCCGTAACAGCGACCAGGCATTACTATCTTCAAAg                                                                                            >  8:121964/1‑85 (MQ=255)                                                                                                                  aaaCCGCCACAGTGAAGGACGGGACCACAGTACAAAAGCACGTCGAATAGCCCGTAACAGCGACCAGGCATTACTATCTTCAAAg                                                                                            <  7:121964/85‑1 (MQ=255)                                                                                                                   aaCCGCCACAGTGAAGGACGGGACCACAGTACAAAAGCACGTCGAATAGCCCGTAACAGCGACCAGGCATTACTATCTTCAAAGACAAACCCACTGGCGACGCCATCTGCAAGGTTCTCAAGAGAACAGTCAGAAACCGTATCAGc                              >  8:197028/1‑146 (MQ=255)                                                                                                                   aaCCCCCACAGTGAAGGACGGGACCACAGTACAAAAGCACGTCGAATAGCCCGTAACAGCGACCAGGCATTACTATCTTCAAAGACAAACCCACTGGCGACGCCATCTGCAAGGTTCTCAAGAGCACAGTCAGAAACCGTATCAGc                              >  2:352692/1‑146 (MQ=255)                                                                                                                    aCCGCCACAGTGAAGGACGGGACCACAGTACAAAAGCACGTCGAATAGCCCGTAACAGCGACCAGGCATTACTATCTTCAAAGACAAACCCACTGGCGACGCCATCTGCAAGGTTCTCAAGAGAACAGTCAGAAACCGTATCAGCAAg                           <  8:210638/148‑1 (MQ=255)                                                                                                                    aCCGCCACAGTGAAGGACGGGACCACAGTACAAAAGCACGTCGAATAGCCCGTAACAGCGACCAGGCATTACTATCTTCAAAGACAAACCCACTGGCGACGCCATCTGCAAGGTTCTCAAGAGAACAGTCAGAAACCGTATCAGCAAg                           <  8:216008/148‑1 (MQ=255)                                                                                                                                    aCGGGACCACAGTACAAAAGCACGTCGAATAGCCCGTAACAGCGACCAGGCATTACTATCTTCAAAGACAAACCCACTGGCGACGCCATCTGCAAGGTTCTCAAGAGAACAGTCAGAAACCGTATCAGCAAGCCCACCGGTGcgccgc           <  2:139233/148‑1 (MQ=255)                                                                                                                                        gACCACAGTACAAAAGCACGTCGAATAGCCCGTAACAGCGCCCAGGCATTACTATCTTCAAAGACAAACCCACTGGCGACGCCATCTGCAAGGTTCTCAAGAGAACAGTCAGAAACCGTATCAGCAAGCCCACCGGTGCGCCGCACt        >  3:250332/1‑147 (MQ=255)                                                                                                                                         aCCACAGTACAAAAGCACGTCGAATAGCCCGTAACAGCGACCAGGCATTACTATCTTCAAAGACAAACCCACTGGCGACGCCATCTGCAAGGTTCTCAAGAGAACAGTCAGAAACCGTATCAGCAAGCCCACCGGTGCGCCGCACTaa      <  6:307604/148‑1 (MQ=255)                                                                                                                                            acaGTACAAAAGCACGTCGAATAGCCCGAAACAGCGACCAGGCATTACTATCATCAAAGACAAACCCACTGGCGACGCCATCTGCAAGGTTCTCAAGAGAACAGTCAGAAACCGTATCAGCAAGCCCACCCGTGCCCCGCACTAAccg   >  5:199272/1‑146 (MQ=255)                                                                                                                                             caGTACAAAAGCACGTCGAATAGCCCGTAACAGCGACCAGGCATTACTATCTTCAAAGACAAACCCACTGGCGACGCCATCTGCAAGGTTCTCAAGAGAACAGTCAGAAACCGTATCAGCAAGCCCACCGGTGCGCCGCACTAACGGc  <  7:197028/148‑1 (MQ=255)                                                                                                                                             |                                                                                                                                                    GGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACCACAGTACAAAAGCACGTCGAATAGCCCGTAACAGCGACCAGGCATTACTATCTTCAAAGACAAACCCACTGGCGACGCCATCTGCAAGGTTCTCAAGAGAACAGTCAGAAACCGTATCAGCAAGCCCACCGGTGCGCCGCACTAACGGC  >  NC_000913/3566571‑3566858

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: ‑