Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,566,786	G→A	6.0%	G416G (GGC→GGT)	glgA ←	glycogen synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,566,786	0	G	A	6.0%	211.9 / 2.5	67	G416G (GGC→GGT)	glgA	glycogen synthase
Reads supporting (aligned to +/- strand):  ref base G (31/32);  new base A (2/2);  total (33/34)
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.86e-01

CGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACCACAGTACAAAAGCACGTCGAATAGCCCGTAACAGCGACCAGGCATTACTATCTTCAAAGACAAACCCACTGGCGACGCCATCTGCAAGGTTCTCAAGAGAACAGTCAGAAACCGTATCAGCAAGCCCACCGGTGCGCCGCACTAACGGCAGCGTACCGTACTTCAATCCATAAAGTTGCGTTAAGCCGCACGGTTCAAAACGGCTGGGCACCAGAATGACGT  >  NC_000913/3566639‑3566931                                                                                                                                                    |                                                                                                                                                  cgcgACCTGCCAGCTAAAATCCATTGCCATAGCCTGCCGTTGCACAAACCGCCACAGTGAAGGACGGGACCACAGTACAAAAGCACGTCGAATAGCCCGTAACAGCGACCAGGCATTACTATCTTCAAAGACAAACCCACTGGCGACGc                                                                                                                                                  >  1:56318/1‑149 (MQ=255)    gACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACCACAGTACAAAAGCACGTCGAATAGCCCGTAACAGCGACCAGGCATTACTATCTTCAAAGACAAACCCACTGGCGACGCCAt                                                                                                                                               <  3:372691/149‑1 (MQ=255)          ccAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACCACAGTACAAAAGCACGTCGAATAGCCCGTAACAGCGACCAGGCATTACTATCTTCAAAGACAAACCCACTGGCGACGCCATCTGc                                                                                                                                           >  7:224072/1‑147 (MQ=255)            aGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACCACAGTACAAAAGCACGTCGAATAGCCCGTAACAGCGACCAGGCATTACTATCTTCAAAGACAAACCCACTGGCGACGCCATCTGCAAgg                                                                                                                                       <  7:357875/149‑1 (MQ=255)            aGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACCACAGTACAAAAGCACGTCGAATAGCCCGTAACAGCGACCAGGCATTACTATCTTCAAAGACAAACCCACTGGCGACACCATCTGCAAgg                                                                                                                                       >  1:347277/1‑149 (MQ=255)            aGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACCACAGTACAAAAGCACGTCGAATAGCCCGTAACAGCGACCAGGCATTACTATCTTCAAAGACAAACCCACTGGCGACACCATCTGCAAgc                                                                                                                                       >  1:341787/1‑148 (MQ=255)             gCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACCACAGTACAAAAGCACGTCGAATAGCCCGTAACAGCGACCAGGCATTACTATCTTCAAAGACAAACCCACTGGCGACGCCATCTGCAAGGt                                                                                                                                      >  8:129874/1‑149 (MQ=255)                 aaaTCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACCACAGTACAAAAGAACGTCGAATAGCACGTAACAGCGACCAGGCATTACTATCTTCAAAGACAAACCCACTGGCGACGCCATCTGCAAGGTtct                                                                                                                                   >  1:269924/1‑148 (MQ=255)                  aaTCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACCACAGTACAAAAGCACGTCGAATAGCCCGTAACAGCGACCAGGCATTACTATCTTCAAAGACAAACCCACTGGCGACGCCAt                                                                                                                                               >  5:45916/1‑135 (MQ=255)                  aaTCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACCACAGTACAAAAGCACGTCGAATAGCCCGTAACAGCGACCAGGCATTACTATCTTCAAAGACAAACCCACTGGCGACGCCAt                                                                                                                                               <  6:45916/135‑1 (MQ=255)                  aaTCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACCACAGTACAAAAGCACGTCGAATAGCCCGTAACAGCGACCAGGCATTACTATCTTCAAAGACAAACCCACTGGCGACGCCATCTGCAAGGTtctc                                                                                                                                  >  4:177224/1‑148 (MQ=255)                      cATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACCACAGTACAAAAGCACGTCGAATAGCCCGTAACAGCGACCAGGCATTACTATCTTCAAAGACAAACCCACTGGCGACGCCATCTGCAAGGTTCTCAaga                                                                                                                              <  3:154210/148‑1 (MQ=255)                        ttGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACCACAGTACAAAAGCACGTCGAATAGCCCGTAACAGCGACCAGGCATTACTATCTTCAAAGACAAACCCACTGGCGACGCCATCTGCAAGGTTCTCAAGAGaa                                                                                                                           <  4:250110/149‑1 (MQ=255)                        ttGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACCACAGTACAAAAGCACGTCGAATAGCCCGTAACAGCGACCAGGCATTACTATCTTCAAAGACAAACCCACTGGCGACGCCATCTGCAAGGTTCTCAAGAGaa                                                                                                                           <  7:372790/149‑1 (MQ=255)                         tggcATAGCCTGTCGTTGAACAAACCGCCACAGTGAAGGAGGGCACCACAGTACAACAGCACGTCGATTAGCCCGTAAGATCGACCAGGCATTACTATCTTCAAAGACAAACCCACTGGCGACGCGATCTGGAAGTTTCTCAAGAGaa                                                                                                                           <  6:273454/145‑1 (MQ=255)                            cATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACCACAGTACAAAAGCACGTCGAATAGCCCGTAACAGCGACCAGGCATTACTATCTTCAAAGACAAACCCACTGGCGACGCCATCTGCAAGGTTCTCAAGAGAACAGt                                                                                                                       <  2:5466/149‑1 (MQ=255)                                     aCGTTGCACAAACCGCCACAGTGAAGGACGGGACCACAGTACAAAAGCACGTCGAATAGCCCGTAACAGCGACCAGGCATTACTATCTTCAAAGACAAACCCACTGGCGACGCCATCTGCAAGGTTCTCAAGAGAACAGTCAGAAAcc                                                                                                               <  1:188475/148‑1 (MQ=255)                                        ttGCACAAACCGCCACAGTGAAGGACGGGACCACAGTACAAAAGCACGTCGAATAGCCCGTAACAGCGACCAGGCATTACTATCTTCAAAGACAAACCCACTGGCGACGCCATCTGCAAGGTTCTCAAGAGAACAGTCAGAAACCGTAt                                                                                                           <  4:2383/149‑1 (MQ=255)                                          gCACAAACCGCCACAGTGAAGGACGGGACCACTGTACAAAAGCACGTCGAATAGCCCGTAACAGCGACCAGGCATTACTCTCTTCAAAGACAAACCCACTGGCGACGCCATCTGCAAGGTTCTCAAGAGAACAGTCAGAAACCGTATc                                                                                                          <  4:65279/148‑1 (MQ=255)                                            acaAACCGCCACAGTGAAGGACGGGACCACAGTACAAAAGCACGTCGAATAGCCCGTAACAGCGACCAGGCATTACTATCTTCAAAGACAAACCCACTGGCGACGCCATCTGCAAGGTTCTCAAGAGAACAGTCAGAAACCGTATCAg                                                                                                        >  8:136490/1‑148 (MQ=255)                                            acaAACCGCCACAGTGAAGGACGGGACCACAGTACAAAAGCAAGTAGAATAGCCAGTAACAGCGAAAAGGAATAACTATCTTCAAAGACAAAACCACTGGCGACGCCATCTGCAAGGTTCTCAAGAGAACAGTCAGAAACCGTATCAg                                                                                                        >  2:279378/1‑148 (MQ=255)                                             caAACCGCCACAGTGAAGGACGGGACCACAGTACAAAAGCACGTCGAATAGCCCGTAACAGCGACCAGGCATTACTATCTTCAAAGACAAACCCACTGGCGACGCCATCTGCAAGGTTCTCAAGAGAACAGTCAGAAACCGTATCAGc                                                                                                       >  8:277451/1‑148 (MQ=255)                                              aaaCCGCCACAGTGAAGGACGGGACCACAGTACAAAAGCACGTCGAATAGCCCGTAACAGCGACCAGGCATTACTATCTTCAAAGACAAACCCACTGGCGACGCCATCTGCAAGGTTCTCAAGAGAACAGTCAGAAACCGTATCAGc                                                                                                       >  7:385089/1‑147 (MQ=255)                                                aCCGCCACAGTGAAGGACGGGACCACAGTACAAAAGCACGTCGAATAGCCCGTAACAGCGACCAGGCATTACTATCTTCAAAGACAAACCCACTGGCGACGCCATCTGCAAGGTTCTCAAGAGAACAGTCAGAAACCGTATCAGCAAGc                                                                                                   <  7:58760/149‑1 (MQ=255)                                                   gCCACAGTGAAGGACGGGACCACAGTACAAAAGCACGTCGAATAGCCCGTAACAGCGACCAGGCATTACTATCTTCAAAGACAAACCCACTGGCGACGCCATCTGCAAGGTTCTCAAGAGAACAGTCAGAAACCGTATCAGCAAGCCCa                                                                                                >  1:300427/1‑149 (MQ=255)                                                      acaGTGAAGGACGGGACCACAGTACAAAAGCACGTCGAATAGCCCGTAACAGCGACCAGGCATTACTATCTTCAAAGACAAACCCACTGGCGACGCCATCTGCAAGGTTCTCAAGAGAACAGTCAGAAAcc                                                                                                               <  2:106108/131‑1 (MQ=255)                                                      acaGTGAAGGACGGGACCACAGTACAAAAGCACGTCGAATAGCCCGTAACAGCGACCAGGCATTACTATCTTCAAAGACAAACCCACTGGCGACGCCATCTGCAAGGTTCTCAAGAGAACAGTCAGAAAcc                                                                                                               >  1:106108/1‑131 (MQ=255)                                                        aGTGAAGGACGGGACCACAGTACAAAAGCACGTCGAATAGCCCGTAACAGCGACCAGGCATTACTATCTTCAAAGACAAACCCACTGGCGACGCCATCTGCAAGGTTCTCAAGAGAACAGTCAGAAACCGTATCAGCAAGCCCACCGGt                                                                                           >  1:78360/1‑149 (MQ=255)                                                             aGGACGGGACCACAGTACAAAAGCACGTCGAATAGCCCGTAACAGCGACCAGGCATTACTATCTTCAAAGACAAACCCACTGGCGACGCCATCTGCAAGGTTCTCAAGAGAACAGTCAGAAACCGTATCAGCAAGCCCACCGGTGcgcc                                                                                      <  8:224072/149‑1 (MQ=255)                                                                  gggACCACAGTACAAAAGCACGGCGAATAGCCCGTAACAGCGACCAGGCATTACTCTCATCAAAGACAAACCCACTGGCGACGCCATCTGCAAGGTTCTCAAGAGAACAGTCAGAAACCGTATCAGCAAGCCCACCGGTGCGCCGCACt                                                                                 <  8:350454/149‑1 (MQ=255)                                                                       cacaGTACAAAAGCACGTCGAATAGCCCGTAACAGCGACCAGGCATTACTATCTTCAAAGACAAACCCACTGGCGACGCCATCTGCAAGGTTCTCAAGAGAACAGTCAGAAACCGTATCAGCAAGCCCACCGGTGCGCCGCACTAACgg                                                                            <  2:56318/149‑1 (MQ=255)                                                                        acaGTACAAAAGCACGTCGAATAGCCCGTAACAGCGACCAGGCATTACTATCTTCAAAGACAAACCCACTGGCGACGCCATCTGCAAGGTTCTCAAGAGAACAGTCAGAAACCGTATCAGCAAGCCCACCGGTGCGCCGCACTAACGGc                                                                           <  3:323507/149‑1 (MQ=255)                                                                         caGTACAAAAGCACGTCGAATAGCCCGTAACAGCGACCAGGCATTACTATCTTCAAAGACAAACCCACTGGCGACGCCATCTGCAAGGTTCTCAAGAGAACAGTCAGAAACCGTATCAGCAAGCCCACCGGTGCGCCGCACTAACgg                                                                            >  8:395580/1‑147 (MQ=255)                                                                              cAAAAGCACGTCGAATAGCCCGTAACAGCGACCAGGCATTACTATCTTCAAAGACAAACCCACTGGCGACGCCATCTGCAAGGTTCTCAAGAGAACAGTCAGAAACCGTATCAGCAAGCCCACCGGTGCGCCGCACTAACGGCAGcgta                                                                     <  6:365670/149‑1 (MQ=255)                                                                               aaaaGCACGTCGAATAGCCCGTAACAGCGACCAGGCATTACTATCTTCAAAGACAAACCCACTGGCGCCGCCATCTGCAAGGTTCTCAAGAGAACAGTCAGAAACCGTATCAGCAAGCCCACCGGTGCGCCGCACTAACGGCAGcgtac                                                                    >  5:86605/1‑149 (MQ=255)                                                                                aaaGCACGTCGAATAGCCCGTAACAGCGACCAGGCATTACTATCTTCAAAGACAAACCCACTGGCGACGCCATCTGCAAGGTTCTCAAGAGAACAGTCAGAAACCGTATCAGCAAGCCCACCGGTGCGCCGCACTAACGGCAGcgtacc                                                                   <  7:339203/149‑1 (MQ=255)                                                                                 aaGCACGTCGAATAGCCCGTAACAGCGACCAGGCATTACTATCTTCAAAGACAAACCCACTGGCGACGCCATCTGCAAGGTTCTCAAGAGAACAGTCAGAAACCGTATCAGCAAGCCCACCGGTGCGCCGCACTAACGGCAGcgtaccg                                                                  >  8:70245/1‑149 (MQ=255)                                                                                     aCGTCGAATAGCCCGTAACAGCGACCAGGCATTACTATCTTCAAAGACAAACCCACTGGCGACGCCATCTGCAAGGTTCTCAAGAGAACAGTCAGAAACCGTATCAGCAAGCCCACCGGTGCGCCGCACTAACGGCAGCGTACCGTACt                                                              >  3:145695/1‑149 (MQ=255)                                                                                       gTCGAATAGCCCGTAACAGCGACCAGGCATTACTATCTTCAAAGACAAACCCACTGGCGACGCCATCTGCAAGGTTCTCAAGAGAACAGTCAGAAACCGTATCAGCAAGCCCACCGGTGCGCCGCACTAACGGCAGCGTACCGTACtt                                                             >  1:73652/1‑148 (MQ=255)                                                                                          gAATAGCCCGTAACAGCGACCAGGCATTACTATCTTCAAAGACAAACCCACTGGCGACACCATCTGCAAGGTTCTCAAGAGAACAGTCAGAAACCGTATCAGCAAGCCCACCGGTGCGCCGCACTAACGGCAGCGTACCGTACTTCAAt                                                         <  2:341787/149‑1 (MQ=255)                                                                                          gAATAGCCCGTAACAGCGACCAGGCATTACTATCTTCAAAGACAAACCCACTGGCGACACCATCTGCAAGGTTCTCAAGAGAACAGTCAGAAACCGTATCAGCAAGCCCACCGGTGCGCCGCACTAACGGCAGCGTACCGTACTTCAAt                                                         <  2:347277/149‑1 (MQ=255)                                                                                                              cAGGCATTACTATCTTCAAAGACAAACCCACTGGCGACGCCATCTGCAAGGTTCTCAAGAGAACAGTCAGAAACCGTATCAGCAAGCCCACCGGTGCGCCGCACTAACGGCAGCGTACCGTACTTCAATCCATAAAGTTGCGTTAAGcc                                     >  8:232316/1‑149 (MQ=255)                                                                                                                  cATTACTATCTTCAAAGACAAACCCACTGGCGACGCCATCTGCAAGGTTCTCAAGAGAACAGTCAGAAACCGTATCAGCAAGCCCACCGGTGCGCCGCACTAACGGCAGCGTACCGTACTTCAATCCATAAAGTTGCGTTAAGCCGCAc                                 >  4:4694/1‑149 (MQ=255)                                                                                                                  cATTACTATCTTCAAAGACAAACCCACTGGCGACGCCATCTGCAAGGTTCTCAAGAGAACAGTCAGAAACCGTATCAGCAAGCCCACCGGTGCGCCGCACTAACGGCAGCGTACCGTACTTCAATCCATAAAGTTGCGTTAAGCCGCAc                                 <  2:78360/149‑1 (MQ=255)                                                                                                                   aTTACTATCTTCAAAGACAAACCCACTGGCGACGCCATCTGCAAGGTTCTCAAGAGAACAGTCAGAAACCGTATCAGCAAGCCCACCGGTGCGCCGCACTAACGGCAGCGTACCGTACTTCAATCCATAAAGTTGCGTTAAGCCGCACg                                >  4:339227/1‑149 (MQ=255)                                                                                                                      aCTATCTTCAAAGACAAACCCACTGGCGACGCCATCTGCAAGGTTCTCAAGAGAACAGTCAGAAACCGTATCAGCAAGCCCACCGGTGCGCCGCACTAACGGCAGCGTACCGTACTTCAATCCAt                                                     >  1:170512/1‑125 (MQ=255)                                                                                                                      aCTATCTTCAAAGACAAACCCACTGGCGACGCCATCTGCAAGGTTCTCAAGAGAACAGTCAGAAACCGTATCAGCAAGCCCACCGGTGCGCCGCACTAACGGCAGCGTACCGTACTTCAATCCAt                                                     <  2:170512/125‑1 (MQ=255)                                                                                                                        tATCTTCAAAGACAAACCCACTGGCGACGCCATCTGCAAGGTTCTCAAGAGAACAGTCAGAAACCGTATCAGCAAGCCCACCGGTGCGCCGCACTAACGGCAGCGTACCGTACTTCAATCCATaaa                                                  <  5:130153/126‑1 (MQ=255)                                                                                                                        tATCTTCAAAGACAAACCCACTGGCGACGCCATCTGCAAGGTTCTCAAGAGAACAGTCAGAAACCGTATCAGCAAGCCCACCGGTGCGCCGCACTAACGGCAGCGTACCGTACTTCAATCCATaaa                                                  >  6:130153/1‑126 (MQ=255)                                                                                                                         aTCTTCAAAGACAAACCCACTGGCGACGCCATCTGCAAGGTTCTCAAGAGAACAGTCAGAAACCGTATCAGCAAGCCCACCGGTGCGCCGCACTAACGGCAGCGTACCGTACTTCAATCCAt                                                     >  2:166152/1‑122 (MQ=255)                                                                                                                            ttCAAAGACAAACCCACTGGCGACGCCATCTGCAAGGTTCTCAAGAGAACAGTCAGAAACCGTATCAGCAAGCCCACCGGTGCGCCGCACTAACGGCAGCGTACCGTACTTCCATCCATAAAGTTGCGTTAAGCCGCACGGTTCAAAAc                       >  5:131245/1‑149 (MQ=255)                                                                                                                            ttCAAAGACAAACCCACTGGCGACGCCATCTGCAAGGTTCTCAAGAGAACAGTCAGAAACCGTATCAGCAAGCCCACCGGTGCGCCGCACTAACGGCAGCGTACCGTACTTCAATCCATAAAGTTGCGTTAAGCCGCACGGTTCAAAAc                       >  1:317708/1‑149 (MQ=255)                                                                                                                              cAAAGACAAACCCACTGGCGACGCCATCTGCAAGGTTCTCAAGAGAACAGTCAGAAACCGTATCAGCAAGCCCACCGGTGCGCCGCACTAACGTCAGCGTACCGTACTTCAATCCATAAAGTTGCGTTAAGCCGCACGGTTCaaaa                        <  6:126225/146‑1 (MQ=255)                                                                                                                                 aGACAAACCCACTGGCGACGCCATCTGCAAGGTTCTCAAGAGAACAGTCAGAAACCGTATCAGCAAGCCCACCGGTGCGCCGCACTAACGGCAGCGTACCGTACTTCAATCCATAAAGTTGCGTTAAGCCGCACGGTTCAAAACGGCTg                  >  2:160593/1‑149 (MQ=255)                                                                                                                                  gACAAACCCACTGGCGACGCCATCTGCAAGGTTCTCAAGAGAACAGTCAGAAACCGTATCAGCAAGCCCACCGGTGCGCCGCACTAACGGCAGCGTACCGTACTTCAATCCATAAAGTTGCGTTAAGCCGCACGGGTCAAAACGGCTgg                 >  6:21891/1‑149 (MQ=255)                                                                                                                                   aCAAACCCACTGGCGACGCCATCTGCAAGGTTCTCAAGAGAACAGTCAGAAACCGTATCAGCAAGCCCACCGGTGCGCCGCACTAACGGCAGCGTACCGTACTTCAATCCATAAAGTTGCGTTAAGCCGCACGGTTCAAAACGGCTggg                <  7:301481/149‑1 (MQ=255)                                                                                                                                      aaCCCACTGGCGACGCCATCTGCAAGGTTCTCAAGAGAACAGTCAGAAACCGTATCAGCAAGCCCACCGGTGCGCCGCACTAACGGCAGCGTACCGTACTTCAATCCATAAAGTTGCGTTAAGCCGCACGGTTCAAAACGGCTGGGCAc             <  2:317708/149‑1 (MQ=255)                                                                                                                                        cccACTGGCGACGCCATCTGCAAGGTTCTCAAGAGAACAGTCAGAAACCGTATCAGCAAGCCCACCGGTGCGCCGCACTAACGGCAGc                                                                        <  2:370374/88‑1 (MQ=255)                                                                                                                                        cccACTGGCGACGCCATCTGCAAGGTTCTCAAGAGAACAGTCAGAAACCGTATCAGCAAGCCCACCGGTGCGCCGCACTAACGGCAGc                                                                        >  1:370374/1‑88 (MQ=255)                                                                                                                                        cccACTGGCGACGCCATCTGCAAGGTTCTCAAGAGAACAGTCAGAAACCGTATCAGCAAGCCCACCGGTGCGCCGCACTAACGGCAGCGTACCGTACTTCAATCCATAAAGTTGCGTTAAGCCGCACGGTTCAAAACGGCTGGGCACCa           <  3:4694/149‑1 (MQ=255)                                                                                                                                          cACTGGCGACGCCATCTGCAAGGTTCTCAAGAGAACAGTCAGAAACCGTATCAGCAAGCCCACCGGTGCGCCGCACTAACGGCAGCGTACCGTACTTCAATCCATAAAGTTGCGTTAAGCCGCACGGTTCAAAACGGCTGGGCACCAg          >  1:177261/1‑148 (MQ=255)                                                                                                                                            cTGGCGACGCCATCTGCAAGGTTCTCAAGAGAACAGTCAGAAACCGTATCAGCAAGCCCACCGGTGCGCCGCACTAACGGCAGCGTACCGTACTTCAATCCATAAAGTTGCGTTAAGCCGCACGGTTCAAAACGGCTGGGCACCAGAAt       <  1:67833/149‑1 (MQ=255)                                                                                                                                            cTGGCGACGCCATCTGCAAGGTTCTCAAGAGAACAGTCAGAAACCGTATCAGCAAGCCCACCGGTGCGCCGCACTAACGGCAGCGTACCGTACTTCAATCCATAAAGTTGCGTTAAGCCGCACGGTTCAAAACGGCTGGGCACCAGAAt       <  3:177224/149‑1 (MQ=255)                                                                                                                                               gCGACGCCATCTGCAAGGTTCTCAAGAGAACAGTCAGAAACCGTATCAGCAAGCCCACCGGTGCGCCGCACTAACGGCAGCGTACCGTACTTCAATCCATAAAGTTGCGTTAAGCCGCACGGTTCAAAACGGCTGGGCACCAGAATGAc    >  7:310360/1‑149 (MQ=255)                                                                                                                                               gCGACGCCATCTGCAAGGTTCTCAAGAGAACAGTCAGAAACCGTATCAGCAAGCCCACCGGTGCGCCGCACTAACGGCAGCGTACCGTACTTCAATCCATAAAGTTGCGTTAAGCCACACGGTTCAAAACGGCTGGGCACCAGAATGAc    <  2:300427/149‑1 (MQ=255)                                                                                                                                                cGACGCCATCTGCAAGGTTCTCAAGAGAACAGTCAGAAACCGTATCAGCAAGCCCACCGGTGCGCCGCACTAACGGCAGCGTACCGTACTTCAATCCATAAAGTTGCGTTAAGCCGCACGGTTCAAAACGGCTGGGCACCAGAATGACg   <  7:70245/149‑1 (MQ=255)                                                                                                                                                 gACGCCATCTGCAAGGTTCTCAAGAGAACAGTCAGAAACCGTATCAGCAAGCCCACCGGTGCGCCGCACTAACGGCAGCGTACCGTACTTCAATCCATAAAGTTGCGTTAAGCCGCACGGTTCAAAACGGCTGGGCACCAGAATGACGt  <  7:277451/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                  CGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACCACAGTACAAAAGCACGTCGAATAGCCCGTAACAGCGACCAGGCATTACTATCTTCAAAGACAAACCCACTGGCGACGCCATCTGCAAGGTTCTCAAGAGAACAGTCAGAAACCGTATCAGCAAGCCCACCGGTGCGCCGCACTAACGGCAGCGTACCGTACTTCAATCCATAAAGTTGCGTTAAGCCGCACGGTTCAAAACGGCTGGGCACCAGAATGACGT  >  NC_000913/3566639‑3566931

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