Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_012971	767,828	A→C	K40Q (AAG→CAA)	ybcW →	protein YbcW

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_012971	767,828	0	A	C	100.0%	138.1 / NA	38	K40Q (AAG→CAA)	ybcW	protein YbcW
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (20/18);  total (20/18)

AATCTGCTGATGATCCCTCCGTGGATCTGATTCGTGTAAAAAATATGCTTAATAGCACCATTTCTATGAGTTACCCTGATGTTGTAATTGCATGTATAGAACATAAGGTGTCTCTGGAAGCATTCAGAGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTATTTAGTCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAACGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCT  >  NC_012971/767724‑768119                                                                                                         |                                                                                                                                                                                                                                                                                                    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                                                             <  1:304887‑M1/300‑122 (MQ=255)                                      aaaaacTATGCTTAATAGCACCATTTCGATGGGTTACCCGGATGTTGTAATTGCATGTATAGAACATCAAGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctacccacgtaatat                                                             <  2:72176‑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)                                                                                                         |                                                                                                                                                                                                                                                                                                    AATCTGCTGATGATCCCTCCGTGGATCTGATTCGTGTAAAAAATATGCTTAATAGCACCATTTCTATGAGTTACCCTGATGTTGTAATTGCATGTATAGAACATAAGGTGTCTCTGGAAGCATTCAGAGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTATTTAGTCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAACGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCT  >  NC_012971/767724‑768119

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