Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,950,784	A→G	5.2%	pseudogene (225/984 nt)	ilvG →	pseudogene, acetolactate synthase 2 large subunit, valine‑insensitive; acetolactate synthase II, large subunit, cryptic, interrupted

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,950,784	0	A	G	5.2%	244.2 / 3.3	78	pseudogene (225/984 nt)	ilvG	pseudogene, acetolactate synthase 2 large subunit, valine‑insensitive; acetolactate synthase II, large subunit, cryptic, interrupted
Reads supporting (aligned to +/- strand):  ref base A (38/35);  new base G (2/2);  total (40/38)
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.71e-01

TTATGCCGGTTTACGATGCATTGTATGACGGCGGCGTGGAGCACTTGCTATGCCGACATGAGCAGGGTGCGGCAATGGCGGCTATCGGTTATGCTCGTGCTACCGGCAAAACTGGCGTATGTATCGCCACGTCTGGTCCGGGCGCAACCAACCTGATAACCGGGCTTGCGGACGCACTGTTAGATTCCATCCCTGTTGTTGCCATCACCGGTCAAGTGTCCGCACCGTTTATCGGCACTGACGCATTTCAGGAAGTGGATGTCCTGGGATTGTCGTTAGCCTGTACCAAGCA  >  NC_000913/3950639‑3950930                                                                                                                                                  |                                                                                                                                                   ttATGCCGGTTTACGATGCATTGTATGACGGCGGCGTGGAGCACTTGCTATGCCGACATGAGCAGGGTGCGGCAATGGCGGCTATCGGTTATGCTCGTGCTACCGGCAAAACTGGCGTATGTATCGCCACGTCTGGTCCGGGCGcaac                                                                                                                                                  >  1:303760/1‑148 (MQ=255)         gTTTACGATGCATTGTATGACGGCGGCGTGGAGCACTTGCTATGCCGACATGAGCAGGGTGCGGCAATGGCGGCTATCGGTTATGCTCGTGCTACCGGCAAAACTGGCGTATGTATCGCCACGTCTGGTCCGGGCGCAACCAACCTGa                                                                                                                                          <  4:69450/148‑1 (MQ=255)          tttACGATGCATTGTATGACGGCGGCGTGGAGCACTTGCTATGCCGACATGAGCAGGGTGCGGCAATGGCGGCTATCGGTTATGCTCGTGCTACCGGCAAAACTGGCGTATGTATCGCCACGTCTGGTCCGGGCGCAACCAACCTGa                                                                                                                                          >  1:365272/1‑147 (MQ=255)             aCGATGCATTGTATGACGGCGGCGTGGAGCACTTGCTATGCCGACATGAGCAGGGTGCGGCAATGGCGGCTATCGGTTATGCTCGTGCTACCGGCAAAACTGGCGTATGTATCGCCACGTCTGGTCCGGGCGCAACCAACCTGATaa                                                                                                                                       <  4:263096/147‑1 (MQ=255)               gATGCATTGTATGACGGCGGCGTGGAGCACTTGCTATGCCGACATGAGCAGGGTGCGGCAATGGCGGCTATCGGTTATGCTCGTGCTACCGGCAAAACTGGCGTATGTATCGCCACGTCTGGTCCGGGCGCAACCAACCTGATAACCg                                                                                                                                    >  5:319269/1‑148 (MQ=255)                   cATTGTATGATGGTGGCTTGGAGCACTTGCTATGCCGACATGAGCAGGGTTCGGTAATGGCGGCTATCGGTTATGCTCGTGCTACCGGCAAAACTGGCGTATGTATCGCCACGTCTGGTCCGGGCGCAACCAACCTGATAACCTGGc                                                                                                                                 <  2:303760/147‑1 (MQ=255)                   cATTGTATGACGGCGGCGTGGAGCACTTGCTATGCCGACATGAGCAGGGTGCGGCAATGGCGGCTATCGGTTATGCTCGTGCTACCGGCAAAACTGGCGTATGTATCGCCACGTCTGGTCCGGGCGCAACCAACCTGATAACCGGGc                                                                                                                                 <  4:91600/147‑1 (MQ=255)                                     tGGAGCACTTGCTATGCCGACATGAGCAGGGTGCGGCAATGGCGGCTATCGGTTATGCTCGTGCTACCGGCAAAACTGGCGTATGTATCGCCACGTCTGGTCCGGGCGCAACCAACCTGATAACCGGGCTTGCGGACGCACTGt                                                                                                                  >  3:265527/1‑144 (MQ=255)                                      ggAGCACTTGCTATTCCGACATGAGCAGGGTGCGGCAATGGCGGCTATCGGTTATGCTCGTGCTACCGGCAAAACTGGCGTATGTATCGCCACGTCTGGTCCGGGCGCAACCAACCTGATAACCGGGCTTGCGGACGCACTGt                                                                                                                  <  2:388040/143‑1 (MQ=255)                                      ggAGCACTTGCTATGCCGACATGAGCAGGGTGCGGCAATGGCGGCTATCGGTTATGCTCGTGCTACCGGCAAAACTGGCGTATGTATCGCCACGTCTGGTCCGGGCGCAACCAACCTGATAACCGGGCTTGCGGACGCACTGt                                                                                                                  >  1:388040/1‑143 (MQ=255)                                      ggAGCACTTGCTATGCCGACATGAGCAGGGTGCGGCAATGGCGGCTATCGGTTATGCTCGTGCTACCGGCAAAACTGGCGTATGTATCGCCACGTCTGGTCCGGGCGCAACCAACCTGATAACCGGGCTTGCGGACGCACTGTTAGAt                                                                                                             >  5:156917/1‑148 (MQ=255)                                      ggAGCACTTGCTATGCCGACATGAGCAGGGTGCGGCAATGGCGGCTATCGGTTATGCTCGTGCTACCGGCAAAACTGGCGTATGTATCGCCACGTCTGGTCCGGGCGCAACCAACCTGATAACCGGGCTTGCGGACGCACTGTTAGAt                                                                                                             >  5:319246/1‑148 (MQ=255)                                        aGCACTTGCTATGCCGACATGAGCAGGGTGCGGCAATGGCGGCTATCGGTTATGCTCGTGCTACCGGCAAAACTGGCGTATGTATCGCCACGTCTGGTCCGGGCGCAACCAACCTGATAACCGGGCTTGCGGACGCAc                                                                                                                     <  2:249890/138‑1 (MQ=255)                                        aGCACTTGCTATGCCGACATGAGCAGGGTGCGGCAATGGCGGCTATCGGTTATGCTCGTGCTACCGGCAAAACTGGCGTATGTATCGCCACGTCTGGTCCGGGCGCAACCAACCTGATAACCGGGCTTGCGGACGCAc                                                                                                                     >  1:249890/1‑138 (MQ=255)                                          cACTTGCTATGCCGACATGAGCAGGGTGCGGCAATGGCGGCTATCGGTTATGCTCGTGCTACCGGCAAAACTGGCGTATGTATCGCCACGTCTGGTCCGGGCGCAACCAACCTGATAACCGGGCTTGCGGACGCACTGTTAGATTCCa                                                                                                         >  8:248287/1‑148 (MQ=255)                                              tGCTATGCCGACATGAGCAGGGTGCGGCAATGGCGGCTATCGGTTATGCTCGTGCTACCGGCAAAACTGGCGTATGTATCGCCACGTCTGGTCCGGGCGCGACCAACCTGATAACCGGGCTTGCGGACGCACTGTTAGATTCCATccc                                                                                                     >  5:13139/1‑148 (MQ=255)                                              tGCTATGCCGACATGAGCAGGGTGCGGCAATGGCGGCTATCGGTTATGCTCGTGCTACCGGCAAAACTGGCGTATGTATCGCCACGTCTGGTCCGGGCGCGACCAACCTGATAACCGGGCTTGCGGACGCACTGTTAGATTCCATccc                                                                                                     >  5:7730/1‑148 (MQ=255)                                              tGCTATGCCGACATGAGCAGGGTGCGGCAATGGCGGCTATCGGTTATGCTCGTGCTACCGGCAAAACTGGCGTATGTATCGCCACGTCTGGTCCGGGCGCAACCAACCTGATAACCGGGCTTGCGGACGc                                                                                                                       <  6:359318/130‑1 (MQ=255)                                              tGCTATGCCGACATGAGCAGGGTGCGGCAATGGCGGCTATCGGTTATGCTCGTGCTACCGGCAAAACTGGCGTATGTATCGCCACGTCTGGTCCGGGCGCAACCAACCTGATAACCGGGCTTGCGGACGc                                                                                                                       >  5:359318/1‑130 (MQ=255)                                              tGCTATGCCGACATGAGCAGGGTGCGGCAATGGCGGCTATCGGTTATGCTCGTGCTACCGGCAAAACTGGCGTATGTATCGCCACGTCTGGTCCGGGCGCAACCAACCTGATAACCGGGCTTGCGGACGCACTGTTAGATTCCATccc                                                                                                     >  5:226677/1‑148 (MQ=255)                                                            gAGCAGGGTGCGGCAATGGCGGCTATCGGTTATGCTCGTGCTACCGGCAAAACTGGCGTATGTATCGCCACGTCTGGTCCGGGCGCAACCAACCTGATAACCGGGCTTGCgg                                                                                                                           >  2:295257/1‑112 (MQ=255)                                                            gAGCAGGGTGCGGCAATGGCGGCTATCGGTTATGCTCGTGCTACCGGCAAAACTGGCGTATGTATCGCCACGTCTGGTCCGGGCGCAACCAACCTGATAACCGGGCTTGCgg                                                                                                                           <  1:295257/112‑1 (MQ=255)                                                            gAGCAGGGTGCGGCAATGGCGGCTATCGGTTATGCTCGTGCTACCGGCAAAACTGGCGTATGTATCGCCACGTCTGGTCCGGGCGCAACCAACCTGATAACCGGGCTTGCGGACGCACTGTTAGATTCCATCCCTGTTGTTGCCATCa                                                                                       >  8:113078/1‑148 (MQ=255)                                                             aGCAGGGTGCGGCAATGGCGGCTATCGGTTATGCTCGTGCTACCGGCAAAACTGGCGTATGTATCGCCACGTCTGGTCCGGGCGCAACCAACCTGATAACCGGGCTTGCGGACGCACTGTTAGATTCCATCCCTGTTGTTGCCATCAc                                                                                      <  2:365272/148‑1 (MQ=255)                                                                   gTGCGGCAATGGCGGCTATCGGTTATGCTCGTGCTACCGGCAAAACTGGCGTATGTATCGCCACGTCTGGTCCGGGCGCAACCAACCTGATAACCGGGCTTGCGGACGCACTGTTAGATTCCATCCCTGTTGTTGCCATCACCGGTc                                                                                 >  5:288058/1‑147 (MQ=255)                                                                   gTGCGGCAATGGCGGCTATCGGTTATGCTCGTGCTACCGGCAAAACTGGCGTATGTATCGCCACGTCTGGTCCGGGCGCAACCAACCTGATAACCGGGCTTGCGGACGCACTGTTAGATTCCATCCCTGTTGTTGCCATCACCGGTc                                                                                 >  5:214186/1‑147 (MQ=255)                                                                   gTGCGGCAATGGCGGCTATCGGTTATGCTCGTGCTACCGGCAAAACTGGCGTATGTATCGCCACGTCTGGTCCGGGCGCAACCAACCTGATAACCGGGCTTGCGGACGCACTGTTAGATTCCATCCCTGTTGTTGCCATCACCGGTc                                                                                 >  7:380202/1‑147 (MQ=255)                                                                   gTGCGGCAATGGCGGCTATCGGTTATGCTCGTGCTACCGGCAAAACTGGCGTATGTATCGCCACGTCTGGTCCGGGCGCAACCAACCTGATAACCGGGCTTGCGGACGCACTGTTAGATTCCATCCCTGTTGTTGCCATCACCGGTCa                                                                                <  8:380202/148‑1 (MQ=255)                                                                         cAATGGCGGCTATCGGTTATGCTCGTGCTACCGGCAAAACTGGCGTATGTATCGCCACGTCTGGTCCGGGCGCAACCAACCTGATAACCGGGCTTGCGGACGCACTGTTAGATTCCATCCCTGTTGTTGCCATCACCGGTCAAGTGTc                                                                          <  6:226677/148‑1 (MQ=255)                                                                          aaTGGCGGCTATCGGTTATGCTCGTGCTACCGGCAAAACTGGCGTATGTATCGCCACGTCTGGTCCGGGCGCAACCAACCTGATAACCGGGCTTGCGGACGCACTGTTAGATTcc                                                                                                          <  4:210820/115‑1 (MQ=255)                                                                          aaTGGCGGCTATCGGTTATGCTCGTGCTACCGGCAAAACTGGCGTATGTATCGCCACGTCTGGTCCGGGCGCAACCAACCTGATAACCGGGCTTGCGGACGCACTGTTAGATTcc                                                                                                          >  3:210820/1‑115 (MQ=255)                                                                          aaTGGCGGCTATCGGTTATGCTCGTGCTACCGGCAAAACTGGCGTATGTATCGCCACGTCTGGTCCGGGCGCAACCAACCTGATAACCGGGCTTGCGGACGCACTGTTAGATTCCATCCCTGTTGTTGCCATCACCgg                                                                                   <  4:188354/138‑1 (MQ=255)                                                                          aaTGGCGGCTATCGGTTATGCTCGTGCTACCGGCAAAACTGGCGTATGTATCGCCACGTCTGGTCCGGGCGCAACCAACCTGATAACCGGGCTTGCGGACGCACTGTTAGATTCCATCCCTGTTGTTGCCATCACCgg                                                                                   >  3:188354/1‑138 (MQ=255)                                                                           aTGGCGGCTATCGGTTATGCTCGTGCTACCGGCAAAACTGGCGTATGTATCGCCACGTCTGGTCCGGGCGCAACCAACCTGATAACCGGGCTTGCGGACGCACTGTTAGATTCCATCCCTGTTGTTGcc                                                                                           <  4:222234/129‑1 (MQ=255)                                                                           aTGGCGGCTATCGGTTATGCTCGTGCTACCGGCAAAACTGGCGTATGTATCGCCACGTCTGGTCCGGGCGCAACCAACCTGATAACCGGGCTTGCGGACGCACTGTTAGATTCCATCCCTGTTGTTGcc                                                                                           >  3:222234/1‑129 (MQ=255)                                                                             ggcggcTATCGGTTATGCTCGTGCTACCGGCAAAACTGGCGTATGTATCGCCACGTCTGGTCCGGGCGCAACCAACCTGATAACCGGGCTTGCGGACGCACTGTTAGATTCCATCCCTGTTGTTGCCATCACCGGTCAAGTGTCCGCa                                                                      >  6:215592/1‑148 (MQ=255)                                                                             ggcggcTATCGGTTATGCTCGTGCTACCGGCAAAACTGGCGTATGTATCGCCACGTCTGGTCCGGGCGCAACCAACCTGATAACCGGGCTTGCGGACGCACTGTTAGATTCCAACCCTGTTGTTGCCATCACCGCTCAAGTGTCCCCa                                                                      >  7:209121/1‑148 (MQ=255)                                                                               cggcTATCGGTTATGCTCGTGCTACCGGCAAAACTGGCGTATGTATCGCCACGTCTGGTCCGGGCGCAACCAACCTGATAACCGGGCTTGCGGACGCACTGTTAGATTCCATCCCTGTTGTTGCCATCACCGGTCAAGTGTCCGCAcc                                                                    >  5:33564/1‑148 (MQ=255)                                                                                 gcTATCGGTTATGCTCGTGCTACCGGCAAAACTGGCGTATGTATCGCCACGTCTGGTCCGGGCGCAACCAACCTGATAACCGGGCTTGCGGACGCACTGTTAGATTCCATcc                                                                                                      >  6:293800/1‑112 (MQ=255)                                                                                 gcTATCGGTTATGCTCGTGCTACCGGCAAAACTGGCGTATGTATCGCCACGTCTGGTCCGGGCGCAACCAACCTGATAACCGGGCTTGCGGACGCACTGTTAGATTCCATcc                                                                                                      <  5:293800/112‑1 (MQ=255)                                                                                     tCGGTTATGCTCGTGCTACCGGCAAAACTGGCGTATGTATCGCCACGTCTGGTCCGGGCGCAACCAACCTGATAACCGGGCTTGCGGACGCACTGTTAGATTCCATCCCTGTTGTTGCCATCACCGGTCAAGTGTCCGCACCGttt                                                                <  5:215592/146‑1 (MQ=255)                                                                                         ttATGCTCGTGCTACCGGCAAAACTGGCGTATGTATCGCCACGTCTGGTCCGGGCGCGACCAACCTGATAACCGGGCTTGCGGACGCACTGTTAGATTCCATCCCTGTTGTTGCCATCACCGGTCAAGTGTCCGCACCGTTTATCGGc                                                          <  6:7730/148‑1 (MQ=255)                                                                                         ttATGCTCGTGCTACCGGCAAAACTGGCGTATGTATCGCCACGTCTGGTCCGGGCGCGACCAACCTGATAACCGGGCTTGCGGACGCACTGTTAGATTCCATCCCTGTTGTTGCCATCACCGGTCAAGTGTCCGCACCGTTTATCGGc                                                          <  6:13139/148‑1 (MQ=255)                                                                                         ttATGCTCGTGCTACCGGCAAAACTGGCGTATGTATCGCCACGTCTGGTCCGGGCGCAACCAACCTGATAACCGGGCTTGCGGACGCACTGTTAGATTCCATCCCTGTTGTTGCCATCACCGGTCAAGTGTCCGCACCGTTTATCGGc                                                          <  6:319246/148‑1 (MQ=255)                                                                                              cTCGTGCTACCGGCAAAACTGGCGTATGTATCGCCACGTCTGTTCCGGGCGCAACCAACCTGATAACCGGGCTTGCGGACGCACTGTTAGATTCCATCCCTGTTGTTGCCATCACCGGTCAAGTGTCCGCACCGTTTATCGGCACTGa                                                     <  6:156917/148‑1 (MQ=255)                                                                                                  tGCTACCGGCAAAACTGGCGTATGTATCGCCACGTCTGGTCCGGGCGCAACCAACCTGATAACCGGGCTTGCGGACGCACTGTTAg                                                                                                               <  4:124975/86‑1 (MQ=255)                                                                                                  tGCTACCGGCAAAACTGGCGTATGTATCGCCACGTCTGGTCCGGGCGCAACCAACCTGATAACCGGGCTTGCGGACGCACTGTTAg                                                                                                               >  3:124975/1‑86 (MQ=255)                                                                                                   gCTACCGGCAAAACTGGCGTATGTATCGCCACGTCTGGTCCGGGCGCAACCAACCTGATAACCGGGTTTGCGGACGCACTGTTAGATTCCATCCCTGTTGTTGCCATCACCGGTCAAGTGTCCGCACCGTTTATCGGCACTGACGCAt                                                >  7:276337/1‑148 (MQ=255)                                                                                                    cTACCGGCAAAACTGGCGTATGTATCGCCACGTCTGGTCCGGGCGCAACCAACCTGATAACCGGGTTTGCGGACGCACTGTTAGATTCCATCCCTGTTGTTGCCATCACCTGTCAAGTGTCCGCACCGTTTATCGGCACTGACGCAtt                                               <  8:276337/148‑1 (MQ=255)                                                                                                      aCCGGCAAAACTGGCGTATGAATCGCCACGTCTGGTCCGGGCGCAACCAACCTGATAACCGGGCTTGCGGACGCACTGTTAGATTCCATCCCTGTTGTTGCCATCACCGGTCAAGTGTCCGCACCGTTTATCGGCACTGACGCAttt                                              >  3:235798/1‑147 (MQ=255)                                                                                                          gCAAAACTGGCGTATGTATCGCCACGTCTGGTCCGGGCGCAACCAACCTGATAACCGGGCTTGCGGACGCACTGTTAGATTCCATCCCTGTTGTTGCCATCACCGGTCAAGTGTCCGCACCGTTTATCGGCACTGACGCATTTCAGGa                                         <  7:248287/148‑1 (MQ=255)                                                                                                           cAAAACTGGCGTATGTATCGCCACGTCTGGTCCGGGCGCAACCAACCTGATAACCGGGCTTGCGGACGCACTGTTAGATTCCATCCCTGTTGTTGCCATCACCGGTCAAGTGTCCGCACCGTTTATCGGCACTGACGCATTTCAGGaa                                        <  7:113078/148‑1 (MQ=255)                                                                                                           cAAAACTGGCGTATGTATCGCCACGTCTGGTCCGGGCGCAACCAACCTGATAACCGGGCTTGCGGACGCACTGTTAGATTCCATCCCTGTTGTTGCCATCACCGGTCAAGTGTCCGCACCGTTTATCGGCACTGACGCATTTCAGGaa                                        <  4:232931/148‑1 (MQ=255)                                                                                                              aaCTGGCGTATGTATCGCCACGTCTGGTCCGGGCGCAACCAACCTGATAACCGGGCTTGCGGACGCACTGTTAGATTCCATCCCTGTTGTTGCCATCACCGGTCAAGTGTCCGCACCGTTTATCGGCACTGACGCATTTCAGGAAGTg                                     <  6:319269/148‑1 (MQ=255)                                                                                                              aaCTGGCGTATGTATCGCCACGTCTGGTCCGGGCGCAACCAACCTGATAACCGGGCTTGCGGACGCACTGTTAGATTCCATCCCTGTTGTTGCCATCACCGGTCAAGTGTCCGCACCGTTTATCGGCACTGACGCATTTCAGGAAGTg                                     <  6:33564/148‑1 (MQ=255)                                                                                                               aCTGGCGTATGTATCGCCACGTCTGGTCCGGGCGCAACCAACCTGATAACCGGGCTTGCGGACGCACTGTTAGATTCCATCCCTGTTGTTGCCATCACCGGTCAAGTGTCCGCACCGTTTATCGGCACTGACGCATTTCAGGAAGTgg                                    <  6:288058/148‑1 (MQ=255)                                                                                                                      tatgtatCGCCACGTCTGGTCCGGGCGCAACCAACCTGATAACCGGGCTTGCGGACGCACTGTTAGATTCCATCCCTGTTGTTGCCATCACCGGTCAAGTGTCCGCACCGTTTATCGGCACTGACGCATTTCAGGAAGTGGATGTCCt                             <  6:214186/148‑1 (MQ=255)                                                                                                                      tatggatCGCCACGTCTGGTCCGGGCGNANNCAANCTGATAACCGGGCTTGCGGACGCACTGTTAGATTCCATCCCTGTTGTTGCCATCACCGGTCAAGTGTCCGCACCGTTTATCGGCACTGACGCATTTCAGGAAGTGGATGTCCt                             <  8:209121/148‑1 (MQ=255)                                                                                                                       atgtatCGCCACGTCTGGTCCGGGCGCAACCAACCTGATAACC‑GGCTTGCGGACGCCCGGTTAGATTCCATCCCGGGTGGTGCCATCACCCGTCAAGTGTCCGCACCGTTTATCGGCa                                                         >  4:100056/1‑118 (MQ=255)                                                                                                                         gtatCGCCACGTCTGGTCCGGGCGCAACCAACCTGATAACCGTGCTTGCGGACGCACTGTTAGATTCCATCCCTGTTGTTGCCATCACCGGTCAAGTGTCCGCACCGTTTATCGGCACTTACGCATTTCAGGAAGTGGATGTCCTgt                           >  8:359378/1‑146 (MQ=255)                                                                                                                         gtatCGCCACGTCTGGTCCGGGCGCAACCAACCTGATAACCGGGCTTGCGGACGCACTGTTAGATTCCATCCCTGTTGTTGCCCTCACCGGTCAAGTGTCCGCACCGTTTATCGGCACTGACGCATTTCAGGAAGTGGATGTCCTggg                          >  3:55442/1‑148 (MQ=255)                                                                                                                         gtatCGCCACGTCTGGTCCGGGCGCAACCAACCTGATAACCGGGCTTGCGGACGCACTGTTAGATTCCATCCCTGTTGTTGCCATCACCGGTCAAGTGTCCGCACCGTTTATCGGCACTGACGCATTTCAGGAAGTGGATGTCCTgg                           >  8:162937/1‑147 (MQ=255)                                                                                                                         gtatCGCCACGTCTGGTCCGGGCGCAACCAACCTGATAACCGGGCTTGCGGACGCACTGTTAGATTCCATCCCTGTTGTTGCCATCACCGGTCAAGTGTCCGCACCGTTTATCGGCACTGACGCATTTCAGGAAGTGGATGTCCTgg                           >  8:133209/1‑147 (MQ=255)                                                                                                                             cGCCACGTCTGGTCCGGGCGCAACCAACCTGATAACCGGGCTTGCGGACGCACTGTTAGATTCCATCCCTGTTGTTGCCATCACCGGTCAAGTGTCCGCACCGTTTATCGGCACTGACGCAt                                                >  6:352180/1‑122 (MQ=255)                                                                                                                             cGCCACGTCTGGTCCGGGCGCAACCAACCTGATAACCGGGCTTGCGGACGCACTGTTAGATTCCATCCCTGTTGTTGCCATCACCGGTCAAGTGTCCGCACCGTTTATCGGCACTGACGCAt                                                <  5:352180/122‑1 (MQ=255)                                                                                                                               ccACTTCTGGTCCGGGAGCTACCAACCTGATAACCGGGCTTGCGTACGCACTGTTAGATTCCATCCCTGTTGTTGCCATCACCGGTCAAGTGTCGGCACCGTTTATCGTCACTGACGCATTTCAGGAAGTGGATGTCCTGGGATTGTc                    <  7:359378/148‑1 (MQ=255)                                                                                                                               ccACGTCTGGTCCGGGCGCAACCAACCTGATAACCGGGCTTGCGGACGCACTGTTAGATTCCATCCCTGTTGTTGCCATCACCGGTCAAGTGTCCGCACCGTTTATCGGCACTGACGCATTTCAGGAAGTGGATGTCCTGGGATTGTc                    >  5:150252/1‑148 (MQ=255)                                                                                                                                 aCGTCTGGTCCGGGCGCAACCAACCTGATAACCGGGCTTGCGGACGCACTGTTAGATTCCATCCCTGTTGTTGCCATCACCGGTCAAGTGTCCGCACCGTTTATCGGCACTGACGCATTTCAGGAAGTGGATGTCCTGGGATTGTCGt                  <  4:235798/148‑1 (MQ=255)                                                                                                                                   gTCTGGTCCGGGCGCAACCAACCTGATAACCGGGCTTGCGGACGCACTGTTAGATTCCATCCCTGTTGTTGCCATCACCGGTCAAGTGTCCGCACCGTTTATCGGc                                                          <  7:188043/106‑1 (MQ=255)                                                                                                                                   gTCTGGTCCGGGCGCAACCAACCTGATAACCGGGCTTGCGGACGCACTGTTAGATTCCATCCCTGTTGTTGCCATCACCGGTCAAGTGTCCGCACCGTTTATCGGc                                                          >  8:188043/1‑106 (MQ=255)                                                                                                                                   gTCTGGTCCGGGCGCAACCAACCTGATAACCGGGCTTGCGGACGCACTGTTAGATTCCATCCCTGTTGTTGCCATCACCGGTCAAGTGTCCGCACCGTTTATCGGCACTGACGc                                                  >  7:230385/1‑114 (MQ=255)                                                                                                                                   gTCTGGTCCGGGCGCAACCAACCTGATAACCGGGCTTGCGGACGCACTGTTAGATTCCATCCCTGTTGTTGCCATCACCGGTCAAGTGTCCGCACCGTTTATCGGCACTGACGc                                                  <  8:230385/114‑1 (MQ=255)                                                                                                                                   gTCTGGTCCGGGCGCAACCAACCTGATAACCGGGCTTGCGGACGCACTGTTAGATTCCATCCCTGTTGTTGCCATCACCGGTCAAGTGTCCGCACCGTTTATCGGCACTGACGCATTTCAGGAAGTGGATGTCCTGGGATTGTCGt                  >  7:25944/1‑146 (MQ=255)                                                                                                                                   gTCTGGTCCGGGCGCAACCAACCTGATAACCGGGCTTGCGGACGCACTGTTAGATTCCATCCCTGTTGTTGCCATCACCGGTCAAGTGTCCGCACCGTTTATCGGCACTGACGCATTTCAGGAAGTGGATGTCCTGGGATTGTCGt                  >  7:32063/1‑146 (MQ=255)                                                                                                                                    tCTGGTCCGGGCGCAACCAACCTGATAACCGGGCTTGCGGACGCACTGTTAGATTCCATCCCTGTTGTTGCCATCACCGGTCAAGTGTCCGCACCGTTTATCGGCACTGACGCATTTCAGGAAGTGGATGTCCTGGGATTGTCGTTAg               <  7:162937/148‑1 (MQ=255)                                                                                                                                    tCTGGTCCGGGCGCAACCAACCTGATAACCGGGCTTGCGGACGCACTGTTAGATTCCATCCCTGTTGTTGCCATCACCGGTCAAGTGTCCGCACCGTTTATCGGCACTGACGCATTTCAGGAAGTGGATGTCCTGGGATTGTCGTTAg               >  3:143463/1‑148 (MQ=255)                                                                                                                                         tCCGGGCGCAACCAACCTGATAACCGGGCTTGCGGACGCACTGTTAGATTCCATCCCTGTTGTTGCCATCACCGGTCAAGTGTCCGCACCGTTTATCGGCACTGACGCATTTCAGGAAGTGGATGTCCTGGGATTGTCGTTAGCCTGt          <  7:133209/148‑1 (MQ=255)                                                                                                                                          ccGGGCGCAACCAACCTGATAACCGGGCTTGCGGACGCACTGTTAGATTCCATCCCTGTTGTTGCCATCACCGGTCAAGTGTCCGCACCGTTTATCGGCACTg                                                      <  6:90514/103‑1 (MQ=255)                                                                                                                                          ccGGGCGCAACCAACCTGATAACCGGGCTTGCGGACGCACTGTTAGATTCCATCCCTGTTGTTGCCATCACCGGTCAAGTGTCCGCACCGTTTATCGGCACTg                                                      >  5:90514/1‑103 (MQ=255)                                                                                                                                                  aaccaaccTGATAACCGGGCTTGCGGACGCACTGTTAGATTCCATCCCTGTTGTTGCCATCACCGGTCAAGTGTCCGCACCGTTTATCGGCACTGACGCATTTCAGGAAGTGGATGTCCTGGGATTGTCGTTAGCCTGTACCAAGca  >  4:63136/1‑147 (MQ=255)                                                                                                                                                  |                                                                                                                                                   TTATGCCGGTTTACGATGCATTGTATGACGGCGGCGTGGAGCACTTGCTATGCCGACATGAGCAGGGTGCGGCAATGGCGGCTATCGGTTATGCTCGTGCTACCGGCAAAACTGGCGTATGTATCGCCACGTCTGGTCCGGGCGCAACCAACCTGATAACCGGGCTTGCGGACGCACTGTTAGATTCCATCCCTGTTGTTGCCATCACCGGTCAAGTGTCCGCACCGTTTATCGGCACTGACGCATTTCAGGAAGTGGATGTCCTGGGATTGTCGTTAGCCTGTACCAAGCA  >  NC_000913/3950639‑3950930

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