Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	128,803	G→A	5.9%	D298N (GAC→AAC)	lpd →	dihydrolipoyl dehydrogenase; E3 component of pyruvate and 2‑oxoglutarate dehydrogenases complexes; glycine cleavage system L protein; dihydrolipoamide dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	128,803	0	G	A	5.9%	208.8 / 2.1	68	D298N (GAC→AAC)	lpd	dihydrolipoyl dehydrogenase; E3 component of pyruvate and 2‑oxoglutarate dehydrogenases complexes; glycine cleavage system L protein; dihydrolipoamide dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (34/29);  new base A (2/2);  total (36/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.05e-01

TGGAAGGCAAAAAAGCACCCGCTGAACCGCAGCGTTACGACGCCGTGCTGGTAGCGATTGGTCGTGTGCCGAACGGTAAAAACCTCGACGCAGGCAAAGCAGGCGTGGAAGTTGACGACCGTGGTTTCATCCGCGTTGACAAACAGCTGCGTACCAACGTACCGCACATCTTTGCTATCGGCGATATCGTCGGTCAACCGATGCTGGCACACAAAGGTGTTCACGAAGGTCACGTTGCCGCTGAAGTTATCGCCGGTAAGAAACACTACTTCGATCCGA  >  NC_000913/128666‑128944                                                                                                                                          |                                                                                                                                              tGGAAGGCAAAAAAGCACCCGCTGCACCGCAGCGTTACGACGCCGTGCTGGTAGCGATTGGTCGGGTGCCGAACGGTAAAAACCTCGCCGCAGGCAAAGCAGGCGTGGAAGTTGACGACCGTGGTTTCATCCGCGTTGACACACAGCTg                                                                                                                                    >  6:97230/1‑149 (MQ=255)   gAAGGCAAAAAAGCACCCGCTGAACCGCAGCGTTACGACGCCGTGCTGGTAGCGATTGGTCGTGTGCCGAACGGTAAAAACCTCGACGCAGGCAAAGCAGGCGTGGAAGTTGACGACCGTGGTTTCATCCGCGTTGACAAACAGCTGc                                                                                                                                   <  1:22512/148‑1 (MQ=255)   gAAGGCAAAAAAGCACCCGCTGAACCGCAGCGTTACGACGCCGTGCTGGTAGCGATTGGTCGTGTGCCGAACGGTAAAAACCTCGACGCAGGCAAAGCAGGCGTGGAAGTTGACGACCGTGGTTTCATCCGCGCTGACAAACAGCTGCg                                                                                                                                  >  4:341310/1‑149 (MQ=255)     aGGCAAAAAAGCACCCGCTGAACCGCAGCGTTACGACGCCGTGCTGGTAGCGATTGGTCGTGTGCCGAACGGTAAAAACCTCGACGCAGGCAAAGCAGGCGTGGAAGTTGACGACCGTGGTTTCATCCGCGTTGACAAACAGCTGCGTa                                                                                                                                >  5:198011/1‑149 (MQ=255)         aaaaaaGCACCCGCTGAACCGCAGCGTTACGACGCCGTGCTGGTAGCGATTGGTCGTGTGCCGAACGGTAAAAACCTCGACGCAGGCAAAGCAGGCGTGGAAGTTGACGACCGTGGTTTCATCCGCGTTGACAAACAGCTGCGTACCaa                                                                                                                            >  2:366043/1‑149 (MQ=255)             aaGCACCCGCTGAACCGCAGCGTTACGACGCCGTGCTGGTAGCGATTGGTCGTGTGCCGAACGGTAAAAACCTCGACGCAGGCAAAGCAGGCGTGGAAGTTGACGACCGTGGTTTCATCCGCGTTGACAAACAGCTGCGTACCAACGTa                                                                                                                        <  1:15276/149‑1 (MQ=255)              aGCACCCGCTGAACCGCAGCGTTACGACGCCGTGCTGGTAGCGATTGGTCGTGTGCCGAACGGTAAAAACCTCGACGCAGGCAAAGCAGGCGTGGAAGTTGACGACCGTGGTTTCATCCGCGTTGACAAACAGCTGCGTACCAACGTAc                                                                                                                       <  4:182428/149‑1 (MQ=255)                 aCCCGCTGAACCGCAGCGTTACGACGCCGTGCTGGTAGCGATTGGTCGTGTGCCGAACGGTAAAAACCTCGACGCAGGCAAAGCAGGAGTGGAAGTTGACGACCGTGGTTTCATCCGCGTTGACAAACAGCTGCGTACCAACGTACCGc                                                                                                                    <  7:268230/149‑1 (MQ=255)                  cccGCTGAACCGCAGCGTTACGACGCCGTGCTGGTAGCGATTGGTCGTGTGCCGAACGGTAAAAACCTCGACGCAGGCAAAGCAGGCGTGGAAGTTGACGACCGTGGTTTCATCCGCGTTGACAAACAGCTGCGTACCAACGTACCGca                                                                                                                   >  5:334838/1‑149 (MQ=255)                    cGCTGAACCGCAGCGTTACGACGCCGTGCTGGTAGCGATTGGTCGTGTGCCGAACGGTAAAAACCTCGACGCAGGCAAAGCAGGCGTGGAAGTTGACGACCGTGGTTTCATCCGCGTTGACAAACAGCTGCGTACCAACGTACCGcaca                                                                                                                 >  6:58378/1‑149 (MQ=255)                       tGAACCGCAGCGTTACGACGCCGTGCTGGTAGCGATTGGTCGTGTGCCGAACGGTAAAAACCTCGACGCAGGCAAAGCAGGCGTGGAAGTTGACGACCGTGGTTTCATCCGCGTTGACAAACAGCTGCGTACCAACGTACCGCACATCt                                                                                                              >  3:280449/1‑149 (MQ=255)                         aaCCGCAGCGTTACGACGCCGTGCTGGTAGCGATTGGTCGTGTGCCGAACGGTAAAAACCTCGACGCAGGCAAAGCAGGCGTGGAAGTTGACGACCGTGGTTTCATCCGCGTTGACAAACAGCTGCGTACCAACGTACCGCACATCttt                                                                                                            >  4:161452/1‑149 (MQ=255)                           ccGCAGCGTTACGACGCCGTGCTGGTAGCGATTGGTCGTGTGCCGAACGGTAAAAACCTCGACGCAGGCAAAGCAGGCGTGGAAGTTGACGACCGTGGTTTCATCCGCGTTGACAAACAGCTGCGTACCAACGTACCGCACATCTTTGc                                                                                                          <  4:188205/149‑1 (MQ=255)                            cGCAGCGTTACGACGCCGTGCTGGTAGCGATTGGTCGTGTGCCGAACGGTAAAAACCTCGACGCAGGCAAAGCAGGCGTGGAAGTTGACGACCGTGGTTTCATCCGCGTTGACAAACAGCTGCGTACCa                                                                                                                             <  4:286130/129‑1 (MQ=255)                            cGCAGCGTTACGACGCCGTGCTGGTAGCGATTGGTCGTGTGCCGAACGGTAAAAACCTCGACGCAGGCAAAGCAGGCGTGGAAGTTGACGACCGTGGTTTCATCCGCGTTGACAAACAGCTGCGTACCAACGTACCGCACATCTTTGc                                                                                                          >  8:362738/1‑148 (MQ=255)                             gCAGCGTTACGACGCCGTGCTGGTAGCGATTGGTCGTGTGCCGAACGGTAAAAACCTCGACGCAGGCAAAGCAGGCGTGGAAGTTGACGACCGTGGTTTCATCCGCGTTGACAAACAGCTGCGTACCAACGTACCGCACATCTTTGCTa                                                                                                        <  1:366043/149‑1 (MQ=255)                              cAGCGTTACGACGCCGTGCTGGTAGCGATTGGTCGTGTGCCGAACGGTAAAAACCTCGACGCAGGCAAAGCAGGCGTGGAAGTTGACGACCGTGGTTTCATCCGCGTTGACAAACAGCTGCGTACCAACGTa                                                                                                                        <  4:328444/132‑1 (MQ=255)                              cAGCGTTACGACGCCGTGCTGGTAGCGATTGGTCGTGTGCCGAACGGTAAAAACCTCGACGCAGGCAAAGCAGGCGTGGAAGTTGACGACCGTGGTTTCATCCGCGTTGACAAACAGCTGCGTACCAACGTa                                                                                                                        >  3:328444/1‑132 (MQ=255)                                     acgacgCCGTGCTGGTAGCGATTGGTCGTGTGCCGAACGGTAAAAACCTCGACGCAGGCAAAGCAGGCGTGGAAGTTGACGACCGTGGTTTCATCCGCGTTGACAAACAGCTGCGTACCAACGTACCGCACATCTTTGCTATCGGCGa                                                                                                 <  1:71706/148‑1 (MQ=255)                                            cGTGCTGGTAGCGATTGGTCGTGTGCCGAACGGTAAAAACCTCGACGCAGGCAAAGCAGGCGTGGAAGTTGACGACCGTGGTTTCATCCGCGTTGACAAACAGCTGCGTACCAACGTACCGCACATCTTTGCTATCGGCGATAtcgtcg                                                                                         <  3:55968/149‑1 (MQ=255)                                             gTGCTTGTAGCGATTGGTCGTGTGCCGAACGGTAAAAACCTCGACGCAGGCAAAGCAGGCGTGGAAGTTGACGACCGTGGTTTCATCCGCGTTGACAAACAGCTGCGTACCAACGTACCGCACATCTTTGCTATCGGCGATATCGTCgg                                                                                        <  1:285197/149‑1 (MQ=255)                                                   gTAGCGATTGGTCGTGTGCCGAACGGTAAAAACCTCGACGCAGGCAAAGCAGGCGTGGAAGTTGACGACCGTGGTTTCATCCGCGTTGACAAACAGCTGCGTACCAACGTACCGCACATCTTTGCTATCGGCGATATCGTCGGTCAAcc                                                                                  >  4:331454/1‑149 (MQ=255)                                                      gCGATTGGTCGTGTGCCGAACGGTAAAAACCTCGACGCAGGCAAAGCAGGCGTGGAAGTTGACGACCGTGGTTTCATCCGCGTTGACAAACAGCTGCGTACCAACGTACCGCACATCTTTGCTATCGGCGATATCGTCGGTCAACCGAt                                                                               >  5:320112/1‑149 (MQ=255)                                                       cGATTGGTCGGGTGCCGAACGGTAAAAACCTCCACGCAGGCAAAGCAGGCGTGGAAGTTGACGACCGTGGTTTCATCCGCGTTGACAAACAGCTGCGTACCAACGTACCGCACATCTTTGCTATCGGCGATATCGTCGGTCAACCGATg                                                                              <  2:149794/149‑1 (MQ=255)                                                             gTCGTGTGCCGAACGGTAAAAACCTCGACGCAGGCAAAGCAGGCGTGGAAGTTGACGACCGTGGTTTCATCCGCGTTGACAAACAGCTGCGTACCAACGTACCGCACATCTTTGCTATCGGCGATATCGTCGGTCAGCCTCTGCTGGca                                                                        >  2:273968/1‑149 (MQ=255)                                                                    gCCGAACGGTAAAAACCTCGACGCAGGCAAAGCAGGCGTGGAAGTTGCCGACCGTGGTTTCATCCGCGTTAACAAACAGCTGCGTACCAACGTACCGCACATCGTTGCTATCGGCGATATCGTCGGTCAACCGATGCTGGCACACAAAg                                                                 >  5:18885/1‑149 (MQ=255)                                                                    gCCGAACGGTAAAAACCTCGACGCAGGCAAAGCAGGCGTGGAAGTTGACGCCCGTGGTTTCATCCGCGTTAACAAACAGCTGCGTACCACCGTACCGCACATCTTTGCTATCGGCGATAGCGTCGGTCCACCGATGCTGGCACACAAAg                                                                 >  5:24119/1‑149 (MQ=255)                                                                     ccGAACGGTAAAAACCTCGACGCAGGCAAAGCAGGCGTGGAAGTTGACGACCGTGGTTTCATCCGCGTTGACAAACAGCTGCGTACCAACGTACCGCACATCTTTGCTATCGGCGATATCGTCGGTCAACCGATGCTGGCACACAAAg                                                                 <  7:309888/148‑1 (MQ=255)                                                                     ccGAACGGTAAAAACCTCGACGCAGGCAAAGCAGGCGTGGAAGTTGACGACCGTGGTTTCATCCGCGTTGACAAACAGCTGCGTACCAACGTACCGCACATCTTTGCTATCGGCGATATCGTCGGTCAACCGATGCTGGCACACAAAg                                                                 <  5:58378/148‑1 (MQ=255)                                                                     ccGAACGGTAAAAACCTCGACGCAGGCAAAGCAGGCGTGGAAGTTGACGACCGTGGTTTCATCAGCGTTGACAAACAGCTGCGTACCAACGTACCGCACATCTTTGCTATCGGCGATATCGTCGGTCAACCGATGCTGGCACACAAAgg                                                                >  5:90046/1‑149 (MQ=255)                                                                      cGAACGGTAAAAACCTCGACGCAGGCAAAGCAGGCGTGGAAGTTGACGACCGTGGTTTCATCCGCGTTGACAAACAGCTGCGTACCAACGTACCGCACATCTTTGCTATCGGCGATATCGTCGGTCAACCGATGCTGGCACACAAAGgt                                                               >  3:146935/1‑149 (MQ=255)                                                                         aCGGTAAAACCCTCGACGCAGGCAAAGCCGGCGTGGAAGTTGACGACCGTGGTTTCATCCGCGTTAACAAACAGCTGCGTACCAACGTACCGCACATCTTTGCTAGCGGCGATATCGTCGGTCAACCGATGCTGGCACACAAAGGTGtt                                                            <  6:24119/149‑1 (MQ=255)                                                                         aCGGTAAAAACCTCGACGCCGGCAAAGCAGGAGTGGAAGTTGACGACCGTGGTTTCATCCGCGTTAACAAACAGCTGCGTACCAACGTACCGCACATCTTTGCTATCGGCGATATCGTCGGTCAACCGATGCTGGCACACAAAGGTGtt                                                            <  6:18885/149‑1 (MQ=255)                                                                          cGGTAAAAACCTCGACTCAGGCAAAGCAGGCGTGGAAGTTGACGACCGTGGTTTCATCCGCGTTGACAAACAGCTGCGTACCAACGTACCGCACATCTTTGCTATCGGCGATATCGTCGGTCAACCGATGCTGGCACACAAAGGTGTTc                                                           <  7:362738/149‑1 (MQ=255)                                                                           ggTAAAAACCTCGACGCAGGCAAAGCAGGCGTGGAAGTTGACGACCGTGGTTTCATCCGCGTTGACAAACAGCTGCGTACCAACGTCCCGCACATCTTTGCTATCGGCGATATCGTCGGTCAACCGATGCTGGCACACAAAGGTGTTCa                                                          >  5:319469/1‑149 (MQ=255)                                                                           ggTAAAAACCTCGACGCAGGCAAAGCAGGCGTGGAAGTTGACGACCGTGGTTTCATCCGCGTTGACAAACAGCTGCGTACCAACGTACCGCACATCTTTGCTATCGGCGATATCGTCGGTCAACCGATGCTGGCACACAAAGGTGTTCa                                                          <  4:64415/149‑1 (MQ=255)                                                                            gTAAAAACCTCGACGCAGGCAAAGCAGGCGTGGAAGTTGACGACCGTGGTTTCATCCGCGCTGACAAACAGCTGCGTACCAACGTACCGCACATCTTTGCTATCGGCGATATCGTCGGTCAACCGATGCTGGCACACAAAGGTGTTCAc                                                         <  3:341310/149‑1 (MQ=255)                                                                              aaaaaCCTCGACGCAGGCAAAGCAGGCGTGGAAGTTGACGACCGTGGTTTCATCCGCGTTGACAAACAGCTGCGTACCAACGTACCGCACATCTTTGCTATCGGCGATATCGTCGGTCAACCGATGCTGGCACACAAAGGTGTTCACg                                                        >  5:242799/1‑148 (MQ=255)                                                                               aaaaCCTCGACGCAGGCAAAGCAGGCGTGGAAGTTGACGACCGTGGTTTCATCCGCGTTGACAAAc                                                                                                                                         >  6:216018/1‑66 (MQ=255)                                                                               aaaaCCTCGACGCAGGCAAAGCAGGCGTGGAAGTTGACGACCGTGGTTTCATCCGCGTTGACAAAc                                                                                                                                         <  5:216018/66‑1 (MQ=255)                                                                                aatccTCGACGCAGGCAAAGCAGGCGTGGAAGTTGACGACCGTGGTTTCATCCGCGTTGACAAACAGCTGCGTACCATCGTACCGCACATCttt                                                                                                            <  2:249525/91‑1 (MQ=255)                                                                                aaaaCTCGACGCAGGCAAAGCAGGATTGGAAGTTGACGACAGTGGTTTCATCCGCGTTGACAAACAGCTGCGTACCAACGTACCGCACATCttt                                                                                                            >  1:249525/1‑94 (MQ=255)                                                                                aaaCCTCGACGCAGGCAAAGCAGGCGTGGAAGTTGACGACCGTGGTTTCATCCGCGTTGACAAACAGCTGCGTACCAACGTACCGCACATCttt                                                                                                            >  1:355688/1‑94 (MQ=255)                                                                                aaaCCTCGACGCAGGCAAAGCAGGCGTGGAAGTTGACGACCGTGGTTTCATCCGCGTTGACAAACAGCTGCGTACCAACGTACCGCACATCttt                                                                                                            <  2:355688/94‑1 (MQ=255)                                                                                 aaCCTCGACGCAGGCAAAGCAGGCGTGGAAGTTGACGACCGTGGTTTCATCCGCGTTGACAAACAGCTGCGTACCAACGTACCGCACATCTTTGCTATCGGCGATATCGTCGGTCAACCGATGCTGGCACACAAAGGTGTTCACGAAg                                                     <  5:236954/148‑1 (MQ=255)                                                                                 aaCCTCGACGCAGGCAAAGCAGGCGTGGAAGTTGACGACCGTGGTTTCATCCGCGTTGACAAACAGCTGCGTACCAACGTACCGCACATCTTTGCTATCGGCGATATCGTCGGTCAACCGATGCTGGCACACAAAGGTGTTCACGAAg                                                     >  6:236954/1‑148 (MQ=255)                                                                                   ccTCGACGCAGGCAAAGCAGGCGTGGAAGTTGACGACCGTGGTTTCATCCGCGTTGACAAACAGCTGCGTACCAACGTACCGCACATCTTTGCTATCGGCGATATCGTCGGTCAACCGATGCTGGCACACAAAGGTGTTCACGAAGGt                                                   >  5:194895/1‑148 (MQ=255)                                                                                    cTCGACGCAGGCAAAGCAGGCGTGGAAGTTGACGACCGTGGTTTCATCCGCGTTGACAAACAGCTGCGTACCAACGTACCGCACATCTTTGCTATCGGCGatat                                                                                              <  8:22727/104‑1 (MQ=255)                                                                                    cTCGACGCAGGCAAAGCAGGCGTGGAAGTTGACGACCGTGGTTTCATCCGCGTTGACAAACAGCTGCGTACCAACGTACCGCACATCTTTGCTATCGGCGatat                                                                                              >  7:22727/1‑104 (MQ=255)                                                                                    cTCGACGCAGGCAAAGCAGGCGTGGAAGTTGACGACCGTGGTTTCATCCGCGTTGACAAACAGCTGCGTACCAACGTACCGCACATCTTTGCTATCGGCGATATCGTCGGTCAACCGATGCTGGCACACAAAGGTGTTCACGAAGGTCa                                                 >  1:266724/1‑149 (MQ=255)                                                                                        aCGCAGGCAAAGCAGGCGTGGAAGTTGACGACCGTGGTTTCATCCGCGTTGACAAACAGCTGCGTACCAACGTACCGCACATCTTTGCTATCGGCGATATCGTCGGTCAACCGATGCTGGCACACAAAGGTGTTCACGAAGGTCACGtt                                             >  2:43812/1‑149 (MQ=255)                                                                                              gCAAAGCAGGCGTGGAAGTTGACGACCGTGGTTTCATCCGCGTTGACAAACAGCTGCGTACCAACGTACCGCACATCTTTGCTATCGGCGATATCGTCGGTCAACCGATGCTGGCACACAAAGGTGTTCACGAAGGTCAc                                                <  1:153914/140‑1 (MQ=255)                                                                                              gCAAAGCAGGCGTGGAAGTTGACGACCGTGGTTTCATCCGCGTTGACAAACAGCTGCGTACCAACGTACCGCACATCTTTGCTATCGGCGATATCGTCGGTCAACCGATGCTGGCACACAAAGGTGTTCACGAAGGTCAc                                                >  2:153914/1‑140 (MQ=255)                                                                                                aaaGCAGGCGTGGAAGTTGACGACCGTGGTTTCATCCGCGTTGACCAACAGCTGCGTACCAACGTCCCGCACATCTTTGCTATCGGCGATATCGTCGGTCAACCGCTGCTGGCCCACAACGGTGTTCACGAAGGTCACGTGGCCGCTg                                      >  5:117195/1‑148 (MQ=255)                                                                                                  aGCAGGCGTGGAAGTTGACGACCGTGGTTTCATCCGCGTTGACAAACAGCTGCGTACCAACGTACCGCACATCTTTGCTATCGGCGATATCGTCGGTCAACCGATGCTGGCACACAAAGGTGTTCACGAAGGTCACGTTGCCGCTGAAg                                   >  5:212945/1‑149 (MQ=255)                                                                                                   gaaGGCTTGGAAGTTGACGACCGTCGTTTCAACCGCGTTGACAAACAGCTGCGTACCAACGTACCGCACATCTTTGCAATCGTCGATATCGTCGGTCAACCGATGCTGGCACACAAAGGTGTTc                                                           >  2:64533/3‑124 (MQ=255)                                                                                                   gCAGGCGTTGAATTTGACGACCGTGGTTTCATCCGAGTTTACAATCAGCTGCGTACCAACGTACCGCACATCTTTGCTATCGGCGATATCGTCGGTCAAGAGTTGCTGGCAGACATTGGTGTTc                                                           <  1:64533/124‑1 (MQ=255)                                                                                                         gTGGAAGTTGACGACCGTGGTTTCATCCGCGTTGACAAACAGCTGCGTACCAACGTACCGCACATCTTTGCTATCGGCGATATCGTCGGTCAACCGATGCTGGCACACAAAGGTGTTCACGAAGGTCACGTTGCCGCTGAAGTTATCGc                            <  3:161452/149‑1 (MQ=255)                                                                                                            gAAGTTGACGACCGTGGTTTCATCCGCGTTGACAAACAGCTGCGTACCAACGTACCGCACATCTTTGCTATCGGCGATATCGTCGGTCAACCGATGCGGGCACACAAAGGTGTTCACGAAGGTCACGTTGCCGCTGAAGTTATCGCCgg                         <  1:273968/149‑1 (MQ=255)                                                                                                                  gacgacCGTGGTTTCATCCGCGTTGACAAACAGCTGCGTACCAACGTACCGCACATCTTTGCTATCGGCGATATCGTCGGTCAACCGATGCTGGCACACAAAGGTGTTCACGAAGGTCACGTTGCCGCTGa                                     >  7:301360/1‑131 (MQ=255)                                                                                                                  gacgacCGTGGTTTCATCCGCGTTGACAAACAGCTGCGTACCAACGTACCGCACATCTTTGCTATCGGCGATATCGTCGGTCAACCGATGCTGGCACACAAAGGTGTTCACGAAGGTCACGTTGCCGCTGa                                     <  8:301360/131‑1 (MQ=255)                                                                                                                     gacCGTGGTTTCATCCGCGTTGACAAACAGCTGCGTACCAACGTACCGCACATCTTTGCTATCGGCGATATCGTCGGTCAACCGATGCTGGCACACAAAGGTGTTCACGAAGGTCACTTTGCCGCTGAAGTTATCGCCGGTAAGAAaca                >  4:250844/1‑149 (MQ=255)                                                                                                                      acCGTGGTTTCATCCGCGTTGACAAACAGCTGCGTACCAACGTACCGCACATCTTTGCTATCGGCGATATCGTCGGTCAACCGATGCTGGCACACAAAGGTGTTCACGAAGGTCACGTTGCCGCTGAAGTTATCGCCGGTAAGAAaca                >  8:56759/1‑148 (MQ=255)                                                                                                                        cGTGGTTTCATCCGCGTTGACAAACAGCTGCGTACCAACGTACCGCACATCTTTGCTATCGGCGATATCGTCGGTCAACCGATGCTGGCACACAAAGGTGTTCACGAAGGTCACGTTGCCGCTGAAGTTATCGCCGGTAAGAAACacta             >  6:168684/1‑149 (MQ=255)                                                                                                                        cGTGGTTTCATCCGCGTTGACAAACAGCTGCGTACCAACGTACCGCACATCTTTGCTATCGGCGATATCGTCGGTCAACCGATGCTGGCACACAAAGGTGTTCACGAAGGTCACGTTGCCGCTGAAGTTATCGCCGGTAAGAAACacta             <  6:242799/149‑1 (MQ=255)                                                                                                                            gTTTCATCCGCGTTGACAAACAGCTGCGTACCAACGTACCGCACATCTTTGCTATCGGCGATATCGTCGGTCAACCGATGCTGGCACACAAAGGTGTTCACGAAGGTCACGTTGCCGCTGAAGTTATCGCCGGTAAGAAACACTACtt          <  4:146935/148‑1 (MQ=255)                                                                                                                                    cgcgTTGACAAACAGCTGCGTACCAACGTACCGCACATCTTTGCTATCGGCGATATCGTCGGTCAACCGATGCTGGCACACAAAGGTGTTCACGAAGGTCACGTTGCCGCTGAAGTTATCGCCGGTAAGAAACACTACTTCGATCCGa  <  6:334838/148‑1 (MQ=255)                                                                                                                                     gcgTTGACAAACAGCTGCGTACCAACGTACCGCACATCTTTGCTATCGGCGATATCGTCGGTCAACCGATGCTGGCACACAAAGGTGTTCACGAAGGTCACGTTGCCGCTGAAGTTATCGCCGGTAAGAAACACTACTTCGATCCGa  >  2:250564/1‑147 (MQ=255)                                                                                                                                          |                                                                                                                                              TGGAAGGCAAAAAAGCACCCGCTGAACCGCAGCGTTACGACGCCGTGCTGGTAGCGATTGGTCGTGTGCCGAACGGTAAAAACCTCGACGCAGGCAAAGCAGGCGTGGAAGTTGACGACCGTGGTTTCATCCGCGTTGACAAACAGCTGCGTACCAACGTACCGCACATCTTTGCTATCGGCGATATCGTCGGTCAACCGATGCTGGCACACAAAGGTGTTCACGAAGGTCACGTTGCCGCTGAAGTTATCGCCGGTAAGAAACACTACTTCGATCCGA  >  NC_000913/128666‑128944

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