Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	139,246	A→G	5.6%	F660F (TTT→TTC)	gcd ←	glucose dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	139,246	0	A	G	5.6%	222.2 / 2.9	71	F660F (TTT→TTC)	gcd	glucose dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (38/29);  new base G (2/2);  total (40/31)
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.98e-01

CATACCCATATTGAACGGCACCGGAACCGGCATCGGGAACGGCATACTGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCATTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCCTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAAC  >  NC_000913/139099‑139390                                                                                                                                                    |                                                                                                                                                 cATACCCATATTGAACGGCACCGGAACCGGCATCGGGAACGGCATACTGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCATTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCa                                                                                                                                                  >  5:10376/1‑148 (MQ=255)    aCCCATATTGAACGGCACCGGAACCGGCATCGGGAACGGCATACTGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCATTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTAGCTGTTTACATGGCAGACCAAATg                                                                                                                                              >  3:171443/1‑149 (MQ=255)      ccATATTGAACGGCACCGGAACCGGCATCGGGAACGGCATACTGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCATTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGt                                                                                                                                            >  6:319814/1‑149 (MQ=255)         tatTGAACGGCACCGGAACCGGCATCGGGAACGGCATACTGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCATTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGt                                                                                                                                            >  8:297747/1‑146 (MQ=255)          atTGAACGGCACCGGAACCGGCATCGGGAACGGCATACTGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCATTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTga                                                                                                                                          >  6:355573/1‑147 (MQ=255)           ttGAACGGCACCGGAACCGGCATCGGGAACGGCATACTGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCATTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGagag                                                                                                                                       <  7:374862/149‑1 (MQ=255)             gAACGGCACCGGAACCGGCATCGGGAACGGCATACTGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACGTCATTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAgg                                                                                                                                      >  4:324103/1‑148 (MQ=255)              aaCGGCACCGGAACCGGCATCGGGAACGGCATACTGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCATTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTAGCTGTTTACATGGCAGACCAAATGGTGAGAGGaa                                                                                                                                    <  4:171443/149‑1 (MQ=255)              aaCGGCACCGGAACCGGCATCGGGAACGGCATACTGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCATTAGTTTTCAGATCAAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGCTCAGAGGaa                                                                                                                                    >  8:313559/1‑149 (MQ=255)                   cACCGGAACCGGCATCGGGAACGGCATACTGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCATTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGt                                                                                                                               >  1:263619/1‑149 (MQ=255)                    aCCGGAACCGGCATCGGGAACGGCATACTGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCATTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGtt                                                                                                                              <  3:76242/149‑1 (MQ=255)                    aCCGGAACCGGCATCGGGAACGGCATACTGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCATTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGtt                                                                                                                              >  2:86902/1‑149 (MQ=255)                       ggAACCGGCATCGGGAACGGCATACTGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCATTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTCAGAGGAACGGGTTGAg                                                                                                                           >  8:384708/1‑149 (MQ=255)                              gCATCGGGAACGGCATACTGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCATTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGaca                                                                                                                    <  6:10376/149‑1 (MQ=255)                                aTCGGGAACGGCATACTGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCATTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACAcc                                                                                                                  <  1:86902/149‑1 (MQ=255)                                   gggAACGGCATACTGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCATTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATa                                                                                                               <  6:238028/149‑1 (MQ=255)                                   gggAACGGCATACTGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCATTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATa                                                                                                               <  5:319814/149‑1 (MQ=255)                                    ggAACGGCATACTGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCATTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATAc                                                                                                              >  4:170278/1‑149 (MQ=255)                                    ggAACGGCATACTGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCATTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATAc                                                                                                              >  7:183146/1‑149 (MQ=255)                                     gAACGGCATACTGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCATTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACg                                                                                                             <  8:183146/149‑1 (MQ=255)                                         ggCATACTGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCATTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATTCTGGCTGTTTACATGGCAGACAAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTa                                                                                                          >  2:379867/1‑148 (MQ=255)                                             tACTGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCATTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAAc                                                                                                                                   <  6:373702/119‑1 (MQ=255)                                             tACTGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCATTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAAc                                                                                                                                   >  5:373702/1‑119 (MQ=255)                                                tGTCCTGCGGCGTACCTATACGTTTCTTCCACACCACTTCATTAGTTTTCTGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAAACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTAc                                                                                                  >  4:277879/1‑149 (MQ=255)                                                  tCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCATTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCGAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTa                                                                                                   <  6:211888/146‑1 (MQ=255)                                                  tCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCATTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCGAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTa                                                                                                   <  6:217245/146‑1 (MQ=255)                                                  tCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCATTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCGAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTa                                                                                                   >  5:211888/1‑146 (MQ=255)                                                  tCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCATTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCGAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTa                                                                                                   >  5:217245/1‑146 (MQ=255)                                                           gTACCAATACGTTTCTTCCACACCACTTCATTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAAt                                                                                       <  5:355573/149‑1 (MQ=255)                                                              ccAATACGTTTCTTCCACACCACTTCATTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGcc                                                                                    <  7:297747/149‑1 (MQ=255)                                                               cAATACGTTTCTTCCACACCACTTCATTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGaa                                                                                        <  1:210692/144‑1 (MQ=255)                                                               cAATACGTTTCTTCCACACCACTTCATTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGaa                                                                                        >  2:210692/1‑144 (MQ=255)                                                                aaTACGTTTCTTCCACACCACTTCATTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCg                                                                                   >  1:133751/1‑148 (MQ=255)                                                                 aTACGTTTCTTCCACACCACTTCATTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGa                                                                                 >  7:232711/1‑149 (MQ=255)                                                                  tACGTTTCGTCGACACCACTTGATTAGTTTTCAGATCCAGCGCCGAGATATAGCCCGAGGGTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGACGGGGTTGAGCGGGACGCCATACGGTACGCCGTGGGGTGGCGGAGTGCCgggt                                                                                >  6:323842/1‑147 (MQ=255)                                                                      tttCTTCCACACCACTTCATTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGcc                                                                                    <  1:309909/141‑1 (MQ=255)                                                                      tttCTTCCACACCACTTCATTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGcc                                                                                    >  2:309909/1‑141 (MQ=255)                                                                            ccACACCACTTCATTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTAc                                                                         <  3:324103/146‑1 (MQ=255)                                                                              acacCACTTCATTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTg                                                                    >  3:183152/1‑149 (MQ=255)                                                                                  cacTTCATTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCCtt                                                                >  1:74685/1‑149 (MQ=255)                                                                                       cATTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGTGGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCCTTTGGc                                                            <  7:149081/148‑1 (MQ=255)                                                                                       cATTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCCTTTGGc                                                            <  2:133751/148‑1 (MQ=255)                                                                                         ttAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCCTTTGGCATc                                                         >  4:277019/1‑149 (MQ=255)                                                                                           aGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTcc                                                                             >  1:383669/1‑127 (MQ=255)                                                                                           aGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTcc                                                                             <  2:383669/127‑1 (MQ=255)                                                                                            gTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCCTTTGGCATCttt                                                      <  2:263619/149‑1 (MQ=255)                                                                                             ttttCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTcc                                                                             >  1:282454/1‑125 (MQ=255)                                                                                             ttttCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTcc                                                                             >  1:288123/1‑125 (MQ=255)                                                                                             ttttCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTcc                                                                             <  2:282454/125‑1 (MQ=255)                                                                                             ttttCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTcc                                                                             <  2:288123/125‑1 (MQ=255)                                                                                             ttttCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCCTTTGGCATCTTTc                                                     <  8:232711/149‑1 (MQ=255)                                                                                              tttCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTgg                                                                                             <  2:108268/108‑1 (MQ=255)                                                                                              tttCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTgg                                                                                             >  1:108268/1‑108 (MQ=255)                                                                                               ttCAGCTCCTGCGCCGAGATATAACCCCATGCTGGCTGTTTTCATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCg                                                                      <  6:128797/130‑1 (MQ=255)                                                                                               ttCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCg                                                                      >  5:128797/1‑130 (MQ=255)                                                                                               ttCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCCTTTGGCATCTTTCgg                                                   >  7:265636/1‑149 (MQ=255)                                                                                                  aGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCATACCCGTGCCTTTGGCATCTTTCGgctg                                                >  2:47197/1‑149 (MQ=255)                                                                                                   gATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCCTTTGGCATCTTTCGgctg                                                >  1:247016/1‑148 (MQ=255)                                                                                                   gATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCCTTTGGCATCTTTCGgctg                                                <  2:247016/148‑1 (MQ=255)                                                                                                           gcCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCCTTTGGCATCTTTCGGCTGCTCCATCgg                                       >  1:133471/1‑149 (MQ=255)                                                                                                                       ccccATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCCTTTGGCAt                                                          >  3:104541/1‑118 (MQ=255)                                                                                                                       ccccATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCCTTTGGCAt                                                          <  4:104541/118‑1 (MQ=255)                                                                                                                       ccccATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCCTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGAcc                           >  8:227230/1‑149 (MQ=255)                                                                                                                               gggctgTTTACCTGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGGGACACCATACGGTACACCGTACTGTGGCGGAATGCCGGATTCCGTACCCGTGCCTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGGCCACGCGGGa                   <  5:323842/148‑1 (MQ=255)                                                                                                                                        aCATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCCTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTc           >  3:126390/1‑148 (MQ=255)                                                                                                                                         cATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCCTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGa         <  8:265636/149‑1 (MQ=255)                                                                                                                                          aTGGCAGCCCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCCTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGa         >  7:362357/1‑148 (MQ=255)                                                                                                                                          aTGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCCTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGa         >  1:107195/1‑148 (MQ=255)                                                                                                                                               aGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCCTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGaaacaaa   >  8:212946/1‑149 (MQ=255)                                                                                                                                                aaCCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCCTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGaaacaaac  <  3:277879/148‑1 (MQ=255)                                                                                                                                                   cAAATGTTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCCTTTGGc                                                            >  8:368250/1‑88 (MQ=255)                                                                                                                                                   cAAATGGTGATAGGAACGGGTTGAGCGTGACACCATACGTTACACCGTACTGTGTCTTAATGCCGGATTCCGTACCCGTGCCTTTGGc                                                            <  7:368250/88‑1 (MQ=255)                                                                                                                                                    aaaTGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCCTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGAc                            >  7:263491/1‑119 (MQ=255)                                                                                                                                                    aaaTGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCCTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGAc                            <  8:263491/119‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                 CATACCCATATTGAACGGCACCGGAACCGGCATCGGGAACGGCATACTGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCATTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCCTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAAC  >  NC_000913/139099‑139390

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