Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_012971	767,851	A→G	R47R (AGA→AGG)	ybcW →	protein YbcW

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_012971	767,851	0	A	G	97.9%	166.3 / NA	47	R47R (AGA→AGG)	ybcW	protein YbcW
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (23/23);  minor base T (0/1);  total (23/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATCTGCTGATGATCCCTCCGTGGATCTGATTCGTGTAAAAAATATGCTTAATAGCACCATTTCTATGAGTTACCCTGATGTTGTAATTGCATGTATAGAACATAAGGTGTCTCTGGAAGCATTCAGAGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTATTTAGTCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAACGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGTC  >  NC_012971/767724‑768134                                                                                                                                |                                                                                                                                                                                                                                                                                            aatacactgatgccccctcccaggagcaggttcgtgtaaaaactatgcttaatagcaccattTCGATGGGTTACCCGGATGTTGTAATTGCATGTATAGAACATCAAGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgc                                                                                                                 <  1:23631‑M1/270‑85 (MQ=255)       ctgctgccccctcccaggagcaggttcgtgtaaaaactatgcttaatagcACCATTTCGATGGGTTACCCGGATGTTGTAATTGCATGTATAGAACATCAAGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatac                                                                                                                                            >  2:313320‑M1/25‑210 (MQ=255)       ctgatgccccctcccaggagcaggttcgtgtaaaaactatgcttaatagcACCATTTCGATGGGTTACCCGGATGTTGTAATTGCATGTATAGAACATCAAGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatac                                                                                                                                            <  1:313320‑M1/243‑58 (MQ=255)                cctcccaggagcaggttcgtgtaaaaactatgCTTAATAGCACCATTTCGATGGGTTACCCGGATGTTGTAATTGCATGTATAGAACATCAAGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtga                                                                                                  >  1:65392‑M1/16‑201 (MQ=255)                cctaccaggagcaggttcgtgtaaaaactatgCTTAATAGCACCATTTCGATGGGTTACCCGGATGTTGTAATTGCATGTATAGAACATCAAGTGTCTATGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtga                                                                                                  >  1:65418‑M1/16‑201 (MQ=255)                  tcccaggagcaggttcgtgtaaaaactaTGCTTAATAGCACCATTTCGATGGGTTACCCGGATGTTGTAATTGCATGTATAGAACATCAAGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaat                                                                                                >  2:465998‑M1/14‑199 (MQ=255)                   cccaggagcaggttcgtgtaaaaactATGCTTAATAGCACCATTTCGATGGGTTACCAGGATGTTGTAATTGCATGTATAGAACATCAAGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTCACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctagtatttaggtgaatt                                                                                               >  2:13822‑M1/13‑198 (MQ=255)                     caggagcaggttcgtgtaaaaaCTATGCTTAATAGCACCATTTCGATGGGTTACCCGGATGTTGTAATTGCATGTATAGAACATCAAGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattga                                                                                             >  1:604525‑M1/11‑196 (MQ=255)                      aggagtaggttcgtgtaaaaATTATGCTTAATAGCACCATTTCGATGGGTTACCCGGATGTTGTAATTGCATGTATAGAACATCAAGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtattaatac                                                                                                    <  2:182119‑M1/283‑98 (MQ=255)                        gagcaggttcgtgtaaAAACTATGCTTAATAGCACCATTTCGATGGGTTACCCGGATGTTGTAATTGCATGTATAGAACATCAAGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcc                                                                                          >  1:506272‑M1/8‑193 (MQ=255)                        gagcaggttcgtgtaaAAACTATGCTTAATAGCACCATTTCGATGGGTTACCCGGATGTTGTAATTGCATGTATAGAACATCAAGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatc                                                                                           <  2:13043‑M1/292‑107 (MQ=255)                                      aaaaacTATGCTTAATAGCACCATTTCGATGGGTTACCCGGATGTTGTAATTGCATGTATAGAACATCAAGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctacccacgtaatat                                                                            <  2:72176‑M1/300‑122 (MQ=255)                                      aaaaacTATGCTTAATAGCACCATTTCGATGGGTTACCCGGATGTTGTAATTGCATGTATAGAACATCAAGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctacccacgtaatat                                                                            <  1:304887‑M1/300‑122 (MQ=255)                                          actatgcTTAATAGCACCATTTCGATGGGTTACCCGGATGTTGTAATTGCATGTATAGAACATCAAGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctacccacgtaatatggac                                                                        <  1:254657‑M1/298‑126 (MQ=255)                                            tatGCTTAATAGCACCATTTCGATTGGTTACCCGGATGTTGTAATTGCATGTATAGAACATCAATTGTCTCTGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaattgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctacccacgtaatatg                                                                           <  2:347316‑M1/295‑123 (MQ=255)                                            tatGCTTAATAGCACCATTTCGATGGGTTACCCGGATGTTGTAATTGCATGTATAGAACATCAAGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctacccacgtaatatggacac                                                                      <  1:552435‑M1/300‑128 (MQ=255)                                            tatGCTTAATAGCACCATTTCGATGGGTTACCCGGATGTTGTAATTGCATGTATAGAACATCAAGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCNAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctacccacgtaatatggacac                                                                      <  2:580525‑M1/300‑128 (MQ=255)                                               gCTTAATAGCACCATTTCGATGGGTTACCCGGATGTTGTAATTGCATGTATAGAACATCAAGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCCCGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttaccaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattccactacagcaatgcctcgtacttaggtgaattgatcctactcacgtagta                                                                             >  2:7519‑M1/1‑170 (MQ=255)                                                  tAATAGCACCATTTCGATGGGTTACCCGGATGTTGTAATTGCATGTATAGAACATCAAGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctccccacgtaatatggacacaggcct                                                                >  2:558333‑M1/1‑167 (MQ=255)                                                      aGCACCATTTCGATGGGTTACCCGGATGTTGTAATTGCATGTATAGAACATCAAGTGTCTCTGTAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctacccacgtaatatggacacaggcctaagc                                                            <  1:194524‑M1/300‑138 (MQ=255)                                                      aGCACCATTTCGATGGGTTACCCGGATGTTGTAATTGCATGTATAGAACATCAAGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaattgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctacccacgtaatatggacacaggcctaagc                                                            <  1:194923‑M1/300‑138 (MQ=255)                                                      aGCACCATTTCGATGGGTTACCCGGATGTTGTAATTGCATGTATAGAACATCAAGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctacccacgtaatatggacacaggcctaagc                                                            <  2:794728‑M1/300‑138 (MQ=255)                                                       gCACCATTTCGATGGGTTACCCGGATGTTGTAATTGCATGTATAGAACATCAAGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctacccacgtaatatggacacaggcctaagcg                                                           <  2:567229‑M1/300‑139 (MQ=255)                                                            gtttcGATGGGTTACCCGGATGTTGTAATTGCATGTATAGAACATCAAGTGTCTATGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAAATGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcactgcaattcaactacagccatgcctagtatttaggtgaattgatcctacccacgtcatatggacaccggcctaagcgaggtt                                                      >  1:263121‑M1/2‑157 (MQ=255)                                                            aTTTCGATGGGTTACCCGGATGTTGTAATTGCATGTATAGAACATCAAGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatgcagcaatgcaattcaactaaagcaatgaatagtatttaggtgaattgatcctaaccacgtaatatggatacagacctaagtgaggtt                                                      >  2:21721‑M1/1‑157 (MQ=255)                                                             tttCGATGGGTTACCCGGATGTTGTAATTGCATGTATAGAACATCAAGTGTCTCTGGCAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctacacacgtaatatggacacaggcctaagcgaggtt                                                      >  2:432040‑M1/1‑156 (MQ=255)                                                                    gggttACCCGGATGTTGTAATTGCATGTATAGAACATCAAGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctagtatttaggtgaattgatcctacccacgtaatatggacacaggcctaagcgaggttcttgtttt                                              >  2:692873‑M1/3‑149 (MQ=255)                                                                      gTTACCCGGATGTTGTAATTGCATGTATAGAACATCAAGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctacccacgtaatatggacacaggcctaagcgaggttcttgttttca                                            >  1:45770‑M1/1‑147 (MQ=255)                                                                         aCCCGGATGTTGTAATTGCATGTATAGAACATCAAGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctacccacgtaatatggacacaggcctaagcgaggttcttgttttcaaat                                         <  2:804905‑M1/300‑157 (MQ=255)                                                                                  ttgtAATTGCATGTATAGAACATCAAGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactncagcaatgcctcgtatttaggtgaattgatcctacccacgtaatatggacacaggcctaagcgaggttcttgttttcaaattgttccgga                                >  2:634531‑M1/1‑135 (MQ=255)                                                                                         tGCATGTATAGAACATCAAGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACNAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctacccacgtaatatggacacaggcctaagcgaggttcttgttttcaaattgttccggactgaggc                         <  1:580659‑M1/300‑173 (MQ=255)                                                                                         tGCATGTATAGAACATCAAGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctacccacgtaatatggacacaggcctaagcgaggttcttgttttcaaattgttccggactgaggc                         >  1:539792‑M1/1‑128 (MQ=255)                                                                                              gTATAGAACATCAAGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATNACCATAACTGCTAACCATTCtgaccatttcacnagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctacccacgtaatatggacacaggcctaagcgaggttcttgttttcaaattgttccggactgaggccgcca                    >  2:605588‑M1/1‑123 (MQ=255)                                                                                              gTATAGAACATCAAGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctacccacgtaatatggacacaggcctaagcgaggttcttgttttcaaattgttccggactgaggccgcca                    >  1:467675‑M1/1‑123 (MQ=255)                                                                                              gTATAGAACATCAAGTGTCTATGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctacccaagtaatgtggacacaggcctaagtgaggtttttgttttccaattgttccggactgaggccgcca                    >  2:223532‑M1/1‑123 (MQ=255)                                                                                              gTATAGAACATCAAGTGTCTATGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctaaccacgtaatatggacacaggcctaagcggggttcttgttttcaaattgttccggactgaggccgcca                    >  2:109919‑M1/1‑123 (MQ=255)                                                                                               tagagaaCATCAAGTGTCTCTGGAAGCATTCAGGGCAACTGAGGCCGCGTTAGTAAAGCACGAGAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcanaagtgacagagtcaatatcgtattctgtcactggcaggctaatacagcaatgcaatgcaactacagcaatgcctcgtatgtaggtgaatgggtcctacccacgtaatatggacacaggcctaagcgaggttcttgttttcaaattgttccggactgacgccaccac                   <  2:263121‑M1/297‑179 (MQ=255)                                                                                                 taGAACATCAAGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctacccacgtaatatggacacaggcctaagcgaggttcttgttttcaaattgttccggactgaggccgccacac                 <  1:303931‑M1/300‑181 (MQ=255)                                                                                                     acatcaagtgtctctggaAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctacccacgtaatatggacacaggcctaagcgaggttcttgttttcaaattgttccggactgaggccgtcacaccaac             <  2:748714‑M1/293‑185 (MQ=255)                                                                                                      catcaagtgtctctggAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctacccacgtaatatggacacaggcctaagcgaggttcttgttttcaaattgttccggactgaggccgccacaccaact            >  2:351512‑M1/7‑115 (MQ=255)                                                                                                       atcaagtgtctctgTAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctacccacgtaatatggacacaggcctaagcgaggttcttgttttcaaattgttccggactgaggccgccacaccaactg           <  1:353709‑M1/295‑187 (MQ=255)                                                                                                       atcaagtgtctctgGAAGCATTCAGTGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctacccacgtaatatggacacaggcctaagcgaggttcttgttttcaaattgttccggactgaggccgccacaccaactg           <  2:350354‑M1/295‑187 (MQ=255)                                                                                                          aagtgtctCTGTAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctacccacgtaatatggacacaggcctaagcgaggttcttgttttcaaattgttccggactgaggccgccacaccaactgtgc        <  1:388047‑M1/298‑190 (MQ=255)                                                                                                           agtgtcTCTGGAAGCATTCAGGGCAATTGAGGCCGTTTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaattgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctacccacgtaatatggacacaggcctaagcgaggttcttgttttcaaattgttccggactgaggccgccacaccaactgtgcc       <  2:320864‑M1/299‑191 (MQ=255)                                                                                                           agtgtcTCTGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaacgacagcaatgcctcgtatttaggtgaattgatcctacccacttaatatggacacaggcctaagcgaggttcttgttttcaaattgttccggactgaggccgccacaccaactgtgcc       <  2:20081‑M1/299‑191 (MQ=255)                                                                                                              gtCTCTGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctagtatttaggtgaattgatcctacccacgtaatatggccacaggcctaagcgaggttcttgtcttcaaattgttccggactgaggccgccacaccaactgtgccgcc    >  2:654902‑M1/1‑107 (MQ=255)                                                                                                                ctctGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctacccacgtaatatggacacaggcctaagcgaggttcttgttttcaaattgttccggactgaggccgccacaccaactgtgccgccgc  >  1:126053‑M1/1‑105 (MQ=255)                                                                                                                                |                                                                                                                                                                                                                                                                                            AATCTGCTGATGATCCCTCCGTGGATCTGATTCGTGTAAAAAATATGCTTAATAGCACCATTTCTATGAGTTACCCTGATGTTGTAATTGCATGTATAGAACATAAGGTGTCTCTGGAAGCATTCAGAGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTATTTAGTCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAACGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGTC  >  NC_012971/767724‑768134

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑