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% 254.4 / NA 68M60T (ATG→ACG) ybcWprotein YbcW
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (25/43);  total (25/43)

TGATGATCCCTCCGTGGATCTGATTCGTGTAAAAAATATGCTTAATAGCACCATTTCTATGAGTTACCCTGATGTTGTAATTGCATGTATAGAACATAAGGTGTCTCTGGAAGCATTCAGAGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTATTTAGTCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAACGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGTCTCCGACGGCAGGCTTCAATGACCCAGGCTGAGAAATTCCCGGACCCTTT  >  NC_012971/767731‑768183
                                                                                                                                                              |                                                                                                                                                                                                                                                                                                      
tgatgccccctcccaggagcaggttcgtgtaaaaactatgcttaatagCACCATTTCGATGGGTTACCCGGATGTTGTAATTGCATGTATAGAACATCAAGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACTAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtat                                                                                                                                                           <  2:404036‑M1/277‑92 (MQ=255)
  atgccccctcccaggagcaggttcgtgtaaaaactatgcttaatAGCACCATTTCGATGGGTTACCCGGATGTTGTAATTGCATGTATAGAACATCAAGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgcattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtattt                                                                                                                                                         <  1:501038‑M1/279‑94 (MQ=255)
  atgccccctcccaggagcaggttcgtgtaaaaacgatgcttaatAGCACCATTTCGATGGGTTACCCGGATGTTGTAATTGCATGTATAGAACATCAAGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaaggcctcgtattt                                                                                                                                                         <  1:499836‑M1/279‑94 (MQ=255)
          tcccatgagcaggttcgtgtaaaaactaTGCTTAATAGCACCATTTCGATGGGTTACCCGGATGTTGTAATTGCATGTATAGAACATCAAGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGTTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttagg                                                                                                                                                      <  2:148296‑M1/282‑97 (MQ=255)
          tcccatgagcaggttcgtgtaaaaactaTGCTTAATAGCACCATTTCGATGGGTTACCCGGATGTTGTAATTGCATGTATAGAACATCAAGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttagg                                                                                                                                                      <  1:21447‑M1/282‑97 (MQ=255)
          tcccaggagcaggttcgtgtaaaaactaTGCTTAATAGCACCATTTCGATGGGTTACCCGGATGTTGTAATTGCATGTATAGAACATCAAGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttagg                                                                                                                                                      <  1:18761‑M1/282‑97 (MQ=255)
           cccaggagcaggttcgtgtaaaaactATGCTTAATAGCACCATTTCGATGGGTTACCCGGATGTTGTAATTGCATGTATAGAACATCAAGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaatt                                                                                                                                                >  2:460623‑M1/13‑198 (MQ=255)
              aggagcaggttcgtgtaaaaACTATGCTTAATAGCACCATTTCGATGGGTTACCCGGATGTTGTAATTGCATGTATAGAACATCAAGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACTAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttagg                                                                                                                                                      <  1:205602‑M1/282‑97 (MQ=255)
              aggagcaagttcgtgtaaaaACTATGCTTAATAGCACCACTTCGATGTGTTACCCGGATGTTGTAATTGCATGTATAGAACATCAAGTGTCTCTGGAAGCATTCAGGGCAATTGATGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttagg                                                                                                                                                      <  2:148005‑M1/282‑97 (MQ=255)
                gagcaggttcgtgtaaAAACTATGCTTAATAGCACCATTTCGATGGGTTACCCGGATGTTGTAATTGCATGTATAGAACATCAAGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcc                                                                                                                                           <  1:596828‑M1/293‑108 (MQ=255)
                     ggttcgTTTAAAAACTATTCTTAATAGCACCATTTCGATGGGTTACCCGGATGTTGTAATTGCATGTATAGAACATCNAGTGTCTCTGGAAGCATTCAGTGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACTAAGGATTATTCCCTGGTGGTTGACTGATCACCAGAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggggaattgatcctacc                                                                                                                                       <  1:568971‑M1/297‑112 (MQ=255)
                     ggttcgTTTAAAAACTATGCTTAATAGCACCATTTCGATGGGTTACCCGGATGTTGTAATTGCATGTATAGAACATCAAGTGTCTCTGGAAGCATTCAGGGCAATTTAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctaccc                                                                                                                                      <  2:437921‑M1/298‑113 (MQ=255)
                               aaaacTATGCTTAATAGCACCATTTCGATGGGTTACCCGGATGTTGTAATTGCATGTATAGAACATCAAGTGTCTCTGNAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctacccacgtaatatg                                                                                                                            <  1:174434‑M1/300‑123 (MQ=255)
                                    tatGCTTAATAGCACCATTTCGATGGGTTACCCGGATGTTGTAATTGCATGTATAGAACATCAAGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctacccacttaatatggacac                                                                                                                       <  2:209808‑M1/300‑128 (MQ=255)
                                    tatGCTTAATAGCACCATCTCGATGGGTTACCCGGGTGTTGTAATTGCATGTATAGAACATCAAGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagttacagagtcaatatcgtattctgtcactggcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctacccacgtaatatggacac                                                                                                                       <  2:207919‑M1/300‑128 (MQ=255)
                                    tatGCTGAGTAGCACCATTTCGATGGGTTACCCGGAGTTTGTAATTGCATGTATAAAACATCAAGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctacccacgtaatatggacac                                                                                                                       <  2:208052‑M1/300‑128 (MQ=255)
                                       gCTTAATAGCACCATTTCGATGGGTTACCCGGATGTTGTAATTGCATGTATAGAACATCAAGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctacccacgtaatatggacacagg                                                                                                                    <  1:279725‑M1/300‑131 (MQ=255)
                                         ttAATAGCACCATTTCGATGGGTTACCCGGATGTTGTAATTGCATGTATAGAACATCAAGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctacccacgtaatatggacacaggcc                                                                                                                  <  2:544866‑M1/300‑133 (MQ=255)
                                         ttAATAGCACCATTTCGATGGGTTACCCGGATGTTGTAATTGCATGTATAGAACATCAAGTGTCTCTGGAAGCATTAAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctacccacgtaatatggacacaggc                                                                                                                   <  1:479984‑M1/299‑132 (MQ=255)
                                          tAATAGCACCATTTCGATGGGTTACCCGGATGTTGTAATTGCATGTATAGAACATCAAGTGTCTATGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCATGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctacccacgtaatatggacataggc                                                                                                                   >  2:421832‑M1/1‑167 (MQ=255)
                                              agcgccatttCGATGGGTTACCCGGATGGTGTAATTGCATGTATAGAACATCAAGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGAGAAGAATACGAAGGATTATTCCCTGGTGGGTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctacccacgtaatatggacacaggcctaagc                                                                                                             <  2:645857‑M1/296‑138 (MQ=255)
                                              aGCACCATTTTGATGTGTTACCCGGATGTTGTAATTGCATGTATAGAACATCAAGTGTCTCGGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAACACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctacccacgtaatatggacacaggcctaagc                                                                                                             <  2:148873‑M1/300‑138 (MQ=255)
                                              aGCACCATTTCGATGGGTTACCCGGATGTTGTAATTGCATGTATAGAACATCAAGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACTAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctacccacgtaatatggacacaggcctaagc                                                                                                             <  2:383109‑M1/300‑138 (MQ=255)
                                              aGCACCATTTCGATGGGTTACCCGGATGTTGTAATTGCATGTATAGAACATCAAGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACTAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctacccacgtaatatggacacaggcctaagc                                                                                                             <  1:624132‑M1/300‑138 (MQ=255)
                                              aGCACCATTTCGATGGGTTACCCGGATGTTGTAATTGCATGTATAGAACATCAAGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctacccacgtaatatggacacaggcctaagc                                                                                                             <  2:298042‑M1/300‑138 (MQ=255)
                                              aGCACCATTTCGATGGGTTACCCGGATGTTGTAATTGCATGTATAGAACATCAAGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAAAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtnaccatttcacaagtgacagagtcaatatcgtantctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctacccacgtaatatggacacaggcctaagc                                                                                                             <  2:592078‑M1/300‑138 (MQ=255)
                                              aGCACCATTTCGATGGGTTACCCGGAGTTTGTAATTGCATGTATAGAACATCAAGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctacccacgtaatatggacacaggcctaagc                                                                                                             <  2:756355‑M1/300‑138 (MQ=255)
                                               gCACCATTTCGATGGGTTACCCGGATGTTGTAATTGCATGTATAGAACATCAAGTGTCTCTGGAAGCATTCAGGGCAANTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctacccacgtaatatggacacaggcctaagcg                                                                                                            <  1:171360‑M1/300‑139 (MQ=255)
                                                            gggttACCCGGATGTTGTAATTGCATGTATAGAACATCAAGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctacccacgtaatatggacacaggcctaagcgaggttcttgtttt                                                                                               >  1:494028‑M1/3‑149 (MQ=255)
                                                              gTTACCCGGATGTTGTAATTGCATGTATAGAACATCAAGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACNGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctacccacgtaatatggacacaggcataagcgaggttcttgttttca                                                                                             >  2:572251‑M1/1‑147 (MQ=255)
                                                                        tgttgtAATTGCATGTATAGAACATCAAGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctacccacgtaatatggacacaggcataagcgaggttcttgttttcaaattgttcc                                                                                    >  2:456343‑M1/1‑137 (MQ=255)
                                                                        tgttgtAACTGCATGTATAGAACACCAAGTGTCTATGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaacatttcacaagtgacagagtcaatatcgtattctgtcactgccaggataatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctacccacgtaatatggacacaggcctaagcgaggttcttgttttcaaattgttc                                                                                     >  2:9357‑M1/1‑137 (MQ=255)
                                                                         gttgtAATTGCATGTATAGAACATCAAGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctacccacgtaatatggacacaggcctaagcgaggttcttgttttcaaattgttcc                                                                                    >  1:600610‑M1/1‑136 (MQ=255)
                                                                          ttgtAATTGCATGTATAGAACATCAAGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctagtatttaggtgaattgatcctacccacgtaatatggacacaggcctaagcgaggttcttgttttcaaattgttccgga                                                                                 >  1:144966‑M1/1‑135 (MQ=255)
                                                                          ttgtAATTGCATGTATAGAACATCAAGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCngaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctacccacgtaatatggacacaggcctaagcgaggttcttgttttcaaattgttccgga                                                                                 >  2:568542‑M1/1‑135 (MQ=255)
                                                                          ttgtAATTGCATGTATAGAACATCAAGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCATGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactaaagcaatgcctagtatttaggtgaattgatcctacccacgtaatatggacacaggcctaagcgaggttcttgttttcaaattgttccgga                                                                                 >  2:145225‑M1/1‑135 (MQ=255)
                                                                           tgtAATTGCATGTATAGAACATCAAGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTNtgaccacttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctacccacgtaatatggacacaggcctaagcgaggttcttgttttcaaattgttccggac                                                                                <  2:259080‑M1/300‑167 (MQ=255)
                                                                                ttGCATGTATAGAACATCAAGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATANGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctacccacgtaatatggacacaggcctaagcgaggttcttgttttcaaattgttccggactgagg                                                                           <  1:149703‑M1/300‑172 (MQ=255)
                                                                                   cATGTATAGAACATCAAGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaattcaattcaactacagcaatgcctagtatttaggtgaattgatcctacccacgtaatatggacacaggcctaagcgaggttcttgttttcaaattgttccggactgaggc                                                                          >  1:317393‑M1/1‑126 (MQ=255)
                                                                                   cATGTATAGAACATCAAGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaattcaattcaactaaagcaatgcctcgtatttaggtgaattgatcctaaccacgtaatatggacacaggcctaagcgaggttcttgttttcaaattgttcaggac                                                                                >  1:317398‑M1/1‑126 (MQ=255)
                                                                                      gTATAGAACATCAAGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctagtatttaggtgaattgatcctaaccacgtaatatggacacaggcctaagcgaggtttttgttttcaaattgttccggactgaggacgcca                                                                     >  2:177427‑M1/1‑123 (MQ=255)
                                                                                      gTATAGAACATCAAGTGTCTATGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAAATGCTAAACATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatctta                                                                                                                                         >  2:34168‑M1/1‑123 (MQ=255)
                                                                                       tataGAACATCAAGTGTCTTTGGAAGCATTCAGGGCAATTGAGGCCGCATTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctacccacgtaatatggacacaggcctaagcgaggttcttgttttcaaattgttccggactgaggccgccac                                                                    <  2:88922‑M1/300‑179 (MQ=255)
                                                                                       tataGAACATCAAGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctacccacgtaatatggacacaggcctaagcgaggttcttgttttcaaattgttccggactgaggccgccac                                                                    <  1:541273‑M1/300‑179 (MQ=255)
                                                                                       tataGAACATCAAGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctacccacgtaatatggacacaggcctaagcgaggttcttgttttcaaattgttccggactgaggccgccac                                                                    <  2:88292‑M1/300‑179 (MQ=255)
                                                                                         taGAACATCAAGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACTAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctacccacgtaatatggacacaggcctaagcgaggttcttgttttcaaattgttccggactgaggccgccacac                                                                  <  1:381305‑M1/300‑181 (MQ=255)
                                                                                         taGAACATCAAGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaattgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctacccacgtaatatggacacaggcctaagcgaggttcttgttttcaaattgttccggactgaggccgccacac                                                                  <  1:382937‑M1/300‑181 (MQ=255)
                                                                                         taGAACATCAAGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctacccacgtaatatggacacaggcctaagcgaggttcttgttttcaaattgttccggactgaggccgccacac                                                                  <  1:381221‑M1/300‑181 (MQ=255)
                                                                                           gAACATCAAGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctacccacgtaatatggacacaggcctaagcgaggttcttgttttcaaattgttccggactgaggccgccacacca                                                                >  2:689352‑M1/1‑118 (MQ=255)
                                                                                             acatcaagtgtctctgtaAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgncagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctacccacgtaatatggacacaggcctaagcgaggttcttgttttcaaattgttccggactgaggccgccacaccaac                                                              <  2:260792‑M1/293‑185 (MQ=255)
                                                                                             acatcaagtgtctctggnAGCATTCAGGGCAATTGANGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctacccacgtaatatggacacaggcctaagcgaggttcttgttttcaaattgttccggactnaggccgccacaccaac                                                              <  1:678106‑M1/293‑185 (MQ=255)
                                                                                               atcaggtGTCTCTGGAAGCAGTTAGGGCAATTTAGGCCGCGTTAGTAAAGCACGATAAGAATACTAAGGATTATTCCTTGGTGGTTGACTGATCACTATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctacccacgtaatatggacacaggcctaagcgaggttcttgttttcaaattgttccggactgaggccgccacaccaactg                                                            <  2:521313‑M1/297‑187 (MQ=255)
                                                                                                        ctctGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctacccacgtaatatggacacaggcctaagcgaggttcttgttttcaaattgttccggactgaggccgccacaccaactgtgccgccgc                                                   >  1:24949‑M1/1‑105 (MQ=255)
                                                                                                                gCATTCAGAGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctagtatttaggtgaattgatcctacccacgtaatatggacacaggactaagcgagggtcttgttttcaagttgttcaggaatgaggccgccacaccaaatgtgccgccgcaaccgat                                            >  2:221390‑M1/1‑97 (MQ=255)
                                                                                                                       gggcaATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaatttaactacagcaatgcatcgtatttaggtggattgatcctacccacgtaatatggacacaggcctacgcgaggttcttgttttcaaattgttccggactgaggccgcccaaccaactgtgccgacgccacagattgtaacc                                     >  2:75070‑M1/3‑90 (MQ=255)
                                                                                                                       gggcaATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagccatgcctcgtatttaggtgaattgatcctacccacgtaatatggacacaggcctaagcgaggttcttgttttcacattgttccggactgaggccgacgcacccactgtgccgacgccaccgattgtaatca                                    >  2:150274‑M1/3‑90 (MQ=255)
                                                                                                                                aggccacgttattaaagcacgataagaagaCGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctacccacgtaatatggacacaggcctaagcgaggttcttgttttcaaattgttccggactgaggccgccacaccaactgtgccgccgccaccgattgtaatcacattcgata                           <  2:394090‑M1/286‑220 (MQ=255)
                                                                                                                                     gCGTTAGTAAAGCACGAGAAGACTACGAAGGATTATTCCCTGGTGGTGGACTGATCACCAGAACTGCTAATCATCCtgaccatttcacaagtgacagattcaatatcgtagtctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctacccacggaatatggacacaggcctaagcgaggttcttgttttcaaattgttccggactgaggccgccacaccaactgtgccgccgccaccgattgtaatcacattcgatataat                       <  2:81359‑M1/299‑224 (MQ=255)
                                                                                                                                      cgttagtaaagcacgataAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctacccacgtaatatggacacaggcctaagcgaggttcttnttttcaaattgttccggactgaggccgccacaccaactgtgccgccgccaccgattgtaatcacattcgatataatta                     >  1:576944‑M1/9‑75 (MQ=255)
                                                                                                                                      cgttagtaaagcacgataAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctagtatttaggtgaattgatcctacccacgtaatatggacacaggcctaagcgaggttcttgttttcaaattgttccggactgaggccgccacaccaactgtgccgccgccaccgattgtactcacattcgatataatta                     >  2:139365‑M1/9‑75 (MQ=255)
                                                                                                                                      cgttagtaaagcacgataAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcatagtatttaggtgaattgatcctacccacgtaatatggacacaggcctaagcgaggttcttgttttcaaattgttccggactgaagccgccacaccaactgtgccgccgccaccgattgtaatcacattcgatataatt                      >  2:764713‑M1/9‑75 (MQ=255)
                                                                                                                                      cgttagtaaagcacgataAGAATACGAAGGATTATTCCCTGGTGGTTGAATGATCACCATAACTGCTAAACATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctagtatttaggtgaattgatcctaaccacgtaatatggacccaggcctaagcgaggtttttgttttcaacttgttccggactgaggccgcaacaccaactgtgccgccgccaccgattgtaatcacattcgatataatta                     >  2:727084‑M1/9‑75 (MQ=255)
                                                                                                                                       gttagtaaagcacgatAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctacccacgtaatatggacacaggcctaagcgaggttcttgttttcaaattgttccggactgaggccgccacaccaactgtgccgccgccaccgattgtaatcacattcgatataatta                     <  1:188774‑M1/292‑226 (MQ=255)
                                                                                                                                       gttagtaaagcacgatAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctacccacgtaatatggacacaggcctaagcgaggttcttgttttcaaattgttccggactgaggccgccacaccaactgtgccgccgccaccgattgtaatcacattcgatataat                       <  2:82151‑M1/290‑224 (MQ=255)
                                                                                                                                          agtaaagcacGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctacccacgtaatatggacacaggcctaagcgaggttcttgttttcaaattgttccggactgaggccgccacaccaactgtgccgccgccaccgattgtaatcacattcgatataattaaac                  >  2:791175‑M1/5‑71 (MQ=255)
                                                                                                                                           gtaaagcaCGATAAGAATACGAAGGATTATTCCCTGTTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctacccacgtaatatggacacaggcctaagcgaggttcttgttttcaaattgttccggactgaggccgccacaccaactgtgccgccgccaccgattgtaatcacattcgatataattaaacac                <  1:27018‑M1/297‑231 (MQ=255)
                                                                                                                                             aaagCACGATAAGAATACTAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctacccacgtaatatggacacaggcctaagcgaggttcttgttttcaaattgttccggactgaggccgccacatcaactgtgccgccgccaccgattgtaatcacattcgatataattaaacaccg              <  1:103834‑M1/299‑233 (MQ=255)
                                                                                                                                             aaagCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctacccacgtaatatggacacaggcctaagcgaggttcttgttttcaaattgttccggactgaggccgccacaccaactgtgccgccgccaccgattgtaatcacattcgatataattaaacaccg              <  1:528025‑M1/299‑233 (MQ=255)
                                                                                                                                                   cGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaataaagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctacccacgtaatatggacacaggcctaagcgaggttcttgttttcaaattgttccggactgaggccgccacaccaactgtgccgccgccaacgattgtaatcacattcgatataattaaacaccgttgccc        >  2:185849‑M1/1‑62 (MQ=255)
                                                                                                                                                         gaatacGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctacccacgtaatatggacacaggcctaagcgaggttcttgttttcaaattgttccggactgaggccgccacaccaactgtgccgccgccaccgattgtaatcacattcgatataattaaacaccgttgcccgcatta  >  1:232659‑M1/2‑56 (MQ=255)
                                                                                                                                                              |                                                                                                                                                                                                                                                                                                      
TGATGATCCCTCCGTGGATCTGATTCGTGTAAAAAATATGCTTAATAGCACCATTTCTATGAGTTACCCTGATGTTGTAATTGCATGTATAGAACATAAGGTGTCTCTGGAAGCATTCAGAGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTATTTAGTCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAACGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGTCTCCGACGGCAGGCTTCAATGACCCAGGCTGAGAAATTCCCGGACCCTTT  >  NC_012971/767731‑768183

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: