Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_012971 767,889 T→C M60T (ATG→ACG)  ybcW → protein YbcW

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_012971767,8890TC100.0% 224.8 / NA 60M60T (ATG→ACG) ybcWprotein YbcW
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (28/32);  total (28/32)

AATCTGCTGATGATCCCTCCGTGGATCTGATTCGTGTAAAAAATATGCTTAATAGCACCATTTCTATGAGTTACCCTGATGTTGTAATTGCATGTATAGAACATAAGGTGTCTCTGGAAGCATTCAGAGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTATTTAGTCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAACGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGTCTCCGACGGCAGGCTTCAATGACCCAGGCTGAGAAATTCCCGGACC  >  NC_012971/767724‑768179
                                                                                                                                                                     |                                                                                                                                                                                                                                                                                                  
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)
                                                                                                    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)
                                                                                                                             atggcaaGTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctacccacgtaatatggacacaggcctaagcgaggttcttgttttcaaattgttccggactgaggccgccacaccaactgtgccgccgccaccgattgtaatc                                 <  1:764687‑M1/297‑210 (MQ=255)
                                                                                                                             agggcaaTTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctacccacgtaatatggacacaggcctaagcgaggttcttgttttcaaattgttccggactgaggccgccacaccaactgtgccgccgccaccgattgtaatc                                 <  1:764681‑M1/297‑210 (MQ=255)
                                                                                                                              gggcaATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctagtatttaggtgaattgatcctacccacgtaatatggacacaggcctaagcgaggttcttgttttcaaattgttccggactgaggccgccaaaccaactgtgccgccgccaccgattgtaacc                                 >  2:489135‑M1/3‑90 (MQ=255)
                                                                                                                              gggcaATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgccaggctaatacagcaatgcaattcaactacagcaatgcctagtatttaggtgaattgatcctacccacgtaatatggacacaggcctaagcgaggttcttgttttcaaattgttccggactgaggccgccacaccaactgtgccgcagccaccgattgtaatca                                >  2:319209‑M1/3‑90 (MQ=255)
                                                                                                                                  aaTTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaattgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctacccacgtaatatggacacaggcctaagcgaggttcttgttttcaaattgttccggactgaggccgccacaccaactgtgccgccgccaccgattgtaatcaca                              <  2:297970‑M1/298‑213 (MQ=255)
                                                                                                                                   aTTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctacccacgtaatatggacacaggcctaagcgaggttcttgttttcaaattgttccggactgaggccgccacaccaactgtgccgccgccaccgattgtaatcacattc                           <  1:394485‑M1/300‑216 (MQ=255)
                                                                                                                                                tagtaaagcacgATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctacccacgtaatatggacacaggcctaagcgaggttcttgttttcaaattgttccggactgaggccgccacaccaactgtgccgccgccaccgattgtaatcacattcgatataattaaac              <  1:498956‑M1/295‑229 (MQ=255)
                                                                                                                                                 agtaaagcacGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctacccacgtaatatggacacaggcctaagcgaggttcttgttttcaaattgttccggactgaggccgccacaccaactgtgccgccgccaccgattgtaatcacattcgatataattaaaca             >  1:129672‑M1/5‑71 (MQ=255)
                                                                                                                                                 agtaaagcacGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaaggcaattcaactacagcaatgcctcgtatttaggtgaattgatcctacccacgtaatatggacacaggcctaagcgaggttcttgttttcaaattgttccggactgaggccgccacaccaactgtgccgccgccaccgattgtaatcacattcgatataattaaac              <  2:727630‑M1/295‑229 (MQ=255)
                                                                                                                                                  gtaaagcaCGATAAGAATACGAAGGATTATTCCTTGGTGGTTGACTGATCACTAGCACTGCTAACCATTCtgaccatttcacaattgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctacccacgtaatatggacacaggcctaagcgaggttcttgttttcaaattgttccggaccgagaccgccacaccaactgtgccgccaccaccaattgtaatcacattcgatataattacacac            <  2:513042‑M1/297‑231 (MQ=255)
                                                                                                                                                  gtaaagcaCGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTGACCATTCtgaccatttcacaagtgacagagtcaatatcggattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctagtatttaggtgaattgatcctacccacgtaatatggacacaggcctaagcgaggttcttgttttcaaattgttccggactgaggccgccacaccaactgtgccgccgccaccgattgtaatcacattcgatataattaaacac            >  2:704316‑M1/4‑70 (MQ=255)
                                                                                                                                                  gtaaagcaCGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctacccacgtaatatggacacaggcctaagcgaggttcttgttttcaaattgttccggactgaggccgccacaccaactgtgccgccgccaccgattgtaatcacattcgatataattaaacac            <  1:788631‑M1/297‑231 (MQ=255)
                                                                                                                                                            ataagaataCGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctacccacgtaatatggacacaggcctaagcgaggttcttgttttcaaattgttccggactgaggccgccacaccaactgtgccgccgccaccgattgtaatcacattcgatataattaaacaccgttgcccgc  >  2:795923‑M1/1‑60 (MQ=255)
                                                                                                                                                                     |                                                                                                                                                                                                                                                                                                  
AATCTGCTGATGATCCCTCCGTGGATCTGATTCGTGTAAAAAATATGCTTAATAGCACCATTTCTATGAGTTACCCTGATGTTGTAATTGCATGTATAGAACATAAGGTGTCTCTGGAAGCATTCAGAGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTATTTAGTCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAACGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGTCTCCGACGGCAGGCTTCAATGACCCAGGCTGAGAAATTCCCGGACC  >  NC_012971/767724‑768179

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: